[从头学数学] 第276节 [计算几何] 中心点与顶点的凹凸性

最新推荐文章于 2022-11-14 13:29:45 发布
最新推荐文章于 2022-11-14 13:29:45 发布
阅读量712 收藏
点赞数
分类专栏： 从头学数学 从头学数学
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/mwsister/article/details/52659072
版权
剧情提要：
阿伟看到了一本比较有趣的书，是关于《计算几何》的，2008年由北清派出版。很好奇
它里面讲了些什么，就来看看啦。


正剧开始：
星历2016年09月25日 11:26:01, 银河系厄尔斯星球中华帝国江南行省。

[工程师阿伟]正在和[机器小伟]一起研究[计算几何]]。




<span style="font-size:18px;">#
>>> 
顶点数量： 8
顶点集： [Point([6, 9]), Point([2.4, 5.4]), Point([0.8, 3.8]), Point([1.33, 3]), Point([2, 3]), Point([6, 3]), Point([6, 3.86]), Point([3.33, 5])]
边集： [SegLine([2.4, 5.4], [6, 9]), SegLine([0.8, 3.8], [2.4, 5.4]), SegLine([1.33, 3], [0.8, 3.8]), SegLine([1.33, 3], [2, 3]), SegLine([2, 3], [6, 3]), SegLine([6, 3], [6, 3.86]), SegLine([6, 3.86], [3.33, 5]), SegLine([3.33, 5], [6, 9])] 

边集
[1.33, 3]: [Point([2, 3]), Point([0.8, 3.8])], len =2
[2, 3]: [Point([1.33, 3]), Point([6, 3])], len =2
[6, 3]: [Point([2, 3]), Point([6, 3.86])], len =2
[0.8, 3.8]: [Point([1.33, 3]), Point([2.4, 5.4])], len =2
[6, 3.86]: [Point([6, 3]), Point([3.33, 5])], len =2
[3.33, 5]: [Point([6, 3.86]), Point([6, 9])], len =2
[2.4, 5.4]: [Point([0.8, 3.8]), Point([6, 9])], len =2
[6, 9]: [Point([3.33, 5]), Point([2.4, 5.4])], len =2
顶点集
[Point([1.33, 3]), Point([2, 3]), Point([6, 3]), Point([0.8, 3.8]), Point([6, 3.86]), Point([3.33, 5]), Point([2.4, 5.4]), Point([6, 9])]


>>> 

>>> 
操作用时：3.159e-02 s
[[3.48, 4.51], [3.73, 4.51], [2.23, 4.13], [2.21, 4.32], [3.07, 4.26], [2.73, 4.51], [5.22, 2.3], [5.72, 2.41], [1.09, 4.99], [4.13, 3.26], [5.62, 1.19], [5.67, 2.11], [5.48, 1.62], [6.12, 1.29], [5.98, 1.73], [6.1, 1.16], [6.4, 0.25], [6.38, 0.17], [6.27, 0.67], [6.06, 0.98]]
>>> 

def tmp3():
    path = [[[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [6, 3], [6, 3.86], [3.33, 5], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [3.33, 5], [6, 3.86], [8, 3], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.5, 1.5], [1.78, 2.33], [2, 3], [3.33, 5], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [1.78, 2.33], [2, 3], [3.33, 5], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [2, 3], [6, 3], [6, 3.86], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [2, 3], [3.33, 5], [6, 3.86], [6, 9]], [[6, 9], [2.4, 5.4], [3.33, 5], [2, 3], [6, 3], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [2.4, 5.4], [3.33, 5], [6, 3.86], [6, 3], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [2.4, 5.4], [-6, 9], [-0.37, 2.56], [0, 3], [1.33, 3], [2, 3], [3.33, 5], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [1.78, 2.33], [3.67, -0.5], [4.29, 0.43], [6, 3], [6, 3.86], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [6, 3], [6, -3], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [6, 3], [6, -3], [8, -7], [8, 3], [6, 3.86], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [6, 3], [4.5, 0], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [3.33, 5], [6, 3.86], [6, 3], [6, -3], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [3.33, 5], [6, 3.86], [6, 3], [4.5, 0], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [4.29, 0.43], [4.5, 0], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [6, -3], [5.56, -3.33], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [6, -3], [5.67, -3.5], [5.35, -3.97], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [4.5, 0], [6, -3], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [4.5, 0], [4, -1], [6, -3], [8, -7], [8, 3], [6, 9]]]
    len_path = len(path);

    centerArray = [];

    #计时开始
    startTime = time.clock();
    
    for i in range(len_path):
        poly = Polygon(path[i]);
        poly.genVertexTree();

        centerArray.append(poly.centerPoint());

    #计时结束
    endTime = time.clock();

    #打印结果
    print('操作用时：{0:.3e} s'.format(endTime-startTime));   

    print(centerArray);

