关于椭球区域面积计算问题的讨论.docVIP

关于椭球区域面积计算问题的讨论 　　摘 要：目前关于椭球区域面积的计算都是采用边界点高斯平面坐标进行，由于高斯投影为等角投影存在面积变形，因而所计算的面积和实际面积有一定差距。椭球面上的不规则面积计算则显得尤为复杂困难。文章讨论顾及地球曲率的面积计算方法，转换为采用等面积投影的方法，采用平面坐标面积计算，对不规则椭球区域面积计算方法进行讨论。 　　关键词：椭球规则梯形计算；等面积投影；高斯投影；等角投影；高斯正反算 　　1 概述 　　目前椭球面面积的计算都是采用边界线的高斯平面坐标进行，没有考虑地球曲率。由于高斯投影存在面积变形，虽然单宗土地面积变形不大，但是全省、全国大面积统计，则影响不可忽视，因此精确计算椭球面区域面积，是国土部门亟待解决的问题。目前提出的方法多种多样，如利用傅立叶级数快速转化实现面积计算、利用freeman链码矢量分析对边界进行综合处理获取边界像素坐标加权求和，求得面积等方法。上述的方法，最终并没有成为解决椭球面区域计算的方法。 　　椭球梯形是椭球面上唯一能直接计算出准确面积的图形，它是由两条子午线和两条平行圈围成的梯形表面。但事实生产工作中并不会简单的计算梯形面积，而是需要对不规则的图形进行计算。本文希望能通过对简单投影方法的运用，得到区域面积的计算简便方法，并利用椭球梯形作为实际面积进行检验。利用等面积投影特有的投影后面积不变的特点和高斯投影直接利用坐标计算面积的方式，将不规则椭球区域坐标转化为等面积投影和高斯投影坐标，再用平面面积计算公式计算不规则区域。通过具体数据比较高斯投影与等面积的投影转化方法能否解决不规则区域面积计算，并用椭球梯形面积进行检验。 　　2 椭球面积计算方法 　　2.1 规则梯形面积计算 　　地球为一个不规则的球体，广大区域上遍布着江河，湖海，高山，盆地，低洼，峡谷等等，因此地球上区域的面积计算就变得很困难，此时需要引入一个类地球的椭球，这个椭球的目的主要是为了方便地球表面上的测量计算工作。这个数学模型为规则的，它非常接近大地体并用来替代大地体。我国的参考椭球经历了几代更迭，从1952年前的海福特椭球，至1953年开始使用克拉索夫斯基椭球，1978年则使用75国际椭球，通过75椭球的基础数据建立了我国自己的80坐标系，现在国家正在大力推行CGS2000坐标。 　　参考椭球虽然替代了大地体，但对于它上面的不规则区域面积计算仍然没有简单有效的解决办法。虽说不规则的区域无法准确计算，但对于椭球梯形，通过微分、积分的方式能够计算其面积，并总结得出以下公式： 　　公式中的xy，需要先经过（x，y）高斯投影→（L，B）大地坐标→（x，y）等面积投影来进行转换，得到等面积投影的横纵坐标。其中的S是一弧度对应的纬差B所构成的梯形面积。 　　3 不规则区域面积计算方法的讨论 　　3.1 规则区域面积计算 　　3.1.1 高斯投影面积计算分析 　　由于高斯投影需要控制投影变形，就必须控制投影区域的范围。面积计算受到分带、形变等的影响，但是其方法简便，直接利用平面坐标计算的方法，依然有可取之处。通过以下数据说明高斯投影与实际面积的差距。 　　表1中的数据取自不同经度差，不同大小的梯形。是为了能得到随机数据，通过上表可以看出，在利用高斯投影直接进行梯形计算时绝对误差较大，得到的相对误差较等面积投影的相对误差明显偏大。根据上表中当经度差等于一度时数据可以看出，实际面积和高斯面积相差很大达到425440平方米，这个误差相当于59个足球场的面积。这个面积差相当明显不能忽视。当经度差减小时，实际面积和高斯面积之间的差值在不断减少，但是这个数据仅仅是对梯形面积进行的比较，还未对不规则图形进行分析。由于高斯投影面积变形与计算区域的中央子午线距离有关，而针对论文中讨论的用梯形作为实际面积的参考物，经度差也就决定了实际面积与投影面积的关系，利用高斯投影做标准梯形的替代情况下，在经度差很小（如表面积相当于121公顷）甚至面积更小的情况下，获得相对误差数量级在E-07这个级别。 　　3.1.2 等面积投影分析 　　由方案一可知单独考虑用高斯投影来解决不规则区域面积的计算，不太可行。考虑到等面积投影虽角度长度都有变形，但整体面积投影后没有发生变化，这正是解决区域面积计算的理想方法，利用等面积投影的方法，理论上能够避免不规则区域面积计算的转换损失。将椭球梯形面积与等面积投影面积进行比较可以得到如下数据（Dc距中央子午线距离，De距赤道距离，δ相对误差，△B经度差）。 　　通过与椭球梯形对比，由上表可以看出： 　　随离中央子午线距离越远，椭球梯形面积与等面积投影面积并无明显差异，距离中央子午线远近，不会决定面积的大小变化。 　　对于同样的椭球梯形，面积不会随经度变化而变换，当纬度变化时，