331 double a1, a2, b1, b2, c1, c2, d1, d2;
332 double A1, A2, B1, B2, C1, C2, D1, D2;
347 A1 = cMo[0][0] * a1 + cMo[0][1] * b1 + cMo[0][2] * c1;
348 B1 = cMo[1][0] * a1 + cMo[1][1] * b1 + cMo[1][2] * c1;
349 C1 = cMo[2][0] * a1 + cMo[2][1] * b1 + cMo[2][2] * c1;
350 D1 = d1 - (cMo[0][3] * A1 + cMo[1][3] * B1 + cMo[2][3] * C1);
352 A2 = cMo[0][0] * a2 + cMo[0][1] * b2 + cMo[0][2] * c2;
353 B2 = cMo[1][0] * a2 + cMo[1][1] * b2 + cMo[1][2] * c2;
354 C2 = cMo[2][0] * a2 + cMo[2][1] * b2 + cMo[2][2] * c2;
355 D2 = d2 - (cMo[0][3] * A2 + cMo[1][3] * B2 + cMo[2][3] * C2);
367 a2 = B1 * C2 - C1 * B2;
368 b2 = C1 * A2 - A1 * C2;
369 c2 = A1 * B2 - B1 * A2;
372 a1 = A2 * D1 - A1 * D2;
373 b1 = B2 * D1 - B1 * D2;
374 c1 = C2 * D1 - C1 * D2;
376 if (fabs(D2) < fabs(D1))
385 d1 = 1.0 / sqrt(a1 * a1 + b1 * b1 + c1 * c1);
386 cP_[0] = A1 = a1 * d1;
387 cP_[1] = B1 = b1 * d1;
388 cP_[2] = C1 = c1 * d1;
393 a1 = b2 * C1 - c2 * B1;
394 b1 = c2 * A1 - a2 * C1;
395 c1 = a2 * B1 - b2 * A1;
398 d1 = 1.0 / sqrt(a1 * a1 + b1 * b1 + c1 * c1);
404 D2 /= (A2 * a1 + B2 * b1 + C2 * c1);
static void displayLine(double rho, double theta, const vpCameraParameters &cam, const vpImage< unsigned char > &I, const vpColor &color=vpColor::green, unsigned int thickness=1)
void setWorldCoordinates(const double &oA1, const double &oB1, const double &oC1, const double &oD1, const double &oA2, const double &oB2, const double &oC2, const double &oD2)