#</span>

 
<span style="font-size:18px;">//
	if (1) {
		var r = 20;      
        config.setSector(1,1,1,1);        
        config.graphPaper2D(0, 0, r);      
        config.axis2D(0, 0, 250, 1.2);   
		
		//坐标轴设定        
        var scaleX = 2*r, scaleY = 2*r;          
        var spaceX = 2, spaceY = 2;           
        var xS = -10, xE = 10;          
        var yS = -10, yE = 10;          
        config.axisSpacing(xS, xE, spaceX, scaleX, 'X');            
        config.axisSpacing(yS, yE, spaceY, scaleY, 'Y');            
                    
        var transform = new Transform();    
		
		//顶点
		var a = [];
		
		for (var i = 0; i < $vertex.length; i++) {
			a.push($vertex[i][0]);
		}
		
		//显示变换        
        if (a.length > 0) {        
            a = transform.scale(transform.translate(a, 0, 0), scaleX/spaceX, scaleY/spaceY);
		}
		
		var lable = [];
		for (var i = 0; i < 100; i++) {
			lable.push(i.toFixed(0));
		}
		
		/*
		//边集
		var b = [];
		for (var i = 0; i < $seg.length; i++) {
			b.push([a[$seg[i][0]], a[$seg[i][1]]]);
		}
		
				
		var edges = b.length;
		for (var i = 0; i < edges; i++) {
			shape.multiLineDraw([].concat(b[i]), 'red');
		}*/
		
		var colorArray = ['red', 'orange', 'yellow', 'green', 'cyan', 'blue', 'purple'];
		var seg = [];		
		
		var idx = xGlobal;
		var len = $path2.length;

		plot.setLineWidth(3);

		seg = transform.scale(transform.translate($path2[idx%len], 0, 0), scaleX/spaceX, scaleY/spaceY);
		shape.multiLineDraw([].concat(seg), 'pink');
		
		cPoint = transform.scale(transform.translate($center, 0, 0), scaleX/spaceX, scaleY/spaceY);
		
		shape.pointDraw([].concat([cPoint[idx%len]]), 'red');
		

		/*
		//次要顶点
		var b = [];
		
		for (var i = 0; i < 6; i++) {	
			b = [].concat($pointType[n][i]);
			//显示变换        
			if (b.length > 0) {        
				b = transform.scale(transform.translate(b, 0, 0), scaleX/spaceX, scaleY/spaceY);
				shape.pointDraw([].concat(b), colorArray[(i-1)%7]);
			}

			
			
		}*/
		
		//主要顶点
		shape.pointDraw([].concat(a), 'blue', 1, 1, lable);		

	}

//</span>


下面各图中红色的为中心点，橙色的为凹顶点，事实证明，中心点的位置与多边形的凹凸性，与顶点的凹凸性，没有任何关系。


路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [6, 3], [6, 3.86], [3.33, 5]]
凹点： [[3.33, 5]]
凸点： [[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [6, 3], [6, 3.86]]


路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [3.33, 5], [6, 3.86], [8, 3]]
凹点： [[3.33, 5], [6, 3.86]]
凸点： [[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [8, 3]]



路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.5, 1.5], [1.78, 2.33], [2, 3], [3.33, 5]]
凹点： [[0.8, 3.8], [2, 3], [3.33, 5]]
凸点： [[6, 9], [2.4, 5.4], [0, 3], [1.5, 1.5], [1.78, 2.33]]



路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [1.78, 2.33], [2, 3], [3.33, 5]]
凹点： [[0.8, 3.8], [1.33, 3], [2, 3], [3.33, 5]]
凸点： [[6, 9], [2.4, 5.4], [0, 3], [1.78, 2.33]]



路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [2, 3], [6, 3], [6, 3.86]]
凹点： [[0.8, 3.8]]
凸点： [[6, 9], [2.4, 5.4], [0, 3], [6, 3]]



路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [2, 3], [3.33, 5], [6, 3.86]]
凹点： [[0.8, 3.8], [3.33, 5]]
凸点： [[6, 9], [2.4, 5.4], [0, 3], [2, 3], [6, 3.86]]



路径： [[6, 9], [2.4, 5.4], [3.33, 5], [2, 3], [6, 3], [6, -3], [8, -7], [8, 3]]
凹点： [[3.33, 5], [6, 3]]
凸点： [[6, 9], [2.4, 5.4], [2, 3], [6, -3], [8, -7], [8, 3]]



路径： [[6, 9], [2.4, 5.4], [3.33, 5], [6, 3.86], [6, 3], [6, -3], [8, -7], [8, 3]]
凹点： [[6, 3.86]]
凸点： [[6, 9], [2.4, 5.4], [3.33, 5], [6, -3], [8, -7], [8, 3]]



路径： [[6, 9], [2.4, 5.4], [-6, 9], [-0.37, 2.56], [0, 3], [1.33, 3], [2, 3], [3.33, 5]]
凹点： [[2.4, 5.4], [0, 3], [3.33, 5]]
凸点： [[6, 9], [-6, 9], [-0.37, 2.56], [2, 3]]



