Rotating​ ​rasters​ ​and​ ​age-based​ ​masking​ ​of Raster​ ​data

Authors:​ ​Christian​ ​Heine​ ​&​ ​Kara​ ​J.​ ​Matthews Edited​ ​by:​ ​Julia​ ​Sheehan

EarthByte​ ​Research​ ​Group,​ ​School​ ​of​ ​Geosciences,​ ​The​ ​University​ ​of​ ​Sydney,​ ​Australia

Rotating​ ​rasters​ ​and​ ​age-based​ ​masking​ ​of​ ​Raster​ ​data

​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​ ​Background Exercise​ ​1:​ ​Rotating​ ​and​ ​cookie​ ​cutting​ ​raster​ ​data Exercise​ ​2:​ ​Age​ ​-based​ ​masking​ ​or​ ​raster​ ​data Appendix

Background With​ ​the​ ​release​ ​of​ ​version​ ​0.9.10​ ​of​ ​GPlates​ ​in​ ​2010,​ ​functionality​ ​to​ ​do age-based​ ​masking​ ​of​ ​raster​ ​data​ ​was​ ​included.​ ​This​ ​means​ ​any​ ​age-grid can​ ​be​ ​used​ ​to​ ​mask​ ​underlying​ ​rasters​ ​which​ ​in​ ​turn​ ​can​ ​be​ ​cookie-cut​ ​by polygons​ ​and​ ​rotated​ ​to​ ​their​ ​position​ ​in​ ​the​ ​past. In​ ​this​ ​tutorial​ ​we​ ​will​ ​be​ ​working​ ​on​ ​importing​ ​and​ ​visualising​ ​raster​ ​data​ ​in GPlates​ ​and​ ​rotating​ ​and​ ​masking​ ​raster​ ​data​ ​back​ ​through​ ​time.​ ​The tutorial​ ​will​ ​use​ ​the​ ​data​ ​included​ ​in​ ​the​ ​GPlates​ ​distribution​ ​in​ ​the​ ​Sample data​ ​folder​ ​(see​ ​the​ ​“Sample​ ​data”​ ​section​ ​under​ ​Appdx.) Today​ ​we​ ​will​ ​be​ ​working​ ​with​ ​raster​ ​files.​ ​For​ ​all​ ​those​ ​computer​ ​illiterate folk​ ​out​ ​there,​ ​a​ ​raster​ ​is​ ​simply​ ​a​ ​file​ ​which​ ​is​ ​made​ ​of​ ​2-dimensional​ ​grid​ ​of pixels​ ​and​ ​is​ ​stored​ ​as​ ​JPEGS​ ​or​ ​grid​ ​files​ ​like​ ​netCDF.​ ​This​ ​is​ ​different​ ​to Vector​ ​data​ ​we​ ​have​ ​used​ ​in​ ​previous​ ​tutorial,​ ​which​ ​is​ ​composed​ ​of​ ​points, lines​ ​and​ ​polygons.

Exercise​ ​1:​ ​ ​Rotating​ ​and​ ​cookie​ ​cutting​ ​raster​ ​data This​ ​tutorial​ ​will​ ​show​ ​how​ ​to​ ​cookie-cut​ ​polygons​ ​from​ ​rasters​ ​and​ ​rotate them​ ​to​ ​paleopositions. In​ ​order​ ​to​ ​split​ ​a​ ​global​ ​raster​ ​file​ ​into​ ​different​ ​polygons,​ ​load​ ​the​ ​sample data​ ​into​ ​GPlates.​ ​Specifically,​ ​load​ ​the​ ​following​ ​files​ ​that​ ​have​ ​been​ ​used​ ​in tutorial​ ​10a. 1.​ ​The​ ​global​ ​rotation​ ​file (Global_EarthByte_GPlates_Rotation_20100927.rot) 2.​ ​The​ ​global​ ​static​ ​polygon​ ​file (Global_EarthByte_GPlates_PresentDay_StaticPlatePolygons_20100927.gpml ) 3.​ ​The​ ​global​ ​topography/bathymetry​ ​image​ ​(color​ ​etopo1_ice_low.jpg supplied​ ​with​ ​this​ ​tutorial)​ ​or​ ​the​ ​global​ ​gravity​ ​image​ ​supplied​ ​with​ ​GPlates (DNSC08GRA_6m.jpg). Once​ ​this​ ​has​ ​been​ ​done,​ y ​ ou​ ​should​ ​have​ ​a​ ​something​ ​on​ ​your​ ​GPlates main​ ​window​ ​which​ ​looks​ l​ ike​ ​in​ ​Fig.1.

