[ncl-talk] Shapefile Problem with 6.2.1

Scott Capps scott at allvertum.com
Tue Nov 11 16:38:38 MST 2014


This issue has been solved.  I had NCL v6.1.2 still in the PATH 
environment variable.

Thank you,

Scott
On 11/11/2014 03:10 PM, Rick Brownrigg wrote:
> Hi,
>
> Is there a typo in this line:
>
>    plres at gsSegments  = f->segments ;(:,0)
>
> I don't think the semicolon should be there.  f->segments should be an 
> N x 2 array.
>
> HTH...
> Rick
>
>
> On Tue, Nov 11, 2014 at 4:00 PM, Scott Capps <scott at allvertum.com 
> <mailto:scott at allvertum.com>> wrote:
>
>     Greetings,
>
>     I have just upgraded to version 6.2.1
>     (ncl_ncarg-6.2.1.Linux_Debian7.6_x86_64_nodap_gcc472.tar.gz) and
>     am happy with the many enhancements.  I am especially interested
>     in the more efficient shapefile plotting abilities.  However, I am
>     getting strange output when trying the new shapefile methods. 
>     Attached is a ps file generated using the new method (Which has an
>     issue).  I can provide the shapefile (Major U.S. highways) as well.
>
>          ; THIS IS THE NEW METHOD WHICH PRODUCES STRANGE OUTPUT:
>          plres             = True           ; resources for polylines
>          plres at gsLineColor = "gray25"
>          plres at gsLineThicknessF     = 3
>          f =
>     addfile(ncl_dir+"shapefiles/intrastate_shapefiles/intrstat.shp",
>     "r") ; Open shapefile
>          plres at gsSegments  = f->segments ;(:,0)
>          poly0   = gsn_add_polyline(wks, plot, f->x, f->y, plres)
>
>          ;
>          ; THIS OLD METHOD WORKS BUT IS VERY SLOW:
>          ; f =
>     addfile(ncl_dir+"shapefiles/intrastate_shapefiles/intrstat.shp",
>     "r") ; Open shapefile
>          ; segments = f->segments
>          ; geometry = f->geometry
>          ; segsDims = dimsizes(segments)
>          ; geomDims = dimsizes(geometry)
>          ; geom_segIndex = f at geom_segIndex
>          ; geom_numSegs  = f at geom_numSegs
>          ; segs_xyzIndex = f at segs_xyzIndex
>          ; segs_numPnts  = f at segs_numPnts
>          ; numFeatures   = geomDims(0)
>          ; lines = new(segsDims(0),graphic)   ; array to hold polylines
>
>          ; plres             = True           ; resources for polylines
>          ; plres at gsLineColor = "gray25"
>          ; plres at gsLineThicknessF     = 3
>
>          ; lon    = f->x
>          ; lat    = f->y
>          ; segNum = 0       ; Counter for adding polylines
>          ; do i=0, numFeatures-1
>          ;   startSegment = geometry(i, geom_segIndex)
>          ;   numSegments  = geometry(i, geom_numSegs)
>          ;   do seg=startSegment, startSegment+numSegments-1
>          ;     startPT = segments(seg, segs_xyzIndex)
>          ;     endPT   = startPT + segments(seg, segs_numPnts) - 1
>          ;     lines(segNum) = gsn_add_polyline(wks, plot,
>     lon(startPT:endPT),  \
>          ;     lat(startPT:endPT), plres)
>          ;     segNum = segNum + 1
>          ;   end do
>          ; end do
>          ;
>
>     The full NCL script is attached.  I can provide sample data upon
>     request (invariant_d01.nc <http://invariant_d01.nc>, wrftest.nc
>     <http://wrftest.nc> and the shapefile).
>
>     Thank you,
>
>     Scott
>
>
>     _______________________________________________
>     ncl-talk mailing list
>     List instructions, subscriber options, unsubscribe:
>     http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20141111/00197adb/attachment.html 


More information about the ncl-talk mailing list