[建议]：缓和曲线类

Lagier

:9
曲线放样是有数学模型的，根据数学模型来计算坐标，下面提供一个类的定义，供参考。
/***
*Engin.hpp - definitions and declarations for additional engineering library
*
*       Copyright (c) 2000-2001, By Liurenzhao. All rights reserved.
*
*Purpose:
*       This file contains constant definitions and external subroutine
*       declarations for the additional math subroutine library.
*
*       [Public]
*
****/

#ifndef _INC_ENGIN_HPP
#define _INC_ENGIN_HPP

///////////////////////////////////
class CircleRoute      //圆曲线测设
///////////////////////////////////
{  
private:      
     double radius;  //半径(米)
     double angle;   //偏角(弧度)
         char *JDstake;  //交点桩号 例如: JDstake[]="K87+441.41"
         double step;
         bool  flag;     //左偏为true,右偏为false

public:
     CircleRoute(double r,double a, char *stake)
         {   radius=r; angle=fabs(a);
             if(a>=0) flag=true;
                 else     flag=false;
             JDstake=new char[strlen(stake)+1];
             strcpy(JDstake,_strupr(strdup(stake)));
         };
     CircleRoute(double r,double a, double l=10.0)
         {   radius=r; angle=fabs(a); step=l;
             if(a>=0) flag=true;
                 else     flag=false;             
         };
     CircleRoute(){};
     ~CircleRoute() { delete JDstake; };

         void setRadius(double);
         void setAngle(double);
         void setJDstake(char *);
         void setStep(double);

         double getRadius(void);
         double getAngle(void);
         double getStep(void);
         bool   getFlag(void);

         double getT(void); //切线长
         double getL(void); //曲线长
         double getE(void); //外矢距
         double getq(void); //切曲差

         char* getJDstake(void);
     char* getZYstake(void);
     char* getYZstake(void);
     char* getQZstake(void);

     bool getCoordinate(AcGePoint3dArray& point3dArray);
} ;

inline void CircleRoute::setRadius(double r) { radius=r; }
inline void CircleRoute::setAngle(double a)
{         angle=fabs(a);
         if(a>=0) flag=true;
         else     flag=false;
}
inline void CircleRoute::setJDstake(char *stake)
{   JDstake=new char[strlen(stake)+1];
    strcpy(JDstake,_strupr(strdup(stake)));
};
inline void CircleRoute::setStep(double l)   { step=l; }

inline double CircleRoute::getRadius(void) { return(radius); }
inline double CircleRoute::getAngle(void)  { return(angle);  }
inline double CircleRoute::getStep(void)   { return(step);   }
inline bool   CircleRoute::getFlag(void)   {  return(flag);  }
///////////////////////////////////////////

inline double CircleRoute::getT(void) {        return(radius*tan(angle*0.5));  }
inline double CircleRoute::getL(void) {        return(angle*radius);  }
inline double CircleRoute::getE(void) {        return(radius/cos(angle*0.5)-radius); }
inline double CircleRoute::getq(void) {        return(2*this->getT()-this->getL()); }
////////////////////////////////////////////

inline char* CircleRoute::getJDstake(void) { return(JDstake); }

inline char* CircleRoute::getZYstake(void)
{        char *ZYstake,strkm[6],strm[10],sk[2][32],ss[2][32];
    double zy_km,zy_m,zy,jd;
        ZYstake=new char[strlen(JDstake)+1]; ZYstake[0]='\0';
    rzrxstrlist(JDstake,'K',sk); rzrxstrlist(sk[1],'+',ss);
    jd=atof(ss[0])*1000+atof(ss[1]);
    zy=jd-this->getT();
    zy_km=(int)(zy/1000);  zy_m=zy-zy_km*1000;
        sprintf(strkm,"%.0lf",zy_km); sprintf(strm, "%.3lf",zy_m );
        strcat(ZYstake,"K"); strcat(ZYstake,strkm);
        strcat(ZYstake,"+"); strcat(ZYstake,strm);
        return(ZYstake);
}

inline char* CircleRoute::getYZstake(void)
{        char *YZstake,strkm[5],strm[10],sk[2][32],ss[2][32];
    double yz_km,yz_m,yz,jd;
        YZstake=new char[strlen(JDstake)+1]; YZstake[0]='\0';
    rzrxstrlist(JDstake,'K',sk); rzrxstrlist(sk[1],'+',ss);
    jd=atof(ss[0])*1000+atof(ss[1]);
    yz=jd-this->getT()+this->getL();
    yz_km=(int)(yz/1000);  yz_m=yz-yz_km*1000;
        sprintf(strkm,"%.0lf",yz_km); sprintf(strm, "%.3lf",yz_m );
        strcat(YZstake,"K"); strcat(YZstake,strkm);
        strcat(YZstake,"+"); strcat(YZstake,strm);
        return(YZstake);
}