路径： [[6, 9], [3.33, 5], [2, 3], [1.78, 2.33], [3.67, -0.5], [4.29, 0.43], [6, 3], [6, 3.86]]
凹点： []
凸点： [[6, 9], [3.33, 5], [2, 3], [1.78, 2.33], [3.67, -0.5], [4.29, 0.43], [6, 3]]



路径： [[6, 9], [3.33, 5], [2, 3], [6, 3], [6, -3], [5.67, -3.5], [8, -7], [8, 3]]
凹点： [[6, 3], [6, -3]]
凸点： [[6, 9], [3.33, 5], [2, 3], [5.67, -3.5], [8, -7], [8, 3]]



路径： [[6, 9], [3.33, 5], [2, 3], [6, 3], [6, -3], [8, -7], [8, 3], [6, 3.86]]
凹点： [[6, 3], [6, 3.86]]
凸点： [[6, 9], [3.33, 5], [2, 3], [6, -3], [8, -7], [8, 3]]



路径： [[6, 9], [3.33, 5], [2, 3], [6, 3], [4.5, 0], [6, -3], [8, -7], [8, 3]]
凹点： [[6, 3]]
凸点： [[6, 9], [3.33, 5], [2, 3], [4.5, 0], [8, -7], [8, 3]]



路径： [[6, 9], [3.33, 5], [6, 3.86], [6, 3], [6, -3], [5.67, -3.5], [8, -7], [8, 3]]
凹点： [[6, 3.86], [6, -3]]
凸点： [[6, 9], [3.33, 5], [5.67, -3.5], [8, -7], [8, 3]]



路径： [[6, 9], [3.33, 5], [6, 3.86], [6, 3], [4.5, 0], [6, -3], [8, -7], [8, 3]]
凹点： [[6, 3.86], [6, 3]]
凸点： [[6, 9], [3.33, 5], [4.5, 0], [8, -7], [8, 3]]



路径： [[6, 9], [6, 3.86], [6, 3], [4.29, 0.43], [4.5, 0], [6, -3], [8, -7], [8, 3]]
凹点： [[6, 3]]
凸点： [[6, 9], [4.29, 0.43], [4.5, 0], [8, -7], [8, 3]]



路径： [[6, 9], [6, 3.86], [6, 3], [6, -3], [5.56, -3.33], [5.67, -3.5], [8, -7], [8, 3]]
凹点： [[6, -3]]
凸点： [[6, 9], [5.56, -3.33], [5.67, -3.5], [8, -7], [8, 3]]



路径： [[6, 9], [6, 3.86], [6, 3], [6, -3], [5.67, -3.5], [5.35, -3.97], [8, -7], [8, 3]]
凹点： [[6, -3], [5.67, -3.5]]
凸点： [[6, 9], [5.35, -3.97], [8, -7], [8, 3]]



路径： [[6, 9], [6, 3.86], [6, 3], [4.5, 0], [6, -3], [5.67, -3.5], [8, -7], [8, 3]]
凹点： [[6, 3], [6, -3]]
凸点： [[6, 9], [4.5, 0], [5.67, -3.5], [8, -7], [8, 3]]



路径： [[6, 9], [6, 3.86], [6, 3], [4.5, 0], [4, -1], [6, -3], [8, -7], [8, 3]]
凹点： [[6, 3], [6, -3]]
凸点： [[6, 9], [4, -1], [8, -7], [8, 3]]


<span style="font-size:18px;">#
>>> 
路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [6, 3], [6, 3.86], [3.33, 5]]
凹点： [[3.33, 5]]
凸点： [[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [6, 3], [6, 3.86]]


路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [3.33, 5], [6, 3.86], [8, 3]]
凹点： [[3.33, 5], [6, 3.86]]
凸点： [[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [8, 3]]


路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.5, 1.5], [1.78, 2.33], [2, 3], [3.33, 5]]
凹点： [[0.8, 3.8], [2, 3], [3.33, 5]]
凸点： [[6, 9], [2.4, 5.4], [0, 3], [1.5, 1.5], [1.78, 2.33]]


路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [1.78, 2.33], [2, 3], [3.33, 5]]
凹点： [[0.8, 3.8], [1.33, 3], [2, 3], [3.33, 5]]
凸点： [[6, 9], [2.4, 5.4], [0, 3], [1.78, 2.33]]


路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [2, 3], [6, 3], [6, 3.86]]
凹点： [[0.8, 3.8]]
凸点： [[6, 9], [2.4, 5.4], [0, 3], [6, 3]]


路径： [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [2, 3], [3.33, 5], [6, 3.86]]
凹点： [[0.8, 3.8], [3.33, 5]]
凸点： [[6, 9], [2.4, 5.4], [0, 3], [2, 3], [6, 3.86]]


路径： [[6, 9], [2.4, 5.4], [3.33, 5], [2, 3], [6, 3], [6, -3], [8, -7], [8, 3]]
凹点： [[3.33, 5], [6, 3]]
凸点： [[6, 9], [2.4, 5.4], [2, 3], [6, -3], [8, -7], [8, 3]]


