Mapinfo在森林资源二类调查面积计算中的应用研究
2010-06-15张士朋
张士朋
(丰宁满族自治县林业局,河北 丰宁 068350)
Mapinfo在森林资源二类调查面积计算中的应用研究
张士朋
(丰宁满族自治县林业局,河北 丰宁 068350)
该文通过对mapinfo软件应用功能的深刻剖析,了解其面积计算环境及各种不同的面积计算方法。结合丰宁县二类调查数据对上述方法的适用性进行检验性分析研究,最终得出正确的面积计算环境和方法的结论。
Mapinfo软件;二类调查;面积计算
1 森林资源二类调查面积计算环境及基础计算方法
基础环境与面积计算原理:我国的地形图采用高斯-克吕格平面直角坐标系。量算面积要求为小班的平面面积即闭合小班线范围内所勾绘的地形图区域平面面积。小班区域为多边形区域,精确的面积计算方法为解析法。
1.1 小班调绘要求
根据实际情况,可分别采用以下方法进行小班调绘
(1)采用由测绘部门绘制的当地最新的比例尺为1∶10000~1:25000的地形图到现地进行勾绘。对于没有上述比例尺的地区可采用由1∶50000放大到1∶25000的地形图。
(2)使用近期拍摄的(以不超过2年为宜)、比例尺不小于 1:25000 或由 1∶50000 放大到 1∶25000 的航片、1∶100000放大到1∶25000的侧视雷达图片在室内进行小班勾绘,然后到现地核对,或直接到现地调绘。
(3)使用近期(以不超过1年为宜)经计算机几何校正及影像增强的比例尺1∶25000的卫片(空间分辨率10m以内)在室内进行小班勾绘,然后到现地核对。
(4)空间分辨率10m以上的卫片只能作为调绘辅助用图,不能直接用于小班勾绘。
(5)现地小班调绘、小班核对以及为林分因子调查或总体蓄积量精度控制调查而布设样地时,可用GPS确定小班界线和样地位置。
1.2 面积量算要求
(1)按照“层层控制,分级量算,按比例平差”的原则进行面积量算。即先量算林业局(县、保护区、森林公园)的面积,再量算林场(乡、管理站)、林班(村)面积,最后量算小班面积。如无特殊情况,县、乡各级行政单位的面积应与民政部门公布的面积一致。各级面积经准确量算后,复查时除非界线发生变化,否则不准变动。
(2)国有林业局(县、保护区、森林公园)、林场(乡、管理站)的面积用理论图幅面积计算,即将分布在各图幅上的部分累加求得。一个图幅上的各部分面积,要分别量测进行平差。
(3)用地理信息系统(GIS)绘制成果图时,可直接用地理信息系统量算林班和小班面积。手工绘制成果图时,可用几何法、网点网格法或求积仪等量算林班和小班面积。
(4)林场(乡)内各林班面积之和与林场面积相差不到1%,林班内各小班面积之和与林班面积相差不到2%时,可进行平差,超出时应重新量算。
(5)面积量算以公顷为单位,精确到0.1hm2。
2 Mapinfo地理信息系统面积计算环境及基础计算方法
2.1 面积计算环境
丰宁县所处的地理位置东经115.90°~117.39°位于高斯-克吕格投影6度投影带的第20带。因此所选择投影类别为 Gauss-kruger(Pulkovo 1942),类别项为 Gk Zone 20(Pulkovo 1942)。不同地域根据其地理位置经度选择适宜的类别项。距离和面积使用平面坐标。
2.2 面积计算方法
(1)在上述面积计算环境下,使用选择工具双击需要计算的单个多边形区域,在弹出的对话框中即可显示该多边形区域的平面面积。
(2)在上述面积计算环境下,当需要计算多个区域的多边形面积时,为了计算方便,应使用面积计算函数Area()、Cartesian area()和 Spherical area()对面积字段进行更新计算。
3 利用Mapinfo地理信息系统对取样区域进行面积计算分析
(1)在上述计算环境下,对地形图上一平方公里网格(四角坐标分别为:1(20470000,4563000),2(20470000,4564000),3(20471000,4564000),4(20471000,4563000))区域的面积计算,结果见表1。
(2)在上述计算环境下,对丰宁县二类调查面积的统计计算结果见表2。
(3)数据分析结果。①对以上表格数据进行分析,得知在二类调查中在投影选择正确的情况下应使用cartesian area()函数计算小班面积。以确保面积数据的相对正确性,虽然有误差,但误差很小,符合二类调查误差要求。②Area()函数和Spherical area()函数显然不能作为用来求算面积的函数。虽然误差不大,在二类调查误差允许范围内,但这两个函数不适用于二类调查面积计算环境。使用这两个函数计算本身就说明选择方法是错误。
表1 区域的面积计算结果
表2 丰宁县二类调查面积的统计计算结果
4 结论
经过以上分析论证得出使用mapinfo计算森林资源二类调查面积时应该选择的面积计算环境为:投影类别为Gauss-kruger(Pulkovo 1942),类别项为 Gk Zone NO(Pulkovo 1942)(注:NO为所在区域的高斯-克吕格6度投影带的带号:NO=[当地经度/6]+1);距离和面积使用平面坐标。正确的面积计算方法是运用Cartesian area()函数对面积进行数据更新。
Tp31
A
1002-3356(2010)01-0027-01
2009-11-25