[pyngl-talk] polyline draw order
haley at ucar.edu
Sat Apr 18 12:10:41 MDT 2015
Sorry for the delay in response. Ironically, the delay was caused by us
being at a Python conference this week at NCAR.
The problem with using Ngl.polyline, is that the lines are not "attached"
to the plot, and hence when you do the maximize_plot call later, the lines
are getting ignored and probably discarded.
There are a couple of things you can try:
- Move the #10 part of your code to *after* the #11 code, so that the
lines get drawn after you maximize the plot and draw the map.
- Updated to PyNGL 1.5.0-beta, and use the special "gsSegments" resource
which allows you to add shapefile outlines with a single call. So instead
of this code:
# 10. Add Natura2000 polylines to map
for i in range(0,numFeatures):
startSegment = geom[i, geom_segIndex]
numSegments = geom[i, geom_numSegs]
for seg in range(startSegment, startSegment+numSegments):
startPT = segm[seg, segs_xyzIndex]
endPT = startPT + segm[seg, segs_numPnts]# - 1
Ngl.polyline(wks, mp, LON[startPT:endPT], LAT[startPT:endPT],
It would look something like this::
LON1D = numpy.ravel(LON)
LAT1D = numpy.ravel(LAT)
natres.gsSegments = segm
id = Ngl.add_polyline(wks, mp, LON1D, LAT1D, natres)
On Mon, Apr 13, 2015 at 4:40 AM, Diego Bruciaferri <
diego.bruciaferri at ingv.it> wrote:
> Dear all,
> I'm trying to plot several polylines from a shapefile over a map which
> contains both scalar (oil concentration) and vector fields (current and
> stokes' drift fields).
> As I've seen in yours example scripts to plot shapefiles, I want to use
> the Ngl.polyline function over the Ngl.add_polyline because the first is
> faster, but it doesn't work (I don't see the lines) for me.
> I've tried to "play" with the DrawOrder resources, without any effect.
> Can you please suggest me some possible solution?
> Attached, you can find the code.
> Thanks in advance!
> Diego Bruciaferri
> Istituto Nazionale di Geofisica e Vulcanologia (INGV)
> National Institute of Geophysics and Volcanology
> Viale Aldo Moro 44
> 40127 Bologna
> Tel: +39 051 3782665
> pyngl-talk mailing list
> List instructions, subscriber options, unsubscribe:
-------------- next part --------------
An HTML attachment was scrubbed...
More information about the pyngl-talk