inline char* CircleRoute::getQZstake(void)
{        char *QZstake,strkm[5],strm[10],sk[2][32],ss[2][32];
    double qz_km,qz_m,qz,jd;
        QZstake=new char[strlen(JDstake)+1]; QZstake[0]='\0';
    rzrxstrlist(JDstake,'K',sk); rzrxstrlist(sk[1],'+',ss);
    jd=atof(ss[0])*1000+atof(ss[1]);   
        qz=jd-this->getT()+this->getL()*0.5;
    qz_km=(int)(qz/1000);  qz_m=qz-qz_km*1000;
        sprintf(strkm,"%.0lf",qz_km); sprintf(strm, "%.3lf",qz_m );
        strcat(QZstake,"K"); strcat(QZstake,strkm);
        strcat(QZstake,"+"); strcat(QZstake,strm);
        return(QZstake);
}


inline bool CircleRoute::getCoordinate(AcGePoint3dArray& point3dArray)
{  int num,i;
   double xi,yi,arc;
   AcGePoint3d point3d;

   arc=step/radius; num=(int)(angle/arc);
   for(i=0;i<num;i++) {
             xi=radius*sin(i*arc);
      if(flag==true) yi=radius*(1-cos(i*arc));
          else           yi=-radius*(1-cos(i*arc));
      point3d.set(xi,yi,0);
          point3dArray.append(point3d);
   }
   xi=radius*sin(angle);
   if(flag==true) yi=radius*(1-cos(angle));
   else           yi=-radius*(1-cos(angle));
   point3d.set(xi,yi,0); point3dArray.append(point3d);
   if(point3dArray.length() < 3) return false;
   else                          return true;
}



///////////////////////////////////////
class AllayCircleRoute   //缓和曲线测设
///////////////////////////////////////
{  
private:      
     double radius;  //半径(米)
     double angle;   //偏角(弧度)
         char *JDstake;  //交点桩号 例如: JDstake[]="K87+441.41"
         double step;
         double allay;   //缓和曲线长度
         bool  flag;     //左偏为true,右偏为false

public:
     AllayCircleRoute(double r, double a, double L0, char *stake)
         {   radius=r; angle=fabs(a); allay=L0;
             if(a>=0) flag=true;
                 else     flag=false;
             JDstake=new char[strlen(stake)+1];
             strcpy(JDstake,_strupr(strdup(stake)));
         };
     AllayCircleRoute(double r, double a, double L0, double l=10.0)
         {   radius=r; angle=fabs(a); allay=L0; step=l;
             if(a>=0) flag=true;
                 else     flag=false;
         };
     AllayCircleRoute(){};
     ~AllayCircleRoute() { delete JDstake; };

         void setRadius(double);
         void setAngle(double);
         void setAllay(double);
         void setJDstake(char *);
         void setStep(double);

         double getRadius(void);
         double getAngle(void);
         double getAllay(void);
         double getStep(void);
         bool   getFlag(void);

         double getT(void); //切线长
         double getL(void); //曲线长
         double getE(void); //外矢距
         double getq(void); //切曲差

         char* getJDstake(void);
     char* getZHstake(void);
     char* getHYstake(void);
     char* getQZstake(void);
     char* getZYstake(void);
     char* getYHstake(void);
     char* getHZstake(void);

     bool getCoordinate(AcGePoint3dArray& point3dArray);
} ;

inline void AllayCircleRoute::setRadius(double r) { radius=r; }
inline void AllayCircleRoute::setAngle(double a)
{    angle=a;
         if(a>=0) flag=true;
         else     flag=false;
}
inline void AllayCircleRoute::setAllay(double L0) { allay=L0; }
inline void AllayCircleRoute::setJDstake(char *stake)
{   JDstake=new char[strlen(stake)+1];
    strcpy(JDstake,_strupr(strdup(stake)));
};
inline void AllayCircleRoute::setStep(double l) { step=l; }

