投影转换公式

常用地图投影转换公式青岛海洋地质研究所戴勤奋(Email:qddginfenecgs.gov.cn)最近几乎天天都有Emai1跟我要这样、那样的坐标系转换或投影转换公式，或问我编的投影程序公式是哪来的，有没有专门介绍投影公式的书等等，让我越来越觉得有必要就此方面写点东西，一来我自己总结一下，二来对那些我没有回Emai1的同行也有个交代，因为那些公式实在太难敲了。我在“海洋地质制图常用地图投影系列小程序”(http:/www.gissky..net)中用的公式来自我原来的积累，同时参考了POSC(http:/www.posc.org，国际石油技术软件开放公司)的文献“CoordinateConversions and Transformation including Formulas”,该文献由EPSG(http:/www.epsg.org,欧洲石油勘探组)编写，比较全面地介绍了各种地图投影与坐标系的转换方法及计算公式，而且最新更新到了2004年，是我目前看到的最全面、最新的相关文档了，只不过是英文的，我正在打算将它们翻成中文，到时与大家共享。投影计算公式往往表达方式不止一种，有时很难分辨谁对谁错，我只把“墨卡托投影”、“高斯-克吕格投影”、“UTM投影”、“兰勃特等角投影”(1：100万地形图规范中称作正轴等角圆锥投影，GB/T14512-93)的正反转换公式列出，因为我基本能保证这些公式的正确性。“海洋地质制图常用地图投影系列小程序”(http:/www.gissky.net)己升级，原下载者请注意下载更新版本。1.约定本文中所列的转换公式都基于椭球体a-椭球体长半轴b一椭球体短半轴f-扁率(a-b)/ae第一偏心率e=√h-(b/a)2e'一第二偏心率e'=Va/b2-1(a2/b)N一卯酉圈曲率半径N=v1+e2*cos2 Ba(1-e2)R-子午圈曲率半径R=(1-e2*sin2B)2B--纬度，L-经度，单位弧度(RAD)Xw-一纵直角坐标，Y一横直角坐标，单位米()