blob: 848b9e89c6b4d4bc183ade6ef4e7a8bb9e006ede (
plain)
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
|
===============================================================================
Test for the deformation operation - Kinematic Gridshifting
For all the deformation tests the alaska and egm96_15.gtx grids are used even
though they are not parts of a deformation model, they are in the proper format
and for testing purposes it doesn't really matter all that much...
The input coordinate is located at lon=60, lam=-160 - somewhere in Alaska.
===============================================================================
<gie>
-------------------------------------------------------------------------------
Test with an extract of nkgrf03vel_realigned with ctable2+gtx
-------------------------------------------------------------------------------
operation +proj=pipeline
+step +proj=cart +ellps=GRS80
+step +proj=deformation
+xy_grids=tests/nkgrf03vel_realigned_xy_extract.ct2
+z_grids=tests/nkgrf03vel_realigned_z_extract.gtx +ellps=GRS80 +dt=1
+step +proj=cart +ellps=GRS80 +inv
-------------------------------------------------------------------------------
tolerance 0.05 mm
accept 21.5 63 0
expect 21.5000000049 62.9999999937 0.0083
roundtrip 5
-------------------------------------------------------------------------------
Test with an extract of nkgrf03vel_realigned with GeoTIFF
-------------------------------------------------------------------------------
operation +proj=pipeline
+step +proj=cart +ellps=GRS80
+step +proj=deformation
+grids=tests/nkgrf03vel_realigned_extract.tif +ellps=GRS80 +dt=1
+step +proj=cart +ellps=GRS80 +inv
-------------------------------------------------------------------------------
tolerance 0.05 mm
accept 21.5 63 0
expect 21.5000000049 62.9999999937 0.0083
roundtrip 5
-------------------------------------------------------------------------------
Test the +dt parameter
-------------------------------------------------------------------------------
operation +proj=deformation +xy_grids=alaska +z_grids=egm96_15.gtx
+ellps=GRS80 +dt=16.0 # 2016.0 - 2000.0
-------------------------------------------------------------------------------
tolerance 0.1 mm
ignore pjd_err_failed_to_load_grid
accept -3004295.5882503074 -1093474.1690603832 5500477.1338251457
expect -3004295.7025 -1093474.2106 5500477.3444
roundtrip 5
# Test that errors are reported for coordinates outside the grid.
# Here we test 120W 40N which is well outside the alaska grid.
accept -2446353.8001 -4237209.0750 4077985.572
expect failure errno grid_area
accept -2446353.8001 -4237209.0750 4077985.572
expect failure errno grid_area
-------------------------------------------------------------------------------
Test using both horizontal and vertical grids
-------------------------------------------------------------------------------
operation +proj=deformation
+xy_grids=alaska +z_grids=egm96_15.gtx +t_epoch=2016.0 +ellps=GRS80
-------------------------------------------------------------------------------
tolerance 0.1 mm
direction inverse
ignore pjd_err_failed_to_load_grid
accept -3004295.5882503074 -1093474.1690603832 5500477.1338251457 2000.0
expect -3004295.7025 -1093474.2106 5500477.3444 2000.0
roundtrip 5
-------------------------------------------------------------------------------
operation proj=deformation xy_grids=alaska +dt=1.0
ellps=GRS80
expect failure pjd_err_no_args
operation proj=deformation z_grids=egm96_15.gtx +dt=1.0
ellps=GRS80
expect failure pjd_err_no_args
operation proj=deformation xy_grids=nonexisting z_grids=egm96_15.gtx
+dt=1.0 ellps=GRS80
expect failure pjd_err_failed_to_load_grid
operation proj=deformation xy_grids=alaska z_grids=nonexisting
+dt=1.0 ellps=GRS80
expect failure pjd_err_failed_to_load_grid
operation proj=deformation xy_grids=alaska z_grids=nonexisting ellps=GRS80
expect failure pjd_err_missing_args
-------------------------------------------------------------------------------
operation +proj=vgridshift +grids=egm96_15.gtx +t_epoch=2010.0 +t_final=2018.0
-------------------------------------------------------------------------------
tolerance 0.1 mm
ignore pjd_err_failed_to_load_grid
accept 12 56 0.0 2000.0
expect 12 56 -36.5966 2000.0
roundtrip 100
accept 12 56 0.0 2011.0
expect 12 56 0.0 2011.0
roundtrip 100
accept 12 56 0.0 2019.0
expect 12 56 0.0 2019.0
roundtrip 100
accept 12 56 0.0
expect 12 56 -36.5966
roundtrip 100
-------------------------------------------------------------------------------
operation +proj=vgridshift +grids=egm96_15.gtx +t_epoch=2010.0 +t_final=now
-------------------------------------------------------------------------------
tolerance 0.1 mm
ignore pjd_err_failed_to_load_grid
accept 12 56 0.0 2000.0
expect 12 56 -36.5966 2000.0
roundtrip 100
accept 12 56 0.0 2011.0
expect 12 56 0.0 2011.0
roundtrip 1000
accept 12 56 0.0 3011.0
expect 12 56 0.0 3011.0
roundtrip 100
-------------------------------------------------------------------------------
operation +proj=hgridshift +grids=alaska +t_epoch=2010.0 +t_final=2018.0
-------------------------------------------------------------------------------
tolerance 0.1 mm
ignore pjd_err_failed_to_load_grid
accept -147.0 64.0 0.0 2000.0
expect -147.0023233121 63.9995792119 0.0 2000.0
roundtrip 100
accept -147.0 64.0 0.0 2011.0
expect -147.0 64.0 0.0 2011.0
roundtrip 100
accept -147.0 64.0 0.0 2011.0
expect -147.0 64.0 0.0 2020.0
roundtrip 100
-------------------------------------------------------------------------------
operation +proj=hgridshift +grids=alaska +t_epoch=2010.0 +t_final=now
-------------------------------------------------------------------------------
tolerance 0.1 mm
ignore pjd_err_failed_to_load_grid
accept -147.0 64.0 0.0 2000.0
expect -147.0023233121 63.9995792119 0.0 2000.0
roundtrip 100
accept -147.0 64.0 0.0 2011.0
expect -147.0 64.0 0.0 2011.0
roundtrip 100
accept -147.0 64.0 0.0 3011.0
expect -147.0 64.0 0.0 3011.0
roundtrip 100
</gie>
|
