<<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

- 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 cutted 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 136 2006-07-10 15:20:19Z pinsard $

History

S. Masson (smasson@lodyc.jussieu.fr) July 5th, 2002

Other attributes

Uses routines

cutsegment.pro
Produced by IDLdoc 2.0.