地图投影变换开源包PROJ.pptVIP

PROJ.4 PROJ.4 是开源的投影变换包 /proj/ 基本函数 pj_init：初始化投影参数 pj_init_plus：初始化投影参数 pj_transform：投影变换 pj_fwd：投影正变换 pj_inv：投影逆变换 pj_free：释放内存 projPJ pj_init(int agrc, char ** argv) 常用于接收命令行参数 projPJ pj_init_plus(const char *definition); 初始化投影参数，返回projPJ结构 int pj_transform( projPJ srcdefn, projPJ dstdefn, long point_count, int point_offset, double *x, double *y, double *z ); 进行投影变换 x,y,z是输入、输出参数 z可为NULL projXY pj_fwd(projLP lp, projPJ pj); 将经纬度转换为平面直角坐标 projLP pj_inv(projXY xy, projPJ pj); 将平面直接坐标转换为经纬度 void pj_free( projPJ pj ); 释放projPJ结构内存 #include proj_api.h main(int argc, char **argv) { projPJ pj_merc, pj_latlong; double x, y; if (!(pj_merc = pj_init_plus(+proj=merc +ellps=clrk66 +lat_ts=33)) ) exit(1); if (!(pj_latlong = pj_init_plus(+proj=latlong +ellps=clrk66)) ) exit(1); while (scanf(%lf %lf, x, y) == 2) { x *= DEG_TO_RAD; y *= DEG_TO_RAD; p = pj_transform(pj_latlong, pj_merc, 1, 1, x, y, NULL ); printf(%.2f\t%.2f\n, x, y); } exit(0); } #include proj_api.h int main(int argc, char **argv) { char *args[] = { proj=merc, ellps=clrk66, lat_ts=33 }; projUV p; projPJ pj; if (!(pj = pj_init(3, args))) exit(1); while (scanf(%lf %lf, p.v, p.u) == 2) { p.u *= DEG_TO_RAD; p.v *= DEG_TO_RAD; p = pj_fwd(p, pj); printf(%.2f\t%.2f\n, p.u, p.v); } exit(0); } 其他重要函数 pj_datum_transform pj_geocentric_to_geodetic pj_geodetic_to_geocentric pj_compare_datums pj_is_latlong pj_is_geocent int pj_datum_transform( projPJ src, projPJ dst, long point_count, int point_offset, double *x, double *y, double *z ) 基准变换（BLH=BLH） x,y,z输入和输出均为BLH 注意： src、dst中必须有towgs84参数 3参数：dx,dy,dz +towgs84=-199.87,74.79,246.62