路径： [[6, 9], [2.4, 5.4], [3.33, 5], [6, 3.86], [6, 3], [6, -3], [8, -7], [8, 3]]
凹点： [[6, 3.86]]
凸点： [[6, 9], [2.4, 5.4], [3.33, 5], [6, -3], [8, -7], [8, 3]]


路径： [[6, 9], [2.4, 5.4], [-6, 9], [-0.37, 2.56], [0, 3], [1.33, 3], [2, 3], [3.33, 5]]
凹点： [[2.4, 5.4], [0, 3], [3.33, 5]]
凸点： [[6, 9], [-6, 9], [-0.37, 2.56], [2, 3]]


路径： [[6, 9], [3.33, 5], [2, 3], [1.78, 2.33], [3.67, -0.5], [4.29, 0.43], [6, 3], [6, 3.86]]
凹点： []
凸点： [[6, 9], [3.33, 5], [2, 3], [1.78, 2.33], [3.67, -0.5], [4.29, 0.43], [6, 3]]


路径： [[6, 9], [3.33, 5], [2, 3], [6, 3], [6, -3], [5.67, -3.5], [8, -7], [8, 3]]
凹点： [[6, 3], [6, -3]]
凸点： [[6, 9], [3.33, 5], [2, 3], [5.67, -3.5], [8, -7], [8, 3]]


路径： [[6, 9], [3.33, 5], [2, 3], [6, 3], [6, -3], [8, -7], [8, 3], [6, 3.86]]
凹点： [[6, 3], [6, 3.86]]
凸点： [[6, 9], [3.33, 5], [2, 3], [6, -3], [8, -7], [8, 3]]


路径： [[6, 9], [3.33, 5], [2, 3], [6, 3], [4.5, 0], [6, -3], [8, -7], [8, 3]]
凹点： [[6, 3]]
凸点： [[6, 9], [3.33, 5], [2, 3], [4.5, 0], [8, -7], [8, 3]]


路径： [[6, 9], [3.33, 5], [6, 3.86], [6, 3], [6, -3], [5.67, -3.5], [8, -7], [8, 3]]
凹点： [[6, 3.86], [6, -3]]
凸点： [[6, 9], [3.33, 5], [5.67, -3.5], [8, -7], [8, 3]]


路径： [[6, 9], [3.33, 5], [6, 3.86], [6, 3], [4.5, 0], [6, -3], [8, -7], [8, 3]]
凹点： [[6, 3.86], [6, 3]]
凸点： [[6, 9], [3.33, 5], [4.5, 0], [8, -7], [8, 3]]


路径： [[6, 9], [6, 3.86], [6, 3], [4.29, 0.43], [4.5, 0], [6, -3], [8, -7], [8, 3]]
凹点： [[6, 3]]
凸点： [[6, 9], [4.29, 0.43], [4.5, 0], [8, -7], [8, 3]]


路径： [[6, 9], [6, 3.86], [6, 3], [6, -3], [5.56, -3.33], [5.67, -3.5], [8, -7], [8, 3]]
凹点： [[6, -3]]
凸点： [[6, 9], [5.56, -3.33], [5.67, -3.5], [8, -7], [8, 3]]


路径： [[6, 9], [6, 3.86], [6, 3], [6, -3], [5.67, -3.5], [5.35, -3.97], [8, -7], [8, 3]]
凹点： [[6, -3], [5.67, -3.5]]
凸点： [[6, 9], [5.35, -3.97], [8, -7], [8, 3]]


路径： [[6, 9], [6, 3.86], [6, 3], [4.5, 0], [6, -3], [5.67, -3.5], [8, -7], [8, 3]]
凹点： [[6, 3], [6, -3]]
凸点： [[6, 9], [4.5, 0], [5.67, -3.5], [8, -7], [8, 3]]


路径： [[6, 9], [6, 3.86], [6, 3], [4.5, 0], [4, -1], [6, -3], [8, -7], [8, 3]]
凹点： [[6, 3], [6, -3]]
凸点： [[6, 9], [4, -1], [8, -7], [8, 3]]


