<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>