当前位置：首页 > news >正文

计算斜着椭圆内某个点到边距离（验证ok）

news 来源：原创 2025/5/9 22:03:06

       sPoint center;
       center.x = sLRpInfo.rdOut.x; center.y = sLRpInfo.rdOut.y;
       float a = sLRpInfo.rdOut.a; // 长轴半径
       float b = sLRpInfo.rdOut.b; // 短轴半径
       float theta = (360.0 - sLRpInfo.rdOut.fAngle) * CV_PI / 180.0f; // 旋转角度，转换为弧度
       float fddtoCen = 0;
       for (i = 0;i < nDftNum;i++)
       {
           rdFit = A_RoundFitting(v_FitRound.at(i), 0, v_FitRound.at(i).size());
           rdFit.x += roi.l; rdFit.y += roi.t;
           sLRpInfo.rdFeaFit[i] = rdFit;
           if (2 * rdFit.r > sLRpInfo.fGetMaxMinRd[1])sLRpInfo.fGetMaxMinRd[1] = 2 * rdFit.r;
           if (2 * rdFit.r < sLRpInfo.fGetMaxMinRd[0])sLRpInfo.fGetMaxMinRd[0] = 2 * rdFit.r;

           if (sLRpInfo.bArith[5])
           {
               float x_prime = (rdFit.x - center.x) * cos(theta) + (rdFit.y - center.y) * sin(theta);
               float y_prime = -(rdFit.x - center.x) * sin(theta) + (rdFit.y - center.y) * cos(theta);

               // 2. 计算参数角 alpha
               float alpha = atan2(y_prime, x_prime);
               if (alpha < 0) alpha += 2 * CV_PI; // 确保 alpha 在 0 到 2π 之间

               // 3. 计算椭圆边界点
               float x_ellipse = a * cos(alpha);
               float y_ellipse = b * sin(alpha);

               // 4. 旋转回全局坐标系
               sLRpInfo.ptOffArcEdge[i].x = center.x + x_ellipse * cos(theta) - y_ellipse * sin(theta);
               sLRpInfo.ptOffArcEdge[i].y = center.y + x_ellipse * sin(theta) + y_ellipse * cos(theta);

// 计算边界点到中心的距离
fddtoCen = A_GetPointDistence(sLRpInfo.ptOffArcEdge[i], center);

               if (fddtoCen > sLRpInfo.fGetMaxMinToArc[1])
                   sLRpInfo.fGetMaxMinToArc[1] = fddtoCen;
               if (fddtoCen < sLRpInfo.fGetMaxMinToArc[0])
                   sLRpInfo.fGetMaxMinToArc[0] = fddtoCen;
           }
       }