<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;">Hi Mary<div> Thank you very much for your useful information ! I will look at it.</div><div><br></div><div>Best wishes</div><div><br></div><div>Xiaoni</div><div><br><div><div>Le 7 déc. 2017 à 22:58, Mary Haley <<a href="mailto:haley@ucar.edu">haley@ucar.edu</a>> a écrit :</div><br class="Apple-interchange-newline"><blockquote type="cite"><div dir="ltr"><div class="gmail_default" style="font-size:small">Hi Xianoi,</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">I'm CC-ing <a href="mailto:pyngl-talk@ucar.edu" target="_blank">pyngl-talk@ucar.edu</a> on this question, since we are encouraging users to go through our email lists in order to get continued help from us. You need to subscribe to pyngl-talk before you can send questions to it:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://mailman.ucar.edu/mailman/listinfo/pyngl-talk" target="_blank">http://mailman.ucar.edu/mailma<wbr>n/listinfo/pyngl-talk</a><br></div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">There's no way to use NCL or PyNGL to query the high-res database for specific lat/lon locations of lakes or other borders for that matter.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">What I would recommend is using a shapefile instead, which you can read using PyNIO, NCL, or other Python tools. There are many free shapefiles out there, and it's easy to google for them.</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default" style="font-size:small">As an example, I googled "shapefiles for lakes", and came across this link:</div><div class="gmail_default" style="font-size:small"><br></div><div class="gmail_default"><a href="http://www.naturalearthdata.com/downloads/10m-physical-vectors/10m-lakes/" target="_blank">http://www.naturalearthdata.co<wbr>m/downloads/10m-physical-vecto<wbr>rs/10m-lakes/</a><br></div><div class="gmail_default"><br></div>I clicked on the "download lakes" button and it download the shapefiles and ancillary data to a directory called<div class="gmail_default" style="font-size:small;display:inline"> </div>"ne_10m_lakes".<div class="gmail_default" style="font-size:small;display:inline"> I was then able to use the attached PyNIO script to query the file and print out the min/max lat/lon for each segment for each "Lake" on the file. I created a separate script to plot the different types of features.</div><div><div class="gmail_default" style="font-size:small;display:inline"><br></div></div><div><div class="gmail_default" style="font-size:small;display:inline">--Mary</div></div><div><div class="gmail_default" style="font-size:small;display:inline"><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Dec 7, 2017 at 12:22 PM, Xiaoni Wang <span dir="ltr"><<a href="mailto:xiaoni.wang@obspm.fr" target="_blank">xiaoni.wang@obspm.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Mary<br>
Thanks for your replay. I changed to another projection method and it makes the figure bigger. It is OK now.<br>
I would like to find the border of some lakes (i.e., lat and lon). I am wondering if I can read the high-resolution coast database myself (which is used in Pyngl) in order to get such information ? If you could give some clue about how to read this database in python, I would appreciate !<br>
<br>
Best wishes<br>
<span class="HOEnZb"><font color="#888888"><br>
Xiaoni</font></span></blockquote></div><br></div>
<span><ne_shp_file.png></span><span><query_and_plot_ne_shp_file.py></span><span><query_ne_shp_file.py></span></blockquote></div><br></div></body></html>