inline double AllayCircleRoute::getRadius(void) { return(radius); }
inline double AllayCircleRoute::getAngle(void)  { return(angle);  }
inline double AllayCircleRoute::getAllay(void)  { return(allay);  }
inline double AllayCircleRoute::getStep(void)   { return(step);  }
inline bool   AllayCircleRoute::getFlag(void)   { return(flag);  }
///////////////////////////////////////////

inline double AllayCircleRoute::getT(void)
{        double m, p;
    m=allay/2.0-pow(allay,3)/(240.0*radius*radius);
        p=allay*allay/(24.0*radius);
    return(m+(radius+p)*tan(angle*0.5));
}

inline double AllayCircleRoute::getL(void)
{   double b, ang;
    b=allay/(2.0*radius); ang=angle-2*b;
        return(ang*radius+2*allay);
}

inline double AllayCircleRoute::getE(void)
{        double p=allay*allay/(24.0*radius);
    return((radius+p)/cos(angle*0.5)-radius);
}

inline double AllayCircleRoute::getq(void) { return(2*this->getT()-this->getL()); }
////////////////////////////////////////////

inline char* AllayCircleRoute::getJDstake(void) { return(JDstake); }

inline char* AllayCircleRoute::getZHstake(void)//计算直缓点的里程桩号
{        char *ZHstake,strkm[6],strm[10],sk[2][32],ss[2][32];
    double zh_km,zh_m,zh,jd;
        ZHstake=new char[strlen(JDstake)+1]; ZHstake[0]='\0';
    rzrxstrlist(JDstake,'K',sk); rzrxstrlist(sk[1],'+',ss);
    jd=atof(ss[0])*1000+atof(ss[1]);
    zh=jd-this->getT();
    zh_km=(int)(zh/1000);  zh_m=zh-zh_km*1000;
        sprintf(strkm,"%.0lf",zh_km); sprintf(strm, "%.3lf",zh_m );
        strcat(ZHstake,"K"); strcat(ZHstake,strkm);
        strcat(ZHstake,"+"); strcat(ZHstake,strm);
        return(ZHstake);
}

inline char* AllayCircleRoute::getHYstake(void)//计算缓圆点的里程桩号
{        char *HYstake,strkm[6],strm[10],sk[2][32],ss[2][32];
    double hy_km,hy_m,hy,zh;
        HYstake=new char[strlen(this->getZHstake())+1]; HYstake[0]='\0';
    rzrxstrlist(this->getZHstake(),'K',sk); rzrxstrlist(sk[1],'+',ss);
    zh=atof(ss[0])*1000+atof(ss[1]);
    hy=zh+allay;
    hy_km=(int)(hy/1000);  hy_m=hy-hy_km*1000;
        sprintf(strkm,"%.0lf",hy_km); sprintf(strm, "%.3lf",hy_m );
        strcat(HYstake,"K"); strcat(HYstake,strkm);
        strcat(HYstake,"+"); strcat(HYstake,strm);
        return(HYstake);
}

inline char* AllayCircleRoute::getQZstake(void)//计算曲中点的里程桩号
{        char *QZstake,strkm[5],strm[10],sk[2][32],ss[2][32];
    double qz_km,qz_m,qz,zh;
        QZstake=new char[strlen(this->getZHstake())+1]; QZstake[0]='\0';
    rzrxstrlist(this->getZHstake(),'K',sk); rzrxstrlist(sk[1],'+',ss);
    zh=atof(ss[0])*1000+atof(ss[1]);   
        qz=zh+this->getL()*0.5;
    qz_km=(int)(qz/1000);  qz_m=qz-qz_km*1000;
        sprintf(strkm,"%.0lf",qz_km); sprintf(strm, "%.3lf",qz_m );
        strcat(QZstake,"K"); strcat(QZstake,strkm);
        strcat(QZstake,"+"); strcat(QZstake,strm);
        return(QZstake);
}

inline char* AllayCircleRoute::getHZstake(void)//计算圆缓直点的里程桩号
{        char *HZstake,strkm[5],strm[10],sk[2][32],ss[2][32];
    double hz_km,hz_m,hz,zh;
        HZstake=new char[strlen(this->getZHstake())+1]; HZstake[0]='\0';
    rzrxstrlist(this->getZHstake(),'K',sk); rzrxstrlist(sk[1],'+',ss);
    zh=atof(ss[0])*1000+atof(ss[1]);
        hz=zh+this->getL();
    hz_km=(int)(hz/1000);  hz_m=hz-hz_km*1000;
        sprintf(strkm,"%.0lf",hz_km); sprintf(strm, "%.3lf",hz_m );
        strcat(HZstake,"K"); strcat(HZstake,strkm);
        strcat(HZstake,"+"); strcat(HZstake,strm);
        return(HZstake);
}

