<div dir="ltr"><div dir="ltr">Hiep, please keep the mailing list in all replies.<div><br></div><div>"Convert from netcdf file back to shapefile?" Can you be a little more specific about what kind of data in the netcdf file, and what you would want the shape file to look like?</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Nov 27, 2019 at 2:07 AM Hiep Duc <<a href="mailto:Hiep.Duc@environment.nsw.gov.au">Hiep.Duc@environment.nsw.gov.au</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
<div lang="EN-AU">
<div class="gmail-m_-1918979364039494199WordSection1">
<p class="MsoNormal"><span>Hi all,<u></u><u></u></span></p>
<p class="MsoNormal"><span><u></u> <u></u></span></p>
<p class="MsoNormal"><span>Is it possible to convert from netcdf file back to shapefile ?.
<u></u><u></u></span></p>
<p class="MsoNormal"><span>Or is this transformation one way only ?
<u></u><u></u></span></p>
<p class="MsoNormal"><span>Thanks<u></u><u></u></span></p>
<p class="MsoNormal"><span>Hiep<u></u><u></u></span></p>
<p class="MsoNormal"><span><u></u> <u></u></span></p>
<p class="MsoNormal"><b><span lang="EN-US">From:</span></b><span lang="EN-US"> ncl-talk <<a href="mailto:ncl-talk-bounces@ucar.edu" target="_blank">ncl-talk-bounces@ucar.edu</a>>
<b>On Behalf Of </b>Dave Allured - NOAA Affiliate via ncl-talk<br>
<b>Sent:</b> Wednesday, 27 November 2019 6:27 PM<br>
<b>To:</b> Amy Hendricks <<a href="mailto:ashendricks@alaska.edu" target="_blank">ashendricks@alaska.edu</a>><br>
<b>Cc:</b> <a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
<b>Subject:</b> Re: [ncl-talk] converting shapefile to netcdf<u></u><u></u></span></p>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<p class="MsoNormal">Correction, the name of this example is shapefiles_14_mask, not just shapefiles_14.<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<p class="MsoNormal">On Tue, Nov 26, 2019 at 8:57 PM Dave Allured - NOAA Affiliate <<a href="mailto:dave.allured@noaa.gov" target="_blank">dave.allured@noaa.gov</a>> wrote:<u></u><u></u></p>
</div>
<blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin-left:4.8pt;margin-right:0cm">
<div>
<div>
<p class="MsoNormal">Example shapefiles_14 is close to what you want. Please notice these two lines. Their function is to select a particular subset of all the state outlines in the example shapefile, then apply a binary mask to the base grid "tc_lev0" so
that only grid points over the listed states are retained. All grid points outside this list of states are set to missing. This demonstrates how to access a subset of a shape file, instead of all regions together,<br>
<br>
opt@shape_names = (/"Texas","New Mexico","Colorado","Kansas","Oklahoma"/)<br>
tc_mask = shapefile_mask_data(tc_lev0,shp_filename1,opt)<br>
<br>
For your application, use this masking function to select one region at a time. Start with a master grid with ERA5 coordinates, set to all zero or missing values. Loop over each region value or region number. Call the function shapefile_mask_data inside
the loop, selecting the grid points for only one region each time. U<span style="font-family:Arial,sans-serif">se the resulting binary mask to add the current region number to the master grid. When the loop is complete, write the master grid with all region
numbers to a Netcdf file.</span><u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<p class="MsoNormal"><u></u> <u></u></p>
<div>
<div>
<p class="MsoNormal">On Tue, Nov 26, 2019 at 3:51 PM Amy Hendricks via ncl-talk <<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a>> wrote:<u></u><u></u></p>
</div>
<blockquote style="border-top:none;border-right:none;border-bottom:none;border-left:1pt solid rgb(204,204,204);padding:0cm 0cm 0cm 6pt;margin-left:4.8pt;margin-right:0cm">
<div>
<p class="MsoNormal">Hello NCL geniuses,<u></u><u></u></p>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">I am trying to write a shapefile to netcdf. The shapefile consists of circumpolar Treshnikov regions, and I would like to create a netCDF file that has the region values like tresh(lon,lat) = [region# from shapefile] on the ERA5 grid. All
the examples I'm finding are for creating binary masks, and the ncl_convert2nc does not seem to work, at least it doesn't return anything useable. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Attached are the ERA5 netCDF file, Treshnikov shapefile, as well as the result of ncl_convert2nc. <u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Thank you in advance,<u></u><u></u></p>
</div>
<div>
<p class="MsoNormal"><u></u> <u></u></p>
</div>
<div>
<p class="MsoNormal">Amy</p></div></div></blockquote></div></div></blockquote></div></div></div></div>
</blockquote></div></div>