GPlates / gplately

GPlately is a Python package to interrogate tectonic plate reconstructions.
https://gplates.github.io/gplately/
GNU General Public License v2.0
55 stars 11 forks source link

Use multipoint coverage to save age gridding debug data #208

Closed michaelchin closed 2 days ago

michaelchin commented 4 days ago

Instead of saving as individual points you can save as a multipoint coverage. Basically, each point gets one or more scalar values – such as age, spreading rate, direction, obliquity, etc.

The details is John's email. Find the email in Michael Chin's email inbox.

jcannon-gplates commented 4 days ago

Thanks, here's the email...

Instead of saving as individual points you can save as a multipoint coverage. Basically, each point gets one or more scalar values – such as age, spreading rate, direction, obliquity, etc. But it’s all compactly stored in the GPML file – so it’s much smaller and loads faster. In isopolate, you can see how it’s done here (where this function shows how the age, spreading rate, direction, obliquity, etc, are added). This ultimately hinges around pygplates.Feature.set_geometry() (see the bottom example in that function).

The following screenshot shows some old isopolate output (it’s actually lines instead of multipoints, but it’s the same idea, each point gets a different colour). And it’s easy to switch between the different output types (eg, age, spreading obliquity, etc) and colour each output type with a different palette. This is in grey layers (that aren’t commonly seen in GPlates).

Also, if you could set the feature time-period on each GPML to [time+0.5, time-0.5] that would be great! Because it means we can load all the GPMLs at once and animate the time slider - and only one GPML will be visible at a time. That allows us to more easily animate the seed points to see where they go over time. Currently, when you animate the slider you also have to disable some layers at later times so they don’t overlap the current time.

image