【地理信息系统中的量算】转载来源：
http://space.itpub.net/2052/viewspace-166491
地理信息系统中牵涉到各种各样的量算，例如距离、方位、面积等等。这些量算的方法或者计算繁复，或者较难找到资料，从今天开始逐渐贴出自己整理的计算方法，也作为自己的资料积累。
距离量算：
距离量算可以通过不同的方式进行：
1、将经纬度坐标通过投影方式投影到二维坐标，然后通过二维坐标的取值计算两点之间的距离。这种方法计算量较大，同时投影通常都会引起误差。在需要进行快速计算时不推荐使用；
2、大圆航线距离计算。按照定义，大圆线是球面上半径等于球体半径的圆弧。大圆线是连接球面上两点最短的路径所在的曲线。大圆距离的公式是：
d=acos(sin(lat1)*sin(lat2)+cos(lat1)*cos(lat2)*cos(lon1-lon2))
其中，lat1和lon1分别为第一点纬度和经度，而lat2和lon2分别为第二点的纬度和经度。而最后的结果d的单位是度或者弧度，而根据定义，60海里为一度，1海里等于1852米。
另一个对于短距离计算误差较小的大圆长度公式是：
d=2*asin(sqrt((sin((lat1-lat2)/2))^2 + cos(lat1)*cos(lat2)*(sin((lon1-lon2)/2))^2))
3、恒向线距离。恒向线的定义是：即沿两点间固定方位角航行的航线。这不是航程最短的航线，而是操纵方便的沿单一航向航行的航线。恒向线的距离要大于两点之间的大圆线距离。
恒向线距离d的公式是：
if (abs(lat2-lat1) < sqrt(TOL))
{ q=cos(lat1) }
else
{ q= (lat2-lat1)/log(tan(lat2/2+pi/4)/tan(lat1/2+pi/4)) } d=sqrt((lat2-lat1)^2+ q^2*(lon2-lon1)^2)
其中，TOL是一个很小的数，例如1E-15。pi是圆周率。当恒向线穿过东经180度线时，这个公式会有问题。改进的公式是这样的：
dlon_W=mod(lon2-lon1,2*pi) dlon_E=mod(lon1-lon2,2*pi) dphi=log(tan(lat2/2+pi/4)/tan(lat1/2+pi/4)) if (abs(lat2-lat1) < sqrt(TOL))
{ q=cos(lat1) }
else
{ q= (lat2-lat1)/dphi } if (dlon_W < dlon_E)
{// Westerly rhumb line is the shortest tc=mod(atan2(-dlon_W,dphi),2*pi) d= sqrt(q^2*dlon_W^2 + (lat2-lat1)^2) }
else
{ tc=mod(atan2(dlon_E,dphi),2*pi) d= sqrt(q^2*dlon_E^2 + (lat2-lat1)^2) }
其中tc是方位角，mod的定义是：
mod=y - x * int(y/x) if ( mod < 0) mod = mod + x
