<<prev file | next file >>    view single page | view frames    summary: fields | routine    details: routine

Interpolation/

cutpar.pro

cut p parallelogram(s) into p*n^2 parallelograms

topcutpar 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

Statistics

McCabe cyclic 1
McCabe essential 1
McCabe modular design 1
Produced by IDLdoc 2.0.