操作用时：5.985e-01 s
[[3.48, 4.51], [3.73, 4.51], [2.23, 4.13], [2.21, 4.32], [3.07, 4.26], [2.73, 4.51], [5.22, 2.3], [5.72, 2.41], [1.09, 4.99], [4.13, 3.26], [5.62, 1.19], [5.67, 2.11], [5.48, 1.62], [6.12, 1.29], [5.98, 1.73], [6.1, 1.16], [6.4, 0.25], [6.38, 0.17], [6.27, 0.67], [6.06, 0.98]]
[[[3.33, 5]], [[3.33, 5], [6, 3.86]], [[0.8, 3.8], [2, 3], [3.33, 5]], [[0.8, 3.8], [1.33, 3], [2, 3], [3.33, 5]], [[0.8, 3.8]], [[0.8, 3.8], [3.33, 5]], [[3.33, 5], [6, 3]], [[6, 3.86]], [[2.4, 5.4], [0, 3], [3.33, 5]], [], [[6, 3], [6, -3]], [[6, 3], [6, 3.86]], [[6, 3]], [[6, 3.86], [6, -3]], [[6, 3.86], [6, 3]], [[6, 3]], [[6, -3]], [[6, -3], [5.67, -3.5]], [[6, 3], [6, -3]], [[6, 3], [6, -3]]]
>>> 

    #计算中心点
    def centerPoint(self):
        xCenter = 0;
        yCenter = 0;

        for i in range(self.vertNum):
            p_ = self.vertex[i].value();

            xCenter += p_[0];
            yCenter += p_[1];

        return [round(xCenter/self.vertNum, 2), round(yCenter/self.vertNum, 2)];

    #判断三角形是顺时针还是逆时针
    def judgeClockwise(self, P1, P2, P3):
        #三个点都是[x, y]形式，即投影在平面上
        x1, y1, x2, y2, x3, y3 = P1[0], P1[1], P2[0], P2[1], P3[0], P3[1];
        cross = (x2-x1)*(y3-y1)-(y2-y1)*(x3-x1);

        if (cross > 0):
            #逆时针
            return 1;
        elif (cross < 0):
            #顺时针
            return -1;
        else:
            return 0;

    #判断凹多边形，顶点的凹凸性
    def judgeConcavePolygon(self, array):
        newArray = array;
        vertexNum = len(newArray);
        #多边形上的凹点
        concavePoint = [];
        #多边形上的凸点
        convexPoint = [];
        #返回矩阵
        returnArray = [];

        if (vertexNum >= 3):
            for i in range(vertexNum):
                P1 = newArray[(i-1+vertexNum)%vertexNum];
                P2 = newArray[i%vertexNum];
                P3 = newArray[(i+1)%vertexNum];

                if (self.judgeClockwise(P1, P2, P3) < 0): #顺序不能错
                    #三角形是顺时针，是凹点
                    concavePoint.append(P2);
                elif (self.judgeClockwise(P1, P2, P3) > 0): #顺序不能错
                    convexPoint.append(P2);
            returnArray.append(newArray);
            returnArray.append(concavePoint);
            returnArray.append(convexPoint);
                      
            return returnArray;
        else:
            return [array, [], []];
			
def tmp3():
    path = [[[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [6, 3], [6, 3.86], [3.33, 5], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [3.33, 5], [6, 3.86], [8, 3], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.5, 1.5], [1.78, 2.33], [2, 3], [3.33, 5], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [1.78, 2.33], [2, 3], [3.33, 5], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [2, 3], [6, 3], [6, 3.86], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [2, 3], [3.33, 5], [6, 3.86], [6, 9]], [[6, 9], [2.4, 5.4], [3.33, 5], [2, 3], [6, 3], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [2.4, 5.4], [3.33, 5], [6, 3.86], [6, 3], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [2.4, 5.4], [-6, 9], [-0.37, 2.56], [0, 3], [1.33, 3], [2, 3], [3.33, 5], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [1.78, 2.33], [3.67, -0.5], [4.29, 0.43], [6, 3], [6, 3.86], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [6, 3], [6, -3], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [6, 3], [6, -3], [8, -7], [8, 3], [6, 3.86], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [6, 3], [4.5, 0], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [3.33, 5], [6, 3.86], [6, 3], [6, -3], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [3.33, 5], [6, 3.86], [6, 3], [4.5, 0], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [4.29, 0.43], [4.5, 0], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [6, -3], [5.56, -3.33], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [6, -3], [5.67, -3.5], [5.35, -3.97], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [4.5, 0], [6, -3], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [4.5, 0], [4, -1], [6, -3], [8, -7], [8, 3], [6, 9]]]
    len_path = len(path);

    centerArray = [];
    concaveArray = [];

    #计时开始
    startTime = time.clock();
    
    for i in range(len_path):
        poly = Polygon(path[i]);
        #poly.genVertexTree();
        centerArray.append(poly.centerPoint());

        a = poly.judgeConcavePolygon(poly.getVertexValue());
        concaveArray.append(a[1]);
        print('路径：', a[0]);
        print('凹点：', a[1]);
        print('凸点：', a[2]);
        print('\n');

    #计时结束
    endTime = time.clock();

    #打印结果
    print('操作用时：{0:.3e} s'.format(endTime-startTime));   

    print(centerArray);
    print(concaveArray);

#</span>

 
<span style="font-size:18px;">//
$center = [[3.48, 4.51], [3.73, 4.51], [2.23, 4.13], [2.21, 4.32], [3.07, 4.26], [2.73, 4.51], [5.22, 2.3], [5.72, 2.41], [1.09, 4.99], [4.13, 3.26], [5.62, 1.19], [5.67, 2.11], [5.48, 1.62], [6.12, 1.29], [5.98, 1.73], [6.1, 1.16], [6.4, 0.25], [6.38, 0.17], [6.27, 0.67], [6.06, 0.98]];

