From bruce.eglington at usask.ca Sun Apr 8 02:44:52 2018 From: bruce.eglington at usask.ca (Eglington, Bruce) Date: Sat, 7 Apr 2018 16:44:52 +0000 Subject: [GPlates-discuss] Need to omit white space in rasters Message-ID: <5997d6770c2d4ae1ae504e758bc32a6c@Mail02.usask.ca> Hi Is there a way to turn off or omit parts of raster images which either have no data or a background colour? In general, these appear bright white in GPlates. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From john.cannon at sydney.edu.au Sun Apr 8 13:46:08 2018 From: john.cannon at sydney.edu.au (John Cannon) Date: Sun, 8 Apr 2018 03:46:08 +0000 Subject: [GPlates-discuss] Need to omit white space in rasters In-Reply-To: <5997d6770c2d4ae1ae504e758bc32a6c@Mail02.usask.ca> References: <5997d6770c2d4ae1ae504e758bc32a6c@Mail02.usask.ca> Message-ID: <1E57973145EDFC41B1D69DE11923922801B95E271E@ex-mbx-pro-05> Hi Bruce, Yes that should be possible in the following situations... Numerical rasters (eg, NetCDF, GeoTIFF) contain a single data value per pixel, and so any pixels equal to the no-data value specified for the raster (eg, NaN) should get rendered as transparent (ie, not drawn). An example is the age grid in the sample data (which uses no-data=NaN). Colour rasters that support transparency (eg, PNG, BMP, TIFF) contain Red/Green/Blue/Alpha channels per pixel, and so any pixels with alpha equal to zero should get rendered as transparent (ie, not drawn). Also note that alpha values between 0 and 1.0 (where 1.0 is equivalent to 255 for 8-bit colour channels) get rendered as translucent (ie, underneath is partially visible). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Sunday, 8 April 2018 2:45 AM To: GPlates general discussion mailing list (gplates-discuss at mail.usyd.edu.au) Subject: [GPlates-discuss] Need to omit white space in rasters Hi Is there a way to turn off or omit parts of raster images which either have no data or a background colour? In general, these appear bright white in GPlates. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From bruce.eglington at usask.ca Sun Apr 8 16:12:50 2018 From: bruce.eglington at usask.ca (Eglington, Bruce) Date: Sun, 8 Apr 2018 06:12:50 +0000 Subject: [GPlates-discuss] Need to omit white space in rasters In-Reply-To: <1E57973145EDFC41B1D69DE11923922801B95E271E@ex-mbx-pro-05> References: <5997d6770c2d4ae1ae504e758bc32a6c@Mail02.usask.ca>, <1E57973145EDFC41B1D69DE11923922801B95E271E@ex-mbx-pro-05> Message-ID: <1523167970316.27405@usask.ca> Hi John Any suggestion for how I might set alpha values for specific pixels in a colour raster? In Arcmap I can specify that some pixel value, normally a background value of 255 or 65535, should be shown as transparent but I don?t know of a way to save this as an alpha value Bruce ________________________________ From: GPlates-discuss on behalf of John Cannon Sent: Saturday, April 7, 2018 9:46:08 PM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] Need to omit white space in rasters Hi Bruce, Yes that should be possible in the following situations? Numerical rasters (eg, NetCDF, GeoTIFF) contain a single data value per pixel, and so any pixels equal to the no-data value specified for the raster (eg, NaN) should get rendered as transparent (ie, not drawn). An example is the age grid in the sample data (which uses no-data=NaN). Colour rasters that support transparency (eg, PNG, BMP, TIFF) contain Red/Green/Blue/Alpha channels per pixel, and so any pixels with alpha equal to zero should get rendered as transparent (ie, not drawn). Also note that alpha values between 0 and 1.0 (where 1.0 is equivalent to 255 for 8-bit colour channels) get rendered as translucent (ie, underneath is partially visible). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Sunday, 8 April 2018 2:45 AM To: GPlates general discussion mailing list (gplates-discuss at mail.usyd.edu.au) Subject: [GPlates-discuss] Need to omit white space in rasters Hi Is there a way to turn off or omit parts of raster images which either have no data or a background colour? In general, these appear bright white in GPlates. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From john.cannon at sydney.edu.au Sun Apr 8 20:47:03 2018 From: john.cannon at sydney.edu.au (John Cannon) Date: Sun, 8 Apr 2018 10:47:03 +0000 Subject: [GPlates-discuss] Need to omit white space in rasters In-Reply-To: <1523167970316.27405@usask.ca> References: <5997d6770c2d4ae1ae504e758bc32a6c@Mail02.usask.ca>, <1E57973145EDFC41B1D69DE11923922801B95E271E@ex-mbx-pro-05> <1523167970316.27405@usask.ca> Message-ID: <1E57973145EDFC41B1D69DE11923922801B95E2804@ex-mbx-pro-05> Hi Bruce, I don't have much experience with this, but maybe you can set the NoData value to the colour you want to be transparent (eg, white), and then export as a numerical raster type (such as GeoTIFF). Although you might need to turn compression off to avoid changing the pixel values during export (such that they no longer equal the NoData value) if the compression is lossy. This link talks a bit about NoData (in Arcmap raster datasets)... https://protect-au.mimecast.com/s/XMpsCE8kz9tOo8pZCNAwpd?domain=desktop.arcgis.com ...maybe you could try setting the NoData value in the raster properties dialog (see bottom of that page). Note that Red/Green/Blue data is usually represented as three bands in a numerical (data) raster so you might have to set a NoData value for each band. Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Sunday, 8 April 2018 4:13 PM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] Need to omit white space in rasters Hi John Any suggestion for how I might set alpha values for specific pixels in a colour raster? In Arcmap I can specify that some pixel value, normally a background value of 255 or 65535, should be shown as transparent but I don't know of a way to save this as an alpha value Bruce ________________________________ From: GPlates-discuss > on behalf of John Cannon > Sent: Saturday, April 7, 2018 9:46:08 PM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] Need to omit white space in rasters Hi Bruce, Yes that should be possible in the following situations... Numerical rasters (eg, NetCDF, GeoTIFF) contain a single data value per pixel, and so any pixels equal to the no-data value specified for the raster (eg, NaN) should get rendered as transparent (ie, not drawn). An example is the age grid in the sample data (which uses no-data=NaN). Colour rasters that support transparency (eg, PNG, BMP, TIFF) contain Red/Green/Blue/Alpha channels per pixel, and so any pixels with alpha equal to zero should get rendered as transparent (ie, not drawn). Also note that alpha values between 0 and 1.0 (where 1.0 is equivalent to 255 for 8-bit colour channels) get rendered as translucent (ie, underneath is partially visible). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Sunday, 8 April 2018 2:45 AM To: GPlates general discussion mailing list (gplates-discuss at mail.usyd.edu.au) Subject: [GPlates-discuss] Need to omit white space in rasters Hi Is there a way to turn off or omit parts of raster images which either have no data or a background colour? In general, these appear bright white in GPlates. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From paul.angrand at get.omp.eu Sat Apr 28 00:49:39 2018 From: paul.angrand at get.omp.eu (Paul Angrand) Date: Fri, 27 Apr 2018 16:49:39 +0200 Subject: [GPlates-discuss] Exporting strain accumulation Message-ID: <43990D56-D637-4C7C-9CE2-A1AE31599ADA@get.omp.eu> Dear GPlates Community, I am working with GPlates 2.0 for plate reconstruction and deformation. I would like to export the strain accumulation that GPlates calculates for a resolved topological network (after analyzing the crustal thickness), and I do not know how to do that. There is no strain accumulation within the export tool. The dilatation strain rate and the crustal thickness are successfully exported. If you have any idea about how to do this. Thanks for any help, Paul Angrand - - Paul ANGRAND, Ph.D Post-doc in Structural Geology and Geodynamics at GET-Observatoire Midi-Pyr?n?es GET-OMP UMR 5563, 14 av Edouard Belin 31400 Toulouse, France +33 (0) 5 61 33 30 27 +33 (0) 6 95 33 94 29 -------------- next part -------------- An HTML attachment was scrubbed... URL: From john.cannon at sydney.edu.au Sat Apr 28 01:01:56 2018 From: john.cannon at sydney.edu.au (John Cannon) Date: Fri, 27 Apr 2018 15:01:56 +0000 Subject: [GPlates-discuss] Exporting strain accumulation In-Reply-To: <43990D56-D637-4C7C-9CE2-A1AE31599ADA@get.omp.eu> References: <43990D56-D637-4C7C-9CE2-A1AE31599ADA@get.omp.eu> Message-ID: <1E57973145EDFC41B1D69DE11923922801B95E871F@ex-mbx-pro-05> Hi Paul, That will be available in GPlates 2.1 (soon), which will export the strain ellipse (of the 2D strain tensor) as the two principal strain magnitudes and the angle of the larger principal (ie, alignment relative to North). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Paul Angrand Sent: Saturday, 28 April 2018 12:50 AM To: gplates-discuss at mailman.sydney.edu.au Subject: [GPlates-discuss] Exporting strain accumulation Dear GPlates Community, I am working with GPlates 2.0 for plate reconstruction and deformation. I would like to export the strain accumulation that GPlates calculates for a resolved topological network (after analyzing the crustal thickness), and I do not know how to do that. There is no strain accumulation within the export tool. The dilatation strain rate and the crustal thickness are successfully exported. If you have any idea about how to do this. Thanks for any help, Paul Angrand - - Paul ANGRAND, Ph.D Post-doc in Structural Geology and Geodynamics at GET-Observatoire Midi-Pyr?n?es GET-OMP UMR 5563, 14 av Edouard Belin 31400 Toulouse, France +33 (0) 5 61 33 30 27 +33 (0) 6 95 33 94 29 -------------- next part -------------- An HTML attachment was scrubbed... URL: From bruce.eglington at usask.ca Sat Apr 28 06:57:49 2018 From: bruce.eglington at usask.ca (Eglington, Bruce) Date: Fri, 27 Apr 2018 20:57:49 +0000 Subject: [GPlates-discuss] GPlates and linking to databases In-Reply-To: <4aebc4cf9aea4761b6c8df5e08e2da64@shell.com> References: <56cd4a91fb324bb5b2191b2f73071eca@Mail01.usask.ca> <4aebc4cf9aea4761b6c8df5e08e2da64@shell.com> Message-ID: <0c26e5d403804d5cbfef0981b73f8674@Mail01.usask.ca> Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID's and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you'd count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the "lossy" shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn't like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From sabin.zahirovic at sydney.edu.au Sat Apr 28 07:56:35 2018 From: sabin.zahirovic at sydney.edu.au (Sabin Zahirovic) Date: Fri, 27 Apr 2018 21:56:35 +0000 Subject: [GPlates-discuss] GPlates and linking to databases Message-ID: Hi Bruce, If you digitize a feature to be a MidOceanicRidge, you can specify the Reconstruction Method as ?Half Stage Rotation? with a Left and Right PlateID. The ridge will then move with a half stage rotation, and will be equally spaced. Regards, Sabin -- DR SABIN ZAHIROVIC | Postdoctoral Research Associate School of Geosciences | Faculty of Science THE UNIVERSITY OF SYDNEY Rm 403, Madsen Building F09 | The University of Sydney | NSW | 2006 M +61 416 775 589 P +61 2 9351 3625 E sabin.zahirovic at sydney.edu.au | W https://protect-au.mimecast.com/s/NhyqCnxyErCnQRW7u9lM_1?domain=earthbyte.org | R http://www.researchgate.net/profile/Sabin_Zahirovic F https://protect-au.mimecast.com/s/2Cm5Cp8AJQt5PRwntYFQ8M?domain=facebook.com | T https://protect-au.mimecast.com/s/v310Cq7BKYt9GM28HElESM?domain=twitter.com CRICOS 00026A This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments. From: GPlates-discuss on behalf of "Eglington, Bruce" Reply-To: GPlates general discussion mailing list Date: Saturday, 28 April 2018 at 6:59 am To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID?s and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you?d count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the ?lossy? shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn?t like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From bruce.eglington at usask.ca Sat Apr 28 09:11:16 2018 From: bruce.eglington at usask.ca (Eglington, Bruce) Date: Fri, 27 Apr 2018 23:11:16 +0000 Subject: [GPlates-discuss] GPlates and linking to databases In-Reply-To: References: Message-ID: Hi Sabin Thanks, I am probably not specifying the half-plate part. Cheers Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Sabin Zahirovic Sent: Friday, April 27, 2018 15:57 To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you digitize a feature to be a MidOceanicRidge, you can specify the Reconstruction Method as ?Half Stage Rotation? with a Left and Right PlateID. The ridge will then move with a half stage rotation, and will be equally spaced. Regards, Sabin -- DR SABIN ZAHIROVIC | Postdoctoral Research Associate School of Geosciences | Faculty of Science THE UNIVERSITY OF SYDNEY Rm 403, Madsen Building F09 | The University of Sydney | NSW | 2006 M +61 416 775 589 P +61 2 9351 3625 E sabin.zahirovic at sydney.edu.au | W https://protect-au.mimecast.com/s/4zx1C81Zj6t98XqAsnzpQQ?domain=earthbyte.org | R http://www.researchgate.net/profile/Sabin_Zahirovic F https://protect-au.mimecast.com/s/Ua7yC0YZWVFj0kB9h2ftJd?domain=facebook.com | T https://protect-au.mimecast.com/s/jNuhCgZowLH85qyjI39PRT?domain=twitter.com CRICOS 00026A This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments. From: GPlates-discuss > on behalf of "Eglington, Bruce" > Reply-To: GPlates general discussion mailing list > Date: Saturday, 28 April 2018 at 6:59 am To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID?s and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you?d count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the ?lossy? shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn?t like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From bruce.eglington at usask.ca Sat Apr 28 10:53:12 2018 From: bruce.eglington at usask.ca (Eglington, Bruce) Date: Sat, 28 Apr 2018 00:53:12 +0000 Subject: [GPlates-discuss] GPlates and linking to databases In-Reply-To: References: Message-ID: Hi Sabin I need to be a bit more certain about what is needed so please excuse what will appear very evident as soon as I get things working. As I understand you, for a given OceanRidge, I need to specify the half-stage left plate plateid value and the half-stage right plate plate id. When I do this do I specify the normal PlateID field as none and the conjugate plateid field as none? Presumably I ought also to set the featuretype as having an appropriate GPGIM_TYPE value for ocean ridges (?gpml:MidOceanRidge? as I remember)? Do I need to have fields with something in my shapefile for reconstruction method and for spreading asymmetry? Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Sabin Zahirovic Sent: Friday, April 27, 2018 15:57 To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you digitize a feature to be a MidOceanicRidge, you can specify the Reconstruction Method as ?Half Stage Rotation? with a Left and Right PlateID. The ridge will then move with a half stage rotation, and will be equally spaced. Regards, Sabin -- DR SABIN ZAHIROVIC | Postdoctoral Research Associate School of Geosciences | Faculty of Science THE UNIVERSITY OF SYDNEY Rm 403, Madsen Building F09 | The University of Sydney | NSW | 2006 M +61 416 775 589 P +61 2 9351 3625 E sabin.zahirovic at sydney.edu.au | W https://protect-au.mimecast.com/s/5lbhCr8DLRtx31JDi7YUzO?domain=earthbyte.org | R http://www.researchgate.net/profile/Sabin_Zahirovic F https://protect-au.mimecast.com/s/H0D0CwVLQmioZArRH93Uzt?domain=facebook.com | T https://protect-au.mimecast.com/s/YWFACxnMRvtZq9NQHw8Sgo?domain=twitter.com CRICOS 00026A This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments. From: GPlates-discuss > on behalf of "Eglington, Bruce" > Reply-To: GPlates general discussion mailing list > Date: Saturday, 28 April 2018 at 6:59 am To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID?s and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you?d count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the ?lossy? shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn?t like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From john.cannon at sydney.edu.au Sat Apr 28 16:03:09 2018 From: john.cannon at sydney.edu.au (John Cannon) Date: Sat, 28 Apr 2018 06:03:09 +0000 Subject: [GPlates-discuss] GPlates and linking to databases In-Reply-To: References: Message-ID: <1E57973145EDFC41B1D69DE11923922801B95E8BDF@ex-mbx-pro-05> Hi Bruce, If you?re digitizing a feature in GPlates then you don?t need to specify the Plate ID and Conjugate ID fields. When you switch to Half Stage Rotation (in Create Feature dialog) those fields disappear, and the Left Plate ID and Right Plate ID fields appear. And the feature type is taken care of when you selected MidOceanRidge. Saving to Shapefile will then map the feature type and spreading asymmetry to the GPGIM_TYPE and SPREAD_ASY shapefile attributes by default. However if you?re creating a feature outside GPlates, and without using pyGPlates, then you can add the GPGIM_TYPE shapefile attribute as ?gpml:MidOceanRidge? (as you noted) and the SPREAD_ASY attribute as a value in the range [-1,1] where 0 is symmetric spreading. And add the left/right plate IDs as L_PLATE and R_PLATE. Note that when GPlates loads your shapefile it will also load plate/conjugate IDs of zero by default, but they won?t get used for a half stage rotation. On a related note, the current public pyGPlates (revision 12) will not work with mid-ocean ridges created in GPlates 2.0 (ie, pyGPlates will not reconstruct them properly). To solve this there will be a simultaneous pyGPlates release alongside GPlates 2.1 (soon). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Saturday, 28 April 2018 10:53 AM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Sabin I need to be a bit more certain about what is needed so please excuse what will appear very evident as soon as I get things working. As I understand you, for a given OceanRidge, I need to specify the half-stage left plate plateid value and the half-stage right plate plate id. When I do this do I specify the normal PlateID field as none and the conjugate plateid field as none? Presumably I ought also to set the featuretype as having an appropriate GPGIM_TYPE value for ocean ridges (?gpml:MidOceanRidge? as I remember)? Do I need to have fields with something in my shapefile for reconstruction method and for spreading asymmetry? Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Sabin Zahirovic Sent: Friday, April 27, 2018 15:57 To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you digitize a feature to be a MidOceanicRidge, you can specify the Reconstruction Method as ?Half Stage Rotation? with a Left and Right PlateID. The ridge will then move with a half stage rotation, and will be equally spaced. Regards, Sabin -- DR SABIN ZAHIROVIC | Postdoctoral Research Associate School of Geosciences | Faculty of Science THE UNIVERSITY OF SYDNEY Rm 403, Madsen Building F09 | The University of Sydney | NSW | 2006 M +61 416 775 589 P +61 2 9351 3625 E sabin.zahirovic at sydney.edu.au | W https://protect-au.mimecast.com/s/o0dwCq7BKYt96030cZIoHQ?domain=earthbyte.org | R http://www.researchgate.net/profile/Sabin_Zahirovic F https://protect-au.mimecast.com/s/eilBCvl0PoCv0N5NhzZF1N?domain=facebook.com | T https://protect-au.mimecast.com/s/mQQ5CwVLQmioZ8m8hKsNak?domain=twitter.com CRICOS 00026A This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments. From: GPlates-discuss > on behalf of "Eglington, Bruce" > Reply-To: GPlates general discussion mailing list > Date: Saturday, 28 April 2018 at 6:59 am To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID?s and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you?d count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the ?lossy? shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn?t like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From john.cannon at sydney.edu.au Sat Apr 28 17:54:28 2018 From: john.cannon at sydney.edu.au (John Cannon) Date: Sat, 28 Apr 2018 07:54:28 +0000 Subject: [GPlates-discuss] GPlates and linking to databases In-Reply-To: <1E57973145EDFC41B1D69DE11923922801B95E8BDF@ex-mbx-pro-05> References: <1E57973145EDFC41B1D69DE11923922801B95E8BDF@ex-mbx-pro-05> Message-ID: <1E57973145EDFC41B1D69DE11923922801B95E8C61@ex-mbx-pro-05> ?I forgot to mention the reconstruction method maps to the RECON_METH shapefile attribute, and can be set to HalfStageRotationVersion3 (or HalfStageRotationVersion2 or HalfStageRotation). However note the issues below (better to use GPlates or pyGPlates to do all this)? Version 3 of the half-stage rotation was introduced in GPlates 2.0 (and, as mentioned in my last post, current pyGPlates only supports version 2 and below). With version 3 the spreading starts at the digitization time (ie, the reconstruction time at which the mid-ocean ridge was digitized in GPlates) such that subsequent adjustments to the spreading asymmetry no longer change the position of the ridge at the digitization time. With version 2, the ridge position would offset quite substantially making it unfeasible to adjust the asymmetry. So now the digitization time needed to be recorded, and this introduced the new property ?gpml:geometryImportTime? (also used in the new ability to reconstruct geometries using deforming and non-deforming topologies in GPlates 2.0). Oh, I just realized that gpml:geometryImportTime is not mapped in GPlates ! Which means it gets lost on saving to, or loading from, Shapefile and defaults to 0Ma. And this puts (version 3) mid-ocean ridges in the wrong location in GPlates 2.0 (when loaded from Shapefile). I?ll fix that for GPlates 2.1 (by mapping to a new attribute, say IMPORT_AGE). Actually, even when that is fixed, it will still be quite difficult to create a mid-ocean ridge without using GPlates or pyGPlates. This is because you need to take your mid-ocean ridge geometry (at say 100Ma) and reverse reconstruct it back to present day, and then store that present day geometry in the feature (in GPML or Shapefile) ? this is how all features are stored. And that (reverse) reconstruction differs between the half-stage rotation versions. I suppose you could emulate version 2 or version 3, and reverse reconstruct yourself, but that requires knowing exactly how these versions are implemented inside GPlates/pyGPlates. It?s better to just use pyGPlates to create a mid-ocean ridge (since it takes care of the reverse reconstruction for you). So in general I would recommend creating features (especially mid-ocean ridges) using GPlates or pyGPlates. This sample code shows how to create a mid-ocean ridge in pyGPlates: https://protect-au.mimecast.com/s/PdpcC71ZgLtLjwPnT8MCKI?domain=gplates.org ?and how to save it to a file (eg, GPML or Shapefile): https://protect-au.mimecast.com/s/wlowC81Zj6t98EJBT1OlVb?domain=gplates.org Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of John Cannon Sent: Saturday, 28 April 2018 4:03 PM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you?re digitizing a feature in GPlates then you don?t need to specify the Plate ID and Conjugate ID fields. When you switch to Half Stage Rotation (in Create Feature dialog) those fields disappear, and the Left Plate ID and Right Plate ID fields appear. And the feature type is taken care of when you selected MidOceanRidge. Saving to Shapefile will then map the feature type and spreading asymmetry to the GPGIM_TYPE and SPREAD_ASY shapefile attributes by default. However if you?re creating a feature outside GPlates, and without using pyGPlates, then you can add the GPGIM_TYPE shapefile attribute as ?gpml:MidOceanRidge? (as you noted) and the SPREAD_ASY attribute as a value in the range [-1,1] where 0 is symmetric spreading. And add the left/right plate IDs as L_PLATE and R_PLATE. Note that when GPlates loads your shapefile it will also load plate/conjugate IDs of zero by default, but they won?t get used for a half stage rotation. On a related note, the current public pyGPlates (revision 12) will not work with mid-ocean ridges created in GPlates 2.0 (ie, pyGPlates will not reconstruct them properly). To solve this there will be a simultaneous pyGPlates release alongside GPlates 2.1 (soon). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Saturday, 28 April 2018 10:53 AM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Sabin I need to be a bit more certain about what is needed so please excuse what will appear very evident as soon as I get things working. As I understand you, for a given OceanRidge, I need to specify the half-stage left plate plateid value and the half-stage right plate plate id. When I do this do I specify the normal PlateID field as none and the conjugate plateid field as none? Presumably I ought also to set the featuretype as having an appropriate GPGIM_TYPE value for ocean ridges (?gpml:MidOceanRidge? as I remember)? Do I need to have fields with something in my shapefile for reconstruction method and for spreading asymmetry? Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Sabin Zahirovic Sent: Friday, April 27, 2018 15:57 To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you digitize a feature to be a MidOceanicRidge, you can specify the Reconstruction Method as ?Half Stage Rotation? with a Left and Right PlateID. The ridge will then move with a half stage rotation, and will be equally spaced. Regards, Sabin -- DR SABIN ZAHIROVIC | Postdoctoral Research Associate School of Geosciences | Faculty of Science THE UNIVERSITY OF SYDNEY Rm 403, Madsen Building F09 | The University of Sydney | NSW | 2006 M +61 416 775 589 P +61 2 9351 3625 E sabin.zahirovic at sydney.edu.au | W https://protect-au.mimecast.com/s/PaIQC91ZkQt8VDZrT3Nqd5?domain=earthbyte.org | R http://www.researchgate.net/profile/Sabin_Zahirovic F https://protect-au.mimecast.com/s/40x5CgZowLH85Dg4TEvWWH?domain=facebook.com | T https://protect-au.mimecast.com/s/y_9NCjZrzqHg1V6Zcnkw7O?domain=twitter.com CRICOS 00026A This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments. From: GPlates-discuss > on behalf of "Eglington, Bruce" > Reply-To: GPlates general discussion mailing list > Date: Saturday, 28 April 2018 at 6:59 am To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID?s and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you?d count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the ?lossy? shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn?t like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From bruce.eglington at usask.ca Sun Apr 29 00:39:38 2018 From: bruce.eglington at usask.ca (Eglington, Bruce) Date: Sat, 28 Apr 2018 14:39:38 +0000 Subject: [GPlates-discuss] Defining mid-ocean ridges Message-ID: <7a9554fa3e114299802f70b04c7c3641@Mail01.usask.ca> Hi John Thanks for all the extra information and the warnings for what not to do. A few more questions though: 1) The polylineonsphere procedure expects a series of nodes. For present day mid-ocean ridges there are lots of offsets with ridge segments interspersed with transform faults so one needs lots of separate ridge polylines. Do these need to be named differently or can each ocean ridge polyline segment of what is geologically one ridge have the same name? 2) How does one judge what is left and what is right? Does GPlates have some preferred frame of reference (convention) for this? When a ridge or subduction zone is oriented north-south, it is easy to say what is left (west) and what is right (east) if one follows standard map viewing conventions but what if the ridge or subduction zone is oriented east-west. Does one view the polyline from an eastern or western perspective? More generally, does one view from the perspective of the starting node or from the ending node? 3) Separately, what is meant by subduction zone polarity? What is the convention for deciding what is left or right? Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of John Cannon Sent: Saturday, April 28, 2018 01:54 To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases ?I forgot to mention the reconstruction method maps to the RECON_METH shapefile attribute, and can be set to HalfStageRotationVersion3 (or HalfStageRotationVersion2 or HalfStageRotation). However note the issues below (better to use GPlates or pyGPlates to do all this)? Version 3 of the half-stage rotation was introduced in GPlates 2.0 (and, as mentioned in my last post, current pyGPlates only supports version 2 and below). With version 3 the spreading starts at the digitization time (ie, the reconstruction time at which the mid-ocean ridge was digitized in GPlates) such that subsequent adjustments to the spreading asymmetry no longer change the position of the ridge at the digitization time. With version 2, the ridge position would offset quite substantially making it unfeasible to adjust the asymmetry. So now the digitization time needed to be recorded, and this introduced the new property ?gpml:geometryImportTime? (also used in the new ability to reconstruct geometries using deforming and non-deforming topologies in GPlates 2.0). Oh, I just realized that gpml:geometryImportTime is not mapped in GPlates ! Which means it gets lost on saving to, or loading from, Shapefile and defaults to 0Ma. And this puts (version 3) mid-ocean ridges in the wrong location in GPlates 2.0 (when loaded from Shapefile). I?ll fix that for GPlates 2.1 (by mapping to a new attribute, say IMPORT_AGE). Actually, even when that is fixed, it will still be quite difficult to create a mid-ocean ridge without using GPlates or pyGPlates. This is because you need to take your mid-ocean ridge geometry (at say 100Ma) and reverse reconstruct it back to present day, and then store that present day geometry in the feature (in GPML or Shapefile) ? this is how all features are stored. And that (reverse) reconstruction differs between the half-stage rotation versions. I suppose you could emulate version 2 or version 3, and reverse reconstruct yourself, but that requires knowing exactly how these versions are implemented inside GPlates/pyGPlates. It?s better to just use pyGPlates to create a mid-ocean ridge (since it takes care of the reverse reconstruction for you). So in general I would recommend creating features (especially mid-ocean ridges) using GPlates or pyGPlates. This sample code shows how to create a mid-ocean ridge in pyGPlates: https://protect-au.mimecast.com/s/BZGwCGvmB5iwZ5vLfKnD8V?domain=gplates.org ?and how to save it to a file (eg, GPML or Shapefile): https://protect-au.mimecast.com/s/-2pSCJyp0qhxLENQiGJrOt?domain=gplates.org Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of John Cannon Sent: Saturday, 28 April 2018 4:03 PM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you?re digitizing a feature in GPlates then you don?t need to specify the Plate ID and Conjugate ID fields. When you switch to Half Stage Rotation (in Create Feature dialog) those fields disappear, and the Left Plate ID and Right Plate ID fields appear. And the feature type is taken care of when you selected MidOceanRidge. Saving to Shapefile will then map the feature type and spreading asymmetry to the GPGIM_TYPE and SPREAD_ASY shapefile attributes by default. However if you?re creating a feature outside GPlates, and without using pyGPlates, then you can add the GPGIM_TYPE shapefile attribute as ?gpml:MidOceanRidge? (as you noted) and the SPREAD_ASY attribute as a value in the range [-1,1] where 0 is symmetric spreading. And add the left/right plate IDs as L_PLATE and R_PLATE. Note that when GPlates loads your shapefile it will also load plate/conjugate IDs of zero by default, but they won?t get used for a half stage rotation. On a related note, the current public pyGPlates (revision 12) will not work with mid-ocean ridges created in GPlates 2.0 (ie, pyGPlates will not reconstruct them properly). To solve this there will be a simultaneous pyGPlates release alongside GPlates 2.1 (soon). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Saturday, 28 April 2018 10:53 AM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Sabin I need to be a bit more certain about what is needed so please excuse what will appear very evident as soon as I get things working. As I understand you, for a given OceanRidge, I need to specify the half-stage left plate plateid value and the half-stage right plate plate id. When I do this do I specify the normal PlateID field as none and the conjugate plateid field as none? Presumably I ought also to set the featuretype as having an appropriate GPGIM_TYPE value for ocean ridges (?gpml:MidOceanRidge? as I remember)? Do I need to have fields with something in my shapefile for reconstruction method and for spreading asymmetry? Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Sabin Zahirovic Sent: Friday, April 27, 2018 15:57 To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you digitize a feature to be a MidOceanicRidge, you can specify the Reconstruction Method as ?Half Stage Rotation? with a Left and Right PlateID. The ridge will then move with a half stage rotation, and will be equally spaced. Regards, Sabin -- DR SABIN ZAHIROVIC | Postdoctoral Research Associate School of Geosciences | Faculty of Science THE UNIVERSITY OF SYDNEY Rm 403, Madsen Building F09 | The University of Sydney | NSW | 2006 M +61 416 775 589 P +61 2 9351 3625 E sabin.zahirovic at sydney.edu.au | W https://protect-au.mimecast.com/s/3QV_CK1qJZtLx5Vrs3b0Na?domain=earthbyte.org | R http://www.researchgate.net/profile/Sabin_Zahirovic F https://protect-au.mimecast.com/s/M9uQCMwvLQTjV4N6CQMbUo?domain=facebook.com | T https://protect-au.mimecast.com/s/3y2ACNLwM9io5m4MS0HV_0?domain=twitter.com CRICOS 00026A This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments. From: GPlates-discuss > on behalf of "Eglington, Bruce" > Reply-To: GPlates general discussion mailing list > Date: Saturday, 28 April 2018 at 6:59 am To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID?s and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you?d count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the ?lossy? shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn?t like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From bruce.eglington at usask.ca Sun Apr 29 03:22:24 2018 From: bruce.eglington at usask.ca (Eglington, Bruce) Date: Sat, 28 Apr 2018 17:22:24 +0000 Subject: [GPlates-discuss] GPlates and linking to databases In-Reply-To: <1E57973145EDFC41B1D69DE11923922801B95E8BDF@ex-mbx-pro-05> References: <1E57973145EDFC41B1D69DE11923922801B95E8BDF@ex-mbx-pro-05> Message-ID: <06e41d2a03eb4efabd21b3a742f759ad@Mail01.usask.ca> Hi John I have just had another look at one of the EarthByte group?s published plateboundary gpml files as I try to better understand how to implement mid-ocean ridges and subduction zones. In some cases I note that the PlateID field in the gpml file has a value 0 (zero) when half stage rotation is going to occur (since in this situation left and right plate ID fields are defined and take precedence over PlateID). I just want to check that, within GPlates you don?t assume that PlateID is zero and ALSO assume that zero is the spin axis. To deal more easily with true polar wander, Dave Evans and I now assume that plate id zero is the absolute reference frame and the spin axis is a different number, the value of which is likely to vary in different versions of our models as we add TPW complexity. As an aside, for our models to be visually correct, we have to provide latitude/longitude graticules tied to the spin axis PlateID we use as separate gpml files and switch off the default GPlates graticule but that allows everything to look right, no matter what plate is considered to be fixed during reconstruction and visualization. As noted in a previous email, it would be very useful if the alpha value for the graticule could also be stored in the GPlates project file so that we don?t have to adjust it every time GPlates gets started. Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of John Cannon Sent: Saturday, April 28, 2018 00:03 To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you?re digitizing a feature in GPlates then you don?t need to specify the Plate ID and Conjugate ID fields. When you switch to Half Stage Rotation (in Create Feature dialog) those fields disappear, and the Left Plate ID and Right Plate ID fields appear. And the feature type is taken care of when you selected MidOceanRidge. Saving to Shapefile will then map the feature type and spreading asymmetry to the GPGIM_TYPE and SPREAD_ASY shapefile attributes by default. However if you?re creating a feature outside GPlates, and without using pyGPlates, then you can add the GPGIM_TYPE shapefile attribute as ?gpml:MidOceanRidge? (as you noted) and the SPREAD_ASY attribute as a value in the range [-1,1] where 0 is symmetric spreading. And add the left/right plate IDs as L_PLATE and R_PLATE. Note that when GPlates loads your shapefile it will also load plate/conjugate IDs of zero by default, but they won?t get used for a half stage rotation. On a related note, the current public pyGPlates (revision 12) will not work with mid-ocean ridges created in GPlates 2.0 (ie, pyGPlates will not reconstruct them properly). To solve this there will be a simultaneous pyGPlates release alongside GPlates 2.1 (soon). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Saturday, 28 April 2018 10:53 AM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Sabin I need to be a bit more certain about what is needed so please excuse what will appear very evident as soon as I get things working. As I understand you, for a given OceanRidge, I need to specify the half-stage left plate plateid value and the half-stage right plate plate id. When I do this do I specify the normal PlateID field as none and the conjugate plateid field as none? Presumably I ought also to set the featuretype as having an appropriate GPGIM_TYPE value for ocean ridges (?gpml:MidOceanRidge? as I remember)? Do I need to have fields with something in my shapefile for reconstruction method and for spreading asymmetry? Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Sabin Zahirovic Sent: Friday, April 27, 2018 15:57 To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you digitize a feature to be a MidOceanicRidge, you can specify the Reconstruction Method as ?Half Stage Rotation? with a Left and Right PlateID. The ridge will then move with a half stage rotation, and will be equally spaced. Regards, Sabin -- DR SABIN ZAHIROVIC | Postdoctoral Research Associate School of Geosciences | Faculty of Science THE UNIVERSITY OF SYDNEY Rm 403, Madsen Building F09 | The University of Sydney | NSW | 2006 M +61 416 775 589 P +61 2 9351 3625 E sabin.zahirovic at sydney.edu.au | W https://protect-au.mimecast.com/s/1pf0CYWL1vi5ovEvu0QmN2?domain=earthbyte.org | R http://www.researchgate.net/profile/Sabin_Zahirovic F https://protect-au.mimecast.com/s/5G5GC1WZXri8w020cpIzwb?domain=facebook.com | T https://protect-au.mimecast.com/s/9iZtC2xZYvC7o4P4T2DEVP?domain=twitter.com CRICOS 00026A This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments. From: GPlates-discuss > on behalf of "Eglington, Bruce" > Reply-To: GPlates general discussion mailing list > Date: Saturday, 28 April 2018 at 6:59 am To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID?s and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you?d count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the ?lossy? shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn?t like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From john.cannon at sydney.edu.au Sun Apr 29 16:15:48 2018 From: john.cannon at sydney.edu.au (John Cannon) Date: Sun, 29 Apr 2018 06:15:48 +0000 Subject: [GPlates-discuss] Defining mid-ocean ridges In-Reply-To: <7a9554fa3e114299802f70b04c7c3641@Mail01.usask.ca> References: <7a9554fa3e114299802f70b04c7c3641@Mail01.usask.ca> Message-ID: <1E57973145EDFC41B1D69DE11923922801B95E915E@ex-mbx-pro-05> Hi Bruce, 1) Each ridge polyline segment can be in a different feature with the same name. Or can be a single feature containing multiple polyline segments with the same property name (eg, ?gpml:centerLineOf?). Or it?s acceptable to join the ridge and transform segments into a single polyline. 2) In general left/right indicates which side of the line you are on when following the order of points in the line (eg, if line goes South to North then left is West, and for North to South then left is East). This is how subduction polarity works ? see (3). Ideally this should also be followed for mid-ocean ridges, but here the concept of left/right is also used for reconstruction. And it?s possible that the left/right plate IDs might, strictly speaking, be incorrect (ie, swapped such that the left plate ID corresponds to the right side, and right plate ID to left side) but you could still get the correct reconstruction provided the spreading asymmetry is also swapped (ie, negated). For symmetric spreading (ie, asymmetry is zero) this swapping won?t affect reconstruction. 3) Subduction polarity indicates which side the over-riding plate is on (eg, if subduction zone line goes South to North and polarity is Left then over-riding plate is West and subducting plate is East). There?s also a description of these features/properties at https://protect-au.mimecast.com/s/TfbBCBNZwLivZvQWizcWta?domain=gplates.org , for example the MidOceanRidge at https://protect-au.mimecast.com/s/f2_dCD1jy9tX6X2Ys55Aez?domain=gplates.org and subduction polarity at https://protect-au.mimecast.com/s/T1KYCE8kz9tBABr4cp5MTc?domain=gplates.org Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Sunday, 29 April 2018 12:40 AM To: GPlates general discussion mailing list Subject: [GPlates-discuss] Defining mid-ocean ridges Hi John Thanks for all the extra information and the warnings for what not to do. A few more questions though: 1) The polylineonsphere procedure expects a series of nodes. For present day mid-ocean ridges there are lots of offsets with ridge segments interspersed with transform faults so one needs lots of separate ridge polylines. Do these need to be named differently or can each ocean ridge polyline segment of what is geologically one ridge have the same name? 2) How does one judge what is left and what is right? Does GPlates have some preferred frame of reference (convention) for this? When a ridge or subduction zone is oriented north-south, it is easy to say what is left (west) and what is right (east) if one follows standard map viewing conventions but what if the ridge or subduction zone is oriented east-west. Does one view the polyline from an eastern or western perspective? More generally, does one view from the perspective of the starting node or from the ending node? 3) Separately, what is meant by subduction zone polarity? What is the convention for deciding what is left or right? Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of John Cannon Sent: Saturday, April 28, 2018 01:54 To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases ?I forgot to mention the reconstruction method maps to the RECON_METH shapefile attribute, and can be set to HalfStageRotationVersion3 (or HalfStageRotationVersion2 or HalfStageRotation). However note the issues below (better to use GPlates or pyGPlates to do all this)? Version 3 of the half-stage rotation was introduced in GPlates 2.0 (and, as mentioned in my last post, current pyGPlates only supports version 2 and below). With version 3 the spreading starts at the digitization time (ie, the reconstruction time at which the mid-ocean ridge was digitized in GPlates) such that subsequent adjustments to the spreading asymmetry no longer change the position of the ridge at the digitization time. With version 2, the ridge position would offset quite substantially making it unfeasible to adjust the asymmetry. So now the digitization time needed to be recorded, and this introduced the new property ?gpml:geometryImportTime? (also used in the new ability to reconstruct geometries using deforming and non-deforming topologies in GPlates 2.0). Oh, I just realized that gpml:geometryImportTime is not mapped in GPlates ! Which means it gets lost on saving to, or loading from, Shapefile and defaults to 0Ma. And this puts (version 3) mid-ocean ridges in the wrong location in GPlates 2.0 (when loaded from Shapefile). I?ll fix that for GPlates 2.1 (by mapping to a new attribute, say IMPORT_AGE). Actually, even when that is fixed, it will still be quite difficult to create a mid-ocean ridge without using GPlates or pyGPlates. This is because you need to take your mid-ocean ridge geometry (at say 100Ma) and reverse reconstruct it back to present day, and then store that present day geometry in the feature (in GPML or Shapefile) ? this is how all features are stored. And that (reverse) reconstruction differs between the half-stage rotation versions. I suppose you could emulate version 2 or version 3, and reverse reconstruct yourself, but that requires knowing exactly how these versions are implemented inside GPlates/pyGPlates. It?s better to just use pyGPlates to create a mid-ocean ridge (since it takes care of the reverse reconstruction for you). So in general I would recommend creating features (especially mid-ocean ridges) using GPlates or pyGPlates. This sample code shows how to create a mid-ocean ridge in pyGPlates: https://protect-au.mimecast.com/s/EyZkCGvmB5iwywx3HQes6K?domain=gplates.org ?and how to save it to a file (eg, GPML or Shapefile): https://protect-au.mimecast.com/s/5EvwCJyp0qhx4xG9CvW1SE?domain=gplates.org Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of John Cannon Sent: Saturday, 28 April 2018 4:03 PM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you?re digitizing a feature in GPlates then you don?t need to specify the Plate ID and Conjugate ID fields. When you switch to Half Stage Rotation (in Create Feature dialog) those fields disappear, and the Left Plate ID and Right Plate ID fields appear. And the feature type is taken care of when you selected MidOceanRidge. Saving to Shapefile will then map the feature type and spreading asymmetry to the GPGIM_TYPE and SPREAD_ASY shapefile attributes by default. However if you?re creating a feature outside GPlates, and without using pyGPlates, then you can add the GPGIM_TYPE shapefile attribute as ?gpml:MidOceanRidge? (as you noted) and the SPREAD_ASY attribute as a value in the range [-1,1] where 0 is symmetric spreading. And add the left/right plate IDs as L_PLATE and R_PLATE. Note that when GPlates loads your shapefile it will also load plate/conjugate IDs of zero by default, but they won?t get used for a half stage rotation. On a related note, the current public pyGPlates (revision 12) will not work with mid-ocean ridges created in GPlates 2.0 (ie, pyGPlates will not reconstruct them properly). To solve this there will be a simultaneous pyGPlates release alongside GPlates 2.1 (soon). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Saturday, 28 April 2018 10:53 AM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Sabin I need to be a bit more certain about what is needed so please excuse what will appear very evident as soon as I get things working. As I understand you, for a given OceanRidge, I need to specify the half-stage left plate plateid value and the half-stage right plate plate id. When I do this do I specify the normal PlateID field as none and the conjugate plateid field as none? Presumably I ought also to set the featuretype as having an appropriate GPGIM_TYPE value for ocean ridges (?gpml:MidOceanRidge? as I remember)? Do I need to have fields with something in my shapefile for reconstruction method and for spreading asymmetry? Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Sabin Zahirovic Sent: Friday, April 27, 2018 15:57 To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you digitize a feature to be a MidOceanicRidge, you can specify the Reconstruction Method as ?Half Stage Rotation? with a Left and Right PlateID. The ridge will then move with a half stage rotation, and will be equally spaced. Regards, Sabin -- DR SABIN ZAHIROVIC | Postdoctoral Research Associate School of Geosciences | Faculty of Science THE UNIVERSITY OF SYDNEY Rm 403, Madsen Building F09 | The University of Sydney | NSW | 2006 M +61 416 775 589 P +61 2 9351 3625 E sabin.zahirovic at sydney.edu.au | W https://protect-au.mimecast.com/s/Cvf0CK1qJZtLALg7iGLosE?domain=earthbyte.org | R http://www.researchgate.net/profile/Sabin_Zahirovic F https://protect-au.mimecast.com/s/fgnsCMwvLQTjojKBtPJG8T?domain=facebook.com | T https://protect-au.mimecast.com/s/TydSCNLwM9ionoKkiyHGcs?domain=twitter.com CRICOS 00026A This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments. From: GPlates-discuss > on behalf of "Eglington, Bruce" > Reply-To: GPlates general discussion mailing list > Date: Saturday, 28 April 2018 at 6:59 am To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID?s and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you?d count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the ?lossy? shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn?t like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From john.cannon at sydney.edu.au Sun Apr 29 16:16:11 2018 From: john.cannon at sydney.edu.au (John Cannon) Date: Sun, 29 Apr 2018 06:16:11 +0000 Subject: [GPlates-discuss] Defining mid-ocean ridges In-Reply-To: <06e41d2a03eb4efabd21b3a742f759ad@Mail01.usask.ca> References: <1E57973145EDFC41B1D69DE11923922801B95E8BDF@ex-mbx-pro-05> <06e41d2a03eb4efabd21b3a742f759ad@Mail01.usask.ca> Message-ID: <1E57973145EDFC41B1D69DE11923922801B95E9182@ex-mbx-pro-05> Hi Bruce, For half-stage rotations GPlates ignores the main reconstruction PlateID (and only uses left/right plate IDs). So the final ?half-stage? rotation (that positions/reconstructs the mid-ocean ridge) uses the rotations of the left and right plates through time, and those are calculated relative to the anchor plate ID, which the user can change (ie, is not assumed to be zero). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Sunday, 29 April 2018 3:22 AM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi John I have just had another look at one of the EarthByte group?s published plateboundary gpml files as I try to better understand how to implement mid-ocean ridges and subduction zones. In some cases I note that the PlateID field in the gpml file has a value 0 (zero) when half stage rotation is going to occur (since in this situation left and right plate ID fields are defined and take precedence over PlateID). I just want to check that, within GPlates you don?t assume that PlateID is zero and ALSO assume that zero is the spin axis. To deal more easily with true polar wander, Dave Evans and I now assume that plate id zero is the absolute reference frame and the spin axis is a different number, the value of which is likely to vary in different versions of our models as we add TPW complexity. As an aside, for our models to be visually correct, we have to provide latitude/longitude graticules tied to the spin axis PlateID we use as separate gpml files and switch off the default GPlates graticule but that allows everything to look right, no matter what plate is considered to be fixed during reconstruction and visualization. As noted in a previous email, it would be very useful if the alpha value for the graticule could also be stored in the GPlates project file so that we don?t have to adjust it every time GPlates gets started. Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of John Cannon Sent: Saturday, April 28, 2018 00:03 To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you?re digitizing a feature in GPlates then you don?t need to specify the Plate ID and Conjugate ID fields. When you switch to Half Stage Rotation (in Create Feature dialog) those fields disappear, and the Left Plate ID and Right Plate ID fields appear. And the feature type is taken care of when you selected MidOceanRidge. Saving to Shapefile will then map the feature type and spreading asymmetry to the GPGIM_TYPE and SPREAD_ASY shapefile attributes by default. However if you?re creating a feature outside GPlates, and without using pyGPlates, then you can add the GPGIM_TYPE shapefile attribute as ?gpml:MidOceanRidge? (as you noted) and the SPREAD_ASY attribute as a value in the range [-1,1] where 0 is symmetric spreading. And add the left/right plate IDs as L_PLATE and R_PLATE. Note that when GPlates loads your shapefile it will also load plate/conjugate IDs of zero by default, but they won?t get used for a half stage rotation. On a related note, the current public pyGPlates (revision 12) will not work with mid-ocean ridges created in GPlates 2.0 (ie, pyGPlates will not reconstruct them properly). To solve this there will be a simultaneous pyGPlates release alongside GPlates 2.1 (soon). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Saturday, 28 April 2018 10:53 AM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Sabin I need to be a bit more certain about what is needed so please excuse what will appear very evident as soon as I get things working. As I understand you, for a given OceanRidge, I need to specify the half-stage left plate plateid value and the half-stage right plate plate id. When I do this do I specify the normal PlateID field as none and the conjugate plateid field as none? Presumably I ought also to set the featuretype as having an appropriate GPGIM_TYPE value for ocean ridges (?gpml:MidOceanRidge? as I remember)? Do I need to have fields with something in my shapefile for reconstruction method and for spreading asymmetry? Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Sabin Zahirovic Sent: Friday, April 27, 2018 15:57 To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you digitize a feature to be a MidOceanicRidge, you can specify the Reconstruction Method as ?Half Stage Rotation? with a Left and Right PlateID. The ridge will then move with a half stage rotation, and will be equally spaced. Regards, Sabin -- DR SABIN ZAHIROVIC | Postdoctoral Research Associate School of Geosciences | Faculty of Science THE UNIVERSITY OF SYDNEY Rm 403, Madsen Building F09 | The University of Sydney | NSW | 2006 M +61 416 775 589 P +61 2 9351 3625 E sabin.zahirovic at sydney.edu.au | W https://protect-au.mimecast.com/s/p9DQCyoNVrcYjYxpIZ5UMt?domain=earthbyte.org | R http://www.researchgate.net/profile/Sabin_Zahirovic F https://protect-au.mimecast.com/s/0lBsCANZvPijmjk8H9GlXj?domain=facebook.com | T https://protect-au.mimecast.com/s/hlrXCBNZwLivZvYyHje-Rt?domain=twitter.com CRICOS 00026A This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments. From: GPlates-discuss > on behalf of "Eglington, Bruce" > Reply-To: GPlates general discussion mailing list > Date: Saturday, 28 April 2018 at 6:59 am To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID?s and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you?d count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the ?lossy? shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn?t like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From john.cannon at sydney.edu.au Sun Apr 29 16:16:37 2018 From: john.cannon at sydney.edu.au (John Cannon) Date: Sun, 29 Apr 2018 06:16:37 +0000 Subject: [GPlates-discuss] Saving global settings to project files In-Reply-To: <06e41d2a03eb4efabd21b3a742f759ad@Mail01.usask.ca> References: <1E57973145EDFC41B1D69DE11923922801B95E8BDF@ex-mbx-pro-05> <06e41d2a03eb4efabd21b3a742f759ad@Mail01.usask.ca> Message-ID: <1E57973145EDFC41B1D69DE11923922801B95E91A8@ex-mbx-pro-05> Hi Bruce, > it would be very useful if the alpha value for the graticule could also be stored in the GPlates project file so that we don?t have to adjust it every time GPlates gets started. Adding graticule alpha to project files is an interesting one. I can definitely see why you?d want that in your situation (where you load your own rotating graticules from a gpml file, and want the default graticule to disappear). Normally we only save state relevant to the individual layers, and not ?global? state (for reasons noted below), but I think it?s probably OK in some cases. /detalMode on/ Currently we don?t save global state such anchor plate ID, or camera location/view/zoom, or whether view is currently in 3D globe or 2D map view. This is mainly so that if you load someone else?s project file it doesn?t change your global settings (such that when you?re finished with their project ? maybe you?ve unloaded all their data files - then you have to manually change your global settings back, such as making the graticule reappear). However there is ?File > Clear Session? which resets all project/session state, and this would reset these global settings back to their default (such as making the graticule reappear) as well as unload all files/layers. This is better, and easier, than manually unloading all the data files (which doesn?t make the graticule reappear). And loading a new project/session has essentially the same reset effect (but resets with loaded state instead of default state). So if you use projects (and recent sessions) to swap back and forth between work sessions then I don?t think saving some global state would be a problem. There is probably some other global state that falls under this category too, perhaps visual state such as background colour (I think you noted in the past) and maybe the text overlay. Though I?m not sure about things like camera and projection settings. Anchor plate ID should probably be left out since it?s not really a ?visual? state, although it?s not clear to me. Then there?s other things like the animation settings and current reconstruction time which are also unclear. Where to draw the line ? As a side note, there are already some cases where global state is saved, such as the settings under ?View > Geometry Visibility? such as ?Show lines?, because those settings should have been individual layer settings (rather than applied globally to all layers). /detailMode off/ Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Sunday, 29 April 2018 3:22 AM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi John I have just had another look at one of the EarthByte group?s published plateboundary gpml files as I try to better understand how to implement mid-ocean ridges and subduction zones. In some cases I note that the PlateID field in the gpml file has a value 0 (zero) when half stage rotation is going to occur (since in this situation left and right plate ID fields are defined and take precedence over PlateID). I just want to check that, within GPlates you don?t assume that PlateID is zero and ALSO assume that zero is the spin axis. To deal more easily with true polar wander, Dave Evans and I now assume that plate id zero is the absolute reference frame and the spin axis is a different number, the value of which is likely to vary in different versions of our models as we add TPW complexity. As an aside, for our models to be visually correct, we have to provide latitude/longitude graticules tied to the spin axis PlateID we use as separate gpml files and switch off the default GPlates graticule but that allows everything to look right, no matter what plate is considered to be fixed during reconstruction and visualization. As noted in a previous email, it would be very useful if the alpha value for the graticule could also be stored in the GPlates project file so that we don?t have to adjust it every time GPlates gets started. Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of John Cannon Sent: Saturday, April 28, 2018 00:03 To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you?re digitizing a feature in GPlates then you don?t need to specify the Plate ID and Conjugate ID fields. When you switch to Half Stage Rotation (in Create Feature dialog) those fields disappear, and the Left Plate ID and Right Plate ID fields appear. And the feature type is taken care of when you selected MidOceanRidge. Saving to Shapefile will then map the feature type and spreading asymmetry to the GPGIM_TYPE and SPREAD_ASY shapefile attributes by default. However if you?re creating a feature outside GPlates, and without using pyGPlates, then you can add the GPGIM_TYPE shapefile attribute as ?gpml:MidOceanRidge? (as you noted) and the SPREAD_ASY attribute as a value in the range [-1,1] where 0 is symmetric spreading. And add the left/right plate IDs as L_PLATE and R_PLATE. Note that when GPlates loads your shapefile it will also load plate/conjugate IDs of zero by default, but they won?t get used for a half stage rotation. On a related note, the current public pyGPlates (revision 12) will not work with mid-ocean ridges created in GPlates 2.0 (ie, pyGPlates will not reconstruct them properly). To solve this there will be a simultaneous pyGPlates release alongside GPlates 2.1 (soon). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Saturday, 28 April 2018 10:53 AM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Sabin I need to be a bit more certain about what is needed so please excuse what will appear very evident as soon as I get things working. As I understand you, for a given OceanRidge, I need to specify the half-stage left plate plateid value and the half-stage right plate plate id. When I do this do I specify the normal PlateID field as none and the conjugate plateid field as none? Presumably I ought also to set the featuretype as having an appropriate GPGIM_TYPE value for ocean ridges (?gpml:MidOceanRidge? as I remember)? Do I need to have fields with something in my shapefile for reconstruction method and for spreading asymmetry? Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Sabin Zahirovic Sent: Friday, April 27, 2018 15:57 To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you digitize a feature to be a MidOceanicRidge, you can specify the Reconstruction Method as ?Half Stage Rotation? with a Left and Right PlateID. The ridge will then move with a half stage rotation, and will be equally spaced. Regards, Sabin -- DR SABIN ZAHIROVIC | Postdoctoral Research Associate School of Geosciences | Faculty of Science THE UNIVERSITY OF SYDNEY Rm 403, Madsen Building F09 | The University of Sydney | NSW | 2006 M +61 416 775 589 P +61 2 9351 3625 E sabin.zahirovic at sydney.edu.au | W https://protect-au.mimecast.com/s/gPhFCq7BKYt949ElSZbIyU?domain=earthbyte.org | R http://www.researchgate.net/profile/Sabin_Zahirovic F https://protect-au.mimecast.com/s/mSsdCvl0PoCvBvYKuzNhDQ?domain=facebook.com | T https://protect-au.mimecast.com/s/FKlPCwVLQmioOoq9UKS2lH?domain=twitter.com CRICOS 00026A This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments. From: GPlates-discuss > on behalf of "Eglington, Bruce" > Reply-To: GPlates general discussion mailing list > Date: Saturday, 28 April 2018 at 6:59 am To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID?s and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you?d count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the ?lossy? shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn?t like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: From bruce.eglington at usask.ca Sun Apr 29 18:14:49 2018 From: bruce.eglington at usask.ca (Eglington, Bruce) Date: Sun, 29 Apr 2018 08:14:49 +0000 Subject: [GPlates-discuss] Saving global settings to project files In-Reply-To: <1E57973145EDFC41B1D69DE11923922801B95E91A8@ex-mbx-pro-05> References: <1E57973145EDFC41B1D69DE11923922801B95E8BDF@ex-mbx-pro-05> <06e41d2a03eb4efabd21b3a742f759ad@Mail01.usask.ca>, <1E57973145EDFC41B1D69DE11923922801B95E91A8@ex-mbx-pro-05> Message-ID: <1524989689259.77245@usask.ca> Hi John Thanks Bruce ________________________________ From: GPlates-discuss on behalf of John Cannon Sent: Sunday, April 29, 2018 12:16:37 AM To: GPlates general discussion mailing list Subject: [GPlates-discuss] Saving global settings to project files Hi Bruce, > it would be very useful if the alpha value for the graticule could also be stored in the GPlates project file so that we don?t have to adjust it every time GPlates gets started. Adding graticule alpha to project files is an interesting one. I can definitely see why you?d want that in your situation (where you load your own rotating graticules from a gpml file, and want the default graticule to disappear). Normally we only save state relevant to the individual layers, and not ?global? state (for reasons noted below), but I think it?s probably OK in some cases. /detalMode on/ Currently we don?t save global state such anchor plate ID, or camera location/view/zoom, or whether view is currently in 3D globe or 2D map view. This is mainly so that if you load someone else?s project file it doesn?t change your global settings (such that when you?re finished with their project ? maybe you?ve unloaded all their data files - then you have to manually change your global settings back, such as making the graticule reappear). However there is ?File > Clear Session? which resets all project/session state, and this would reset these global settings back to their default (such as making the graticule reappear) as well as unload all files/layers. This is better, and easier, than manually unloading all the data files (which doesn?t make the graticule reappear). And loading a new project/session has essentially the same reset effect (but resets with loaded state instead of default state). So if you use projects (and recent sessions) to swap back and forth between work sessions then I don?t think saving some global state would be a problem. There is probably some other global state that falls under this category too, perhaps visual state such as background colour (I think you noted in the past) and maybe the text overlay. Though I?m not sure about things like camera and projection settings. Anchor plate ID should probably be left out since it?s not really a ?visual? state, although it?s not clear to me. Then there?s other things like the animation settings and current reconstruction time which are also unclear. Where to draw the line :) As a side note, there are already some cases where global state is saved, such as the settings under ?View > Geometry Visibility? such as ?Show lines?, because those settings should have been individual layer settings (rather than applied globally to all layers). /detailMode off/ Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Sunday, 29 April 2018 3:22 AM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi John I have just had another look at one of the EarthByte group?s published plateboundary gpml files as I try to better understand how to implement mid-ocean ridges and subduction zones. In some cases I note that the PlateID field in the gpml file has a value 0 (zero) when half stage rotation is going to occur (since in this situation left and right plate ID fields are defined and take precedence over PlateID). I just want to check that, within GPlates you don?t assume that PlateID is zero and ALSO assume that zero is the spin axis. To deal more easily with true polar wander, Dave Evans and I now assume that plate id zero is the absolute reference frame and the spin axis is a different number, the value of which is likely to vary in different versions of our models as we add TPW complexity. As an aside, for our models to be visually correct, we have to provide latitude/longitude graticules tied to the spin axis PlateID we use as separate gpml files and switch off the default GPlates graticule but that allows everything to look right, no matter what plate is considered to be fixed during reconstruction and visualization. As noted in a previous email, it would be very useful if the alpha value for the graticule could also be stored in the GPlates project file so that we don?t have to adjust it every time GPlates gets started. Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of John Cannon Sent: Saturday, April 28, 2018 00:03 To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you?re digitizing a feature in GPlates then you don?t need to specify the Plate ID and Conjugate ID fields. When you switch to Half Stage Rotation (in Create Feature dialog) those fields disappear, and the Left Plate ID and Right Plate ID fields appear. And the feature type is taken care of when you selected MidOceanRidge. Saving to Shapefile will then map the feature type and spreading asymmetry to the GPGIM_TYPE and SPREAD_ASY shapefile attributes by default. However if you?re creating a feature outside GPlates, and without using pyGPlates, then you can add the GPGIM_TYPE shapefile attribute as ?gpml:MidOceanRidge? (as you noted) and the SPREAD_ASY attribute as a value in the range [-1,1] where 0 is symmetric spreading. And add the left/right plate IDs as L_PLATE and R_PLATE. Note that when GPlates loads your shapefile it will also load plate/conjugate IDs of zero by default, but they won?t get used for a half stage rotation. On a related note, the current public pyGPlates (revision 12) will not work with mid-ocean ridges created in GPlates 2.0 (ie, pyGPlates will not reconstruct them properly). To solve this there will be a simultaneous pyGPlates release alongside GPlates 2.1 (soon). Regards, John From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Saturday, 28 April 2018 10:53 AM To: GPlates general discussion mailing list Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Sabin I need to be a bit more certain about what is needed so please excuse what will appear very evident as soon as I get things working. As I understand you, for a given OceanRidge, I need to specify the half-stage left plate plateid value and the half-stage right plate plate id. When I do this do I specify the normal PlateID field as none and the conjugate plateid field as none? Presumably I ought also to set the featuretype as having an appropriate GPGIM_TYPE value for ocean ridges (?gpml:MidOceanRidge? as I remember)? Do I need to have fields with something in my shapefile for reconstruction method and for spreading asymmetry? Thanks Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Sabin Zahirovic Sent: Friday, April 27, 2018 15:57 To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, If you digitize a feature to be a MidOceanicRidge, you can specify the Reconstruction Method as ?Half Stage Rotation? with a Left and Right PlateID. The ridge will then move with a half stage rotation, and will be equally spaced. Regards, Sabin -- DR SABIN ZAHIROVIC | Postdoctoral Research Associate School of Geosciences | Faculty of Science THE UNIVERSITY OF SYDNEY Rm 403, Madsen Building F09 | The University of Sydney | NSW | 2006 M +61 416 775 589 P +61 2 9351 3625 E sabin.zahirovic at sydney.edu.au | W https://protect-au.mimecast.com/s/bTdACq7BKYt94mRPFZBB6M?domain=earthbyte.org | R http://www.researchgate.net/profile/Sabin_Zahirovic F https://protect-au.mimecast.com/s/-gghCvl0PoCvBM4jhzk9hH?domain=facebook.com | T https://protect-au.mimecast.com/s/G2-dCwVLQmioOYPzhKGFrP?domain=twitter.com CRICOS 00026A This email plus any attachments to it are confidential. Any unauthorised use is strictly prohibited. If you receive this email in error, please delete it and any attachments. From: GPlates-discuss > on behalf of "Eglington, Bruce" > Reply-To: GPlates general discussion mailing list > Date: Saturday, 28 April 2018 at 6:59 am To: GPlates general discussion mailing list > Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Christian Thanks for the suggestions. Separate from the database issue, can you explain how to make a mid ocean ridge stay equally spaced between two continents or oceanic plates? I assumed that one would link the MOR line to one of the adjacent PlateID?s and specify the other adjacent PlateID as the conjugate plate but that does not seem to work. I have also tried to link the adjacent plates as LeftPlate and RightPlate but also no joy. Any suggestions will be a great help as none of the tutorials seem to explain this, at least not in a way that makes sense to me. Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Christian.Heine at shell.com Sent: Wednesday, March 7, 2018 09:28 To: gplates-discuss at mailman.sydney.edu.au Subject: Re: [GPlates-discuss] GPlates and linking to databases Hi Bruce, AFAIK there is no such capability, unless you?d count hooking up to a WFS server as one. We had discussions a few years back of linking to PostGIS, but also we recently talked about the new OGC Geopackage standard (being sqlite-based I would also classify this as database, see geopackage.org) but so far this has not matured further (I think). In theory the implementation should be relative straightforward as GPlates is already using the GDAL library (which reads PostGIS, sqlite, geopackage etc). However, the complications come in when turning relatively arbitrary database attributes/columns into the GPlates feature model (GPGIM). I suppose some way you could go about this (thinking out loud)could be a crude combination of pyGPlates (or even using the built in Python shell) together with either the GDAL python library or something like psycopg2 to connect to your (geospatial) database (for some older project I used psycopg2+postgis/postgresql with GMT quite successfully). However, this will also require you to map the feature you import to the GPGIM. You could write some script to write out native GPML files from your DB as well to avoid the ?lossy? shapefile format (ie col header truncation if > 12 chars, 255 char limits in attributes etc etc) which is the only way to bridge the ESRI world with GPlates currently. I would very much support any effort to get GPlates and the GPGIM to read from DBs as this would also help a lot to improve the interface to commonly used GIS tools, especially in the corporate environment which somehow doesn?t like open formats. Not sure if John & co will have some surprises in this direction in the upcoming 2.1 release. Cheers, Christian From: GPlates-discuss [mailto:gplates-discuss-bounces at mailman.sydney.edu.au] On Behalf Of Eglington, Bruce Sent: Wednesday, March 07, 2018 3:55 PM To: 'gplates-discuss at mailman.sydney.edu.au' > Subject: [GPlates-discuss] GPlates and linking to databases Good day Is there currently any design capability to link to tables or queries in any form of database already built in to GPlates or are the only data input options ROT, GPML, SHP, NETCDF and various raster formats? If not currently available, are there any plans to add database capabilities in the foreseeable future? Regards Bruce Bruce Eglington (Ph.D.) Murray Pyke Chair Geological Sciences University of Saskatchewan 114 Science Place Saskatoon SK S7N 5E2 Canada bruce.eglington at usask.ca +1-306-966-5732 -------------- next part -------------- An HTML attachment was scrubbed... URL: