NAME: cutsegment
PURPOSE: cut p segments into p*n equal parts
CATEGORY: basic work
CALLING SEQUENCE: res = cutsegment(x0, y0, x1, y1, n)
INPUTS:
x0,y0 and x1,y1, 1d arrays of p elements, the coordinates of
the endpoints of the p segmements
n: the number of pieces we want to cut each segment
KEYWORD PARAMETERS:
/endpoints: see ouputs
/onsphere: to specify that the points are located on a
sphere. In this case, x and y corresponds to longitude and
latitude in degrees.
OUTPUTS:
defaut: a 3d array (2,n,p) that gives the coordinates of the
middle of the cutted segments.
if /endpoints, a 3d array (2,n+1,p) that gives the
coordinates of the endpoints of the cutted segments.
COMMON BLOCKS: no
SIDE EFFECTS: no
RESTRICTIONS: ?
EXAMPLE:
IDL> x0=[2,5]
IDL> y0=[5,1]
IDL> x1=[9,3]
IDL> y1=[1,8]
IDL> res=cutsegment(x0, y0, x1, y1, 10)
IDL> splot, [0,10], [0,10], xstyle = 1, ystyle = 1,/nodata
IDL> oplot, [x0[0], x1[0]], [y0[0], y1[0]]
IDL> oplot, [res[0,*,0]], [res[1,*,0]], color = 20, psym = 1, thick = 3
IDL> oplot, [x0[1], x1[1]], [y0[1], y1[1]]
IDL> oplot, [res[0,*,1]], [res[1,*,1]], color = 40, psym = 1, thick = 3
MODIFICATION HISTORY: