|
3 | 3 | from geodepy.transform import ( |
4 | 4 | conform7, |
5 | 5 | conform14, |
| 6 | + plate_motion_transformation, |
6 | 7 | atrf2014_to_gda2020, |
7 | 8 | transform_gda2020_to_atrf2014, |
8 | 9 | transform_atrf2014_to_gda2020, |
@@ -46,6 +47,22 @@ def test_conform14(self): |
46 | 47 | assert abs(alic_gda2020_comp[1] - alic_gda2020[1]) < 5e-5 |
47 | 48 | assert abs(alic_gda2020_comp[2] - alic_gda2020[2]) < 5e-5 |
48 | 49 |
|
| 50 | + def test_plate_motion_transformation(self): |
| 51 | + alic_atrf14at2018 = (-4052052.6588, 4212835.9938, -2545104.6946) |
| 52 | + alic_atrf14at2025 = (-4052052.9336, 4212835.9577, -2545104.3168) |
| 53 | + alic_atrf14at2025_comp = plate_motion_transformation( |
| 54 | + *alic_atrf14at2018, date(2018,1,1), date(2025,1,1),itrf2014_to_gda2020 |
| 55 | + ) |
| 56 | + alic_atrf14at2018_comp = plate_motion_transformation( |
| 57 | + *alic_atrf14at2025, date(2025,1,1), date(2018,1,1), itrf2014_to_gda2020 |
| 58 | + ) |
| 59 | + assert abs(alic_atrf14at2018_comp[0] - alic_atrf14at2018[0]) < 5e-5 |
| 60 | + assert abs(alic_atrf14at2018_comp[1] - alic_atrf14at2018[1]) < 5e-5 |
| 61 | + assert abs(alic_atrf14at2018_comp[2] - alic_atrf14at2018[2]) < 5e-5 |
| 62 | + assert abs(alic_atrf14at2025_comp[0] - alic_atrf14at2025[0]) < 5e-5 |
| 63 | + assert abs(alic_atrf14at2025_comp[1] - alic_atrf14at2025[1]) < 5e-5 |
| 64 | + assert abs(alic_atrf14at2025_comp[2] - alic_atrf14at2025[2]) < 5e-5 |
| 65 | + |
49 | 66 | def test_transform_atrf2014_to_gda2020(self): |
50 | 67 | alic_gda2020 = (-4052052.7373, 4212835.9835, -2545104.5867) |
51 | 68 | alic_itrf14at2018 = (-4052052.6588, 4212835.9938, -2545104.6946) |
|
0 commit comments