2 This file is part of I am car.
4 I am car is free software: you can redistribute it and/or modify
5 it under the terms of the GNU General Public License as published by
6 the Free Software Foundation, either version 3 of the License, or
7 (at your option) any later version.
9 I am car is distributed in the hope that it will be useful,
10 but WITHOUT ANY WARRANTY; without even the implied warranty of
11 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12 GNU General Public License for more details.
14 You should have received a copy of the GNU General Public License
15 along with I am car. If not, see <http://www.gnu.org/licenses/>.
22 #include "slotplanner.h"
24 ParallelSlot::ParallelSlot()
28 std::vector<std::vector<RRTNode *>> &ParallelSlot::cusp()
33 float ParallelSlot::DH() const
38 PolygonObstacle &ParallelSlot::slot()
43 float ParallelSlot::slotHeading()
45 float y0 = this->slot().bnodes()[0]->y();
46 float x0 = this->slot().bnodes()[0]->x();
47 float y3 = this->slot().bnodes()[3]->y();
48 float x3 = this->slot().bnodes()[3]->x();
54 SlotSide ParallelSlot::slotSide()
56 if (!this->slotSide_) {
57 float y0 = this->slot().bnodes()[0]->y();
58 float x0 = this->slot().bnodes()[0]->x();
59 float y1 = this->slot().bnodes()[1]->y();
60 float x1 = this->slot().bnodes()[1]->x();
61 float y3 = this->slot().bnodes()[3]->y();
62 float x3 = this->slot().bnodes()[3]->x();
63 float h1 = atan2(y0 - y3, x0 - x3);
64 float h2 = atan2(y1 - y3, x1 - x3);
66 this->slotSide_ = LEFT;
68 this->slotSide_ = RIGHT;
70 return this->slotSide_;
73 SlotType ParallelSlot::slotType()
75 if (!this->slotType_) {
77 this->slot().bnodes()[0],
78 this->slot().bnodes()[1]
81 this->slot().bnodes()[1],
82 this->slot().bnodes()[2]
85 this->slotType_ = PERPENDICULAR;
87 this->slotType_ = PARALLEL;
89 return this->slotType_;
93 void ParallelSlot::DH(float dh)
99 void ParallelSlot::fip()
101 // see https://courses.cs.washington.edu/courses/cse326/03su/homework/hw3/bfs.html
102 // RRTNode.s() works as iteration level
103 std::queue<BicycleCar *, std::list<BicycleCar *>> q;
105 // TODO add init nodes
106 // for now just copy fpose()
107 bool left = false; // right parking slot
109 // new pose for parallel parking to right slot
114 // temporary tnx is angle
115 tnx = this->slotHeading() + M_PI;
116 if (this->slotSide() == RIGHT)
120 nx = this->fposecenter()->x() + 0.01 * cos(tnx);
121 ny = this->fposecenter()->y() + 0.01 * sin(tnx);
122 BicycleCar *CC = new BicycleCar(nx, ny, this->slotHeading());
123 // move left by car width / 2
124 tnx = CC->x() + CC->width() / 2 * cos(CC->h() + M_PI / 2);
125 tny = CC->y() + CC->width() / 2 * sin(CC->h() + M_PI / 2);
127 nx = tnx - (CC->length() + CC->wheelbase()) / 2 * cos(CC->h());
128 ny = tny - (CC->length() + CC->wheelbase()) / 2 * sin(CC->h());
129 if (this->slotSide() == LEFT) {
130 std::cerr << "left PS" << std::endl;
133 // move right by car width / 2
134 tnx = CC->x() + CC->width() / 2 * cos(CC->h() - M_PI / 2);
135 tny = CC->y() + CC->width() / 2 * sin(CC->h() - M_PI / 2);
137 nx = tnx - (CC->length() + CC->wheelbase()) / 2 * cos(CC->h());
138 ny = tny - (CC->length() + CC->wheelbase()) / 2 * sin(CC->h());
140 BicycleCar *B = new BicycleCar(nx, ny, CC->h());
141 this->DH(di * 0.01 / CC->out_radi());
144 c = B->move(CC, -i * di * 0.01 / CC->diag_radi());
145 while (!this->slot().collide(c->frame())) {
147 c = B->move(CC, -i * di * 0.01 / CC->diag_radi());
150 delete c; // not in q and collide
155 if (this->DH() > 0 && c->rfx() <= 0 && c->rrx() <= 0) {
156 goto createcuspandfinish;
157 } else if (this->DH() < 0 && c->lfx() >= 0 && c->lrx() >= 0) {
158 goto createcuspandfinish;
159 } else if (c->s() < 9) {
160 BicycleCar *cc = this->flnc(c);
165 delete c; // not in q and collide
169 std::vector<RRTNode *> cusp;
171 cusp.push_back(new RRTNode(c->x(), c->y(), c->h()));
174 std::reverse(cusp.begin(), cusp.end());
175 this->cusp().push_back(cusp);
176 std::queue<BicycleCar *, std::list<BicycleCar *>> empty;
180 void ParallelSlot::fipr(BicycleCar *B)
182 // TODO for right parallel parking also
183 // it's only for lpar scenario now
185 std::vector<RRTNode *> cusp;
186 cusp.push_back(new RRTNode(B->x(), B->y(), B->h()));
187 // just copied from fip()
188 this->DH(-0.01 / B->out_radi());
191 while (c->lfx() < 0) {
192 cusp.push_back(new RRTNode(c->x(), c->y(), c->h()));
193 BicycleCar *cc = this->flncr(c);
198 cusp.push_back(new RRTNode(c->x(), c->y(), c->h()));
199 std::reverse(cusp.begin(), cusp.end());
200 this->cusp().push_back(cusp);
203 BicycleCar *ParallelSlot::flnc(BicycleCar *B)
205 // TODO find last not colliding
206 // for now just copy flast()
208 if (int(B->s()) % 2 == 0)
209 cc = BicycleCar(B->x(), B->y(), B->h()).ccr();
211 cc = BicycleCar(B->x(), B->y(), B->h()).ccl();
214 p = B->move(cc, i * this->DH());
215 while (!this->slot().collide(p->frame())
217 (this->DH() > 0 && p->x() <= 0)
218 || (this->DH() < 0 && p->x() >= 0)
222 p = B->move(cc, i * this->DH());
225 p = B->move(cc, i * this->DH());
226 while (!this->slot().collide(p->frame())
228 (this->DH() > 0 && p->x() <= 0)
229 || (this->DH() < 0 && p->x() >= 0)
231 if (this->DH() > 0 && p->rfx() <= 0 && p->rrx() <= 0) {
235 if (this->DH() < 0 && p->lfx() >= 0 && p->lrx() >= 0) {
241 p = B->move(cc, i * this->DH());
244 return B->move(cc, (i - 1) * this->DH());
247 BicycleCar *ParallelSlot::flncr(BicycleCar *B)
249 // TODO find last not colliding
250 // for now just copy flast()
252 if (int(B->s()) % 2 == 0)
253 cc = BicycleCar(B->x(), B->y(), B->h()).ccr();
255 cc = BicycleCar(B->x(), B->y(), B->h()).ccl();
258 p = B->move(cc, i * this->DH());
259 while (!this->slot().collide(p->frame())
261 (this->DH() > 0 && p->x() >= 0)
262 || (this->DH() < 0 && p->lfx() <= 0)
266 p = B->move(cc, i * this->DH());
269 p = B->move(cc, i * this->DH());
270 while (!this->slot().collide(p->frame())) {
271 if (this->DH() > 0 && p->x() <= 0) {
275 if (this->DH() < 0 && p->lfx() >= 0) {
281 p = B->move(cc, i * this->DH());
284 return B->move(cc, (i - 1) * this->DH());
287 RRTNode *ParallelSlot::fposecenter()
289 return this->slot().bnodes().front();
292 bool ParallelSlot::flast(
296 std::vector<RRTNode *> &cusp
299 BicycleCar *B = new BicycleCar(P->x(), P->y(), P->h());
302 cc = BicycleCar(B->x(), B->y(), B->h()).ccr();
304 cc = BicycleCar(B->x(), B->y(), B->h()).ccl();
307 p = B->move(cc, i * this->DH());
308 while (!this->slot().collide(p->frame())
310 (this->DH() > 0 && p->x() <= 0)
311 || (this->DH() < 0 && p->x() >= 0)
315 p = B->move(cc, i * this->DH());
318 p = B->move(cc, i * this->DH());
319 while (!this->slot().collide(p->frame())
321 (this->DH() > 0 && p->x() <= 0)
322 || (this->DH() < 0 && p->x() >= 0)
324 if (this->DH() > 0 && p->rfx() <= 0 && p->rrx() <= 0) {
328 if (this->DH() < 0 && p->lfx() >= 0 && p->lrx() >= 0) {
334 p = B->move(cc, i * this->DH());
337 p = B->move(cc, (i - 1) * this->DH());
338 if (this->DH() > 0 && p->rfx() <= 0 && p->rrx() <= 0) {
341 } else if (this->DH() < 0 && p->lfx() >= 0 && p->lrx() >= 0) {
346 return this->flast(p, !right, il + 1, cusp);
351 void ParallelSlot::fpose()
353 bool left = false; // right parking slot
355 BicycleCar *CC = new BicycleCar(
356 this->fposecenter()->x(),
357 this->fposecenter()->y() - 0.01,
360 BicycleCar *B = new BicycleCar(
361 CC->x() - CC->width() / 2,
362 CC->y() - (CC->length() + CC->wheelbase()) / 2,
365 if (this->slot().bnodes()[0]->x() > this->slot().bnodes()[1]->x()) {
370 CC->x() + CC->width() / 2,
371 CC->y() - (CC->length() + CC->wheelbase()) / 2,
375 this->DH(di * 0.01 / CC->out_radi());
378 p = B->move(CC, -i * di * 0.01 / CC->diag_radi());
379 while (!this->slot().collide(p->frame())) {
380 std::vector<RRTNode *> tmpcusp;
381 tmpcusp.push_back(new BicycleCar(p->x(), p->y(), p->h()));
382 if (this->flast(p, left, 0, tmpcusp)) {
383 this->cusp().push_back(tmpcusp);
388 p = B->move(CC, -i * di * 0.01 / CC->diag_radi());
392 struct SamplingInfo ParallelSlot::getSamplingInfo()
394 struct SamplingInfo si;
395 BicycleCar *CC = new BicycleCar(
396 this->fposecenter()->x(),
397 this->fposecenter()->y() - 0.01,
400 si.x = this->slot().bnodes()[0]->x();
401 si.y = this->slot().bnodes()[0]->y();
402 si.r = CC->diag_radi();
403 si.h = this->slotHeading() - acos(EDIST( // TODO generalize
404 this->slot().bnodes()[0],
405 this->slot().bnodes()[1]