|
1 | | -from __future__ import annotations |
| 1 | +#!/usr/bin/env python3 |
| 2 | + |
| 3 | +""" |
| 4 | +Geoscience Australia - Python Geodesy Package |
| 5 | +Point in Polygon Module |
| 6 | +""" |
2 | 7 |
|
3 | 8 | import re |
4 | 9 | import math |
|
12 | 17 |
|
13 | 18 | from pyproj import CRS, Transformer |
14 | 19 |
|
15 | | -import geodepy.convert |
16 | | -from geodepy.constants import Transformation, TransformationSD |
17 | | -import geodepy.transform |
| 20 | +import convert |
| 21 | +from constants import Transformation, TransformationSD |
| 22 | +import transform |
18 | 23 |
|
19 | 24 | PLATE_ID_RE = re.compile(r"^[A-Z]{2}$") |
20 | 25 |
|
@@ -260,7 +265,7 @@ def plate_from_xyz(x, y, z, plate_index): |
260 | 265 | :return: plate ID |
261 | 266 | :rtype: Array of strings or None |
262 | 267 | """ |
263 | | - lat, lon, h = geodepy.convert.xyz2llh(x, y, z) |
| 268 | + lat, lon, h = convert.xyz2llh(x, y, z) |
264 | 269 | return plate_from_ll(lat, lon, plate_index) |
265 | 270 |
|
266 | 271 | @dataclass(frozen=True) |
@@ -496,6 +501,6 @@ def universal_plate_motion_transformation(x ,y, z, |
496 | 501 | transformation = plate_transformation(plate_id, poles_file, ref_epoch=ref_epoch) |
497 | 502 |
|
498 | 503 | # Complete transformation between given dates |
499 | | - xtrans, ytrans, ztrans, vcv = geodepy.transform.plate_motion_transformation(x, y, z, from_epoch, to_epoch, transformation, vcv) |
| 504 | + xtrans, ytrans, ztrans, vcv = transform.plate_motion_transformation(x, y, z, from_epoch, to_epoch, transformation, vcv) |
500 | 505 |
|
501 | 506 | return xtrans, ytrans, ztrans, vcv |
0 commit comments