基于TIN设计的地形图精度检测软件

２０ｌ１年６月　第３期　文章编号：１６７２—８２６２（２０１０）０３—１３４—０４　城市勘测　Ｊｕｎ．２０１１　Ｎｏ．３　Ｕｒｂａｎ　Ｇｅｏｔｅｃｈｎｉｃａｌ　Ｉｎｖｅｓｔｉｇａｔｉｏｎ＆Ｓｕｒｖｅｙｉｎｇ　中图分类号：Ｐ２０９　文献标识码：Ｂ　基于ＴＩＮ设计的地形图精度检测软件　吴美容　（江苏兰德数码科技有限公司，江苏南京２１０００７）　摘　要：目前国内行业商用软件如ＣＡＳＳ、广州开思等在地形图高程精度评定方面相对较弱。以往都是采用手工逐点　量距内插计算检测点位的高程并进行高程中误差统计，若碰到几百个检测点，手５－－计算工作量非常大而且容易出错，　因此根据工作需要研制开发基于ＡｕｔｏＣＡＤ平台的地形图精度检测软件就变得很有意义，借助程序计算代替传统的手　工计算使得工作效率大为提高。实践证明，软件设计思路明晰，算法可靠。　关键词：地形图；精度；高程；中误差；软件　１　引　言　工程建设前期均须测绘大比例尺原始地形图，以　作为设计人员规划设计之用。目前，｛９＿１０绘单位常采用　全站仪、ＲＴＫ或近景摄影测量进行外业碎部点数据采　集，然后运用行业成图软件生成地形图。地形图精度　为两大类¨　；即静态三角网和动态三角网。静态算法　中以三角网生长算法较为典型，动态算法中则以数据　逐点插入算法较为普遍。软件基于文献　介绍的构　网思路采用三角网生长算法中的递归生长算法。该算　法首先生成第一个三角形，然后由第一个三角形向外　扩展直到最后联结全部离散点构成Ｄｅｌａｕｎａｙ三角网。　该算法涉及的相关判别法则限于篇幅，在此不作详细　阐述，直接给出构网结果，如图１所示。　是否达到规范要求，需要进行精度检测。检测分为平　面精度检测和高程精度检测，二者可以分开进行也可　以同时进行。ＣＡＳＳ　７．１检查入库菜单里面点位误差　检查和边长误差检查即是针对平面精度检测。实践中　一般采取以下三种方法针对高程精度进行检测：①点　位还原法，即从所有碎部点数据中随机抽取总点数的　１０％一２０％个数据，全站仪利用点的平面坐标采用极　坐标法于实地放样出该点，然后测量点位二次高程并　且与点位的一次高程进行比较，最后根据所抽取的子　样统计得出的中误差来估计总体的精度；②剖面法，即　检测人员于实地进行断面测量，内业根据断面数据文　图１　Ｄｅｌａｕｎａｙ三角网　件生成断面图并与根据地形图所剖得到的断面图进行　比较；③检查点法，即检测人员在测区内采集任意点位　三维数据，点位要求覆盖整个测区要有一定的密度且　分布均匀，然后展点到所生成的地形图中去，采用内插　２．２展点　该部分将检测数据（　．ｄａｔ格式）展点到图形巾　去。相关代码：　Ｃａｌｌ　Ｇｅｔ０ｐｅｎＦｉｌｅＮａｍｅ（Ｄａｔａ，”请输入数据文件”）　Ｏｐｅｎ　Ｄａｔａ　Ｆｏｒ　Ｉｎｐｕｔ　Ａｓ＃３　Ｄｏ　的方法计算点位的图面高程并与点位的抽测高程进行　比较，最后计算中误差。该方法是实践中比较常用的　方法。本文正是基于第３种方法编写地形图精度检测　Ｉｎｐｕｔ＃３，ｎｏ，ＢｉａｎＭａ，ｊｕｚｙ，ｊｕｚｘ，ｊｕｚｈ　Ｓｅｔ　ｂｌｏｃｋＲｅｆＯｂｊ＝ＴｈｉｓＤｒａｗｉｎｇ．ＭｏｄｅｌＳｐａｃｅ．１ｎｓｅｒｔＢｌｏｅｋ　ｆ　ｉｎ—　软件。　ｓＰｏｉｎｔ，”Ｂｌｐｏｉｎｔ”，ｌ＃，１＃，１＃，０）　Ｌｏｏｐ　Ｕｎｔｉｌ　ＥＯＦ（３）Ｏｒ　ｎｏ＝Ｅｍｐｔｙ　２编程思路　２．１构ＴＩＮ　２．３检测　该部分根据上面所生成的三角网及展点进行点位　原始数据为　．ｄａｔ格式文件。三角网构网算法归　高程中误差统计。　算法的基本流程：　收稿日期：２Ｏ１Ｏ一１０—２９　作者简介：吴美容（１９８２～），女，工程师，主要从事航测遥感内业技术　作。　第３期　吴美容．基于ＴＩＮ设计的地形图精度检测软件　１３５　程序首先建立两个选择集合，分别用于存放检测　点及三角形元素，其次建立双循环结构对于检测点集　Ｌ　＝０或Ｌ　＝０或Ｌ，＝０时点Ｐ在△　边上；　的某一　合中的每一个点元素遍历三角网集合中的每一个三角　形元素，并进行点与三角形位置关系的判断，若位于某　三角形则调用该三角形的顶点数据，采用线性插值的　内插方法进行相应的内插计算，将结果放入指定的数　组，最终程序建立与Ｅｘｃｅｌ的通讯将结果输出到Ｅｘｃｅｌ　指定区域。　算法的数据结构：　Ｐｕｂｌｉｃ　Ｔｙｐｅ　ＤｉｎｇＤｉａｎ　Ｘ　Ａｓ　Ｄｏｕｂｌｅ　Ｙ　Ａｓ　Ｄｏｕｂｌｅ　Ｈ　Ａｓ　Ｄｏｕｂｌｅ　Ｅｎｄ　Ｔｙｐｅ　Ｐｕｂｌｉｃ　Ｔｙｐｅ　ＳａｎＪｉａｏＷａｎｇ　Ａ（０　Ｔｏ　２）Ａｓ　Ｄｏｕｂｌｅ　Ｂ（０　Ｔｏ　２）Ａｓ　Ｄｏｕｂｌｅ　Ｃ（０　Ｔｏ　２）Ａｓ　Ｄｏｕｂｌｅ　Ｅｎｄ　Ｔｙｐｅ　Ｐｕｂｌｉｃ　Ｔｙｐｅ　ＸｉａｎＤｕａｎ　ｄｂｌＳｔａｒｔ（０　Ｔｏ　３）Ａｓ　Ｄｏｕｂｌｅ　ｄｂｌＥｎｄ（０　Ｔｏ　３）Ａｓ　Ｄｏｕｂｌｅ　Ｅｎｄ　Ｔｙｐｅ　Ｐｕｂｌｉｃ　Ｔｙｐｅ　ｓｉｄｅ——２——Ｐｏｉｎｔ　ＩＤ　Ａｓ　Ｌｏｎｇ　Ａ（）Ａｓ　Ｌｏｎｇ　ｃｏｕｎｔ　Ａｓ　Ｌｏｎｇ　ＥｎｄＴｙｐｅ　ＤｉａｎＴｒｉａｎｇｌｅ　类模块　Ｐｕｂｌｉｃ　ｙｉ　Ａｓ　Ｌｏｎｇ　Ｐｕｂｌｉｃ　ｅｒ　Ａｓ　Ｌｏｎｇ　Ｐｕｂｌｉｃ　ｓａｎ　Ａｓ　Ｌｏｎｇ　（１）判断检测点与三角网中三角形的关系　］　Ｙ　０　Ｘ　图２点位关系判断　思路：如图２，ＡＶ１　的三个顶点的坐标分别为　Ｖ　（　，Ｙ　），　（　，Ｙ　），　（　，Ｙ　），任一点Ｐ（　，Ｙ）。则　点Ｐ和△　的位置关系如下判定，　、　、　定义　如下：　Ｌ　＝［（　一　：）（ｙ－ｙ：）一（Ｙ。一Ｙ：）（　—　：）］　Ｌ：＝［（　一　３）（ｙ－ｙ３）一（ｙ１－ｙ３）（Ｘ－－Ｘ３）］　Ｌ，＝［（　：一　。）（ｙ－ｙ　）一（Ｙ２－Ｙ　）（Ｘ－－Ｘ　）］　判断条件：　Ｌｌ＞Ｏ且，Ｊ２＞０且　＞０或Ｌ１＜０且　２＜０且Ｌ３＜Ｏ　时点Ｐ在△　内；　除去以上两种情况时点Ｐ在△　外。　代码如下：　Ｐｕｂｌｉｃ　Ｆｕｎｃｔｉｏｎ　Ｆｎ—ｉｎ（Ｄｉａｎｏ（）Ａｓ　Ｄｏｕｂｌｅ，ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ　Ａｓ　ＳａｎＪｉａｏＷａｎｇ）Ａｓ　Ｂｏｏｌｅａｎ　Ｄｉｍ　Ｌ１　Ａｓ　Ｄｏｕｂｌｅ，１２　Ａｓ　Ｄｏｕｂｌｅ，Ｌ３　Ａｓ　Ｄｏｕｂｌｅ　Ｆｕｉｎ＝Ｆａｌｓｅ　＿Ｌ１＝（ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ｃ（０）一ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ｂ（０））￥　（Ｄｉａｎｏ（１）－ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ｂ（１））一（Ｄｉａｎｏ（０）－Ｓａ￣ｉａｏＸｉｎｇ—　ａｂｃ．Ｂ（０））　（ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ｃ（１）－ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ｂ（１））　２１＝（ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ａ（０）一ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ｃ（０））　（Ｄｉａｎｏ（１）－ＳａｎＪｉａｏＸｉｎｇ＿ａｂｃ．Ｃ（１））一（Ｄｉａｎｏ（０）－Ｓａｎ．］ｉａｏＸｉｎｇ—　ａｂｃ．Ｃ（０））　（ＳａｎＪｉａｏＸｉｎｇ＿ａｂｃ．Ａ（１）－ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ｃ（１））　Ｉ２＝（ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ｂ（０）一ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ａ（０））　（Ｄｉａｎｏ（１）－ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ａ（１））一（Ｄｉａｎｏ（０）－ＳａｎＪｉａｏＸｉｎｇ—　ａｂｃ．Ａ（０））　（ＳａｎＪｉａｏＸｉｎｇ＿ａｂｃ．Ｂ（１）－ＳａｎＪｉａｏＸｉｎｇ—ａｂｃ．Ａ（１））　ＩｆＡｂｓ（ｕ）＜０．００００１　Ｏｒ　Ａｂｓ（１２）＜０．００００１　Ｏｒ　Ａｂｓ（ＩＪ３）＜　０．００００１　Ｔｈｅｎ　ＩｆＡｂｓ（Ｌ１）＜０．００００１　Ｔｈｅｎ　Ｉｆ（Ｓｇｎ（１２）＝Ｓｇｎ（Ｌ３））０ｒ（Ａｂｓ（Ｓｇｎ（１２）　Ｓｇｎ　（Ｉ３））＜０．００００００１）Ｔｈｅｎ　Ｆｕｉｎ＝Ｔｒｕｅ　Ｅｘｉｔ　Ｆｕｎｃｔｉｏｎ　Ｅｎｄ　Ｉｆ　Ｅｎｄ　Ｉｆ　Ｉｆ　Ａｂｓ（１２）＜０．００００１　Ｔｈｅｎ　Ⅱ（Ｓｇｎ（ＬＪ３）＝Ｓｇｎ（Ｌ１））０ｒ（Ａｂｓ（Ｓｇｎ（Ｉ３）　Ｓｇｎ（Ｌ１））＜０．００００００１）Ｔｈｅｎ　Ｆｕｉｎ＝Ｔｒｕｅ　Ｅｘｉｔ　Ｆｕｎｃｔｉｏｎ　Ｅｎｄ　Ｉｆ　Ｅｎｄ　Ｉｆ　ＩｆＡｂｓｆＬ３）＜０．００００１　Ｔｈｅｎ　Ｉｆ（Ｓｇｎ（Ｌ１）＝Ｓｇｎ（１２））０ｒ（Ａｂｓ（Ｓｇｎ（Ｌ１）￥Ｓｇｎ　（１２））＜０．００００００１）Ｔｈｅｎ　Ｆｕｉｎ＝Ｔｒｕｅ　—Ｅｘｉｔ　Ｆｕｎｃｔｉｏｎ　Ｅｎｄ　Ｉｆ　Ｅｎｄ　Ｉｆ　Ｅｎｄ　Ｉｆ　Ｉｆ　Ｓｇｎ（Ｌ１）＝Ｓｇｎ（１２）Ａｎｄ　Ｓｇｎ（ＬＩ）＝Ｓｇｎ（Ｉ２）Ｔｈｅｎ　当在三条线段上的值都有一样的符号时，点落入三角形中　Ｆｕ—ｉｎ＝Ｔｒｕｅ　Ｅｌｓｅ　Ｆｕ—ｉｎ＝Ｆａｌｓｅ　Ｅｎｄ　Ｉｆ　Ｅｎｄ　Ｆｕｎｃｔｉｏｎ　（２）计算检测点在三角形中的高程　思路：　使用空间三角形３个顶点数据确定一个平面，继　而计算出内插点的高程。　１３６　城市勘测　２０１１午６　采用平面线性内插公式：ａｘ＋ｂｙ＋ｃｚ＋ｄ＝０，其中：　ａ＝（ｙｂ－ｙａ）ｘ（ｚｃ－ｚａ）一（ｙｃ—ｙａ）×（ｚｂ—ｚａ）　ｂ＝（ｘｂ—ｘａ）×（ｚｃ—ｚａ）一（ＸＣ—ｘａ）×（ｚｂ—ｚａ）×（一１）　ＳｉＧｅＭａＸ　ＢｉａｏＧｅ．Ｃｅｌｌｓ（ｊＪｊ，６）：（（ＳｉＧｅＭａＤｅｈａＬ—　２／　ｎ）／（ｎ一１））　０．５　样本标准差计算公式　ＢｉａｏＧｅ．Ｃｅｌｌｓ（ｊｊｊ＋１，６）＝（（ＳｉＧｅＭａＤｅｈａＬ）　ｃ＝（ｘｂ—ｘａ）×（ｙｃ—ｙａ）一（ｘｏ－ｘａ）×（ｙｂ－ｙａ）　ｄ＝（ａ×ｘａ＋ｂ×ｙａ＋ｃ×ｚａ）×（一１）　真误差计算中误差计算公式¨　Ｉｔ　ｎ＜＞０　Ｔｈｅｎ　代码如下：　Ｐｕｂｌｉｃ　Ｓｕｂ　ＴｏｕＹｉｎｇＨ（ＪｉａｏＪｉＤｉａｎ　Ａｓ　ＤｉｎｇＤｉａｎ．ＳａｎＪＸ　Ａｓ　ＯｕＭｉＧａ＝（（ＳｉＧｅＭａＤｅｌｔａＬ—ＳｉＧｅＭａＸ　２／ｎ）／（ｎ一　１））　０．５　Ｖａｒｉａｎｔ）　ＭｓｇＢｏｘ”需要检查的点数是：”＆（ＩｎｇＴｏｔａｌＰｏｉｎｔＮｕｍ　Ｄｉｍ　Ａ　Ａｓ　Ｄｏｕｂｌｅ，Ｂ　Ａｓ　Ｄｏｕｂｌｅ，Ｃ　Ａｓ　Ｄｏｕｂｌｅ，Ｄ　Ａｓ　Ｄｏｕｂｌｅ　Ａ＝［ＳａｎＪＸ（４）一ＳａｎＪＸ（１）］　［ＳａｎＪＸ（８）一ＳａｎＪＸ（２）］一　【ＳａｎＪＸ（７）一ＳａｎＪＸ（１）］　［ＳａｎＪＸ（５）一ＳａｎＪＸ（２）］　Ｂ＝－Ｉ　［ＳａｎＪＸ（３）－ＳａｎＪＸ（０）］　［ＳａｎＪＸ（８）一ＳａｎＪＸ（２）］　～［ＳａｎＪＸ（６）一ＳａｎＪＸ（０）］　［ＳａｎＪＸ（５）一ＳａｎＪＸ（２）］　Ｃ＝［ＳａｎＪＸ（３）一ＳａｎＪＸ（０）］　［ＳａｎＪＸ（７）一ＳａｎＪＸ（１）］一　［ＳａｎＪＸ（６）一ＳａｎＪＸ（０）］　［ＳａｎＪＸ（４）一ＳａｎＪ￣（１）］　Ｄ＝－１　Ａ　ＳａｎＪＸ（０）＋Ｂ　ＳａｎＪＸ（１）＋ｃ￥ＳａｎＪＸ（２）］　Ｉｆ　Ｃ＝０　Ｔｈｅｎ　ＪｉａｏＪｉＤｉａｎ．Ｈ＝０　Ｅｌｓｅ　ＪｉａｏＪｉＤｉａｎ．Ｈ＝一（Ｄ＋Ａ　ＪｉａｏＪｉＤｉａｎ．Ｘ＋Ｂ　ＪｉａｏＪｉＤｉ—　ａｎ．Ｙ）／Ｃ　Ｅｎｄ　Ｓｕｂ　（３）循环计算并输出　＇　］　Ｆｏｒ　ｉ＝０　Ｔｏ　ｌｎｇＴｏｔａ１ＰｏｉｎｔＮｕｍ　该部分设计每一个抽测点遍历三角网集合中的三角形元　素，若位于某三角形则进行相应的内插计算。　Ｎｅｘｔ　ｉ　Ｄｉｍ　ＡｐｐＥｘｃｅｌ　Ａｓ　Ｏｂｊｅｃｔ　Ｄｉｍ　ＮｅｗＷｏｒｋｂｏｏｋ　Ａｓ　Ｏｂｊｅｃｔ　Ｄｉｍ　ＢｉａｏＧｅ　Ａｓ　Ｏｂｉｅｅｔ　Ｄｉｍ　ＭａｔＲａｇｅ　Ａｓ　Ｏｂｉｅｃｔ　Ｄｉｍ　ＲｅｔｕｒｎＲａｎｇｅ　Ａｓ　Ｏｂｉｅｅｔ　Ｓｅｔ　ＡｐｐＥｘｃｅｌ＝ＣｒｅａｔｅＯｂｊｅｃｔ（”Ｅｘｃｅ１．Ａｐｐｌｉｃａｔｉｏｎ”）　ＡｐｐＥｘｃｅ１．Ｖｉｓｉｂｌｅ＝Ｔｒｕｅ　Ｓｅｔ　ＮｅｗＷｏｒｋｂｏｏｋ＝ＡｐｐＥｘｃｅ１．Ｗｏｒｋｂｏｏｋｓ．Ａｄｄ　Ｓｅｔ　ＢｉａｏＧｅ＝ＮｅｗＷｏｒｋｂｏｏｋ．ＡｃｔｉｖｅＳｈｅｅｔ　ＢｉａｏＧｅ．Ｎａｍｅ＝”地形图高程中误差统计表”　ＢｉａｏＧｅ．Ｃｅｌｌｓ（１，１）＝”地形图高程中误差统计表　ＢｉａｏＧｅ．Ｃｅｌｌｓ（２，１）＝”点号”　ＢｉａｏＧｅ．Ｃｅｌｌｓ（２，２）＝”Ｘ坐标”　ＢｉａｏＧｅ．Ｃｅｉｌｓ（２，３）：”Ｙ坐标”　ＢｉａｏＧｅ．Ｃｅｌｌｓ（２，４）＝”抽测高程”　ＢｉａｏＧｅ．Ｃｅｌｌｓ（２，５）＝”内插高程”　ＢｉａｏＧｅ．Ｃｅｌｌｓ（２，６）＝”误差”　Ｄｉｍ　ｊＪｊ　Ａｓ　Ｉｎｔｅｇｅｒ：ｊＪ　３　ＢｉａｏＧｅ．Ｃｅｌｌｓ（ＪＪＪ，１）＝ｎ　ＢｉａｏＧｅ．Ｃｅｌｌｓ（　，２）＝Ｒｏｕｎｄ（ＡｌｌＰｏｉｎｔ（ｉ，０），４）　ＢｉａｏＧｅ．Ｃｅｌｌｓ（ｊＪｊ，３）＝Ｒｏｕｎｄ（ＡｌｌＰｏｉｎｔ（ｉ，１），４）　ＢｉａｏＧｅ．Ｃｅｌｌｓ（俩，４）＝Ｒｏｕｎｄ（ＡｌｌＰｏｉｎｔ（ｉ，２），４）　ＢｉａｏＧｅ．Ｃｅｌｌｓ（ｊＪｊ，５）＝Ｒｏｕｎｄ（ＪｉａｏＤｉａｎ．Ｈ，４）　ＢｉａｏＧｅ．Ｃｅｌｌｓ（ＪＪＪ，６）＝Ｒｏｕｎｄ（ｄｂｌＤｉｖＥ１，４）　俩＝Ｊ‘Ｊ。ｊ＋１　＋１）＆ｖｂＣｒＬｆ＆一　”内插出高程的点数是：”＆ｎ＆ｖｂＣｒＬｆ＆一　”高程中误差是：　”＆Ｒｏｕｎｄ（ＯｕＭｉＧａ，４）　Ｆｎｄ　Ｔｆ　３　工程实例　以雅砻江某水电站坝区供水泵站原始地形检测为　例。　地形图高程中误差统计表　表１　ｌ　粼　辩褪　¨§§　ｉ｝｝｝　§ｉｌｌ　第３期　吴美容．基于ＴＩＮ设计的地形图精度检测软件　１３７　高程误差折线圈　［检查（＆ｊ）］　Ｃ＾Ｃ一一ｖｂａｒｕｎ”ｄ：／ｒｊｂ／ｄｇｘ／ｗａｎｇｊｇ．ｄｖｂ！ｄｘｊｃ　ｗａｎｇｊｇ．ＧＣＺａｎＧａｏＣｅｎｇ”　（２）测绘人员将测绘原理与计算机语言结合起来进　行二次开发能够解决实际中出现的新问题，也提高了工　作的效率。该软件经过多个工程地形图检测实践，证明　其正确、实用，提高了内业成图检测的作业效率。　点号　（３）程序设计每个点需要遍历三角网并进行判　断、内插计算等操作。因此，点数越多，三角网越复杂，　图３高程误差折线图　程序进行内插计算所花费的时间也越长。本程序的数　结论：　据结构优化将是接下来需要研究的课题。　参考文献　（１）第２７号点的高程偏差较大，分析可能系野外　测量过程中棱镜高输入错误，棱镜高值输小了，需予以　检查。　［１］　李志林，朱庆．数字高程模型（第二版）［Ｍ］．武汉：武汉　大学出版社，２００３　（２）高程中误差０．６　ｍ满足行业规范中关于山　地、高山地高程中误差不大于１　ｍ的要求；但是根据高　程误差折线图统计有２３个点出现负值，占到总点数的　７６．７％，分析认为该原始地形有抬高的趋势。对于工　程建设而言，大面积抬高原始地形意味着增加土石方　［２］　刘友光，黄桂兰，黄全义等．工程中数字地面模型的建立　与应用及大比例尺数字测图［Ｍ］．武汉：武汉测绘科技　大学出版社．１９９７　［３］　吴力新，史文中．地理信息系统原理与算法［Ｍ］．北京：　科学出版社，２００３　开挖量，无形中增加业主的投资成本。因此，即使高程　中误差满足要求，也必须对抬高原始地形的趋势予以　关注控制。　［４］　龚健雅．地理信息系统基础［Ｍ］．北京：科学出版社，　２００１　［５］　张宏，温永宁，刘爱利等．地理信息系统算法基础［Ｍ］．　北京：科学出版社，２００６　４　总　结　（１）程序运用ＶＢＡ语言编写，在Ｗｉｎ　ＸＰ系统Ａｔ—　［６］　姜友谊，黎晓．数字地面模型内插方法的优劣分析［Ｊ］．　西安科技学院学报，２００１，２１（３）　［７］　佟士懋，邢芳芳，夏齐霄．ＡｕｔｏＣＡＤ　ＡｃｔｉｖｅＸ／ＶＢＡ二次开　发技术基础及应用实例［Ｍ］．北京：国防工业出版社，　［８］　佟彪．ＶＢ语言与测量程序设计［Ｍ］．北京：中国电力出　版社，２００７　ｕｏＣＡＤ　２００４平台下运行成功。具体运用时可将宏代　码添加到自定义菜单下面运用十分方便，具体代码为：　［展点（＆ｚ）］　Ｃ　Ｃ一一ｖｂａｒｕｎ”ｄ：／ｒｊｂ／ｄｇｘ／ｗａｎｇｊｇ．ｄｖｂ！　ｚｄｇｗ．ＧＣＺａｎＧａｏＣｅｎｇ”　［构网（＆ｇ）］　Ｃ　Ｃ一一ｖｂａｒｕｎ”ｄ：／ｒｊｂ／ｄｇｘ／ｗａｎｇｊｇ．ｄｖｂ！　ｚｄｇｗ．ＥｘｔｅｎｄＴｒｉａｎｇｌｅＢｙＦｉｌｅ”　［９］　陈佳丽．Ｖｉｓｕａｌ　Ｂａｓｉｃ程序设计基础与实训教程［Ｍ］．北　京：清华大学出版社，２００５　Ｄｅｓｉｇｎ　Ｐｒｅｃｉｓｉｏｎ　Ｅｘａｍｉｎａｔｉｏｎ　Ｓｏｆｔｗａｒｅ　ｏｆ　Ｔｏｐｏｇｒａｐｈｉｃ　Ｍａｐ　Ｂａｓｅｄ　ｏｎ　ＴＩＮ　Ｗｕ　ＭｅｉＲｏｎｇ　（Ｊｉａｎｇｓｕ　Ｄｉｇｉｔａｌａｎｄ　Ｔｅｃｈｎｏｌｏｇｙ　Ｃｏ．，Ｌｔｄ．Ｎａｎｊｉｎｇ　２１０００７，Ｃｈｉｎａ）　Ａｂｓｔｒａｃｔ：Ｃｏｍｍｅｒｃｉａｌ　ｓｏｆｔｗａｒｅ，ｓｕｃｈ　ａｓ　ＣＡＳＳ　ａｎｄ　Ｇｕａｎｇｚｈｏｕ　ＫＡＩＳ１　ｗｈｉｃｈ　ａｒｅ　ｗｉｄｅｌｙ　ｕｓｅｄ　ｉｎ　ｄｏｍｅｓｔｉｃ　Ｓ＆Ｍ　ｉｎｄｕｓ—　ｔｒｙ，ａｒｅ　ｒｅｌａｔｉｖｅｌｙ　ｗｅａｋ　ｉｎ　ｔｅｒｍｓ　ｏｆ　Ｅｌｅｖａｔｉｏｎ　ａｃｃｕｒａｃｙ　ｅｖａｌｕａｔｉｏｎ．Ｉｎ　ｔｈｅ　ｐａｓｔ，ｔｈｅ　ｅｌｅｖａｔｉｏｎ　ｉｓ　ｃａｌｃｕｌａｔｅｄ　ｂｙ　ｍａｎｕａｌ　ｉｎｔｅｒ—　ｐｏｌａｔｉｏｎ　ｔｈｅｎ　ｃａｒｒｉｅｓ　ｏｕｔ　ｓｔａｔｉｓｔｉｃｓ　ｗｏｒｋ　ｏｆ　ｔｈｅ　ｅｌｅｖａｔｉｏｎ　ＲＭＥ．Ｉｆ　ｅｎｃｏｕｎｔｅｒｉｎｇ　ｈｕｎｄｒｅｄｓ　ｏｆ　ｃｈｅｃｋｐｏｉｎｔｓ　ｏｒ　ｅｖｅｎ　ｍｏｒｅ，ｔｈｅ　ｗｏｒｋｌｏａｄ　ｉｓ　ｅｎｏｒｍｏｕｓ　ａｎｄ　ｍａｎｕａｌ　ｃａｌｃｕｌａｔｉｏｎ　ｅｒｒｏｒ　ｉｓ　ｐｒｏｎｅ．Ｓｏ　ｉｔ　ｉｓ　ｖｅｒｙ　ｎｅｃｅｓｓａｒｙ　ａｎｄ　ｍｅａｎｉｎｇｆｕｌ　ｔｏ　ｒｅｓｅａｒｃｈ　ａｎｄ　ｄｅ—　ｖｅｌｏｐ，ｂａｓｅｄ　ｏｎ　ＡｕｔｏＣＡＤ　ｐｌａｔｆｏｒｍ，ｔｏｐｏｇｒａｐｈｉｃ　ｍａｐ　ａｃｃｕｒａｃｙ　ｔｅｓｔｉｎｇ　ｓｏｆｔｗａｒｅ．Ｉｎｓｔｅａｄ　ｏｆ　ｕｓｉｎｇ　ｔｒａｄｉｔｉｏｎａｌ　ｍａｎｕａｌ　ｅａｌｃｕ—　ｌａｔｉｏｎ，ｉｔ　ｍａｋｅｓ　ｅｆｆｉｃｉｅｎｃｙ　ｏｆ　ｃａｌｃｕｌａｔｉｏｎ　ｇｒｅａｔｌｙ　ｉｍｐｒｏｖｅｄ．Ｐｒａｃｔｉｃｅ　ｓｈｏｗｓ　ｔｈａｔ　ｔｈｅ　ｄｅｓｉｇｎ　ｉｄｅａｓ　ｉｓ　ｃｌｅａｒ，ａｌｇｏｒｉｔｈｍ　ｉｓ　ｒｅｌｉａ—　ｂｌｅ．　Ｋｅｙ　ｗｏｒｄｓ：Ｔｏｐｏｇｒａｐｈｉｃａｌ　Ｍａｐ；Ｐｒｅｃｉｓｉｏｎ；Ｅｌｅｖａｔｉｏｎ；ＲＭＳＥ；Ｓｏｆｔｗａｒｅ