如何求最小包围盒？

Lispboy

守仁格竹GM

，这是最小外接面积矩形源码，不知道怎么设置权限，积分超过多少才能看，设置1个豆吧，以示尊重吧[/sell]Lispboy 发表于 2013-5-18 20:48
这话怎么说，咋在你眼里还成习惯了

{:soso_e113:} 试试下面的行不行，有点忘记 这个只适合多线段 别的暂时没考虑到 可以给你做个参考

[sell=1]
[pcode=lisp,true];===================================================================
;;;两矢量的叉积
(defun G:CrossProductor        (vec1 vec2 / aaa bbb ccc ddd eee fff)                 
  (setq aaa (car   vec1))
  (setq bbb (cadr  vec1))
  (setq ccc (caddr vec1))
  (setq ddd (car   vec2))
  (setq eee (cadr  vec2))
  (setq fff (caddr vec2))
  (list
    (- (* bbb fff) (* ccc eee))
    (- (* ccc ddd) (* aaa fff))
    (- (* aaa eee) (* bbb ddd))
  )
)
;;;矩阵的变换与逆变换
(defun GetMatrix (lst org Revflag / mat iii jjj)
  (setq mat (vlax-make-safearray vlax-vbDouble '(0 . 3) '(0 . 3)))
  (setq iii 0)
  (repeat 3
    (vlax-safearray-put-element mat iii 3 (nth iii org))         
    (setq jjj 0)
    (repeat 3                                                     
      (if RevFlag
        (vlax-safearray-put-element mat iii jjj (nth jjj (nth iii lst)))
        (vlax-safearray-put-element mat iii jjj (nth iii (nth jjj lst)))  
      );endif
      (setq jjj (+ jjj 1))                                             
    );endrepeat
    (setq iii (+ iii 1))
  );endrepeat
  (vlax-safearray-put-element mat 3 3 1)                                      
  mat                                                            
)
;;;构造矩形
(defun Make-Rectange (pt1 pt2)
  (entmake
    (list
      '(0 . "LWPOLYLINE")                                    
      '(100 . "AcDbEntity")                              
      '(100 . "AcDbPolyline")
      '(90 . 4)                                             
      '(70 . 1)                                             
      (cons 38 (caddr pt1))                                 
      (cons 10 (list (car pt1) (cadr pt1)))                  
      (cons 10 (list (car pt2) (cadr pt1)))                  
      (cons 10 (list (car pt2) (cadr pt2)))                  
      (cons 10 (list (car pt1) (cadr pt2)))                     
      (cons 210 '(0 0 1))                                       
    )
  )
)
;;;主程序
(vl-load-com)
(defun c:zxj ()
  (setq sel (ssget '((0 . "lwpolyline"))))              
  (setq i 0)
  (repeat (sslength sel)
     (setq num 0)
     (setq pts '() pts1 '())
     ;(setq mianji 1e309)
     (setq mianji 1e308)
     (setq ent0 (ssname sel i))                           
     (setq obj0 (vlax-ename->vla-object ent0))               
     (setq plinedata0 (entget ent0))
     (foreach x plinedata0
        (if (= (car x) 10)
            (progn
               (setq pts (append pts (list (cdr x))))
               (setq num (+ num 1))
            );endprogn   
        );endif
     );endforeach
;------------------
     (setq pts1 (append pts (list (nth 0 pts))))
     (setq j 0)
     (command "RECTANG" '(0 0 0) '(1000 2000 0))
     (repeat num
       (setq p_1 (nth j pts1))
       (setq p_2 (nth (+ j 1) pts1))
       (command "ucs" "z" p_1 p_2)
       (setq UcsFlag (getvar "WORLDUCS"))                           
       (if (= UcsFlag 0)                                               
          (progn                                             
             (setq UcsFlag T)                                            
             (setq xdir (getvar "UCSXDIR"))                              
             (setq ydir (getvar "UCSYDIR"))                              
             (setq zdir (G:CrossProductor xdir ydir))                     
             (setq origin  (getvar "UCSORG"))                        
             (setq WcsOrg  (trans '(0 0 0) 0 1))                           
             (setq matLst  (list xdir ydir zdir))                       
             (setq matrix  (GetMatrix matLst origin nil))              
             (setq revMat  (GetMatrix matLst WcsOrg T))              
          );endprogn
          (setq UcsFlag nil)                                          
       );endif
;;=======================================================
       (setq ent (ssname sel i))                              
       (setq obj (vlax-ename->vla-object ent))                        
       (and UcsFlag (vla-TransformBy obj revMat))               
       (vla-GetBoundingBox obj 'point1 'point2)                    
       (setq point1 (vlax-safearray->list point1))
       (setq point2 (vlax-safearray->list point2))  
       (and UcsFlag (vla-TransformBy obj matrix))              
       (setq point10x (car point1))
       (setq point10y (cadr point1))
       (setq point20x (car point2))
       (setq point20y (cadr point2))
       (setq linex (abs (- point20x point10x)))
       (setq liney (abs (- point20y point10y)))
       (if (<= (* linex liney) mianji)
           (progn
              (setq mianji (* linex liney))
              (entdel (entlast))
              (and
                  (make-Rectange (list point10x point10y 0) (list point20x point20y 0))
                  UcsFlag                                             
                  (vla-TransformBy (vlax-ename->vla-object (entlast)) matrix)
              )
           );endprogn
       );endif
       (setq j (+ j 1))
       (command "ucs" "w")
       );endrepea
    (command "ucs" "w")                             
    (setq i (+ i 1))
  );endrepea
  (command "ucs" "w")
  (princ)
);enddefun
[/pcode]
[/sell]

守仁格竹GM

最小面积矩形 一定过某一边 求出 算出最小值

守仁格竹GM

楼主多少 自己也得学习一下  别要东西养成习惯

Lispboy

守仁格竹GM 发表于 2013-5-18 20:33
最小面积矩形 一定过某一边 求出 算出最小值

:L
怎么求那某一条边？

Lispboy

守仁格竹GM 发表于 2013-5-18 20:34
楼主多少 自己也得学习一下  别要东西养成习惯

这话怎么说，咋在你眼里还成习惯了:dizzy:

wowan1314

以下是个人一点理解！也许根本不正确，供参考！
如果是文字或者规则图形！可以把文字旋转后求包围盒！然后再转回来！
但如楼主图示那样的图，只想到设定旋转角度ang！旋转后求包围盒与之前比较，一直转一圈后得出最小矩阵！其精度取决于ang的值的大小

Free-Lancer

这个 高飞鸟 写过滴

liuyun242

wowan1314 发表于 2013-5-18 21:14
以下是个人一点理解！也许根本不正确，供参考！
如果是文字或者规则图形！可以把文字旋转后求包围盒！然后 ...

我也是怎么认为的，只是感觉效率不高

Lispboy

Free-Lancer 发表于 2013-5-18 21:17
这个 高飞鸟 写过滴

论坛信息量太大，可否提供个链接，谢谢！

守仁格竹GM

这个肯定是比把物体旋转 效率高得多，那种试过批量很慢。

守仁格竹GM

Lispboy 发表于 2013-5-18 20:48
这话怎么说，咋在你眼里还成习惯了

:lol 哈哈 看错人名了 有个跟你比较像 扫瑞  看看我的回复吧 希望对你有帮助

Lispboy

守仁格竹GM 发表于 2013-5-18 22:29
哈哈 看错人名了 有个跟你比较像 扫瑞  看看我的回复吧 希望对你有帮助

谢谢，代码我试试。以后别在认错人了，差点冤枉死我。:dizzy:

对了，看到代码里面用到矩阵，是高飞鸟斑竹的矩阵代码吗？程序是你写的？

守仁格竹GM

Lispboy 发表于 2013-5-19 00:01
谢谢，代码我试试。以后别在认错人了，差点冤枉死我。

对了，看到代码里面用到矩阵，是高飞鸟 ...

不是我写的 是网上常用的 你能看明白，然后会用就可以了

haoryh

路过顶一下！