search

kmaterna / Strain_2D

Strain rate modeling from GNSS velocity fields
MIT License
44 stars 19 forks source link

Add a Python-based VISR implementation #40

Open tobiscode opened 1 month ago

tobiscode commented 1 month ago

Just a meta issue tracking my fork and then subsequent PR to merge a Python-based, simple VISR method to Strain_2D.

tobiscode commented 1 month ago

@kmaterna I had a quick question about coordinate systems while porting my code. Do I understand correctly that the Delaunay and other Python-based method assume that lon/lat is a cartesian grid? Because in my implementation I convert everything to the closest UTM representation and then do the calculations. Since pyproj isn't a requirement (or a similar package), I wouldn't be able to do that.

Also feel free to assign me to this issue.

jlmaurer commented 2 weeks ago

@tobiscode we implement a regular grid in lat/lon space, which is a distorted grid in cartesian space. To calculate grid cell area we use a cosine correction (see e.g. here).