计算线段的交点

/*获得X1,Y1,X2,Y2的交点，无交点时返回FALSE,Flag为1时交点必须位于两个矩形交叉区域内，Flag为2时交点为广义交点，即直线的延长线有效
注意：Flag为2时前四个参数一定为剪切线
参数为3时，求出的交点不被范围限制 }
*/
bool CGeoUtils::getIntersectedPoint(double X1,double Y1,double X2,double Y2,double X3,double Y3,double X4,double Y4, double &XX,double &YY, int nFlag)
{
 double K1,K2,B1,B2;  //二条直线的斜率}
    bool V1,V2;  //有斜率为无穷大的线存在}

 //X1,Y1,X2,Y2 为第一条直线,X3,Y3,X4,Y4 为第二条直线}
 V1 = V2 = false;

    K1 = DBL_MAX;  //
    K2 = DBL_MAX;
    B1 = K1;
    B2 = K2;

 if (abs(X2-X1)>10e-10){
  K1 = (Y2-Y1)/(X2-X1);
  B1 = Y1-K1*X1;
 }
 else {
     XX = X1;
  V1 = true;
 };

 if (abs(X4-X3) > 10e-10){
     K2 = (Y4-Y3)/(X4-X3);
  B2 = Y3-K2*X3;
 }
 else{
     XX = X3;
        V2 = true;
 }
 //如果两直线斜率相等，则退出
 if (abs(K1 - K2) < 10e-10 || (V1 && V2)) return false;

 if (!V1 && ! V2){
         XX = (B1-B2)/(K2-K1);
         YY = K1*XX+B1;
 }
 else if( V2){
  YY = K1*XX+B1;
 }
 else {
     YY = K2*XX+B2;
 }

 switch (nFlag){
  case 1 : {
  //如果计算出来的交点位于两个矩形的公共区域，则为真交点}
   if (isPtInRect(XX,YY,X1,Y1,X2,Y2) && isPtInRect(XX,YY,X3,Y3,X4,Y4))
       return true;
      else 
       return false;
   break;
   }
  case 2 : {
   if (isPtInRect(XX,YY,X3,Y3,X4,Y4))
       return true;
      else
       return false;
   break;
   }
  case 3 : {
      return true;
   break;
   }
  }
 
}