XDGE几何库应用(23)--精度对Ge模型影响

st788796 · 发表于 2014-10-2 11:06:47

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

您需要登录才可以下载或查看，没有账号？立即注册

×

本帖最后由 st788796 于 2014-11-1 07:55 编辑

一个求两个曲线首尾连接后闭合区域的测试，Ge 的 AcGeCompositeCurve3d 和 AcDbRegion 对比，在 xdrxapi 中默认长度精度是 0.01 ，在 ARX 内部模式 Tol 是 1e-10 (Provides a "global" default tolerance object. This is an AcGeTol object that contains the default AcGeTol tolerance values (equalPoint and equalVector tolerances set to 1.e-10). ，这个应用中可以看到面积差在 5% 左右，长度差不大，不过即使用直线段测试也存在误差，搞不懂 Ge 是怎么算的

原帖地址 http://bbs.xdcad.net/forum.php?m ... 4025&extra=page%3D1

(defun c:tt (/ e1 e2 ge1 ge2 sp1 ep1 sp2 ep2 tmp pl ge3d1 ge3d2 gcc gel1
gel2 in endparam lst region old
) ; (setq old (xdrx_document_getprec))
; (xdrx_document_setprec 1e-8 (cadr
; old));_单纯提高长度精度会造成 area
; 为nil
(if (and
(setq e1 (car (xdrx_entsel "\nPicks Pline: " '((0 . "*Polyline,Line,Arc,Ellipse,Spline")))))
(setq e2 (car (xdrx_entsel "\nPicks Pline: " '((0 . "*Polyline,Line,Arc,Ellipse,Spline"))))) ; _两条单一线两端直连形成的闭合区域投影面积
(setq ge1 (xdge::constructor e1))
(not (xdge::getpropertyvalue ge1 "isClosed"))
(setq ge2 (xdge::constructor e2))
(not (xdge::getpropertyvalue ge2 "isClosed"))
(setq gcc (xdge::constructor "kCurveCurveInt3d" ge1 ge2))
(zerop (xdge::getpropertyvalue gcc "numIntPoints"))
)
(progn
(mapcar
'set
'(sp1 ep1)
(xdge::getpropertyvalue ge1 "hasStartPoint" "hasEndPoint")
)
(mapcar
'set
'(sp2 ep2)
(xdge::getpropertyvalue ge2 "hasStartPoint" "hasEndPoint")
)
(if (inters
ep1
sp2
ep1
sp1
) ; _检查方向,交叉时对其中之一换方向
(progn
(xdge::getpropertyvalue ge2 "reverseParam")
(setq tmp sp2
sp2 ep2
ep2 tmp
)
)
)
(grvecs (list -1 ep1 sp2 -1 ep2 sp1))
(setq ln1 (xdge::constructor "kLineSeg3d" ep1 sp2)
ln2 (xdge::constructor "kLineSeg3d" ep2 sp1)
) ; _两端封口
; 构建符合线的 Ge
; 须是简单线，不能复合线与简单线混合
(if (eq (xdge::type ge1) "kCompositeCrv3d")
(setq gel1 (xdge::getpropertyvalue ge1 "getCurveList"))
(setq gel1 (list ge1))
)
(if (eq (xdge::type ge2) "kCompositeCrv3d")
(setq gel2 (xdge::getpropertyvalue ge2 "getCurveList"))
(setq gel2 (list ge2))
)
(setq pl (xdge::constructor "kCompositeCrv3d" (append
gel1
(list ln1)
gel2
(list ln2)
) ; _简单线表
)
in (car (xdge::getpropertyvalue pl "getInterval")) ; _要用
; Interval 查询 endparam
endparam (xdge::getpropertyvalue in "upperBound") ; _上限
lst (xdge::getpropertyvalue pl "area" 0.0 endparam "length" 0.0
endparam
)
) ; _构建连续复合线
(princ "\nCompositeCurve Area = ")
(princ (rtos (car lst) 2 4)) ; _此处和 Region测得的面积有误差
(princ "\nCompositeCurve Length = ")
(princ (rtos (cadr lst) 2 4))
(setq region (xdrx_getpropertyvalue (ssname
(xdrx_curve->region
(xdge::entity:make ge1 ge2 ln1 ln2)
) 0
) "area" "Perimeter"
)
)
(princ "\n---------------------------------")
(princ "\nRegion Area = ")
(princ (rtos (car region) 2 4)) ; _此处和 Region测得的面积有误差
(princ "\nRegion Length = ")
(princ (rtos (cadr region) 2 4))
(princ "\n---------------------------------")
(princ "\n")
(princ (mapcar
'-
region
lst
)
)
(princ "\n---------------------------------")
(princ "\n")
(princ (mapcar
'/
(mapcar
'abs
(mapcar
'-
region
lst
)
)
region
)
)
(princ)
)
) ; (apply 'xdrx_document_setprec old)
(xdge::free)
(princ)
)

(defun c:tt (/ e1 e2 ge1 ge2 sp1 ep1 sp2 ep2 tmp pl ge3d1 ge3d2 gcc gel1 
               gel2 in endparam lst region old 
            )                               ; (setq old (xdrx_document_getprec)) 
                                       ; (xdrx_document_setprec 1e-8 (cadr 
                                       ; old));_单纯提高长度精度会造成 area 
                                       ; 为nil 
  (if (and 
        (setq e1 (car (xdrx_entsel "\nPicks Pline: " '((0 . "*Polyline,Line,Arc,Ellipse,Spline"))))) 
        (setq e2 (car (xdrx_entsel "\nPicks Pline: " '((0 . "*Polyline,Line,Arc,Ellipse,Spline"))))) ; _两条单一线两端直连形成的闭合区域投影面积 
        (setq ge1 (xdge::constructor e1)) 
        (not (xdge::getpropertyvalue ge1 "isClosed")) 
        (setq ge2 (xdge::constructor e2)) 
        (not (xdge::getpropertyvalue ge2 "isClosed")) 
        (setq gcc (xdge::constructor "kCurveCurveInt3d" ge1 ge2)) 
        (zerop (xdge::getpropertyvalue gcc "numIntPoints")) 
      ) 
    (progn 
      (mapcar 
        'set 
        '(sp1 ep1) 
        (xdge::getpropertyvalue ge1 "hasStartPoint" "hasEndPoint") 
      ) 
      (mapcar 
        'set 
        '(sp2 ep2) 
        (xdge::getpropertyvalue ge2 "hasStartPoint" "hasEndPoint") 
      ) 
      (if (inters 
            ep1 
            sp2 
            ep1 
            sp1 
          )                               ; _检查方向,交叉时对其中之一换方向 
        (progn 
          (xdge::getpropertyvalue ge2 "reverseParam") 
          (setq tmp sp2 
                sp2 ep2 
                ep2 tmp 
          ) 
        ) 
      ) 
      (grvecs (list -1 ep1 sp2 -1 ep2 sp1)) 
      (setq ln1 (xdge::constructor "kLineSeg3d" ep1 sp2) 
            ln2 (xdge::constructor "kLineSeg3d" ep2 sp1) 
      )                                       ; _两端封口 
                                       ; 构建符合线的 Ge 
                                       ; 须是简单线，不能复合线与简单线混合 
      (if (eq (xdge::type ge1) "kCompositeCrv3d") 
        (setq gel1 (xdge::getpropertyvalue ge1 "getCurveList")) 
        (setq gel1 (list ge1)) 
      ) 
      (if (eq (xdge::type ge2) "kCompositeCrv3d") 
        (setq gel2 (xdge::getpropertyvalue ge2 "getCurveList")) 
        (setq gel2 (list ge2)) 
      ) 
      (setq pl (xdge::constructor "kCompositeCrv3d" (append 
                                                      gel1 
                                                      (list ln1) 
                                                      gel2 
                                                      (list ln2) 
                                                    ) ; _简单线表 
               ) 
            in (car (xdge::getpropertyvalue pl "getInterval")) ; _要用 
                                       ; Interval 查询 endparam 
            endparam (xdge::getpropertyvalue in "upperBound") ; _上限 
            lst (xdge::getpropertyvalue pl "area" 0.0 endparam "length" 0.0 
                                        endparam 
                ) 
      )                                       ; _构建连续复合线 
      (princ "\nCompositeCurve Area = ") 
      (princ (rtos (car lst) 2 4))     ; _此处和 Region测得的面积有误差 
      (princ "\nCompositeCurve Length = ") 
      (princ (rtos (cadr lst) 2 4)) 
      (setq region (xdrx_getpropertyvalue (ssname 
                                                  (xdrx_curve->region 
                                                                      (xdge::entity:make ge1 ge2 ln1 ln2) 
                                                  ) 0 
                                          ) "area" "Perimeter" 
                   ) 
      ) 
      (princ "\n---------------------------------") 
      (princ "\nRegion Area = ") 
      (princ (rtos (car region) 2 4))  ; _此处和 Region测得的面积有误差 
      (princ "\nRegion Length = ") 
      (princ (rtos (cadr region) 2 4)) 
      (princ "\n---------------------------------") 
      (princ "\n") 
      (princ (mapcar 
               '- 
               region 
               lst 
             ) 
      ) 
      (princ "\n---------------------------------") 
      (princ "\n") 
      (princ (mapcar 
               '/ 
               (mapcar 
                 'abs 
                 (mapcar 
                   '- 
                   region 
                   lst 
                 ) 
               ) 
               region 
             ) 
      ) 
      (princ) 
    ) 
  )                                       ; (apply 'xdrx_document_setprec old) 
  (xdge::free) 
  (princ) 
)

forum.php?mod=attachment&aid=ODU4NXwxY2ZkYTJkM3wxNDEyMjA3Nzg4fDE4NjA4fDY3NDAyNQ%3D%3D.gif

st788796 · 发表于 2014-10-2 11:09:23

命令: tt
Picks Pline:
Picks Pline:
CompositeCurve Area = 3943111.9298
CompositeCurve Length = 8778.0769
---------------------------------
Region Area = 3963976.0667

Region Length = 8778.0667
---------------------------------
(20864.1 -0.010213)
---------------------------------
(0.00526344 1.16346e-006)

st788796 · 发表于 2014-10-2 11:14:37

AcGeCompositeCurve 的结束参数可以这样求得

(setq glst (xdge::getpropertyvalue gcom "getCurveList")
endparam (xdge::getpropertyvalue
gcom
"localToGlobalParam"
1.0
(1- (length glst))
)
)

csharp · 发表于 2014-10-2 12:16:59

本帖最后由 csharp 于 2014-10-2 12:27 编辑

在 Net 中用系统默认精度这样测试了下，误差没有二楼那么大

然后又在 Autocad 中试着将 API 内部 prec 两个值设置一致后 Ge 模型和 Region 误差就很小了，和 net 测试结果基本一致

(xdrx_document_setprec 1e-10 1e-10)

命令: tt
Picks Pline:
Picks Pline:
CompositeCurve Area = 1916591.5121
CompositeCurve Length = 5770.6397
命令:
命令: aa AREA
指定第一个角点或 [对象(O)/加(A)/减(S)]: e
选择对象:
面积 = 1916591.5440，周长 = 6404.5092

构造 CompositeCurve 时省略了最后一段， area 自动将首尾连接计算

Net 简单测试代码，三条首尾相接 Curve

[CommandMethod("test1")]
public void Test01()
{
Document doc = Application.DocumentManager.MdiActiveDocument;
Editor ed = doc.Editor;
Database db = doc.Database;
Transaction tr = db.TransactionManager.StartTransaction();
PromptSelectionResult psr = ed.GetSelection();
if (psr.Status != PromptStatus.OK) return;
using (tr)
{
BlockTableRecord btr = (BlockTableRecord) tr.GetObject(db.CurrentSpaceId, OpenMode.ForWrite);
ObjectId[] ids = psr.Value.GetObjectIds();
Curve3d[] cc3ds = ids.Select(c => ((Curve)c.GetObject(OpenMode.ForRead)).GetGeCurve()).ToArray();
CompositeCurve3d cc = new CompositeCurve3d(cc3ds);
double area = cc.GetArea(cc.GetParameterOf(cc.StartPoint),
cc.GetParameterOf(cc.EndPoint),
Tolerance.Global);
ed.WriteMessage(area.ToString());
tr.Commit();
}
}

[CommandMethod("test1")] 
        public void Test01() 
        { 
            Document doc = Application.DocumentManager.MdiActiveDocument; 
            Editor ed = doc.Editor; 
            Database db = doc.Database; 
            Transaction tr = db.TransactionManager.StartTransaction(); 
            PromptSelectionResult psr = ed.GetSelection(); 
            if (psr.Status != PromptStatus.OK) return; 
            using (tr) 
            { 
                BlockTableRecord btr = (BlockTableRecord) tr.GetObject(db.CurrentSpaceId, OpenMode.ForWrite); 
                ObjectId[] ids = psr.Value.GetObjectIds(); 
                Curve3d[] cc3ds = ids.Select(c => ((Curve)c.GetObject(OpenMode.ForRead)).GetGeCurve()).ToArray(); 
                CompositeCurve3d cc = new CompositeCurve3d(cc3ds); 
                double area = cc.GetArea(cc.GetParameterOf(cc.StartPoint), 
                    cc.GetParameterOf(cc.EndPoint), 
                    Tolerance.Global); 
                ed.WriteMessage(area.ToString()); 
                tr.Commit(); 
            } 
        }

iLisp · 发表于 2014-10-2 12:56:45

哪个准确只有AutoDESK知道，差一点点可以接受，反正在AutoCAD里面也没有别的方法测量

账号		自动登录	找回密码
密码			立即注册

[研讨] XDGE几何库应用(23)--精度对Ge模型影响

马上注册，结交更多好友，享用更多功能，让你轻松玩转社区。

本帖被以下淘专辑推荐: