(d2) Overlay of 3D Coordinates on Image or Contourplot

After you have done the stereoscopic height measurements of a set of curvilinear features (e.g., loops or filaments) you want to generate an overlay plot of the 3D coordinates on the image. Let us suppose you have stored the coaligned images in the array IMAGE_PAIR(2, NX, NY) in a savefile="loop_A.sav" (see tutorial a2) and the 3D coordinates in the output loopfile="loop_A.dat", you can generate an overlay greyscale/color plot with the following commands (for spacecraft A):

IDL>
savefile='loop_A.sav'
restore,savefile ;restore saved data from step (2a)
nsm =3 ;nsm=3, 5, ..., smoothing boxcar for highpass filtering
fov =[640,820,800,980] ;i1,j1,i2,j2 coordinates of field-of-view
hmax =0.1 ;e.g. hmax=0.1 or 0.5 solar radii (maximum height range)
loopfile ='loop_A.dat' ;output filename where 3D coordinates of loops are stored
filter =1 ;0=highpass filter, -1=no filter, linear scale, -2=no filter, logarithmic scale
window =0 ;IDL display window
ct =1 ;IDL color table
sc ='a' ;spacecraft 'a' or 'b'
cont =0 ;0=greyscale plot, 1=contour plot
gamma =0.2 ;slope of color table
inv =0. ;no inversion of color table
euvi_contour,image_pair,para,nsm,fov,hmax,loopfile,filter,window,sc,cont,dcont,ct,gamma,inv

The output on the screen will look like this:



A postscript file of the plot is also generated with the name over_col.ps


If you want a contourplot change to cont=1 and input contour levels (dcont)::
cont =1 ;0=greyscale plot, 1=contour plot
dcont =100. ;lowest (positive) contour level

gamma =1 ;slope of contour levels euvi_contour,image_pair,para,nsm,fov,hmax,loopfile,filter,window,sc,cont,dcont,ct,gamma,inv

A postscript file of the plot is also generated with the name cont_col.ps

The output on the screen will look like this: