<div dir="ltr"><div>You may have encountered a known bug whereby the X and Y coordinates from a shapefile are swapped with each other. We think this was caused by a change in an underlying library. Please see this report for details and a possible easy workaround. Let us know if that works.</div><div><br></div><div> <a href="https://github.com/NCAR/ncl/issues/176">https://github.com/NCAR/ncl/issues/176</a><br></div><div><br></div><div>What version of NCL are you using, and what was the source of your NCL installation?</div><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Apr 5, 2022 at 4:29 PM Vanessa Vincente via ncl-talk <<a href="mailto:ncl-talk@mailman.ucar.edu">ncl-talk@mailman.ucar.edu</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left-width:1px;border-left-style:solid;border-left-color:rgb(204,204,204);padding-left:1ex"><div dir="ltr">Hello,<div><br></div><div>I am having trouble plotting polygon data from a shapefile of a burn scar perimeter over a map of Colorado with counties. I am using the example shapefiles_15.ncl from the <a href="https://www.ncl.ucar.edu/Applications/shapefiles.shtml" target="_blank">ncl page</a>. I am able to successfully plot the state of Colorado and its counties, but not the burn scar perimeter. Burn scar shapefile can be found <a href="https://ftp.wildfire.gov/public/incident_specific_data/rocky_mtn/2018/SpringCreek/IR/20180715/PerimeterFileProvidedByIncident/" target="_blank">here</a>. Can you please help? Thanks in advance!</div><div><br></div><div>Below is a copy of my ncl script:</div><div><br></div><div>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<br>load "./shapefile_utils.ncl"<br><br>;----------------------------------------------------------------------<br>; This function creates a cylindrical equidistant map of Colorado.<br>;----------------------------------------------------------------------<br><br>function create_colorado_map(wks,res,draw_ncl_outlines)<br>local a, mpres<br>begin<br> mpres = res<br><br> mpres@gsnMaximize = True<br> mpres@gsnPaperOrientation = "portrait"<br><br> mpres@gsnDraw = False<br> mpres@gsnFrame = False<br><br> mpres@mpFillOn = False <br><br>;---Turn on fancier tickmark labels.<br> mpres@pmTickMarkDisplayMode = "Always"<br> mpres@tmXBLabelFontHeightF = 0.008 ; smaller tickmark labels<br><br>;---Zoom in on area of interest<br> mpres@mpLimitMode = "LatLon"<br> mpres@mpMinLatF = 37<br> mpres@mpMaxLatF = 41<br> mpres@mpMinLonF = -109.05<br> mpres@mpMaxLonF = -102.05<br> mpres@mpFillOn = False<br><br> if(draw_ncl_outlines) then<br> mpres@mpOutlineOn = True<br> mpres@mpOutlineBoundarySets = "AllBoundaries"<br> mpres@mpDataBaseVersion = "MediumRes"<br> mpres@mpDataSetName = "Earth..4" ; U.S. counties<br> else<br> mpres@mpOutlineOn = False<br> end if<br><br>;---Create map.<br> map = gsn_csm_map(wks,mpres)<br><br> return(map)<br>end<br><br>;--------------------------------------------------<br>; Main code<br>;--------------------------------------------------<br>begin<br> wtype = "png" ; send graphics to PNG file<br> wtype@wkWidth = 2000<br> wtype@wkHeight = 2000<br> wks = gsn_open_wks(wtype,"shapefiles")<br><br> ncl_version = get_ncl_version()<br><br>;---Create two maps of Colorado<br> res = True<br> res@tiMainFontHeightF = 0.015<br> <br> res@tiMainString = "Colorado counties - shapefile"<br> map_shp = create_colorado_map(wks,res,False)<br><br> res@tiMainString = "Colorado counties - NCL (version " + ncl_version + ")"<br> map_ncl = create_colorado_map(wks,res,True)<br><br>;---Add shapefiles to one of the maps<br> lnres = True<br> lnres@gsLineColor = "red"<br><br> plot_line_shp = gsn_add_shapefile_polylines(wks,map_shp,"./usa/gadm40_USA_2.shp",lnres)<br><br>;---Add a slightly transparent marker to both maps to show location of county updates<br> mkres = True<br> mkres@gsMarkerIndex = 16 ; filled dot<br> mkres@gsMarkerOpacityF = 0.5 ; make the marker half transparent<br> mkres@gsMarkerColor = "red"<br> mkres@gsMarkerSizeF = 50.<br><br> counties_lat_center = 39.8<br> counties_lon_center = -104.9<br> <br> plot_marker_ncl = gsn_add_polymarker(wks,map_ncl,counties_lon_center,counties_lat_center,mkres)<br> plot_marker_shp = gsn_add_polymarker(wks,map_shp,counties_lon_center,counties_lat_center,mkres)<br><br>print("Adding polygons to ncl map and shp map...")<br><br>filename = "Perimeter_20180715_SpringCreek_COHUX001313.shp"<br> <br> ;---Set some options for the polygons<br> bscar = True<br> bscar@gsLineColor = "Blue"<br> bscar@gsLineThicknessF = 3.0 ; 3x thickness<br> <br>id1 = gsn_add_shapefile_polygons(wks,map_ncl,filename,bscar)<br>id2 = gsn_add_shapefile_polygons(wks,map_shp,filename,bscar)<br><br>;---Panel both plots. Markers, lines, and polygons will be drawn.<br> pres = True<br> pres@gsnMaximize = True<br> gsn_panel(wks,(/map_ncl,map_shp/),(/1,2/),pres)<br> <br>end</div>-- <br><div dir="ltr"><div dir="ltr"><div><div dir="ltr"><div><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><font color="#000000">Vanessa Vincente<br>Associate Scientist<br>The COMET® Program<br>University Corporation for Atmospheric Research<br></font><br>*My working hours may differ from your own. Please do not feel obligated to reply outside of standard business hours.*<div dir="ltr"></div></div></div></div></div></div></div></div></div></div></div>
_______________________________________________<br>
ncl-talk mailing list<br>
<a href="mailto:ncl-talk@mailman.ucar.edu" target="_blank">ncl-talk@mailman.ucar.edu</a><br>
List instructions, subscriber options, unsubscribe:<br>
<a href="https://mailman.ucar.edu/mailman/listinfo/ncl-talk" rel="noreferrer" target="_blank">https://mailman.ucar.edu/mailman/listinfo/ncl-talk</a></blockquote></div></div>