Re: setting regions for drawing polygons

[Ansley Manke <ansley@xxxxxxxxxxxxx>]

Hi Glenn,
The polymark script uses the K axis on which to generate the polygon
coordinates, and so a confusion arises with the use of regions in K.
We're looking at an improvement to the script to let the user optionally
specify an axis that's "Free" for the script to use for the axis of the polygon
coordinates.

By the way, before you answered your own question, I thought you were
talking about restricting the range of lat and lon -- just while we're on the
topic, they could be restricted by setting a mask for sigma:

   LET/BAD=999 sigma_restricted = IF(lon GE 200 AND lon LE 250 AND lat GE 0) THEN zpts ELSE 999
   GO polymark  POLY  lon, lat, sigma_restricted

Ansley


Glenn Carver wrote:

> Nothing like answering your own question.  8-)
>
> I neglected to mention that sigma is actually sigma[i=1:2562,k=1:40,l=1:400].
>
> If I set the region for k by set reg/k=1 and then call polymark, it doesn't work.
>
> Instead if I don't specify a region for k and then do:
>   polymark poly/key lon lat sigma[k=1] circle
>
> it works!  Quite what difference specifying the region this way makes I don't understand.
>
>  Glenn
>
> On 2002.02.07 16:28 Glenn Carver wrote:
> >
> > Is it possible to use a region with the polymark script (or just the polygon command)?
> >
> > I have a dataset which has three arrays, lon[i=1:2562], lat[i=1:2562] and a value to be plotted at those coords, sigma[i=1:2562]. I want to plot a marker at each (lon,lat) point coloured according to the values in sigma. Polymark.jnl does the job fine for the full lat & long range. But I can't get it or polygon to work if I only want to plot markers over a region.
> >
> > Unfortunately, although the sequence of points is ordered in some manner I can't simply select a subset of the data. Ideally what I need to do is use set region/x=200:250/y=-30:0 but this doesn't work. I get an error that the number of polygons does not match the data.
> >
> > Is there a way around this?
> >
> > Thanks
> >            Glenn
> >
> >
> > --
> >
> > Dr. Glenn Carver, Centre for Atmospheric Science, Univ. of Cambridge,
> > Chemistry Dept., Lensfield Road, Cambridge, CB2 1EW, UK.
> >
> > mailto:Glenn.Carver@atm.ch.cam.ac.uk  http://www.atm.ch.cam.ac.uk/~glenn/
> > Phone: +44 (1223) 763827      ** please note new numbers **
> > Fax:   +44 (1223) 763823
> >
> >
> --
>
> Dr. Glenn Carver, Centre for Atmospheric Science, Univ. of Cambridge,
> Chemistry Dept., Lensfield Road, Cambridge, CB2 1EW, UK.
>
> mailto:Glenn.Carver@atm.ch.cam.ac.uk  http://www.atm.ch.cam.ac.uk/~glenn/
> Phone: +44 (1223) 763827      ** please note new numbers **
> Fax:   +44 (1223) 763823

--
Ansley Manke  Pacific Marine Environmental Laboratory  Seattle WA  (206)526-6246