找回密码
 立即注册

QQ登录

只需一步,快速开始

扫一扫,访问微社区

查看: 3005|回复: 14

[点表] 根据点表求点表构成多边形型心

[复制链接]

已领礼包: 145个

财富等级: 日进斗金

发表于 2013-6-7 14:48:01 | 显示全部楼层 |阅读模式
函数发布
函数名称: XD::Pnts:Centroid
调用格式: (XD::Pnts:Centroid pts)
参数说明: pts ---- 点表
返回值:
函数简介: 根据点表求点表构成多边形型心,不允许自交
函数来源: 原创
函数作者: XDSoft
适用版本: 不限 
最后更新时间: 2013-06-07
备注: 顶点组成的多边形不允许自交。
演示图片: -

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

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

×
  1. ;| 多边形形心
  2. 根据点表求点表构成多边形型心,不允许自交
  3. 参数 pts ---- 点表
  4. 返回值: 点
  5. |;
  6. (defun XD::Pnts:Centroid (pts / ar isum p0 p1 p2 pts1 x0 x1 x2 y0 y1 y2)
  7.   (defun polygon_sumx (pts / pts1 isum p0 p1 p2 x0 x1 x2 y0 y1 y2)
  8.     (setq pts1 pts)
  9.     (setq isum 0
  10.           p0 (car pts1)
  11.           x0 (car p0)
  12.           y0 (cadr p0)
  13.     )
  14.     (while (cdr pts1)
  15.       (setq p1 (car pts1)
  16.             x1 (car p1)
  17.             y1 (cadr p1)
  18.             p2 (cadr pts1)
  19.             x2 (car p2)
  20.             y2 (cadr p2)
  21.             isum (+ isum (* (+ x1 x2) (- (* x1 y2) (* x2 y1))))
  22.             pts1 (cdr pts1)
  23.       )
  24.     )
  25.     (+ isum (* (+ (caar pts1) x0) (- (* (caar pts1) y0) (* x0 (cadar pts1)))))
  26.   )
  27.   (defun polygon_sumy (pts / pts1 isum p0 p1 p2 x0 x1 x2 y0 y1 y2)
  28.     (setq pts1 pts)
  29.     (setq isum 0
  30.           p0 (car pts1)
  31.           x0 (car p0)
  32.           y0 (cadr p0)
  33.     )
  34.     (while (cdr pts1)
  35.       (setq p1 (car pts1)
  36.             x1 (car p1)
  37.             y1 (cadr p1)
  38.             p2 (cadr pts1)
  39.             x2 (car p2)
  40.             y2 (cadr p2)
  41.             isum (+ isum (* (+ y1 y2) (- (* x1 y2) (* x2 y1))))
  42.             pts1 (cdr pts1)
  43.       )
  44.     )
  45.     (+ isum (* (+ (cadar pts1) y0) (- (* (caar pts1) y0) (* x0 (cadar pts1)))))
  46.   )
  47.   (setq ar (XD::Pnts:Area pts))
  48.   (list (/ (polygon_sumx pts) (* 6 ar)) (/ (polygon_sumy pts) (* 6 ar)) 0.0)
  49. )





评分

参与人数 1D豆 +5 收起 理由
xshrimp + 5 技术引导讨论和指点奖!

查看全部评分

论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!

已领礼包: 8121个

财富等级: 富甲天下

发表于 2013-6-7 15:00:12 | 显示全部楼层
本帖最后由 Highflybird 于 2013-6-7 17:25 编辑

这个是我的方式:

