diff options
| author | Kristian Evers <kristianevers@gmail.com> | 2017-10-23 11:52:59 +0200 |
|---|---|---|
| committer | Kristian Evers <kristianevers@gmail.com> | 2017-10-29 14:15:36 +0100 |
| commit | 5877ffac2fd4c2f89deb5b9d962a1969192856c4 (patch) | |
| tree | 1c564825fbee3b4d7ddf3c112bbb353c2a5fdacc /test | |
| parent | 3a2bd267a67d41a461946a6f7b0a99262f47d8a7 (diff) | |
| download | PROJ-5877ffac2fd4c2f89deb5b9d962a1969192856c4.tar.gz PROJ-5877ffac2fd4c2f89deb5b9d962a1969192856c4.zip | |
Addition of 'deformation': Kinematic grid shifting.
Kinematic deformation models are used in some geodetic transformations.
This commit introduces the ability to do transformations involving a
gridded deformation/velocity model. For practical reasons a gridded deformation
model needs to be split into two seperate files, one for the horizontal
components and one for the vertical component. For this we use formats
already known to PROJ.4, namely the CTable/CTable2 and the GTX formats.
Grids are specified in the proj-string with +xy_grids and +z_grids.
Grid values are expected to be in m/year. The kinematic part of the
operation is controlled by the +t_epoch parameter, which is the central
epoch of the transformation. An observation epoch is also needed. It can
be specified either in the PJ_OBS input as the fourth element in the
coordinate, or in the proj-string with +t_obs. If +t_obs is present in
the proj-string it takes presedence over the value in the PJ_OBS
coordinate.
Diffstat (limited to 'test')
| -rw-r--r-- | test/gie/deformation.gie | 73 |
1 files changed, 73 insertions, 0 deletions
diff --git a/test/gie/deformation.gie b/test/gie/deformation.gie new file mode 100644 index 00000000..74a6b25d --- /dev/null +++ b/test/gie/deformation.gie @@ -0,0 +1,73 @@ +BEGIN +=============================================================================== +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. + +=============================================================================== + + +------------------------------------------------------------------------------- +Test using only horizontal grid and +tobs parameter +------------------------------------------------------------------------------- +OPERATION +proj=deformation +xy_grids=alaska +t_epoch=2016.0 +t_obs=2000.0 +------------------------------------------------------------------------------- +TOLERANCE 0.000001 m +ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 +EXPECT -3004295.5888766116 -1093474.1688513425 5500477.1338251457 +ROUNDTRIP 1000 + +------------------------------------------------------------------------------- +Test using only vertical grid and +tobs parameter +------------------------------------------------------------------------------- +OPERATION +proj=deformation +z_grids=egm96_15.gtx +t_epoch=2016.0 +t_obs=2000.0 +------------------------------------------------------------------------------- +TOLERANCE 0.000001 m +ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 +EXPECT -3004295.5882503074 -1093474.1690603832 5500234.008855661 +ROUNDTRIP 1000 + +------------------------------------------------------------------------------- +Test using both horizontal and vertical grids as well as the +tobs parameter +------------------------------------------------------------------------------- +OPERATION +proj=deformation +xy_grids=alaska +z_grids=egm96_15.gtx +t_epoch=2016.0 +t_obs=2000.0 +------------------------------------------------------------------------------- +TOLERANCE 0.000001 m +ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 +EXPECT -3004295.5888766116 -1093474.1688513425 5500234.008855661 +ROUNDTRIP 1000 + +------------------------------------------------------------------------------- +Test using only horizontal grid +------------------------------------------------------------------------------- +OPERATION +proj=deformation +xy_grids=alaska +t_epoch=2016.0 +------------------------------------------------------------------------------- +TOLERANCE 0.000001 m +ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 2000.0 +EXPECT -3004295.5888766116 -1093474.1688513425 5500477.1338251457 2000.0 +ROUNDTRIP 1000 + +------------------------------------------------------------------------------- +Test using only vertical grid +------------------------------------------------------------------------------- +OPERATION +proj=deformation +z_grids=egm96_15.gtx +t_epoch=2016.0 +------------------------------------------------------------------------------- +TOLERANCE 0.000001 m +ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 2000.0 +EXPECT -3004295.5882503074 -1093474.1690603832 5500234.008855661 2000.0 +ROUNDTRIP 1000 + +------------------------------------------------------------------------------- +Test using both horizontal and vertical grids +------------------------------------------------------------------------------- +OPERATION +proj=deformation +xy_grids=alaska +z_grids=egm96_15.gtx +t_epoch=2016.0 +------------------------------------------------------------------------------- +TOLERANCE 0.000001 m +ACCEPT -3004295.5882503074 -1093474.1690603832 5500477.1338251457 2000.0 +EXPECT -3004295.5888766116 -1093474.1688513425 5500234.008855661 2000.0 +ROUNDTRIP 1000 +END |