PAGE

PAGE #/ 7

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

利用网络上开源的资料，可以很容易的实现从经纬度坐标向各种投影坐标 的转换，美国地质调查局开发 USGS勺GCTP就是很好的东西之一，有C语言版 本的支持各种投影类型的源代码，比如 UNIVERSAL TRANSVERSE MERCA，TOR ALBERSCONICAL EQUAL AREELAMBERT CONFORMAL CON等等，我们机房使用 的就是LAMBERTCONFORMALCONQAMBERT圆锥等角投影)，利用 GCTP提 供的源代码详见source文件夹下的lamccfor.c和lamccinv.c,可以轻松实现经纬 度投影坐标和LAMBERTCONFORMALCON坐标的相互转化，其他投影方式使用 的比较少，还要进一步研究，就以后在写了．

lamccforint函数设置LAMBERTCONFORMALCON投影的各个参数，比如长 半轴，短半轴，中心点经纬度坐标，标准纬线等等信息．

long

lamccforint(r_maj,r_min,lat1,lat2,c_lon,c_lat,false_east,false_north)double

r_maj;/* majoraxis */

r_maj;

double r_min;/*

double r_min;

/* minoraxis */

double lat1;/* first standardparallel*/double lat2;/* second standardparallel*/double c_lon; longitude double c_lat;/* center*//*

double lat1;

/* first standard

parallel

*/

double lat2;

/* second standard

parallel

*/

double c_lon; longitude double c_lat;

/* center

*/

/* 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; variable double ms1;1double ms2;2

double con; variable double ms1;

1

double ms2;

2

double temp; variable double ts0;

0 double ts1;

1 double ts2;

2

*/

/* small m

*/

/* small m

*/

/* temporary

*/

/* small t

*/

/* small t

*/

/* small t

*/

r_major = r_maj;

r_minor = r_min;

false_northing = false_north;

false_easting = false_east;

/*StandardParallelscannotbeequalandonoppositesidesoftheequator

*/

if (fabs(lat1 lat2) < EPSLN){p_error("Equal latitudes for St. Parallels on opposite sides ofequator",

"lamcc-for");

return

(41);}temp = r_minor / r_major;

es =

1.0 - SQUARE(temp);

e = sqrt(es);

center_lon = c_lon;

center_lat = c_lat;

sincos(lat1,&sin_po,&cos_po);

con = sin_po;

ms1 = msfnz(e,sin_po,cos_po);

ts1 = tsfnz(e,lat1,sin_po);

sincos(lat2,&sin_po,&cos_po);

ms2 = msfnz(e,sin_po,cos_po);

ts2 = tsfnz(e,lat2,sin_po);

sin_po = sin(center_lat);

ts0 = tsfnz(e,center_lat,sin_po);

if (fabs(lat1 - lat2) > EPSLN)

ns = log (ms1/ms2)/ log (ts1/ts2)