inline char* AllayCircleRoute::getYHstake(void)//计算圆缓点的里程桩号
{        char *YHstake,strkm[5],strm[10],sk[2][32],ss[2][32];
    double yh_km,yh_m,yh,hz;
        YHstake=new char[strlen(this->getHZstake())+1]; YHstake[0]='\0';
    rzrxstrlist(this->getHZstake(),'K',sk); rzrxstrlist(sk[1],'+',ss);
    hz=atof(ss[0])*1000+atof(ss[1]);
    yh=hz-allay;
    yh_km=(int)(yh/1000);  yh_m=yh-yh_km*1000;
        sprintf(strkm,"%.0lf",yh_km); sprintf(strm, "%.3lf",yh_m );
        strcat(YHstake,"K"); strcat(YHstake,strkm);
        strcat(YHstake,"+"); strcat(YHstake,strm);
        return(YHstake);
}


inline bool AllayCircleRoute::getCoordinate(AcGePoint3dArray& point3dArray)
{  int num,i;
   double xi,yi,Li,ai,m,p;
   AcGePoint3d point3d;
   ////缓和曲线坐标计算
   num=(int)(allay/step);
   for(i=0;i<num;i++) {
          Li=i*step;
          xi=Li-pow(Li,5)/(40.0*pow((radius*allay),2));
          if(flag==true) yi=pow(Li,3)/(6.0*radius*allay);
          else           yi=-pow(Li,3)/(6.0*radius*allay);
      point3d.set(xi,yi,0);
          point3dArray.append(point3d);
   }
   xi=allay-pow(allay,3)/(40.0*radius*radius);
   if(flag==true) yi=pow(allay,2)/(6.0*radius);
   else           yi=-pow(allay,2)/(6.0*radius);
   point3d.set(xi,yi,0); point3dArray.append(point3d);

   ////圆曲线坐标计算
   m=allay/2.0-pow(allay,3)/(240.0*radius*radius);
   p=pow(allay,2)/(24.0*radius);
   num=(int)((this->getL()-2.0*allay)/(2.0*step));
   for(i=1;i<num;i++) {
          Li=allay+i*step; ai=(Li-0.5*allay)/radius;
      xi=radius*sin(ai)+m;
          if(flag==true) yi=radius*(1-cos(ai))+p;
          else           yi=-radius*(1-cos(ai))+p;
      point3d.set(xi,yi,0);
          point3dArray.append(point3d);
   }
   ai=(this->getL()-allay)/(2.0*radius);
   xi=radius*sin(ai)+m;
   if(flag==true) yi=radius*(1-cos(ai))+p;
   else           yi=-radius*(1-cos(ai))+p;
   point3d.set(xi,yi,0); point3dArray.append(point3d);
   if(point3dArray.length() < 3) return false;
   else                          return true;
}

////////////////////////////////////////////////////
#endif // _INC_ENGIN_HPP/***
*Engin.hpp - definitions and declarations for additional engineering library
*
*       Copyright (c) 2000-2001, By Liurenzhao. All rights reserved.
*
*Purpose:
*       This file contains constant definitions and external subroutine
*       declarations for the additional math subroutine library.
*
*       [Public]
*
****/

#ifndef _INC_ENGIN_HPP
#define _INC_ENGIN_HPP

///////////////////////////////////
class CircleRoute      //圆曲线测设
///////////////////////////////////
{  
private:      
     double radius;  //半径(米)
     double angle;   //偏角(弧度)
         char *JDstake;  //交点桩号 例如: JDstake[]="K87+441.41"
         double step;
         bool  flag;     //左偏为true,右偏为false

public:
     CircleRoute(double r,double a, char *stake)
         {   radius=r; angle=fabs(a);
             if(a>=0) flag=true;
                 else     flag=false;
             JDstake=new char[strlen(stake)+1];
             strcpy(JDstake,_strupr(strdup(stake)));
         };
     CircleRoute(double r,double a, double l=10.0)
         {   radius=r; angle=fabs(a); step=l;
             if(a>=0) flag=true;
                 else     flag=false;             
         };
     CircleRoute(){};
     ~CircleRoute() { delete JDstake; };

         void setRadius(double);
         void setAngle(double);
         void setJDstake(char *);
         void setStep(double);

         double getRadius(void);
         double getAngle(void);
         double getStep(void);
         bool   getFlag(void);

         double getT(void); //切线长
         double getL(void); //曲线长
         double getE(void); //外矢距
         double getq(void); //切曲差

         char* getJDstake(void);
     char* getZYstake(void);
     char* getYZstake(void);
     char* getQZstake(void);

     bool getCoordinate(AcGePoint3dArray& point3dArray);
} ;

