[ARX程序]：ARX程序中怎样求多义线的长度

bluesolar

ARX程序中怎样求多义线的长度?
十万火急,恳请帮忙.
我用先求各顶点之间的距离,再根据凸度求出曲线段的长度的方法
累加得到其长度,但误差很大.
根据凸度求出曲线段的长度的方法为:
根据两点距离和凸度可求出弓高,再用三点画弧求出半径,
再求出圆心角,最后算出弧长.
各位有何简便的方法,万望不吝赐教!

XDSoft

最初由 bluesolar 发布
[B]ARX程序中怎样求多义线的长度?
十万火急,恳请帮忙.
我用先求各顶点之间的距离,再根据凸度求出曲线段的长度的方法
累加得到其长度,但误差很大.
根据凸度求出曲线段的长度的方法为:
根据两点距离和凸度可求出弓高... [/B]

因为ARX提供了强大的AcGe几何库，所以几何图形学算法都不用我们操心了。

如果是AcDbPolyline类对象，那么可以方便的直接用AcDbCurve方法，首先用AcDbCurve::getEndParam(...)得到终点的参数值，然后用AcDbCurve::getDistAtParam(...)得到终点参数值处到曲线起点的长度，就是你要的了。

如果是AcDb2dPolyline对象，那么简单的方法是用：AcDbPolyline::convertFrom(...)转换成AcDbPolyline对象，然后用上面的方法求。

复杂的方法是遍历AcDb2dPolyline的AcDbVertex子实体对象，然后用：然后根据段类型，构筑AcGe对象，然后用AcDbLineSeg2d,AcGeCircArc2d类的方法，计算长度。

XDSoft

给你一个求“老”的多义线AcDb2dPolyline对象长度的代码

1. 
   
2. [FONT=courier new]
   
3. Acad::ErrorStatus acdcGetAcDb2dPolylineLength(double& length,AcDb2dPolyline*
   
4. p2dPoly){
   
5.    if (p2dPoly==NULL)
   
6.     return Acad::eNullObjectPointer;
   
7.    length =0;
   
8.    Acad::ErrorStatus es;
   
9.    AcDbObjectIterator* pI;
   
10.    if ((pI=p2dPoly->vertexIterator())==NULL)
    
11.     return Acad::eOutOfMemory;
    
12.    AcDb2dVertex* pV1,*pV2;
    
13.    AcGePoint3d v1,v2;
    
14.    Adesk::Boolean bDone = Adesk::kFalse;
    
15.    AcDb::Vertex2dType type;
    
16.    if ((es=p2dPoly->openVertex(pV1,pI->objectId(),AcDb::kForRead))==Acad::eOk){
    
17.     while (Adesk::kTrue){
    
18.         if (!bDone)
    
19.         pI->step();
    
20.         if
    
21. ((es=p2dPoly->openVertex(pV2,pI->objectId(),AcDb::kForRead))==Acad::eOk){
    
22.         switch (type = pV1->vertexType()){
    
23.         case AcDb::k2dVertex:
    
24.         case AcDb::k2dSplineFitVertex:
    
25.         case AcDb::k2dCurveFitVertex:
    
26.             {
    
27.             while (pV2->vertexType()==AcDb::k2dSplineCtlVertex){
    
28.                 pV2->close();
    
29.                 pI->step();
    
30.                 if
    
31. ((es=p2dPoly->openVertex(pV2,pI->objectId(),AcDb::kForRead))!=Acad::eOk)
    
32.                 break;
    
33.             }
    
34.             if (es!=Acad::eOk)
    
35.                 break;
    
36.             v1 = pV1->position();
    
37.             v2 = pV2->position();
    
38.             double A = 4*atan(pV1->bulge()); //angle
    
39.             double D = (v2-v1).length();
    
40.             if (fabs(A)<AcGeContext::gTol.equalPoint()){
    
41.                 length+=D;
    
42.             } else {
    
43.                 double R= fabs(D/(2*sin(A/2)));//radius
    
44.                 length+= R*fabs(A);//length
    
45.             }
    
46.             pV1->close();
    
47.             pV1=pV2;
    
48.             break;
    
49.             }
    
50.         case AcDb::k2dSplineCtlVertex:
    
51.             pV1->close();
    
52.             pV1=pV2;
    
53.             break;
    
54.         default:
    
55.             assert(0);//should never happen
    
56.         }
    
57.         if (bDone)
    
58.             break;
    
59.         } else {
    
60.         if (p2dPoly->isClosed())
    
61.             pI->start();
    
62.         else
    
63.             break;
    
64.         bDone = Adesk::kTrue;
    
65.         }
    
66.     }
    
67.     pV1->close();
    
68.    }
    
69.    delete pI;
    
70.    return Acad::eOk;
    
71. }
    
72. [/FONT]

复制代码

bluesolar

非常感谢晓东，已经作出来了。
再问一下，将弧度化为CSting的角度，该怎么办呢？
谢谢！

XDSoft

最初由 bluesolar 发布
[B]非常感谢晓东，已经作出来了。
再问一下，将弧度化为CSting的角度，该怎么办呢？
谢谢！ [/B]

int
acdbAngToS(
    ads_real v,
    int unit,
    int prec,

    char* str);

在写什么程序，能把代码和大家共享吗？谢谢

bluesolar

但不知int unit,
int prec，各代表什么意思？不好意思，手头的资料全还到北图了，很远的，现在一本资料也没有。
在作一个铁路线路平面设计的小软件，功能很简单的。程序在我的机器上（宿舍内，上传文件很慢），我现在在机房上机，没带。后天毕业答辩，之后我再贴上来吧。
功能太简单了，都有点不好意思贴上来，怕大家笑话。

XDSoft

最初由 bluesolar 发布
[B]但不知int unit,
int prec，各代表什么意思？不好意思，手头的资料全还到北图了，很远的，现在一本资料也没有。
在作一个铁路线路平面设计的小软件，功能很简单的。程序在我的机器上（宿舍内，上传文件很慢），我... [/B]

unit--单位，见系统变量AUNITS
prec--精度，见系统变量AUPREC

bluesolar

所有的文件都在。

shw1311

:):funk::L:o:P