[GPlates-discuss] Generate Crustal Thickness Points in pyGPlates

John Cannon john.cannon at sydney.edu.au
Fri Sep 13 21:03:24 AEST 2019


Hi Thomas,

It should be available in the next pyGPlates release (sometime in the first half of next year). That release is aimed to bring deformation on par with current GPlates, so deformation-related things like generating crustal thickness points should get included.

The crustal point locations are based on subdivision of a Rhombic Triacontahedron, with optional random perturbation of the positions. It's not perfectly uniform on the sphere but quite close. The source code file to look at is "src/maths/GeneratePoints.h": https://protect-au.mimecast.com/s/rhN2CJyp0qhJJLwoIVfcfG?domain=sourceforge.net . There are other algorithms like HEALPix (Simon Williams uses this) which has a library (but not on Windows).

There's two separate aspects here, (1) the deforming topologies and their quantities like strain rate, and (2) reconstructing points in those regions. It's the latter aspect that you are probably interested in, and this is covered briefly in "GPlates: Building a Virtual Earth Through Deep Time" in the section "Reconstructing feature using topologies":
https://protect-au.mimecast.com/s/g6b2CK1qJZtooxWkcv-Elr?domain=doi.org
https://protect-au.mimecast.com/s/Eo8uCL7rK8tWWM92IPkOJu?domain=researchgate.net


But it might be easier to wait until it's exposed via pyGPlates since moving these points through a deforming region is done by spatially interpolating triangulations of the deforming region as it evolves through time, and there's a few implementation gotchas to look out for.

Regards,
John

From: GPlates-discuss <gplates-discuss-bounces at mailman.sydney.edu.au> On Behalf Of Linden, T.J.M. van der (Thomas)
Sent: Friday, 13 September 2019 7:31 PM
To: gplates-discuss at mailman.sydney.edu.au
Subject: [GPlates-discuss] Generate Crustal Thickness Points in pyGPlates

Hi everyone,

Is the "Generate Crustal Thickness Points" from GPlates also available in pyGPlates or will it be in the near future? I would really like to use it in scripts.
Otherwise is there a publication on the methods behind this feature? I looked into coding something similar myself, but ran into problems making a proper grid and with assigning plate-IDs.

Hope you can help me!
Best from Berlin,
Thomas
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.sydney.edu.au/pipermail/gplates-discuss/attachments/20190913/0c00262c/attachment-0001.html>


More information about the GPlates-discuss mailing list