inline void CircleRoute::setRadius(double r) { radius=r; }
inline void CircleRoute::setAngle(double a)
{         angle=fabs(a);
         if(a>=0) flag=true;
         else     flag=false;
}
inline void CircleRoute::setJDstake(char *stake)
{   JDstake=new char[strlen(stake)+1];
    strcpy(JDstake,_strupr(strdup(stake)));
};
inline void CircleRoute::setStep(double l)   { step=l; }

inline double CircleRoute::getRadius(void) { return(radius); }
inline double CircleRoute::getAngle(void)  { return(angle);  }
inline double CircleRoute::getStep(void)   { return(step);   }
inline bool   CircleRoute::getFlag(void)   {  return(flag);  }
///////////////////////////////////////////

inline double CircleRoute::getT(void) {        return(radius*tan(angle*0.5));  }
inline double CircleRoute::getL(void) {        return(angle*radius);  }
inline double CircleRoute::getE(void) {        return(radius/cos(angle*0.5)-radius); }
inline double CircleRoute::getq(void) {        return(2*this->getT()-this->getL()); }
////////////////////////////////////////////

inline char* CircleRoute::getJDstake(void) { return(JDstake); }

inline char* CircleRoute::getZYstake(void)
{        char *ZYstake,strkm[6],strm[10],sk[2][32],ss[2][32];
    double zy_km,zy_m,zy,jd;
        ZYstake=new char[strlen(JDstake)+1]; ZYstake[0]='\0';
    rzrxstrlist(JDstake,'K',sk); rzrxstrlist(sk[1],'+',ss);
    jd=atof(ss[0])*1000+atof(ss[1]);
    zy=jd-this->getT();
    zy_km=(int)(zy/1000);  zy_m=zy-zy_km*1000;
        sprintf(strkm,"%.0lf",zy_km); sprintf(strm, "%.3lf",zy_m );
        strcat(ZYstake,"K"); strcat(ZYstake,strkm);
        strcat(ZYstake,"+"); strcat(ZYstake,strm);
        return(ZYstake);
}

inline char* CircleRoute::getYZstake(void)
{        char *YZstake,strkm[5],strm[10],sk[2][32],ss[2][32];
    double yz_km,yz_m,yz,jd;
        YZstake=new char[strlen(JDstake)+1]; YZstake[0]='\0';
    rzrxstrlist(JDstake,'K',sk); rzrxstrlist(sk[1],'+',ss);
    jd=atof(ss[0])*1000+atof(ss[1]);
    yz=jd-this->getT()+this->getL();
    yz_km=(int)(yz/1000);  yz_m=yz-yz_km*1000;
        sprintf(strkm,"%.0lf",yz_km); sprintf(strm, "%.3lf",yz_m );
        strcat(YZstake,"K"); strcat(YZstake,strkm);
        strcat(YZstake,"+"); strcat(YZstake,strm);
        return(YZstake);
}

inline char* CircleRoute::getQZstake(void)
{        char *QZstake,strkm[5],strm[10],sk[2][32],ss[2][32];
    double qz_km,qz_m,qz,jd;
        QZstake=new char[strlen(JDstake)+1]; QZstake[0]='\0';
    rzrxstrlist(JDstake,'K',sk); rzrxstrlist(sk[1],'+',ss);
    jd=atof(ss[0])*1000+atof(ss[1]);   
        qz=jd-this->getT()+this->getL()*0.5;
    qz_km=(int)(qz/1000);  qz_m=qz-qz_km*1000;
        sprintf(strkm,"%.0lf",qz_km); sprintf(strm, "%.3lf",qz_m );
        strcat(QZstake,"K"); strcat(QZstake,strkm);
        strcat(QZstake,"+"); strcat(QZstake,strm);
        return(QZstake);
}


