[讨论]：如何求不规则封闭pline线内最大圆

xiao_longxin

准确地说，不规则多边形没有与每条边都相切的内圆，但只要求出最大的内圆即可，与多少条这相切没有关系。而我要求的这个也不一定是多边形，准确地说是一个封闭的面域，围合线中有可能使用曲线，但可以把它模拟成近似的不规则多边形（如果有不用模拟的算法则最好了）。有没有人做过这方面的偿试？

或者这个题目可以换个说法：
当有一块不规则钢板在手（其边界可能包含有各种类型的曲线，但可以把它模拟成近似的不规则多边形），如何截出一个最大的圆？

WHGF

问题挺新影，关注中

PFD8667

我也想知道,不知那位高手编个程序就好了,先谢谢了!

xiao_longxin

我是想用这个圆心来作为不规则面域的标志点。我想此点应该是这个面域最空旷的地方，是作为标注诸如：面积、地块号之类属性的最佳位置。

chuanxue

很相知道。

cpu1gmhz

画相切圆出来不就可以了吗?

xiao_longxin

楼上的请详细地说明你的思路

fylinwater

我不动脑筋，说个穷尽法，
多边形所有边逐一搭配取三条边求得内切圆， 所有内切圆里半径最大的就是你要的结果

xiao_longxin

楼上的，这个封闭的面域是一条spline线时怎么办呢？

fylinwater

spl本身特性比较复杂，如果要求精度和速度不是很高，那就先把spl模拟成多线段

xiao_longxin

楼上说的是可行的。可是据我所知，三条直线可以做出四个与之同时相切的圆，在这四个圆中，是否会同时有两个或者两个以上的圆的圆心落在面域之内呢？

f5612140

最初由 fylinwater 发布
[B]我不动脑筋，说个穷尽法，
多边形所有边逐一搭配取三条边求得内切圆， 所有内切圆里半径最大的就是你要的结果 [/B]

NO!,頂點"外擴型"是如此,有內縮時不然

xiao_longxin

或者这个题目可以换个说法：
当有一块不规则钢板在手（其边界可能包含有各种类型的曲线，但可以把它模拟成近似的不规则多边形），如何截出一个最大的圆？

xyp1964

三边形只有唯一的内切圆，在偶的“实用程序汇总”里就有：test121。

其他还是让高人解决吧……

xyp1964

没有使用穷尽法，发挥主观能动性大概看看，几下就可解决。
三角形内切圆程序（只适用于LINE线）：
[php](load "xyp_lib");加载通用函数
;|
如果已经下载xyp_lib并放到搜索路径下可以不再下载!
利用以下任何一种方式即可加载和运行通用函数内的所有子程序：
1.在acad.lsp中增加(load"xyp_lib")
2.在每个程序内增加(load"xyp_lib")
3.在command下，输入(load"xyp_lib")
4.在菜单.mnl中增加(load"xyp_lib")
通用函数下载地址：http://www.xdcad.net/forum/showthread.php?s=&threadid=325268
|;

;;;三角形内切圆并标注圆半径和面积
(defun c:test121 ()
  (cmdla0)
  (setq        s1   (usel1 0 "LINE" "LINE线")
        e-s1 (entget (car s1))
  )
  (grdraw (dxf 10 e-s1) (dxf 11 e-s1) 1)
  (setq        s2   (usel1 0 "LINE" "LINE线")
        e-s2 (entget (car s2))
  )
  (grdraw (dxf 10 e-s2) (dxf 11 e-s2) 1)
  (setq        s3   (usel1 0 "LINE" "LINE线")
        e-s3 (entget (car s3))
  )
  (grdraw (dxf 10 e-s3) (dxf 11 e-s3) 1)
  (setq        pt1   (cadr s1)
        pt2   (cadr s2)
        pt3   (cadr s3)
        pt-12 (inters (dxf 10 e-s1)
                      (dxf 11 e-s1)
                      (dxf 10 e-s2)
                      (dxf 11 e-s2)
                      nil
              )
        pt-23 (inters (dxf 10 e-s2)
                      (dxf 11 e-s2)
                      (dxf 10 e-s3)
                      (dxf 11 e-s3)
                      nil
              )
        pt-31 (inters (dxf 10 e-s3)
                      (dxf 11 e-s3)
                      (dxf 10 e-s1)
                      (dxf 11 e-s1)
                      nil
              )
  )
  (if (and pt-12 pt-23 pt-31)
    (progn
      (setq
        a1   (distance pt-12 pt-23)
        a2   (distance pt-23 pt-31)
        a3   (distance pt-31 pt-12)
        s    (/ (+ a1 a2 a3) 2.0)
        a    (rtos (sqrt (* s (- s a1) (- s a2) (- s a3))) 2 3)
        pt12 (_midp pt1 pt2)
        pt23 (_midp pt2 pt3)
      )
      (jzhz)
      (setvar "osmode" 0)
      (command "dimangular" s1 s2 pt12)
      (setq pt12 (osnap pt12 "mid"))
      (command "erase" "l" "")
      (command "dimangular" s2 s3 pt23)
      (setq pt23 (osnap pt23 "mid"))
      (command "erase" "l" "")
      (setq pt-ci (inters pt-12 pt12 pt-23 pt23 nil)
            pt0          (inters pt-ci
                          (polar pt-ci
                                 (+ (angle (dxf 10 e-s1) (dxf 11 e-s1)) (* pi 0.5))
                                 1000
                          )
                          (dxf 10 e-s1)
                          (dxf 11 e-s1)
                          nil
                  )
            rad          (distance pt-ci pt0)
      )
      (mkla "标注" 3)
      (command "circle" pt-ci pt0)
      (setq s           (entlast)
            Sa           (vla-get-Area (obj-e s))
            txt-bz (strcat "R=" (rtos rad 2 3) " ; " "S=" (rtos Sa 2 3))
      )
      (command "text" "j" "mc" pt-ci (* sc 300) "0" txt-bz)
    )
  )
  (redraw)
  (cmdla1)
)[/php]