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
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
|
comment: >
Test various transformations depending on datum
files that are not always available
exe: cs2cs
env:
PROJ_SKIP_READ_USER_WRITABLE_DIRECTORY: YES
tests:
- comment: 1st through ntv1, 2nd through conus
args: +proj=latlong +ellps=clrk66 +nadgrids=ntv1_can.dat,conus +to +proj=latlong +datum=NAD83
in: |
111d00'00.000"W 44d00'00.000"N 0.0
111d00'00.000"W 39d00'00.000"N 0.0
out: |-
111d0'3.208"W 43d59'59.732"N 0.000
111d0'2.604"W 38d59'59.912"N 0.000
- comment: As above, but without ntv1 everything goes through conus file
# see data/CMakeLists.txt for "dir with space/myconus"
args: +proj=latlong +ellps=clrk66 '+nadgrids="dir with space/myconus"' +to +proj=latlong +datum=NAD83
in: |
111d00'00.000"W 44d00'00.000"N 0.0
111d00'00.000"W 39d00'00.000"N 0.0
out: |-
111d0'2.788"W 43d59'59.725"N 0.000
111d0'2.604"W 38d59'59.912"N 0.000
- comment: Test --bbox -141.01,40.04,-47.74,83.17 NAD27 NAD83 (using NTv2)
args: --bbox -141.01,40.04,-47.74,83.17 NAD27 NAD83
in: 40.5 -60
out: |
40d30'0.368"N 59d59'56.617"W 0.000
- comment: Test --area "USA - CONUS including EEZ" NAD27 NAD83 (using conus)
args: --area "USA - CONUS including EEZ" NAD27 NAD83
in: |
44 -111
40.5 -60
out: |
43d59'59.725"N 111d0'2.788"W 0.000
* * inf
- comment: Test MD used where available
args: +proj=latlong +ellps=clrk66 +nadgrids=MD,conus +to +proj=latlong +datum=NAD83
in: |
79d58'00.000"W 37d02'00.000"N 0.0
79d58'00.000"W 36d58'00.000"N 0.0
out: |-
79d58'0.005"W 37d1'59.998"N 0.000
79d57'59.128"W 36d58'0.501"N 0.000
- comment: Similar to previous test, but using only conus
args: +proj=latlong +ellps=clrk66 +nadgrids=conus +to +proj=latlong +datum=NAD83
in: |
79d58'00.000"W 37d02'00.000"N 0.0
79d58'00.000"W 36d58'00.000"N 0.0
out: |-
79d57'59.126"W 37d2'0.501"N 0.000
79d57'59.128"W 36d58'0.501"N 0.000
- comment: >
Test that we use grid shift files even if we are right on the
edge or even a wee bit outside (#141).
# Our test points are (1) right on mesh corner, (2) outside but within
# epsilon (3) inside a bit (4) outside by more than epsilon
args: +proj=latlong +ellps=WGS84 +nadgrids=ntf_r93.gsb,null +to +proj=latlong +datum=WGS84 -f %.12f
in: |
-5.5 52.0
-5.5000000000001 52.0000000000001
-5.4999 51.9999
-5.5001 52.0
out: |-
-5.501106465528 51.999890470284 0.000000000000
-5.501106465529 51.999890470284 0.000000000000
-5.501006458305 51.999790470257 0.000000000000
-5.500100000000 52.000000000000 0.000000000000
- comment: Inverted version of previous test
args: +proj=latlong +datum=WGS84 +to +proj=latlong +ellps=WGS84 +nadgrids=ntf_r93.gsb,null -f %.12f
in: |
-5.5 52.0
-5.5000000000001 52.0000000000001
-5.4999 51.9999
-5.5001 52.0
out: |-
-5.498893534472 52.000109529716 0.000000000000
-5.498893534472 52.000109529717 0.000000000000
-5.498793593803 52.000009531513 0.000000000000
-5.500100000000 52.000000000000 0.000000000000
- comment: "NAD27 -> NAD83: 1st through ntv1 or ntv2, 2nd through conus"
args: NAD27 NAD83 -f %.4f
in: |
55d00'00.000"N 111d00'00.000"W 0.0
39d00'00.000"N 111d00'00.000"W 0.0
out: |-
55.0001 -111.0009 0.0000
39.0000 -111.0007 0.0000
- comment: WGS84 (2D) -> WGS84+EGM96
args: +init=epsg:4326 +to +init=epsg:4326 +geoidgrids=egm96_15.gtx
in: 2dE 49dN 0
out: |-
2.00 49.00 0.00
- comment: WGS84 (3D) -> WGS84+EGM96
args: +init=epsg:4979 +to +init=epsg:4326 +geoidgrids=egm96_15.gtx
in: 2dE 49dN 0
out: |-
2.00 49.00 -45.06
- comment: WGS84 (2D), promoted to 3D -> WGS84+EGM96
args: --3d +init=epsg:4326 +to +init=epsg:4326 +geoidgrids=egm96_15.gtx
in: 2dE 49dN 0
out: |-
2.00 49.00 -45.06
|