inline bool CircleRoute::getCoordinate(AcGePoint3dArray& point3dArray)
{  int num,i;
   double xi,yi,arc;
   AcGePoint3d point3d;

   arc=step/radius; num=(int)(angle/arc);
   for(i=0;i<num;i++) {
             xi=radius*sin(i*arc);
      if(flag==true) yi=radius*(1-cos(i*arc));
          else           yi=-radius*(1-cos(i*arc));
      point3d.set(xi,yi,0);
          point3dArray.append(point3d);
   }
   xi=radius*sin(angle);
   if(flag==true) yi=radius*(1-cos(angle));
   else           yi=-radius*(1-cos(angle));
   point3d.set(xi,yi,0); point3dArray.append(point3d);
   if(point3dArray.length() < 3) return false;
   else                          return true;
}



///////////////////////////////////////
class AllayCircleRoute   //缓和曲线测设
///////////////////////////////////////
{  
private:      
     double radius;  //半径(米)
     double angle;   //偏角(弧度)
         char *JDstake;  //交点桩号 例如: JDstake[]="K87+441.41"
         double step;
         double allay;   //缓和曲线长度
         bool  flag;     //左偏为true,右偏为false

public:
     AllayCircleRoute(double r, double a, double L0, char *stake)
         {   radius=r; angle=fabs(a); allay=L0;
             if(a>=0) flag=true;
                 else     flag=false;
             JDstake=new char[strlen(stake)+1];
             strcpy(JDstake,_strupr(strdup(stake)));
         };
     AllayCircleRoute(double r, double a, double L0, double l=10.0)
         {   radius=r; angle=fabs(a); allay=L0; step=l;
             if(a>=0) flag=true;
                 else     flag=false;
         };
     AllayCircleRoute(){};
     ~AllayCircleRoute() { delete JDstake; };

         void setRadius(double);
         void setAngle(double);
         void setAllay(double);
         void setJDstake(char *);
         void setStep(double);

         double getRadius(void);
         double getAngle(void);
         double getAllay(void);
         double getStep(void);
         bool   getFlag(void);

         double getT(void); //切线长
         double getL(void); //曲线长
         double getE(void); //外矢距
         double getq(void); //切曲差

         char* getJDstake(void);
     char* getZHstake(void);
     char* getHYstake(void);
     char* getQZstake(void);
     char* getZYstake(void);
     char* getYHstake(void);
     char* getHZstake(void);

     bool getCoordinate(AcGePoint3dArray& point3dArray);
} ;

inline void AllayCircleRoute::setRadius(double r) { radius=r; }
inline void AllayCircleRoute::setAngle(double a)
{    angle=a;
         if(a>=0) flag=true;
         else     flag=false;
}
inline void AllayCircleRoute::setAllay(double L0) { allay=L0; }
inline void AllayCircleRoute::setJDstake(char *stake)
{   JDstake=new char[strlen(stake)+1];
    strcpy(JDstake,_strupr(strdup(stake)));
};
inline void AllayCircleRoute::setStep(double l) { step=l; }

inline double AllayCircleRoute::getRadius(void) { return(radius); }
inline double AllayCircleRoute::getAngle(void)  { return(angle);  }
inline double AllayCircleRoute::getAllay(void)  { return(allay);  }
inline double AllayCircleRoute::getStep(void)   { return(step);  }
inline bool   AllayCircleRoute::getFlag(void)   { return(flag);  }
///////////////////////////////////////////

inline double AllayCircleRoute::getT(void)
{        double m, p;
    m=allay/2.0-pow(allay,3)/(240.0*radius*radius);
        p=allay*allay/(24.0*radius);
    return(m+(radius+p)*tan(angle*0.5));
}

inline double AllayCircleRoute::getL(void)
{   double b, ang;
    b=allay/(2.0*radius); ang=angle-2*b;
        return(ang*radius+2*allay);
}

inline double AllayCircleRoute::getE(void)
{        double p=allay*allay/(24.0*radius);
    return((radius+p)/cos(angle*0.5)-radius);
}

inline double AllayCircleRoute::getq(void) { return(2*this->getT()-this->getL()); }
////////////////////////////////////////////

inline char* AllayCircleRoute::getJDstake(void) { return(JDstake); }

inline char* AllayCircleRoute::getZHstake(void)//计算直缓点的里程桩号
{        char *ZHstake,strkm[6],strm[10],sk[2][32],ss[2][32];
    double zh_km,zh_m,zh,jd;
        ZHstake=new char[strlen(JDstake)+1]; ZHstake[0]='\0';
    rzrxstrlist(JDstake,'K',sk); rzrxstrlist(sk[1],'+',ss);
    jd=atof(ss[0])*1000+atof(ss[1]);
    zh=jd-this->getT();
    zh_km=(int)(zh/1000);  zh_m=zh-zh_km*1000;
        sprintf(strkm,"%.0lf",zh_km); sprintf(strm, "%.3lf",zh_m );
        strcat(ZHstake,"K"); strcat(ZHstake,strkm);
        strcat(ZHstake,"+"); strcat(ZHstake,strm);
        return(ZHstake);
}