Figure​ ​1.​​ ​GPlates​ ​windows​ ​with​ ​sample​ ​data​ ​for​ ​tutorial​ ​1​ ​loaded.Here,​ ​we​ ​have​ ​two​ ​raster images​ ​loaded​ ​(red​ ​rectangle):​ ​the​ ​global​ ​topography​ ​and​ ​the​ ​global​ ​gravity.​ ​Both​ ​are automatically​ ​classified​ ​as​ ​“Reconstructed​ ​raster”.

The​ ​next​ ​step​ ​involves​ ​telling​ ​GPlates​ ​to​ ​cut​ ​the​ ​raster​ ​into​ ​different​ ​pieces by​ ​using​ ​our​ ​global​ ​static​ ​polygon​ ​layer.​ ​It​ ​is​ ​important​ ​to​ ​note​ ​here​ ​that​ ​the polygon​ ​coverage​ ​needs​ ​to​ ​be​ ​global​ ​and​ ​it​ ​needs​ ​to​ ​assign​ ​PlateIDs​ ​to​ ​the individual​ ​pieces​ ​of​ ​the​ ​raster​ ​in​ ​order​ ​to​ ​be​ ​able​ ​to​ ​rotate​ ​them​ ​back through​ ​time.​ ​In​ ​case​ ​you​ ​find​ ​this​ ​confusing,​ ​consult​ ​the​ ​“Rotations tutorial”.​ ​To​ ​cut​ ​the​ ​raster​ ​into​ ​different​ ​pieces​ ​do​ ​the​ ​following: 1.​ ​Make​ ​sure​ ​your​ ​layers​ ​are​ ​in​ ​the​ ​right​ ​order​ ​with​ ​the​ ​raster​ ​images​ ​in​ ​the back​ ​and​ ​the​ ​vector​ ​data​ ​(polygons)​ ​on​ ​top.​ ​If​ ​this​ ​is​ ​not​ ​the​ ​case,​ ​drag​ ​the layers​ ​into​ ​the​ ​proper​ ​order. 2.​ ​Expand​ ​the​ ​image​ ​layer​ ​(either​ ​topography​ ​or​ ​gravity​ ​image)​ ​in​ ​the​ ​Layer window​ ​by​ ​clicking​ ​the​ ​little​ ​black​ ​triangle​ ​to​ ​the​ ​left​ ​in​ ​the​ ​coloured rectangle​ ​of​ ​the​ ​layer.

3.​ ​In​ ​the​ ​“Inputs”​ ​section​ ​of​ ​the​ ​layer,​ ​click​ ​the​ ​“Add​ ​new​ ​connection”​ ​button under​ ​"Reconstructed​ ​Polygons:"​ ​and​ ​select​ ​the​ ​static​ ​polygons​ ​file​ ​from​ ​the list​ ​(Figure​ ​2).

