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

坐标转换从经纬度坐标到大地坐标利用网络上开源的资料，可以很容易的实现从经纬度坐标向各种投影坐标的转换，美国地质调查局开发USGS的GCTP就是很好的东西之一，有Ｃ语言版本的支持各种投影类型的源代码，比如UNIVERSAL TRANSVERSE MERCATOR， ALBERS CONICAL EQUAL AREA ，LAMBERT CONFORMAL CONIC等等，我们机房使用的就是LAMBERT CONFORMAL CONIC（LAMBERT 圆锥等角投影），利用GCTP提供的源代码详见source文件夹下的lamccfor.c　和　lamccinv.c　，可以轻松实现经纬度投影坐标和LAMBERT CONFORMAL CONIC坐标的相互转化，其他投影方式使用的比较少，还要进一步研究，就以后在写了．lamccforint函数设置LAMBERT CONFORMAL CONIC投影的各个参数，比如长半轴，短半轴，中心点经纬度坐标，标准纬线等等信息．long lamccforint(r_maj,r_min,lat1,lat2,c_lon,c_lat,false_east,false_north)double r_maj;???????????????? /* majoraxis */double r_min;???????????????? /* minoraxis */double lat1;??????????????????? /* first standard parallel????????????? */double lat2;??????????????????? /* second standard parallel???????????? */double c_lon;?????????????????? /* center longitude???????????????????? */double c_lat;?????????????????? /* center latitude????????????????????? */double false_east;????????????? /* x offset in meters?????????????????? */double false_north;???????????? /* y offset in meters?????????????????? */{double sin_po;????????????????? /* sin value??????????????????????????? */double cos_po;????????????????? /* cos value??????????????????????????? */double con;???????????????????? /* temporary variable?????????????????? */double ms1;???????????????????? /* small m 1??????????????????????????? */double ms2;???????????????????? /* small m 2??????????????????????????? */double temp;??????????????????? /* temporary variable?????????????????? */double ts0;???????????????????? /* small t 0??????????????????????????? */double ts1;???????????????????? /* small t 1??????????????????????????? */double ts2;???????????????????? /* small t 2??????????????????????????? */r_major = r_maj;r_minor = r_min;false_northing = false_north;false_easting = false_east;/* Standard Parallels cannot be equal and on opposite sides of the equator------------------------------------------------------------------------*/if (fab