<div dir="ltr">Hi Rick,<div><br></div><div>Thank you so much for taking a look at the underlying code for me. I really appreciate it. I also understand the constraints that the pivot to Python might place on resolving this bug. For anyone else who might be seeking a quick open source solution for a "shapefile island" visualization, Spatial Features for R (<a href="https://github.com/r-spatial/sf">https://github.com/r-spatial/sf</a>) is able to properly display the interior feature (see attached figure). R script below:</div><div><br></div><div>****************************************************</div><div>library(sf)</div><div>library(ggplot2)<br></div><div><br></div><div>test <- st_read("***PATH TO SHAPEFILE***/test_shapefile.shp")</div><div>ggplot() + geom_sf(data = test, fill = "lightgreen")</div><div>****************************************************</div><div><br></div><div><div><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div>Best,</div><div>Adam</div></div></div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 29, 2019 at 10:47 PM Rick Brownrigg <<a href="mailto:brownrig@ucar.edu">brownrig@ucar.edu</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 dir="ltr"><div>Hi Adam,</div><div><br></div><div>Firstly, much thanks for providing the concise, nicely illustrated example of the problem!</div><div><br></div><div>After looking through the code, I can believe this is a bug, and unfortunately, not one easily fixed. The underlying code properly knows how to decode shapefile "features" that are comprised of multiple "segments" or closed loops. But in drawing, it treats each sub-polygon independently. This works well, for example, for a coastal nation with a mainland and many islands, but it does not work for the case you provided, one with an outer hull containing interior islands/regions. Alas, the code that gsn_add_shapefile_polygons() is based upon does not provide interface to specify multiple nested interior loops. We can file this as a github issue, but with our current prioritization, I can't promise a timely fix.</div><div><br></div><div>I wish I had a better response...</div><div><br></div><div>Rick</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, May 29, 2019 at 11:15 AM Adam Sisco <<a href="mailto:aws0006@uah.edu" target="_blank">aws0006@uah.edu</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 dir="ltr">Hello NCL users,<div><br></div><div>I have a question regarding the visualization of shapefiles using NCL. For the sake of this help request, I have constructed a sample shapefile with a simplified geometry to illustrate the issue I am having. The shapefile contains a polygon within the opening of larger polygon, but I'm not sure NCL's gsn_add_shapefile_polygons function is capturing the polygon-within-polygon nature of this shapefile. So that my issue can be reproduced, I have attached:</div><div><ol><li style="margin-left:15px">The sample shapefile I am attempting to plot</li><li style="margin-left:15px">The desired visualization of the shapefile</li><li style="margin-left:15px">The visualization that results from the gsn_add_shapefile_polygons function</li><li style="margin-left:15px">The NCL script used to plot the shapefile</li></ol></div><div>I suspect that a more nuanced approach may be needed to produce the proper visualization. If this is within the capabilities of NCL, I would appreciate any assistance you can offer. I am using NCL 6.4.0 on a Windows machine.</div><div><br></div><div>******************** <br></div><div>begin<br> wks = gsn_open_wks("x11","shapefile_testing")<br><br> res = True<br> res@gsnDraw = False<br> res@gsnFrame = False<br> res@gsnMaximize = True<br><br> res@mpMinLonF = -110<br> res@mpMaxLonF = -75<br> res@mpMinLatF = 25<br> res@mpMaxLatF = 55<br><br> plot = gsn_csm_map(wks,res)<br><br> polyRes = True<br> polyRes@gsEdgesOn = True<br><br> poly = gsn_add_shapefile_polygons(wks,plot,"***PATH TO SHAPEFILE***/test_shapefile.shp",polyRes)<br><br> draw(plot)<br> frame(wks)<br>end<br></div><div>********************</div><div><br></div><div>Regards,</div><div><br clear="all"><div><div dir="ltr" class="gmail-m_8221075066115862076gmail-m_9134443247958919851gmail-m_8397596732373598451gmail_signature"><div dir="ltr"><b>Adam Sisco</b><br>Research Associate<br>NASA MSFC IMPACT<br>University of Alabama in Huntsville</div><div dir="ltr"><br></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@ucar.edu" target="_blank">ncl-talk@ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="http://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">http://mailman.ucar.edu/mailman/listinfo/ncl-talk</a><br>
</blockquote></div>
</blockquote></div>