Figure​ ​2.​​ ​Adding​ ​a​ ​polygon​ ​connection​ ​to​ ​the​ ​Gravity​ ​raster​ ​(DNSC08GRA_6m.gpml). 4.​ ​Depending​ ​on​ ​your​ ​graphics​ ​card​ ​power,​ ​you​ ​will​ ​see​ ​that​ ​GPlates​ ​will need​ ​some​ ​little​ ​time​ ​to​ ​think​ ​before​ ​the​ ​main​ ​window​ ​becomes​ ​responsive again. 5.​ ​Now​ ​you​ ​should​ ​be​ ​ready​ ​to​ ​go​ ​and​ ​able​ ​to​ ​drag​ ​the​ ​time​ ​slider​ ​to​ ​a desired​ ​time​ ​(or​ ​punch​ ​in​ ​the​ ​numbers)​ ​to​ ​rotate​ ​your​ ​global​ ​raster​ ​data​ ​to paleo-positions.​ ​See​ ​Figure​ ​3​ ​for​ ​an​ ​example​ ​of​ ​the​ ​ETOPO1​ ​dataset​ ​rotated back​ ​to​ ​50​ ​Ma​ ​in​ ​an​ ​Australia-centric​ ​view.

6.​ ​If​ ​you​ ​would​ ​like​ ​to​ ​see​ ​only​ ​the​ ​raster​ ​data​ ​and​ ​not​ ​have​ ​the​ ​polygons superimposed,simply​ ​toggle​ ​the​ ​polygon​ ​visibility​ ​off​ ​in​ ​the​ ​Layer​ ​manager.

Figure​ ​3.​​ ​Raster​ ​data​ ​cut​ ​to​ ​polygons​ ​and​ ​rotated​ ​back​ ​to​ ​50​ ​Ma.​ ​Notice​ ​that​ ​GPlates​ ​has automatically​ ​removed​ ​polygons​ ​and​ ​raster​ ​data​ ​which​ ​did​ ​not​ ​exist​ ​at​ ​this​ ​time​ ​(using​ ​the FromAge​ ​and​ ​ToAge​ ​feature​ ​attributes).

Exercise​ ​2:​ ​ ​Age​ ​-based​ ​masking​ ​or​ ​raster​ ​data Masking​ ​a​ ​raster​ ​files​ ​simply​ ​means​ ​that​ ​the​ ​sections​ ​of​ ​the​ ​globe​ ​which​ ​you select​ ​will​ ​not​ ​show​ ​the​ ​raster​ ​changing​ ​through​ ​time​ ​eg.Could​ ​mask​ ​the continents,​ ​if​ ​we​ ​had​ ​a​ ​sea​ ​floor​ ​spreading​ ​rotation​ ​data​ ​as​ ​having​ ​seafloor spreading​ ​data​ ​on​ ​top​ ​of​ ​a​ ​continent​ ​would​ ​make​ ​little​ ​sense. In​ ​this​ ​next​ ​exercise​ ​we​ ​will​ ​be​ ​masking​ ​raster​ ​data​ ​based​ ​on​ ​age.​ ​The Appearance​ ​ages​ ​of​ ​the​ ​polygons​ ​in​ ​color_etopol1-ice_low.gpml-​ ​which​ ​in the​ ​last​ ​exercise​ ​was​ ​derived​ ​from​ ​the​ ​static​ ​polygon​ ​ ​set (Global_EarthByte_GPlates_PresentDay_StaticPlatePolygons_20100927.gpml )​ ​will​ ​be​ ​masked​ ​and​ ​the​ ​appearance​ ​age​ ​data​ ​from​ ​the​ ​seafloor​ ​age-grid​ ​( agegrid_6m)​ ​will​ ​be​ ​used​ ​instead.​ ​This​ ​gives​ ​a​ ​smoother​ ​growth​ ​ ​of​ ​the seafloor​ ​through​ ​time. 1.​ ​File​ ​→​ ​Import​ ​Raster​ ​→​ ​select​ ​agegrid_6m.nc​ ​(Fig.4).

Figure​ ​4.​​ ​Import​ ​raster​ ​dialogue.​ ​Chose​ ​the​ ​“age”​ ​as​ ​the​ ​raster​ ​band​ ​when​ ​loading​ ​age grids.