$concave = [[[3.33, 5]], [[3.33, 5], [6, 3.86]], [[0.8, 3.8], [2, 3], [3.33, 5]], [[0.8, 3.8], [1.33, 3], [2, 3], [3.33, 5]], [[0.8, 3.8]], [[0.8, 3.8], [3.33, 5]], [[3.33, 5], [6, 3]], [[6, 3.86]], [[2.4, 5.4], [0, 3], [3.33, 5]], [], [[6, 3], [6, -3]], [[6, 3], [6, 3.86]], [[6, 3]], [[6, 3.86], [6, -3]], [[6, 3.86], [6, 3]], [[6, 3]], [[6, -3]], [[6, -3], [5.67, -3.5]], [[6, 3], [6, -3]], [[6, 3], [6, -3]]];

$path2 = [[[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [6, 3], [6, 3.86], [3.33, 5], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [1.33, 3], [2, 3], [3.33, 5], [6, 3.86], [8, 3], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.5, 1.5], [1.78, 2.33], [2, 3], [3.33, 5], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [1.78, 2.33], [2, 3], [3.33, 5], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [2, 3], [6, 3], [6, 3.86], [6, 9]], [[6, 9], [2.4, 5.4], [0.8, 3.8], [0, 3], [1.33, 3], [2, 3], [3.33, 5], [6, 3.86], [6, 9]], [[6, 9], [2.4, 5.4], [3.33, 5], [2, 3], [6, 3], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [2.4, 5.4], [3.33, 5], [6, 3.86], [6, 3], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [2.4, 5.4], [-6, 9], [-0.37, 2.56], [0, 3], [1.33, 3], [2, 3], [3.33, 5], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [1.78, 2.33], [3.67, -0.5], [4.29, 0.43], [6, 3], [6, 3.86], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [6, 3], [6, -3], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [6, 3], [6, -3], [8, -7], [8, 3], [6, 3.86], [6, 9]], [[6, 9], [3.33, 5], [2, 3], [6, 3], [4.5, 0], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [3.33, 5], [6, 3.86], [6, 3], [6, -3], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [3.33, 5], [6, 3.86], [6, 3], [4.5, 0], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [4.29, 0.43], [4.5, 0], [6, -3], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [6, -3], [5.56, -3.33], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [6, -3], [5.67, -3.5], [5.35, -3.97], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [4.5, 0], [6, -3], [5.67, -3.5], [8, -7], [8, 3], [6, 9]], [[6, 9], [6, 3.86], [6, 3], [4.5, 0], [4, -1], [6, -3], [8, -7], [8, 3], [6, 9]]]

	if (1) {
		var r = 20;      
        config.setSector(1,1,1,1);        
        config.graphPaper2D(0, 0, r);      
        config.axis2D(0, 0, 250, 1.2);   
		
		//坐标轴设定        
        var scaleX = 2*r, scaleY = 2*r;          
        var spaceX = 2, spaceY = 2;           
        var xS = -10, xE = 10;          
        var yS = -10, yE = 10;          
        config.axisSpacing(xS, xE, spaceX, scaleX, 'X');            
        config.axisSpacing(yS, yE, spaceY, scaleY, 'Y');            
                    
        var transform = new Transform();    
		
		//顶点
		var a = [];
		
		for (var i = 0; i < $vertex.length; i++) {
			a.push($vertex[i][0]);
		}
		
		//显示变换        
        if (a.length > 0) {        
            a = transform.scale(transform.translate(a, 0, 0), scaleX/spaceX, scaleY/spaceY);
		}
		
		var lable = [];
		for (var i = 0; i < 100; i++) {
			lable.push(i.toFixed(0));
		}
		
		/*
		//边集
		var b = [];
		for (var i = 0; i < $seg.length; i++) {
			b.push([a[$seg[i][0]], a[$seg[i][1]]]);
		}
		
				
		var edges = b.length;
		for (var i = 0; i < edges; i++) {
			shape.multiLineDraw([].concat(b[i]), 'red');
		}*/
		
		var colorArray = ['red', 'orange', 'yellow', 'green', 'cyan', 'blue', 'purple'];
		var seg = [];		
		
		var idx = xGlobal;
		var len = $path2.length;

		plot.setLineWidth(3);

		seg = transform.scale(transform.translate($path2[idx%len], 0, 0), scaleX/spaceX, scaleY/spaceY);
		shape.multiLineDraw([].concat(seg), 'pink');
		
		
		//主要顶点
		shape.pointDraw([].concat(a), 'blue', 1, 1, lable);
		
		cPoint = transform.scale(transform.translate($center, 0, 0), scaleX/spaceX, scaleY/spaceY);
		
		shape.pointDraw([].concat([cPoint[idx%len]]), 'red');
		
		
		concavePoint = transform.scale(transform.translate($concave[idx%len], 0, 0), scaleX/spaceX, scaleY/spaceY);
		
		shape.pointDraw([].concat(concavePoint), 'orange');
		

		/*
		//次要顶点
		var b = [];
		
		for (var i = 0; i < 6; i++) {	
			b = [].concat($pointType[n][i]);
			//显示变换        
			if (b.length > 0) {        
				b = transform.scale(transform.translate(b, 0, 0), scaleX/spaceX, scaleY/spaceY);
				shape.pointDraw([].concat(b), colorArray[(i-1)%7]);
			}

			
			
		}*/
	}