inline char* AllayCircleRoute::getHYstake(void)//计算缓圆点的里程桩号
{        char *HYstake,strkm[6],strm[10],sk[2][32],ss[2][32];
    double hy_km,hy_m,hy,zh;
        HYstake=new char[strlen(this->getZHstake())+1]; HYstake[0]='\0';
    rzrxstrlist(this->getZHstake(),'K',sk); rzrxstrlist(sk[1],'+',ss);
    zh=atof(ss[0])*1000+atof(ss[1]);
    hy=zh+allay;
    hy_km=(int)(hy/1000);  hy_m=hy-hy_km*1000;
        sprintf(strkm,"%.0lf",hy_km); sprintf(strm, "%.3lf",hy_m );
        strcat(HYstake,"K"); strcat(HYstake,strkm);
        strcat(HYstake,"+"); strcat(HYstake,strm);
        return(HYstake);
}

inline char* AllayCircleRoute::getQZstake(void)//计算曲中点的里程桩号
{        char *QZstake,strkm[5],strm[10],sk[2][32],ss[2][32];
    double qz_km,qz_m,qz,zh;
        QZstake=new char[strlen(this->getZHstake())+1]; QZstake[0]='\0';
    rzrxstrlist(this->getZHstake(),'K',sk); rzrxstrlist(sk[1],'+',ss);
    zh=atof(ss[0])*1000+atof(ss[1]);   
        qz=zh+this->getL()*0.5;
    qz_km=(int)(qz/1000);  qz_m=qz-qz_km*1000;
        sprintf(strkm,"%.0lf",qz_km); sprintf(strm, "%.3lf",qz_m );
        strcat(QZstake,"K"); strcat(QZstake,strkm);
        strcat(QZstake,"+"); strcat(QZstake,strm);
        return(QZstake);
}

inline char* AllayCircleRoute::getHZstake(void)//计算圆缓直点的里程桩号
{        char *HZstake,strkm[5],strm[10],sk[2][32],ss[2][32];
    double hz_km,hz_m,hz,zh;
        HZstake=new char[strlen(this->getZHstake())+1]; HZstake[0]='\0';
    rzrxstrlist(this->getZHstake(),'K',sk); rzrxstrlist(sk[1],'+',ss);
    zh=atof(ss[0])*1000+atof(ss[1]);
        hz=zh+this->getL();
    hz_km=(int)(hz/1000);  hz_m=hz-hz_km*1000;
        sprintf(strkm,"%.0lf",hz_km); sprintf(strm, "%.3lf",hz_m );
        strcat(HZstake,"K"); strcat(HZstake,strkm);
        strcat(HZstake,"+"); strcat(HZstake,strm);
        return(HZstake);
}

inline char* AllayCircleRoute::getYHstake(void)//计算圆缓点的里程桩号
{        char *YHstake,strkm[5],strm[10],sk[2][32],ss[2][32];
    double yh_km,yh_m,yh,hz;
        YHstake=new char[strlen(this->getHZstake())+1]; YHstake[0]='\0';
    rzrxstrlist(this->getHZstake(),'K',sk); rzrxstrlist(sk[1],'+',ss);
    hz=atof(ss[0])*1000+atof(ss[1]);
    yh=hz-allay;
    yh_km=(int)(yh/1000);  yh_m=yh-yh_km*1000;
        sprintf(strkm,"%.0lf",yh_km); sprintf(strm, "%.3lf",yh_m );
        strcat(YHstake,"K"); strcat(YHstake,strkm);
        strcat(YHstake,"+"); strcat(YHstake,strm);
        return(YHstake);
}


