坐标转换从经纬度坐标到大地坐标及源码.doc

《坐标转换从经纬度坐标到大地坐标及源码.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、坐标转换从经纬度坐标到大地坐标及源码利用网络上开源的资料，可以很容易的实现从经纬度坐标向各种投影坐标的转换，美国地质调查局开发USGS的GCTP就是很好的东西之一，有Ｃ语言版本的支持各种投影类型的源代码，比如UNIVERSALTRANSVERSEMERCATOR，ALBERSCONICALEQUALAREA，LAMBERTCONFORMALCONIC等等，我们机房使用的就是LAMBERTCONFORMALCONIC（LAMBERT圆锥等角投影），利用GCTP提供的源代码详见source文件夹下的lamccfor.c　和　lamccinv.c　，可以轻松实现经纬度投影坐标和LAMBERT

2、CONFORMALCONIC坐标的相互转化，其他投影方式使用的比较少，还要进一步研究，就以后在写了．lamccforint函数设置LAMBERTCONFORMALCONIC投影的各个参数，比如长半轴，短半轴，中心点经纬度坐标，标准纬线等等信息．longlamccforint(r_maj,r_min,lat1,lat2,c_lon,c_lat,false_east,false_north)doubler_maj;                /*majoraxis*/doubler_min;                /*minoraxis*/doublelat1;        

3、           /*firststandardparallel             */doublelat2;                   /*secondstandardparallel            */doublec_lon;                  /*centerlongitude                    */doublec_lat;                  /*centerlatitude                     */doublefalse_east;             /*xoffsetin

4、meters                  */doublefalse_north;            /*yoffsetinmeters                  */{doublesin_po;                 /*sinvalue                           */doublecos_po;                 /*cosvalue                           */doublecon;                    /*temporaryvariable               

5、   */doublems1;                    /*smallm1                           */doublems2;                    /*smallm2                           */doubletemp;                   /*temporaryvariable                  */doublets0;                    /*smallt0                           */doublets1;        

6、            /*smallt1                           */doublets2;                    /*smallt2                           */r_major=r_maj;r_minor=r_min;false_northing=false_north;false_easting=false_east;/*StandardParallelscannotbeequalandonoppositesidesoftheequator------------------------------------

7、------------------------------------*/if(fabs(lat1lat2)