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
by 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 returns
ONSPHERE
to specify that the points are located on a
sphere. In this case, x and y correspond 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 372 2008-08-08 12:31:53Z pinsard $
History
S. Masson (smasson@lodyc.jussieu.fr)
- July 5th, 2002
Other attributes
Uses routines
cutsegment
Produced by IDLdoc 2.0.