[GPlates-discuss] Reverse reconstruct a raster in GPlates

Michael Chin michael.chin at sydney.edu.au
Tue Jul 16 09:01:39 AEST 2019


Hi Aleece,

Here is the example code of using pyGPlates to reconstruct a raster. It is quite straightforward.

https://protect-au.mimecast.com/s/aoZICP7yOZtmmBnYszJP8E?domain=github.com

________________________________
From: GPlates-discuss <gplates-discuss-bounces at mailman.sydney.edu.au> on behalf of John Cannon <john.cannon at sydney.edu.au>
Sent: Monday, 15 July 2019 10:20 PM
To: GPlates general discussion mailing list
Cc: Nanfito, Aleece
Subject: Re: [GPlates-discuss] Reverse reconstruct a raster in GPlates


Hi Aleece,



Unfortunately it’s not possible in GPlates. Also the Assign Plate IDs dialog is only for vector data.



However, this post by Simon Williams talks about his pyGPlates code to reconstruct a raster from one time to another:

https://protect-au.mimecast.com/s/1wgBCQnzP0tqqnLGIP3U74?domain=mailman.sydney.edu.au<https://protect-au.mimecast.com/s/1wgBCQnzP0tqqnLGIP3U74?domain=mailman.sydney.edu.au>



My understanding (Simon can correct me if I’m wrong) is he essentially:

  1.  has a grid of points at the ‘to’ time,
  2.  assigns them plate IDs using the static polygons (reconstructed to the ‘to’ time),
  3.  reconstructs the points back to the ‘from’ time (using their newly assigned plate IDs),
  4.  samples the raster (which is assumed to represent raster data at the ‘from’ time) at those reconstructed point locations,
  5.  writes out the sampled values at the reconstructed points as a new raster at the ‘to’ time.



In your case the ‘from’ time is the time of your raster, and the ‘to’ time is any time you want to reconstruct it to. At least that’s my understanding ??  So that’s a good example of how to use pyGPlates to reconstruct a raster by reconstructing points and sampling rasters.



Regards,

John



From: GPlates-discuss <gplates-discuss-bounces at mailman.sydney.edu.au> On Behalf Of Nanfito, Aleece
Sent: Monday, 15 July 2019 6:53 AM
To: gplates-discuss at mailman.sydney.edu.au
Subject: [GPlates-discuss] Reverse reconstruct a raster in GPlates



Hello –

This is likely a novice question, but I am having problems in GPlates trying to reverse reconstruct a raster, in past time/position, to present day. I reconstruct the plate polygons to the past time represented in the raster, import the raster, and then connect the raster to the plate polygons.  The raster is automatically cut by the present day position of the plates. I’ve tried this exercise on rasters (geotiffs) with a global extent and regional extent. The raster and plate polygons are nicely aligned in the past position, before I connect the two. Is there a way to specify a past time/position for when the raster-to-polygon connection occurs?



I’ve also tried using the Assign Plate ID tool on the rasters. I can select all specifications within the assign plate id dialog box, but the process seems to immediately cancel/fail without providing an error.



I appreciate any input on whether this task is possible in GPlates & solutions, or any possible workaround in PyGPlates. Thanks so much in advance!



Best regards,

Aleece





[cid:image002.png at 01D53B59.F5EBD820]



Aleece Nanfito

Geologist

aleece.nanfito at bhpbilliton.com<mailto:aleece.nanfito at bhpbilliton.com>

1500 Post Oak Blvd.

Houston TX 77056 USA



bhpbilliton.com<https://protect-au.mimecast.com/s/xKDBCROAQotllwX9iPePDR?domain=bhpbilliton.com>





________________________________

This message and any attached files may contain information that is confidential and/or subject of legal privilege intended only for use by the intended recipient. If you are not the intended recipient or the person responsible for delivering the message to the intended recipient, be advised that you have received this message in error and that any dissemination, copying or use of this message or attachment is strictly forbidden, as is the disclosure of the information therein. If you have received this message in error please notify the sender immediately and delete the message.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.sydney.edu.au/pipermail/gplates-discuss/attachments/20190715/84f43e6b/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image002.png
Type: image/png
Size: 2555 bytes
Desc: image002.png
URL: <http://mailman.sydney.edu.au/pipermail/gplates-discuss/attachments/20190715/84f43e6b/attachment.png>


More information about the GPlates-discuss mailing list