计算机制图课程实践探究.docxVIP

兰州交通大学《计算机地图制图实践报告》学生姓名殷成钰 学 号 201408312 班 级 地信1401班专 业 地理信息科学学 院 测 地 学 院任课教师刘涛 二O一六年十二月目录1.系统描述12.页面布局13.功能实现13.1点与三角形位置关系判断13.1.1算法原理13.1.2关键代码23.1.3运行结果33.2点与多边形位置关系判断43.2.1算法原理43.2.2关键代码63.2.3运行结果93.3格雷厄姆算法93.3.1算法原理93.3.2关键代码103.3.3运行结果123.4 Douglas-Peuker算法133.4.1算法原理133.4.2关键代码143.4.3运行结果153.5栅格数据处理实验163.5.1地形指标提取163.5.2地形特征信息提取203.6其他算法243.6.1算法原理243.6.2关键代码253.6.3运行结果254.实践总结261.系统描述画图系统设计的画图板主要包括工具栏，画图区，位置关系判断，要素综合和画笔设置四个部分。工具菜单包含清空画板、橡皮擦，求点到线段距离等工具菜单。工具栏可根实际需要选择不同的工具完成画图操作；画图区能够对进行绘制图画和做分析处理。2.页面布局3.功能实现本系统主要实现了，点与线、线与线、点与三角形、点与多边形的位置关系判断；格雷厄姆算法求解平面点集的凸壳，四比特串编码法矩形开窗以及道格拉斯扑克算法综合线要素。另外提供画笔的预设、橡皮擦、清空画板、结构体绘制、求点到线段距离等工具。3.1点与三角形位置关系判断3.1.1算法原理如右图，?ABC及点P（xp，yp），设直线AB对应的方程为：f1=a1x+b1y+c1;直线BC对应的方程为：f2=a2x+b2x+c2;直线CA对应的方程为：f3=a3x+b3x+c3;对于点P（xp，yp），若满足则点P位于三角形的内部，否则点P位于三角形外部。3.1.2关键代码Point p1 = new Point((int)PiX[0], (int)PiY[0]);Point p2 = new Point((int)PiX[1], (int)PiY[1]);Point p3 = new Point((int)PiX[2], (int)PiY[2]);if (LineEquation(p1, p2, pArray[0]) == LineEquation(p1, p2, p3)){if (LineEquation(p1, p3, pArray[0]) == LineEquation(p1, p3, p2)) {if (LineEquation(p2, p3, pArray[0]) == LineEquation(p2, p3, p1)) {MessageBox.Show(点在三角形内！); }else {MessageBox.Show(点在三角形外！); }}else {MessageBox.Show(点在三角形外！); }}else{MessageBox.Show(点在三角形外！);}其中LineEquation(Point p1, Point p2, Point p)方法是判断p点在直线p1p2的一侧代码如下：if ((double)(p.X - p1.X) / (double)(p2.X - p1.X) == (double)(p.Y - p1.Y) / (double)(p2.Y - p1.Y)) {return 0; }else if ((double)(p.X - p1.X) / (double)(p2.X - p1.X) (double)(p.Y - p1.Y) / (double)(p2.Y - p1.Y)) {return -1; }else {return 1; }3.1.3运行结果3.2点与多边形位置关系判断3.2.1算法原理1)夹角和算法：如果多边形各边不相交，且连续，无岛屿，可以用此算法判断简单多边形与点的包含关系。现设有一简单n边形，其顶点可以表示P(x,y)，i=1,2,…,n，另有带判别的点A(x,y)。将A点与各个顶点相连线，计算其夹角和，且规定顺时针方向为正，逆时针为负如图：若有：2)铅垂线算法：从待判别点引铅垂线，由该铅垂线与多边形交点个数的奇偶性来判断点是否在多边形内。若交点个数为奇数，则点在多边形内；若交点个数为偶数，则该点在多边形外。第一步，计算多边形是最小投影矩形，若点在最小投影外，则，点一定在多边形外，算法结束；否则执行第二步。第二步，设置记录交点个数的计数器N