UTM投影正转算法实现

    在Gissky上发现了戴勤奋老师关于坐标转换的程序以及公式，程序的界面和转换精度都值得称赞，只是可惜没有提供实现的源码，最近由于工作需要正好需要实现这部分功能，通过对戴勤奋老师提供的公式以及OGP（欧洲石油勘探组织）最新的《Coordinate Conversions and Transformations including Formulas》文档进行了研究，利用程序将UTM投影正转实现，现发布出来供大家共同学习，也希望大家指出程序中的不足。

 

/**/ /****************************************************************************
 ** 函数名：LLtoUTM
 ** 输入：   a                 椭球体长半轴
 **          f                 椭球体扁率  f=(a-b)/a 其中b代表椭球体的短半轴
 **          Lat               经过UTM投影之前的纬度
 **          Long              经过UTM投影之前的经度
 **          LongOrigin        中央经度线
 **          FN                纬度起始点，北半球为0，南半球为10000000.0m
 **          UTMNorthing 经过UTM投影后的纬度方向的坐标
 **          UTMEasting  经过UTM投影后的经度方向的坐标      
 ** 功能描述：UTM投影正转
 ** 作者：   CiCi
**  单位：   中国地质大学（北京）地球科学与资源学院
 ** 创建日期：2007年3月28日
 ** 版本：1.0
 ***************************************************************************/
LLtoUTM( double  a, double  f, const   double  Lat,  const   double  Long,
           double  LongOrigin, double  FN, double