2.​ ​The​ ​age​ ​grid​ ​is​ ​now​ ​loaded​ ​in​ ​the​ ​Layer​ ​manager.​ ​If​ ​you​ ​expand​ ​the​ ​layer, by​ ​clicking​ ​the​ ​small​ ​black​ ​triangle​ ​to​ ​the​ ​left​ ​of​ ​the​ ​eye,​ ​you​ ​will​ ​see​ ​that GPlates​ ​recognises​ ​this​ ​raster​ ​as​ ​an​ ​age​ ​grid​ ​(Band:age). 3.​ ​You​ ​now​ ​need​ ​load​ ​the​ ​plate​ ​rotation​ ​file​ ​and​ ​the​ ​static​ ​polygon​ ​set​ ​as described​ ​above​ ​in​ ​Tutorial​ ​1​ ​(Sec.​ ​5.1)​ ​into​ ​the​ ​GPlates​ ​application.​ ​This​ ​can also​ ​be​ ​done​ ​by​ ​dragging​ ​and​ ​dropping​ ​the​ ​files​ ​into​ ​the​ ​main​ ​GPlates window. 4.​ ​In​ ​order​ ​to​ ​be​ ​able​ ​to​ ​rotate​ ​the​ ​raster​ ​data,​ ​you​ ​will​ ​again​ ​need​ ​to​ ​assign plate​ ​IDs​ ​to​ ​subset​ ​of​ ​the​ ​raster​ ​by​ ​connecting​ ​the​ ​DNSC08GRA​ ​6m reconstructable​ ​raster​ ​layer​ ​to​ ​the​ ​static​ ​polygon​ ​features,​ ​as​ ​in​ ​Tutorial​ ​1, Step​ ​3​ ​(Sec.​ ​5.1). 5.​ ​In​ ​addition​ ​to​ ​assigning​ ​a​ ​polygon​ ​“connection”,​ ​you​ ​will​ ​now​ ​also​ ​connect an​ ​age​ ​grid​ ​feature​ ​to​ ​the​ ​global​ ​gravity​ ​data.​ ​Click​ ​the​ ​"Add​ ​new​ ​connection button"​ ​below​ ​the​ ​“Age​ ​grid​ ​raster”​ ​heading​ ​in​ ​the​ ​Inputs​ ​subsection​ ​of​ ​the layer​ ​(Fig.​ ​5)​ ​and​ ​select​ ​agegrid_6m.nc​ ​from​ ​the​ ​list. 6.​ ​You​ ​should​ ​now​ ​have​ ​loaded: •​ ​a​ ​rotation​ ​file •​ ​a​ ​global​ ​raster​ ​file​ ​which​ ​has​ ​age​ ​grid​ ​and​ ​polygon​ ​input​ ​channels •​ ​a​ ​static​ ​polygon​ ​file •​ ​an​ ​age​ ​mask/age​ ​grid​ ​file 7.​ ​You​ ​should​ ​now​ ​be​ ​able​ ​to​ ​reconstruct​ ​rasters​ ​again​ ​back​ ​through​ ​time, but​ ​with​ ​the​ ​age-based​ ​masking​ ​functionality​ ​enabled.​ ​If​ ​you​ ​interrogate​ ​for example​ ​the​ ​South​ ​Atlantic,​ ​going​ ​back​ ​in​ ​time,​ ​you​ ​will​ ​see​ ​that​ ​the​ ​seafloor is​ ​succesively​ ​“eaten​ ​up”​ ​during​ ​the​ ​reconstruction​ ​at​ ​any​ ​timestep​ ​(Fig.​ ​6). The​ ​transitions​ ​are​ ​very​ ​smooth​ ​and​ ​not​ ​like​ ​the​ ​polygon-based disappearance​ ​as​ ​you​ ​have​ ​seen​ ​in​ ​Tutorial​ ​1. A​ ​quick​ ​explanation​ ​of​ ​why​ ​it​ ​is​ ​more​ ​smooth​ ​when​ ​attached​ ​to​ ​‘​ ​seafloor age-grid’​ ​is​ ​given​ ​in​ ​the​ ​appendix.a

