Interpolation/
cutpar.pro
cut p parallelogram(s) into p*n^2 parallelograms
cutpar
basic work
result = cutpar(x0, y0, x1, y1, x2, y2, x3, y3, n, ENDPOINTS=ENDPOINTS, ONSPHERE=ONSPHERE)
Return value
- default: a 3d array(2,n^2,p) giving the center position of each
piece of the parallelograms
- if /ENDPOINTS : a 3d array(2,(n+1)^2,p) giving the edge positions
of each piece of the parallelograms
Parameters
x0
in
required
y0
in
required
x1
in
required
y1
in
required
x2
in
required
y2
in
required
x3
in
required
y3
in
required
1d arrays of p elements, giving the edge positions.
The edges must be given as in plot to draw the parallelogram. (see example).
n
in
required
each parallelogram will be cut in n^2 pieces
Keywords
ENDPOINTS
see outputs
ONSPHERE
to specify that the points are located on a
sphere. In this case, x and y corresponds to longitude and
latitude in degrees.
Examples
IDL> x0 = [2,6,2]
IDL> y0 = [0,2,6]
IDL> x1 = [3,8,4]
IDL> y1 = [4,4,6]
IDL> x2 = [1,6,4]
IDL> y2 = [5,6,8]
IDL> x3 = [0,4,2]
IDL> y3 = [1,4,8]
IDL> n = 4
IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata
IDL> for i=0,2 do oplot, [x0[i],x1[i],x2[i],x3[i],x0[i]],[y0[i],y1[i],y2[i],y3[i],y0[i]]
IDL> res=cutpar(x0, y0, x1, y1, x2, y2, x3, y3, n)
IDL> for i=0,2 do oplot, [res[0,*,i]], [res[1,*,i]], color = 20+10*i, psym = 1, thick = 3
Version history
Version
$Id: cutpar.pro 231 2007-03-19 17:15:51Z pinsard $
History
S. Masson (smasson@lodyc.jussieu.fr)
July 5th, 2002
Other attributes
Uses routines
cutsegment
Produced by IDLdoc 2.0.