idl中几何运算.docx

IDL中的几何图形数学运算几何形状，分为点、线、多边形（面）、体几类，利用IDL对这些形状的几何运算，大致分下面几个部分。点集运算点与点之间求距离：DISTANCE_MEASURE（IDL自带）线相关计算点到直线的距离PNT_LINE（IDL自带）或CalDistancePtoLine.pro计算两直线的交点CAL2LINESINTERSECTPOINT.pro线段与坐标轴的夹角（数学坐标系）cal2pointsangle.pro多边形相关（面）求多边形面积：poly_area、(IDLanROI)- ComputeGeometry,area = result（IDL自带）求多边形周长：(IDLanROI)- ComputeGeometry, PERIMETER = result（IDL自带）点坐标是否在多边形范围内：(IDLanROI)-ContainsPoints（IDL自带）三点求通过该系列三个点的圆心坐标和圆半径；CIR_3PNT（IDL自带）多边形与曲面交集：MESH_CLIP（IDL自带）两个多边形合并：MESH_MERGE（IDL自带）多边形是否空间闭合：MESH_ISSOLID（IDL自带）多边形包含的三角形个数：MESH_NUMTRIANGLES（IDL自带）复杂多边形正确显示：IDLgrTessellator（IDL自带）求两个平面的夹角cal2planeangle.pro四面体与平面相交：TETRA_CLIP（IDL自带）体相关体数据任意方向切面：EXTRACT_SLICE（IDL自带）重点推荐MESH_OBJ，可以创建三角形、矩形、圆锥、圆柱、任意多边形旋转等。;Mesh_Obj,2, Vertex_List, Polygon_List, $ ; 调用格式 ; Mesh_Obj, Type, Vertex_List, Polygon_List , $ ; Array1 [, Array2] [, /CLOSED] [, /DEGREES] , $ ; [P1, P2, P3, P4, P5 = value] ;Type = 6 创建旋转后模型的点坐标及连接关系 ;Vertex_List 创建后的点坐标集 ;Polygon_List 旋转模型的点连接关系 ;Array1 创建旋转模型的起始数据集 ;Array2 Type 为6时不支持 ;/Closed 旋转后模型闭合 ;/DEGREES 旋转的输入参数为角度，默认是弧度 ;P1 －旋转后的边数 ;P2 - 旋转的中心点 ;P3 - 旋转的方向[0,0,1] ;P4 - 旋转起始角度;P5 - 旋转结束角度附源码：附源码：;≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌;;计算两个点的距离;Write By DYQ;可用distance_measure函数替代;function CalDistance, point1, point2??? compile_opt idl2??? ;??? Return, SQRT((point1[0]-point2[0])^2+(point1[1]-point2[1])^2+(point1[2]-point2[2])^2)end;≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌; Write By ZML;计算点到直线的距离;s = SQRT(p*(p-a)*(p-b)*(p-c)){p = (a+b+c)/2};Modified By DYQ;function CalDistancePtoLine, point0,linePos1,linePos2??? a = CalDistance(point0,linePos1)??? b = CalDistance(point0,linePos2);SQRT((point[0]-linePos2[0])^2+(point[1]-linePos2[1])^2)??? c = CalDistance(linePos1,linePos2);SQRT((linePos2[0]-linePos1[0])^2+(linePos2[1]-linePos1[1])^2)??? p = (a+b+c)*0.5??? s = SQRT(p*(p-a)*(p-b)*(p-c))??? Return, s*2/cend;≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌≌;;求集合的并集; copy From DAVIDs Code? ^_^;FUNCTION SETUNION, a, b? ;? COMPILE_OPT StrictArr?