Figure​ ​5.​ ​Layer​ ​manager​ ​with​ ​loaded​ ​age​ ​grid​ ​feature​ ​(bottom,​ ​indicated​ ​by​ ​turquoise rectangle).

Appendix In​ ​the​ ​first​ ​exercise​ ​of​ ​this​ ​tutorial​ ​we​ ​used​ ​the​ ​Static​ ​Plate​ ​polygons​ ​to assign​ ​plate​ ​polygons​ ​Id​ ​and​ ​their​ ​respective​ ​age​ ​of​ ​appearance​ ​to​ ​raster data​ ​ ​(​ ​color_etopol1-ice_low.gpml​ ​or​ ​ ​DNSC08GRA​ ​6m).​ ​However​ ​the​ ​plate polygons​ ​are​ ​large​ ​and​ ​age​ ​of​ ​appearance​ ​are​ ​far​ ​apart​ ​so​ ​that​ ​big​ ​grey​ ​gaps were​ ​left​ ​at​ ​the​ ​MOR​ ​(​ ​fig​ ​20a.)

20a.​Color_etop1_ice_low​ ​raster​ ​at​ ​27ma​ ​with​ ​no​ ​age​ ​based​ ​masking​ ​around​ ​MOR.​ ​Note​ ​the large​ ​grey​ ​gaps.

20b.​Color_etop1_ice_low​ ​raster​ ​at​ ​27ma​ ​with​ ​age​ ​based​ ​masking​ ​around​ ​MOR.