inline bool AllayCircleRoute::getCoordinate(AcGePoint3dArray& point3dArray)
{  int num,i;
   double xi,yi,Li,ai,m,p;
   AcGePoint3d point3d;
   ////缓和曲线坐标计算
   num=(int)(allay/step);
   for(i=0;i<num;i++) {
          Li=i*step;
          xi=Li-pow(Li,5)/(40.0*pow((radius*allay),2));
          if(flag==true) yi=pow(Li,3)/(6.0*radius*allay);
          else           yi=-pow(Li,3)/(6.0*radius*allay);
      point3d.set(xi,yi,0);
          point3dArray.append(point3d);
   }
   xi=allay-pow(allay,3)/(40.0*radius*radius);
   if(flag==true) yi=pow(allay,2)/(6.0*radius);
   else           yi=-pow(allay,2)/(6.0*radius);
   point3d.set(xi,yi,0); point3dArray.append(point3d);

   ////圆曲线坐标计算
   m=allay/2.0-pow(allay,3)/(240.0*radius*radius);
   p=pow(allay,2)/(24.0*radius);
   num=(int)((this->getL()-2.0*allay)/(2.0*step));
   for(i=1;i<num;i++) {
          Li=allay+i*step; ai=(Li-0.5*allay)/radius;
      xi=radius*sin(ai)+m;
          if(flag==true) yi=radius*(1-cos(ai))+p;
          else           yi=-radius*(1-cos(ai))+p;
      point3d.set(xi,yi,0);
          point3dArray.append(point3d);
   }
   ai=(this->getL()-allay)/(2.0*radius);
   xi=radius*sin(ai)+m;
   if(flag==true) yi=radius*(1-cos(ai))+p;
   else           yi=-radius*(1-cos(ai))+p;
   point3d.set(xi,yi,0); point3dArray.append(point3d);
   if(point3dArray.length() < 3) return false;
   else                          return true;
}

////////////////////////////////////////////////////
#endif // _INC_ENGIN_HPP

欧特

好样的！
虽然，在野外这段程序，对于使用E500、Fx-4500p、PDA的用处不一定有多大的用处。有空还是研究一下。
鞠躬，谢谢！

xzg_hello

这是我搞施工时编的放样程序，里面有关于边桩放样的算法（即切线方位角和法线方位角的计算），可以作为对楼上的补充。程序是用QB编写的，如果要运行还要QBASIC。程序中关于缓和曲线的算法仍然用的是常规算法，不是另一个帖子中4500计算器中的算法，因为当时还不知道那算法。
希望晓东能将楼上兄弟的函数和我程序中关于切线方位角和法线方位角的部分加到API中，这样以后要编制道路方面的程序就轻松了。

xzg_hello

源程序：

ymd

这里有个4000p计算器测设曲线程序，可置镜任意点，测设曲线（S－R—S型）上任意点，坐标
系统可为假定坐标系统或绝对坐标系统，其他计算器可参照修改！

xzg_hello

文章好像还没有写完，能否补充完？
谢谢！

wweien

我可能用的着，谢谢

xb4270293

不完整的缓和曲线无法使用，复曲线中起点半径不等于∞，不能使用教材上的计算公式。

xzg_hello

最初由 xb4270293 发布
不完整的缓和曲线无法使用，复曲线中起点半径不等于∞，不能使用教材上的计算公式。

确实如此，你知道什么书上有这方面的介绍么？
谢谢！！！！！

ymd

有一本书大概是《铁路缓和曲线》，有各种形式的缓和曲线计算公式
，包括连接不同半径的缓和曲线计算公式，但这种不完整缓和曲线计
算公式复杂，若取的项数多，公式一常串。其实用完整缓和曲线计算
公式加简单的平移旋转即可得到坐标。

xzg_hello

谢谢！！！我去找找看……

下面是在超星找到的几本书：
1、《曲线广义公式及其应用》
   测站设在中线上，综合公式在曲线测设中具有以下功能：
    (1)测站设置在直线上任一点，测设缓和曲线和圆曲线上任
—点。
    (2)测站设置在缓和曲线上任一点，测设缓和曲线、圆曲线和
直线上任一点。
    (3)测站设置在圆曲线上任点，测设圆曲线、缓和曲线和直
线上任一点。
    (4)计算缓和曲线和圆曲线上任一段弧的偏角
    (5)汁其缓和曲线和圆曲线上任一段弧的弦长
    (6)计算曲线综合要素。

2、《铁路曲线》
里面真可谓洋洋大观，铁老大这名字可真不是白叫的，以前看的曲线测设在这仅能算冰山一角

3、《铁路缓和曲线理论和类型（第二版）》
如果上面那本称之为广博，那么这本称为精专好了，当然第一本就只能称为实用了：）

xiaowa

太复杂。。。不懂

forme

不是有个测量的《公路坐标计算系统》吗？

fjbl

只有把它编成便携机的程序这样才具有强大的生命力。

nyfgp

你说，Gps是干什么的？！公路放线没问题能用！