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
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
|
-------------------------------------------------------------------------------
===============================================================================
Test the 4D API handling of cs2cs style transformation options.
These tests are mostly based on the same material as those in
more_builtins.gie, since we are testing the same kinds of things,
but provided through a different interface.
===============================================================================
<gie>
-------------------------------------------------------------------------------
Test the handling of the +towgs84 parameter.
-------------------------------------------------------------------------------
(additional tests of the towgs84 handling can be found in DHDN_ETRS89.gie)
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
This example is from Lotti Jivall: "Simplified transformations from
ITRF2008/IGS08 to ETRS89 for maritime applications" (see also more_builtins.gie)
-------------------------------------------------------------------------------
operation proj=geocent
towgs84 = 0.676780, 0.654950, -0.528270,
-0.022742, 0.012667, 0.022704,
-0.01070
-------------------------------------------------------------------------------
tolerance 1 um
direction inverse
accept 3565285.00000000 855949.00000000 5201383.00000000
expect 3565285.41342351 855948.67986759 5201382.72939791
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
This example is a random point, transformed from ED50 to ETRS89 using KMStrans2.
-------------------------------------------------------------------------------
operation proj=latlong ellps=intl
towgs84 = -081.07030, -089.36030, -115.75260,
000.48488, 000.02436, 000.41321, -0.540645
-------------------------------------------------------------------------------
tolerance 25 mm
accept 16.82 55.17 61.0
expect 16.8210462130 55.1705688946 29.0317
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
operation proj=latlong nadgrids=nzgd2kgrid0005.gsb ellps=GRS80
-------------------------------------------------------------------------------
This functionality is also tested in DHDN_ETRS89.gie
-------------------------------------------------------------------------------
tolerance 1 mm
ignore pjd_err_failed_to_load_grid
accept 173 -45
expect 172.999892181021551 -45.001620431954613
direction inverse
accept 172.999892181021551 -45.001620431954613
expect 173 -45
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
operation proj=latlong geoidgrids=egm96_15.gtx ellps=GRS80
-------------------------------------------------------------------------------
tolerance 15 cm # lax tolerance due to widespread bad egm96 file
ignore pjd_err_failed_to_load_grid
accept 12.5 55.5 0
expect 12.5 55.5 -36.0213
direction inverse
accept 12.5 55.5 -36.0213
expect 12.5 55.5 0
-------------------------------------------------------------------------------
operation proj=merc geoidgrids=egm96_15.gtx ellps=GRS80
-------------------------------------------------------------------------------
tolerance 0.1 mm
ignore pjd_err_failed_to_load_grid
accept 12.5 55.5 0
expect 1391493.63492 7424275.19462 -36.0213
direction inverse
accept 1391493.63492 7424275.19462 -36.0213
expect 12.5 55.5 0
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
Same as the two above, but also do axis swapping.
-------------------------------------------------------------------------------
NOTE: A number of the tests below are commented out. The actually do the
right thing, but the gie distance computation is not yet able to cope
with "unusual" axis orders
-------------------------------------------------------------------------------
operation proj=latlong geoidgrids=egm96_15.gtx axis=neu ellps=GRS80
-------------------------------------------------------------------------------
tolerance 15 cm # lax tolerance due to widely distributed, bad egm96 file
ignore pjd_err_failed_to_load_grid
accept 12.5 55.5 0
expect 55.5 12.5 -36.0213
direction inverse
accept 55.5 12.5 -36.0213
expect 12.5 55.5 0
-------------------------------------------------------------------------------
operation proj=latlong geoidgrids=egm96_15.gtx axis=dne ellps=GRS80
-------------------------------------------------------------------------------
tolerance 15 cm # lax tolerance due to widely distributed, bad egm96 file
ignore pjd_err_failed_to_load_grid
# accept 12.5 55.5 0
# expect 36.0213 55.5 12.5
# direction inverse
# accept 36.0213 55.5 12.5
# expect 12.5 55.5 0
-------------------------------------------------------------------------------
operation proj=merc geoidgrids=egm96_15.gtx ellps=GRS80
-------------------------------------------------------------------------------
tolerance 0.1 mm
ignore pjd_err_failed_to_load_grid
accept 12.5 55.5 0
expect 1391493.63492 7424275.19462 -36.0213
direction inverse
accept 1391493.63492 7424275.19462 -36.0213
expect 12.5 55.5 0
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
Some more complex axis swapping.
-------------------------------------------------------------------------------
operation proj=latlong geoidgrids=egm96_15.gtx axis=nue ellps=GRS80
-------------------------------------------------------------------------------
tolerance 15 cm # lax tolerance due to widely distributed, bad egm96 file
ignore pjd_err_failed_to_load_grid
accept 12.5 55.5 0
expect 55.5 -36.0213 12.5
# direction inverse
# accept 55.5 -36.0213 12.5
# expect 12.5 55.5 0
-------------------------------------------------------------------------------
operation proj=merc geoidgrids=egm96_15.gtx axis=sue ellps=GRS80
-------------------------------------------------------------------------------
tolerance 15 cm
ignore pjd_err_failed_to_load_grid
accept 12.5 55.5 0
expect -7424275.1946 -36.0213 1391493.6349 0.0000
# direction inverse
# accept -7424275.1946 -36.0213 1391493.6349 0.0000
# expect 12.5 55.5 0
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
A test case from a comment by Github user c0nk
-------------------------------------------------------------------------------
operation proj=somerc
lat_0=46.95240555555556 lon_0=7.439583333333333 k_0=1
x_0=2600000 y_0=1200000 ellps=bessel
towgs84=674.374,15.056,405.346 no_defs
-------------------------------------------------------------------------------
tolerance 20 cm
accept 7.438632495 46.951082877
expect 2600000.0 1200000.0
-------------------------------------------------------------------------------
Same test, but now implemented as a pipeline. This is for testing a nasty bug,
where, at the end of pipeline creation, a false warning about missing ellps was
left behind from the creation of the Helmert step (now repaired in pj_init).
-------------------------------------------------------------------------------
operation proj=pipeline
step proj=cart ellps=WGS84 no_defs
step proj=helmert x=674.37400 y=15.05600 z=405.34600 inv no_defs
step proj=cart ellps=bessel no_defs inv
step proj=somerc lat_0=46.95240555555556 lon_0=7.439583333333333
k_0=1 x_0=2600000 y_0=1200000 ellps=bessel units=m no_defs
-------------------------------------------------------------------------------
tolerance 20 cm
accept 7.438632495 46.951082877
expect 2600000.0 1200000.0
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
Make sure that transient errors are returned correctly.
-------------------------------------------------------------------------------
operation +proj=geos +lon_0=0.00 +lat_0=0.00 +a=6378169.00 +b=6356583.80 +h=35785831.0
-------------------------------------------------------------------------------
accept 85.05493299 46.5261074
expect failure
accept 85.05493299 46.5261074 0
expect failure
accept 85.05493299 46.5261074 0 0
expect failure
-------------------------------------------------------------------------------
Test that Google's Web Mercator works as intended (see #834 for details).
-------------------------------------------------------------------------------
operation proj=pipeline step init=epsg:26915 inv step init=epsg:3857
-------------------------------------------------------------------------------
tolerance 20 cm
accept 487147.594520173 4934316.46263998
expect -10370728.80 5552839.74
accept 487147.594520173 4934316.46263998 0
expect -10370728.80 5552839.74 0
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
Test Google's Web Mercator with +proj=webmerc
-------------------------------------------------------------------------------
operation proj=pipeline step init=epsg:26915 inv step proj=webmerc datum=WGS84
-------------------------------------------------------------------------------
tolerance 20 cm
accept 487147.594520173 4934316.46263998
expect -10370728.80 5552839.74
accept 487147.594520173 4934316.46263998 0
expect -10370728.80 5552839.74 0
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
Web Mercator test data from EPSG Guidance Note 7-2, p. 44.
-------------------------------------------------------------------------------
operation proj=webmerc
tolerance 1 cm
accept -100.33333333 24.46358028
expect -11169055.58 2810000.00
accept -100.33333333 24.38178694
expect -11169055.58 2800000.00
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
Test that +datum parameters are handled correctly in pipelines.
See #872 for details.
-------------------------------------------------------------------------------
operation +proj=pipeline
+step +proj=longlat +datum=GGRS87 +inv
+step +proj=longlat +datum=WGS84
-------------------------------------------------------------------------------
tolerance 20 cm
accept 23.7275 37.9838 0
expect 23.729194873180 37.986398897578 31.289740102
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
Test that +towgs84=0,0,0 parameter is handled as still implying cart
transformation
-------------------------------------------------------------------------------
operation +proj=pipeline
+step +proj=utm +zone=11 +ellps=clrk66 +towgs84=0,0,0 +inv
+step +proj=utm +zone=11 +datum=WGS84
-------------------------------------------------------------------------------
tolerance 20 cm
accept 440720 3751320 0
expect 440719.958709357 3751294.2109841 -4.44340920541435
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
Test that pipelines with unit mismatch between steps can't be constructed.
-------------------------------------------------------------------------------
operation +proj=pipeline
+step +proj=merc
+step +proj=merc
expect failure pjd_err_malformed_pipeline
operation +proj=pipeline
+step +proj=latlong
+step +proj=merc
+step +proj=helmert +x=200 +y=100
expect failure pjd_err_malformed_pipeline
operation +proj=pipeline
+step +proj=merc
+step +proj=unitconvert +xy_in=m +xy_out=km
accept 12 56
expect 1335.8339 7522.963
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
Test bugfix of https://github.com/OSGeo/proj.4/issues/1002
(do not interpolate nodata values)
-------------------------------------------------------------------------------
operation +proj=latlong +ellps=WGS84 +geoidgrids=tests/test_nodata.gtx
-------------------------------------------------------------------------------
ignore pjd_err_failed_to_load_grid
accept 4.05 52.1 0
expect 4.05 52.1 -10
-------------------------------------------------------------------------------
-------------------------------------------------------------------------------
Test bug fix of https://github.com/OSGeo/proj.4/issues/1025.
Using geocent in the new API with a custom ellipsoid should return coordinates
that correspond to that particular ellipsoid and not WGS84 as demonstrated in
the bug report.
-------------------------------------------------------------------------------
operation +proj=pipeline +step
+proj=longlat +a=3396190 +b=3376200 +inv +step
+proj=geocent +a=3396190 +b=3376200 +lon_0=0 +units=m
accept 0.0 0.0 0.0
expect 3396190.0 0.0 0.0
roundtrip 1
operation +proj=geocent +a=3396190 +b=3376200 +lon_0=0 +units=m
accept 0.0 0.0 0.00
expect 3396190.0 0.0 0.0
roundtrip 1
-------------------------------------------------------------------------------
Check that geocent and cart take into account to_meter (#1053)
-------------------------------------------------------------------------------
operation +proj=geocent +a=1000 +b=1000 +to_meter=1000
accept 90 0 0
expect 0 1 0
roundtrip 1
operation +proj=cart +a=1000 +b=1000 +to_meter=1000
accept 90 0 0
expect 0 1 0
roundtrip 1
-------------------------------------------------------------------------------
Check that vunits / vto_meter is honored
-------------------------------------------------------------------------------
operation +proj=longlat +a=1 +b=1 +vto_meter=1000
accept 0 0 1000
expect 0 0 1
roundtrip 1
operation +proj=longlat +a=1 +b=1 +vto_meter=1000 +geoc
accept 0 0 1000
expect 0 0 1
roundtrip 1
operation +proj=longlat +a=1 +b=1 +vunits=km
accept 0 0 1000
expect 0 0 1
roundtrip 1
operation +proj=merc +a=1 +b=1 +vto_meter=1000
accept 0 0 1000
expect 0 0 1
roundtrip 1
operation +proj=merc +a=1 +b=1 +vunits=km
accept 0 0 1000
expect 0 0 1
roundtrip 1
</gie>
|