//</span>

本节到此结束，欲知后事如何，请看下回分解。







mwsister
关注
专栏目录
仅仅三行JAVA代码计算多边形的几何中心点
洛阳泰山的博客
10-24 4750
前言：因为工作设计到gis相关的内容，需要计算采煤机工作面的中心点。如果套用数学的计算公式，用java去实现，太多麻烦还费时比较久，于是我找到java几何计算的工具包，几行代码就能求出多变形的中心，简直yyds!如图所示，顺序按照A->B>C>D,且最后闭合点为A点，即开始和最后输入点为A点。除了四边形形，五边形，六边形，不规则多边形，都可以计算。1.坐标点的顺序必须保证是几何图形的相邻点，顺时针相连或者逆时针都可以。2.最后要完成闭合，闭合的方式即最后一个点必须和第一个点相同。
[从头学数学] 第275节 [计算几何] 多边形顶点的五种类型
微笑的彩虹的博客
09-24 1148
剧情提要： 阿伟看到了一本比较有趣的书，是关于《计算几何》的，2008年由北清派出版。很好奇 它里面讲了些什么，就来看看啦。 正剧开始： 星历2016年09月24日 12:50:50, 银河系厄尔斯星球中华帝国江南行省。 [工程师阿伟]正在和[机器小伟]一起研究[计算几何]]。 # >>> [1.33, 3]: [Point([2, 3]), Po
多边形方向及顶点凹凸性的判断
10-13
自己用C#实现的关于多边形方向及顶点凹凸性的判断的方法。里边的函数包括了构造多边形、方向判断、凹凸性判断。适合初学者学习参考之用。
matlab 凸多边形,matlab练习程序（多边形顶点凹凸性）
weixin_39707612的博客
03-19 1489
生成简单多边形后，有时还需要对多边形各顶点的凹凸性做判断。先计算待处理点与相邻点的两个向量，再计算两向量的叉乘，根据求得结果的正负可以判断凹凸性。结果为负则为凹顶点，为正则为凸顶点。凹顶点用o表示，凸顶点用*表示。结果如下：matlab代码如下：clear all;close all;clc;n=;p=rand(n,);p=createSimplyPoly(p); %创建简单多边形hold on...
多边形重心和中心
潘达小新
10-02 1万+
目录 计算任意多边形的中心： 计算多边形重心： （1）划分多边形为三角形： （2）求每个三角形的重心和面积： （3）求原多边形的重心： 总结： 计算任意多边形的中心： 获取多边形外接矩形，获取最大最小XY（也就是外接矩形的四个角点的坐标：Xmax、Xmin、Ymax、Ymin） static double XMax, XMin, YMax, YMin; double CenterX = (XMax + XMin) / 2; dou.
数学科学的完整课程大纲（工科自学必看）
最新发布
weixin_46596227的博客
11-14 4834
为了能够系统性学习数学理论打好更扎实的基础，我特地网上搜索，整理了一份数学科学大纲。按照这份材料去学习，基本上能够满足理论上的需求。
浙江大学软件学院三维动画与交互技术考试概念整理
乞力马扎罗的雪CYF的博客
01-16 8348
第一讲1.    增强现实技术AR：--融合了三维动画、立体视觉和图像处理；--建模、渲染、位置标定、图像融合；2.    OpenGL是一种用于创建实时3D图像的编程接口。3.    三维这个术语表示一个正在描述或显示的物体具有三维维度：宽度、高度、深度；--计算机3D图形实质上也是平面的；--在计算机屏幕上显示的二维图像，提供深度（或第三维）的错觉；2D+透视 = 3D透视使人产生深度的错觉。
Stochastic Light Culling for VPLs on GGX Microsurfaces论文研读
麻豆_matou
04-23 1097
前言 好久没写博客了，今天来共享一下最近研读的一篇论文吧，文章中[***.pdf]为参考文献可自行谷歌学术下载。因为本人还是个图形学菜鸟，所以有什么问题希望大家多多指正。 Stochastic Light Culling for VPLs on GGX Microsurfaces 词汇解析 VPL VPL - virtual point light 虚拟点光源 以虚拟的点光源为间接光提供光照。可以...
"细节重建的多分辨率曲面混合算法在逆向工程中的应用
图形和视觉计算6（2022）200043技术部分用于细节重建的多分辨率曲面混合算法HonoSalval，Andy Keane，David Toal计算工程和设计组，南安普敦大学，劳斯莱斯大学计算工程技术中心，Boldrewood Innovation Campus，...
计算不规则多边形的面积、中心、重心
Keep Moving--Leo Li
04-04 1万+
转载地址：http://blog.csdn.net/shao941122/article/details/53671643 最近项目用到：在不规则多边形的中心点加一个图标。（e.g: xx地区发生暴雪，暴雪区域是多边形，给多边形中心加一个暴雪的图标） 之前的设计是，计算不规则多边形范围矩形bounds的中心点。这个比较简单，对于一些圆，矩形，凸多边形都比较适合。但是遇到凹多边形就会出现
[几何] 计算任意多边形的面积、中心、重心
热门推荐
MapleShao
12-15 3万+
最近项目用到：在不规则任意多边形的中心点加一个图标。（e.g: xx地区发生暴雪，暴雪区域是多边形，给多边形中心加一个暴雪的图标） 之前的设计是，计算不规则多边形范围矩形bounds的中心点。这个比较简单，对于一些圆，矩形，凸多边形都比较适合。但是遇到凹多边形就会出现问题，比如一个月牙型的不规则多边形，bounds的中心点，就落到月牙外了。就有点难以接受了。 经过讨论，决定将中心改为重心。 ...
判断多边形的凹凸性和计算多边形面积：利用向量叉乘
Mikchy的博客
08-08 8955
根据百度百科的讲解： 凸多边形 现在重点讲解顶点凹凸性法（最常用也是较为简单的方法）：计算总结在最后。 利用向量叉乘的相关知识进行计算：假设当前连续的三个顶点分别是P1，P2，P3。计算向量P1P3,P1P2的叉乘，也可以计算三角形P1P2P3的面积，得到的结果如果大于0，则表示P2点在线段P1和P3的右侧，此时P2对应的角度小于180。然后依次计算下一个前后所组成向量的叉乘，如...
【计算几何】点在多边形内部
weixin_34007879的博客
08-01 134
问题描述：已知点P(x,y)和多边形Poly，判断点P(x,y)是否在多边形内部。 基本方法：射线法 以点P为端点，向左方作射线L，由于多边形是有界的，所以射线L的左端一定在多边形外部，考虑沿着L从无究远处开始自左向右移动。 遇到和多边形的第一个交点的时候，进入到了多边形的内部，遇到第二个交点的时候，离开了多边形... 因而当L和多边形的交点数目C是奇数的时候，P在多边形内，是偶数，则P在...
凸包初学整合
K键盘里的青春K
04-13 1759
首先说明，此篇博客是我在初学凸包时将看到过的所有很好很经典博客里简单易懂或者精彩的部分摘取组合而成的。每一篇博客侧重点不同，有些讲的比较粗略但是另一些博客很详细，所以我觉得将这些精华聚合在一起还是很重要的，每篇博客我都会标明出处。另外有些地方我也会增加一些我自己的见解，如果有不对的地方，希望各位大佬指出。如果您觉得好的话欢迎转载，但请说明出处： http://blog.csdn.net/qq
计算不规则多边形的中心、重心 （js 实现）
往事随风的博客
03-10 1万+
在项目中经常会遇到这样一个问题：在不规则多边形的中心点加一个图标。 之前的设计是，计算不规则多边形范围矩形bounds的中心点。这个比较简单，对于一些圆，矩形，凸多边形都比较适合。但是遇到凹多边形就会出现问题，比如一个月牙型的不规则多边形，bounds的中心点，就落到月牙外了。就有点难以接受了。 这个问题之前困扰我很久，今天突然看到一个名词：”重心“ 问题瞬间迎刃而解！！特此记录一下 首...
线性代数-向量空间
levin blog
09-09 5033
线性代数学习笔记向量空间 向量空间 向量空间的严格定义： 设 VVV 为一向量组，如果 VVV 非空，且 VVV 对于向量的加法以及数乘两种运算封闭，那么就称 VVV 为 向量空间。 所谓封闭，是指在 VVV 中向量进行数乘和加减，其结果依然在 VVV 中。具体的说，就是： 若 a∈V,b∈V,a \in V, b \in V,a∈V,b∈V, 则 a+b∈Va + b \in Va+b∈V 若 a∈V,k∈R,a \in V, k \in \mathbb{R},a∈V,k∈R, 则 ka∈Vka \i
matlab练习程序（多边形顶点凹凸性）
weixin_33894992的博客
01-08 1204
生成简单多边形后，有时还需要对多边形各顶点的凹凸性做判断。 先计算待处理点与相邻点的两个向量，再计算两向量的叉乘，根据求得结果的正负可以判断凹凸性。 结果为负则为凹顶点，为正则为凸顶点。 凹顶点用o表示，凸顶点用*表示。 结果如下： matlab代码如下： clear all;close all;clc; n=20; p=rand(n,2); p=createSi...
[计算几何]一些模型的积累
srand(time(0))
02-16 735
计算几何
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值