Proj.4 API 中文参考

ProjAPI

https://github.com/OSGeo/proj.4/wiki/ProjAPI

Tom Kralidis在2015年5月27日编辑此页·修订4

简介

执行pj_init()选择并初始化一个由它的参数决定参数的制图投影坐标系。argc是argv数组中的元素个数，argv数组每个元素是单个投影的关键字参数(+proj arguments)。该列表中必须包含proj=投影、地球半径或椭圆参数。如果初始化的投影是有效的，返回一个projPJ结构体的地址，失败返回NULL。

pj_init_plus()函数的工作方式与pj_init()类似，但将投影参数定义在单个字符串内，且每个参数带有一个+前缀。例如:+proj=utm +zone=11 +ellps=WGS84。

pj_init()的返回值用作投影正算(forward)或者投影逆算(inverse)中的proj参数。投影计算中的参数结构体projUV中的U代表经度或者x，V代表维度或者y，经度和维度都是以弧度单位的。如果操作失败，projUV中的两个值都设置为HUGE_VAL(在math.h中定义http://blog.csdn.net/poechant/article/details/7235594)。

注意：所有投影都有正算(forward)模式，但却不是都有逆算(inverse)模式。如果投影没有逆算模式，那么projPJ中成员inv将设置为NULL。

pj_transform函数用于转换两个坐标系统之间的点。除了投影坐标与地理坐标之间的转换，也可以进行基准面(datum)之间的转换。不同于pj_fwd()和pj_inv()，这里也允许坐标系定义参数projPJ*是地理坐标(使用+proj=latlong定义)。x、y、z数组用于传入计算值，传出计算结果。该函数成功返回0，失败返回错误号(查看inpj_errno)。

由projection(projPJ)关联的内存，可能需要使用pj_free()释放。

示例

The following program reads latitude and longitude values in decimal degress, performs Mercator projection with a Clarke 1866 ellipsoid and a 33° latitude of true scale and prints the projected cartesian values in meters:

下面程序读取十进制的经度和维度值，进行克拉克1866椭球正确范围33°维度的墨卡托投影计算，并打印米单位的笛卡尔坐标值。

北美大地基准维基百科

示例

The following program reads latitude and longitude values in decimal degress, performs Mercator projection with a Clarke 1866 ellipsoid and a 33° latitude of true scale and prints the projected cartesian values in meters:

下面程序读取十进制的经度和维度值，进行克拉克1866椭球正确范围33°维度的墨卡托投影计算，并打印米单位的笛卡尔坐标值。

#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);
}

在此程序中,输入-16 20.25将得到结果-1495284.21 1920596.79.

API函数

基本API

int pj_transform( projPJ srcdefn, projPJ dstdefn, long point_count, int point_offset, double *x, double *y, double *z );

转换点x/y/z，从源坐标系统到目标坐标系统(更多信息).

projPJ pj_init_plus(const char *definition);

从字符串definition创建一个projPJ坐标系统对象(更多信息).

void pj_free( projPJ pj );

释放与pj关联的所有资源.

高级函数

int pj_is_latlong( projPJ pj );

返回TRUE表示是地理坐标系统(proj=latlong).

int pj_is_geocent( projPJ pj );

返回TRUE表示是地心坐标系统(proj=geocent).

char *pj_get_def( projPJ pj, int options);

返回PROJ.4初始化字符串，用于pj_init_plus()可产生此坐标系统。但会尽可能多的扩大定义(例如添加:+init= 和 +datum= definitions).

projPJ pj_latlong_from_proj( projPJ pj_in );

返回一个新的坐标系定义。就是pj_in内部的地理坐标(lat/long)系统.

环境函数

void pj_set_finder( const char *(*new_finder)(const char *) );

Install a custom function for finding init and grid shift files.

安装用于查找init和grid shift(网格转换)文件的自定义函数。

void pj_set_searchpath ( int count, const char **path );

Set a list of directories to search for init and grid shift files.

设置一个目录列表去搜索init和grid shift文件。

void pj_deallocate_grids( void );

Frees all resources associated with loaded and cached datum shift grids.

释放加载和缓存的基准转换网格关联的所有资源。

char *pj_strerrno( int );

Returns the error text associated with the passed in error code.

返回错误代码关联的错误文本。

int *pj_get_errno_ref( void );

Returns a pointer to the global pj_errno error variable.

返回全局的pj_errno错误变量地址。

const char *pj_get_release( void );

Returns an internal string describing the release version.

返回描述发行版本的内部字符串。

废弃(过时)函数

XY pj_fwd( LP lp, PJ *P );

LP pj_inv( XY xy, PJ *P );

projPJ pj_init(int argc, char **argv);