However​ ​in​ ​the​ ​second​ ​exercise​ ​the​ ​time​ ​of​ ​appearance​ ​was​ ​gathered​ ​from the​ ​age_grid6m​ ​file.​ ​Now​ ​each​ ​individual​ ​pixel​ ​in​ ​this​ ​file,​ ​instead​ ​of​ ​a​ ​large polygon,​ ​has​ ​a​ ​respective​ ​appearance​ ​time​ ​according​ ​to​ ​the​ ​colour​ ​ ​palette​ ​( Dark​ ​Red=​ ​150ma,​ ​Dark​ ​Blue=0ma.​ ​Thus​ ​the​ ​growth​ ​of​ ​the​ ​ocean​ ​floor​ ​is smooth(​ ​Fig.20b). Sample​ ​data​ ​When​ ​you​ ​download​ ​GPlates​ ​from​ ​http://www.gplates.org, some​ ​sample​ ​data​ ​is​ ​included​ ​in​ ​your​ ​download.​ ​On​ ​Windows,​ ​this​ ​will​ ​be available​ ​after​ ​the​ ​installation​ ​in​ ​the​ ​GPlates​ ​folder​ ​at​ ​C:\Program Files\GPlates\GPlates​ ​[version]\Sample​ ​data.​ ​For​ ​the​ ​Mac,​ ​the​ ​download​ ​will leave​ ​you​ ​with​ ​a​ ​disk​ ​image​ ​(*.dmg)​ ​file.​ ​Mount​ ​the​ ​file​ ​by​ ​double-clicking, drag​ ​the​ ​GPlates​ ​application​ ​bundle​ ​into​ ​the​ ​Applications​ ​folder.​ ​The​ ​sample data​ ​is​ ​included​ ​as​ ​directory​ ​(“sample-data”)​ ​in​ ​the​ ​top​ ​level​ ​of​ ​the​ ​disk image.

Rotating rasters and age-based masking of Raster data

tutorial will use the data included in the GPlates distribution in the Sample .... polygons are large and age of appearance are far apart so that big grey gaps.

1MB Sizes 1 Downloads 160 Views

Recommend Documents

Introduction to rasters and Time-dependent rasters
reconstructable datasets in order to analyse and investigate features in the geological .... with plate reconstruction software (GPlates) can help geoscientists to.

data masking pdf
There was a problem loading more pages. data masking pdf. data masking pdf. Open. Extract. Open with. Sign In. Main menu. Displaying data masking pdf.

pdf raster
File: Pdf raster. Download now. Click here if your download doesn't start automatically. Page 1 of 1. pdf raster. pdf raster. Open. Extract. Open with. Sign In.

Blinding or masking
investigators, care providers, outcome assessors, data collectors, data analysts, and any other trial staff. The term “single blind” indicates that only patients or investigators are unaware ... the outcomes are measured on dental casts, scraping

raster pdf to vector pdf
File: Raster pdf to vector pdf. Download now. Click here if your download doesn't start automatically. Page 1 of 1. raster pdf to vector pdf. raster pdf to vector pdf.

summary of Masking .pdf
Page 1 of 1. summary of Masking .pdf. summary of Masking .pdf. Open. Extract. Open with. Sign In. Main menu. Displaying summary of Masking .pdf.

Stability and equilibrium of charged rotating drops
Abstract. We study by means of the Boundary Element Method (BEM) the effect that rotation at constant angular momentum L has on the evolution of a conducting and viscous drop when it holds an amount of charge Q on its surface or is immersed in an ext

Asymmetric pedestal masking of S-cone increments ...
The polarity of the sawtooth makes little difference even for strong pedestal conditions. - Contrast polarity seems to determine which pathway (ON or OFF) is processing the stimulus. - The temporal polarity has little effect, even at high contrast le

Design of Rotating Electrical Machines- Wiley.pdf
Design of Rotating Electrical Machines- Wiley.pdf. Design of Rotating Electrical Machines- Wiley.pdf. Open. Extract. Open with. Sign In. Main menu. Displaying ...

pdf to raster converter
Sign in. Page. 1. /. 1. Loading… Page 1 of 1. File: Pdf to raster converter. Download now. Click here if your download doesn't start automatically. Page 1 of 1.

Instabilities and vortex-lattice formation in rotating ... - Springer Link
namic studies of condensate solutions in the rotating frame [4] and ...... 0.5 0.6 0.7 0.8 0.9 1.0. 0.72. 0.68. 0.64. 0.60. 0.58. 100. 101. 102. 10–1. 10–2 εdd. Ωb/ω⊥.

Hopf bifurcation in wakes behind a rotating and ...
Guo-Hui Hu, De-Jun Sun, Xie-Yuan Yin, and Bing-Gang Tong. Department of Modern Mechanics, University of Science and Technology of China, Hefei, Anhui 230026 and Laboratory for Nonlinear Mechanics, Institute of Mechanics, Academia Sinica, Beijing 1000

Dynamo model with thermal convection and free-rotating ... - CiteSeerX
Canuto et al., 1988) for the solution of the Navier–Stokes equation ... tial feature of all these solutions was that the rotation of the ..... London A 456, 1669–1683.

A study of rotating globular clusters - Astronomy & Astrophysics
important for shaping the HB, as He-enriched stars reach hotter temperatures at ..... Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525. Searle, L.

Dynamo model with thermal convection and free-rotating ... - CiteSeerX
dromagnetic model with a free-rotating inner core powered by thermal convection. ..... high degree of uncertainty in the observations of the angular velocity.

EBNP,EBNB,EBNC Rotating Beacon.pdf
For IEC.60079-0 to 1 Standard Zone 1 and Zone 2 A QUALITY PRODUCT. Technical data. Type of Protection: Ex d IIA T6. Conforming to : Standards IEC ...

Breathing modes of a fast rotating Fermi gas
May 23, 2007 - 3Quantum Information Group, Institute for Scientific Interchange (ISI), Viale ... associated with the increase of the radial size of the cloud.

STABILITY OF FERROFLUID FLOW IN ROTATING ...
STABILITY OF FERROFLUID FLOW IN ROTATING. POROUS CYLINDERS WITH RADIAL FLOW. J. Singh, R. Bajaj. Centre for Advanced Study in Mathematics, ...