[pcode=lisp,true];;;----------------------------------------------------;
;;;功能: 获取多边形信息(质心,面积,周长)              ;
;;;输入: Pts---多边形顶点列表                          ;
;;;输出: 列表:第一个为多边形的面积中心(质心),用2d点表示;
;;;      第二个为数值,正数表示多边形方向是CCW(逆时针) ;
;;;      负数表示顺时针;第三个为周长.                 ;
;;;参考: http://en.wikipedia.org/wiki/Centroid        ;
;;;----------------------------------------------------;
(defun POLY:Infomation (Pts / Pts1 Ai S lst cen)
  (setq Pts1 (MISC:1st->Last Pts))                      ;another point of every side
  (setq Ai   (mapcar 'MAT:Det2V Pts Pts1))              ;area of every side
  (setq S    (* (apply '+ Ai) 0.5))                     ;Total area
  (Setq Cen  (MAT:SxVs (mapcar 'MAT:v+v Pts Pts1) Ai))
  (setq Cen  (MAT:vxs Cen (/ 0.166666666666666667 S)))  ;base on the formula
  (list Cen S (apply '+ (mapcar 'distance pts pts1)))   ;Return Centroid,Total area and Perimeter
)[/pcode]

一些函数参见我的矩阵和关于几何计算的帖子:
http://bbs.xdcad.net/forum.php?mod=viewthread&tid=667495&reltid=668705&pre_thread_id=0&pre_pos=1&ext=

点评

这个要等这些矩阵和向量函数整理出来才能用,里面一些矩阵的习惯我认为需要调整下。  详情 回复 发表于 2013-6-7 15:19

评分

参与人数 2D豆 +17 收起 理由
XDSoft + 12 很给力!经验;技术要点;资料分享奖!
xshrimp + 5 技术引导讨论和指点奖!

查看全部评分

论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 145个

财富等级: 日进斗金

 楼主| 发表于 2013-6-7 15:19:52 | 显示全部楼层
Highflybird 发表于 2013-6-7 15:00
这个是我的方式:

;;;----------------------------------------------------;

这个要等这些矩阵和向量函数整理出来才能用,里面一些矩阵的习惯我认为需要调整下。
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 1个

财富等级: 恭喜发财

发表于 2013-6-7 16:40:12 | 显示全部楼层
(XD::Pnts:Area pts)  
---where

点评

就在今天发的帖子里面,找找,更快的方式是点击 “函数发布” 链接。  详情 回复 发表于 2013-6-7 17:49
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

发表于 2013-6-7 17:14:38 | 显示全部楼层
高飞版主的没有测试,LZ的程序获取的型心和用Region获取的有出入,不知是误差所致还是计算方法问题

点评

代码是按照数学方法求的,等验证下,大家也可以帮测试下。  详情 回复 发表于 2013-6-7 17:52
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 51个

财富等级: 招财进宝

发表于 2013-6-7 17:49:32 | 显示全部楼层
dwg001 发表于 2013-6-7 16:40
(XD:nts:Area pts)  
---where

就在今天帖子里面,找找,更快的方式是点击 “函数发布” 链接。
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 145个

财富等级: 日进斗金

 楼主| 发表于 2013-6-7 17:52:06 | 显示全部楼层
Free-Lancer 发表于 2013-6-7 17:14
高飞版主的没有测试,LZ的程序获取的型心和用Region获取的有出入,不知是误差所致还是计算方法问题

代码是按照数学方法求的,等验证下,大家也可以帮测试下。
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 1632个

财富等级: 堆金积玉

发表于 2013-6-7 21:24:52 | 显示全部楼层
以前很多人都是借鉴质心函数讨论,数学计算的方法另有一凡意义
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 2476个

财富等级: 金玉满堂

发表于 2013-6-16 18:01:02 | 显示全部楼层
发布函数是不是也给个范例,新手的理解能力很差,完全不会用呀!

点评

告诉你了函数名,也告诉你了参数是 点表 ,还理解不了啊?  详情 回复 发表于 2013-6-16 19:13
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 19个

财富等级: 恭喜发财

发表于 2013-6-16 19:13:44 | 显示全部楼层
q3_2006 发表于 2013-6-16 18:01
发布函数是不是也给个范例,新手的理解能力很差,完全不会用呀!

告诉你了函数名,也告诉你了参数是 点表 ,还理解不了啊?

点评

的确没理解呀,i不如楼上的来段范例!  详情 回复 发表于 2013-6-16 19:48
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 2476个

财富等级: 金玉满堂

发表于 2013-6-16 19:47:37 | 显示全部楼层
我还没入门,当然看不懂。意思是说能求随意几个点构成多边形的形心吗?的确没理解呀,i不如楼上的来段范例!
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 2476个

财富等级: 金玉满堂

发表于 2013-6-16 19:48:52 | 显示全部楼层
Lisphk 发表于 2013-6-16 19:13
告诉你了函数名,也告诉你了参数是 点表 ,还理解不了啊?

的确没理解呀,i不如楼上的来段范例!

点评

你这个还真不好教,自己多看看论坛帖子吧。  详情 回复 发表于 2013-6-16 20:17
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 19个

财富等级: 恭喜发财

发表于 2013-6-16 20:17:01 | 显示全部楼层
q3_2006 发表于 2013-6-16 19:48
的确没理解呀,i不如楼上的来段范例!

你这个还真不好教,自己多看看论坛帖子吧。

点评

我的表达有问题,这个函数怎么用已经明白了,问题在之前“顶点组成的多边形不允许自交 ”,也就是我获取的点表是有顺序的,我不知道如何给点表排序,还请指教一下,谢谢了!  详情 回复 发表于 2013-6-17 06:02
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

已领礼包: 2476个

财富等级: 金玉满堂

发表于 2013-6-17 06:02:01 | 显示全部楼层
Lisphk 发表于 2013-6-16 20:17
你这个还真不好教,自己多看看论坛帖子吧。

我的表达有问题,这个函数怎么用已经明白了,问题在之前“顶点组成的多边形不允许自交 ”,也就是我获取的点表是有顺序的,我不知道如何给点表排序,还请指教一下,谢谢了!
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

发表于 2015-9-5 15:28:10 | 显示全部楼层
{:soso_e160:}{:soso_e161:}{:soso_e156:}
论坛插件加载方法
发帖求助前要善用【论坛搜索】功能,那里可能会有你要找的答案;
如果你在论坛求助问题,并且已经从坛友或者管理的回复中解决了问题,请把帖子标题加上【已解决】;
如何回报帮助你解决问题的坛友,一个好办法就是给对方加【D豆】,加分不会扣除自己的积分,做一个热心并受欢迎的人!
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

QQ|申请友链|Archiver|手机版|小黑屋|辽公网安备|晓东CAD家园 ( 辽ICP备15016793号 )

GMT+8, 2024-3-29 04:31 , Processed in 0.452145 second(s), 59 queries , Gzip On.

Powered by Discuz! X3.5

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表