1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
|
#include <iostream>
#include <proj.h>
int test_transform() {
PJ *P;
PJ_COORD a, b;
P = proj_create_crs_to_crs(
PJ_DEFAULT_CTX,
"EPSG:4326",
"+proj=utm +zone=32 +datum=WGS84", // or EPSG:32632
NULL);
if (0 == P) {
std::cerr << "Oops" << std::endl;
return 1;
}
// Copenhagen: 55d N, 12d E
a = proj_coord(55, 12, 0, 0);
b = proj_trans(P, PJ_FWD, a);
std::cout.precision(2);
std::cout.setf( std::ios::fixed, std::ios::floatfield );
std::cout << "easting: " << b.enu.e << ", northing: " << b.enu.n;
b = proj_trans(P, PJ_INV, b);
std::cout << ", latitude: " << b.lp.lam << ", longitude: " << b.lp.phi << std::endl;
proj_destroy(P);
return 0;
}
int main(int argc, char *argv[]) {
PJ_INFO info;
info = proj_info();
if(argc == 2 && argv[1][0] == '-') {
switch(argv[1][1]) {
case 't':
return(test_transform());
case 's':
std::cout << info.searchpath << std::endl;
return(0);
case 'v':
std::cout << info.major << '.' << info.minor << '.' << info.patch << std::endl;
return(0);
}
}
std::cerr << "Use option -t, -s or -v" << std::endl;
return(1);
}
|
