View source for Plotting Spectra with pplot

{{DavinciWiki_Procedure|
Run through a variety of [[pplot|pplot()]] cases.  The commands can be run automatically by using [[pplot_tutorial|pplot_tutorial()]].
|
Setup all the data we need...lots of output will be printed (but has been omitted here)
 dv> xlab1 &#61; make_band(lab1&#61;1)
 dv> xtes &#61; make_band(tes&#61;1)
 dv> xmtes &#61; make_band(mtes&#61;1)
 dv> xthemis &#61; make_band(themis&#61;1)
 dv> xthemiscm &#61; make_band(themis&#61;1, cm&#61;1)
 dv> xastercm &#61; make_band(aster&#61;1, cm&#61;1)
 dv> xtimscm &#61; make_band(tims&#61;1, cm&#61;1)
 
 # mini-TES data
 dv> minites_emiss &#61; read($DV_EX+"/minites_cube.vicar")
 dv> array&#61;minites_emiss[8:11, 5:6]
 dv> xx&#61;array[1:4,1]
 dv> yy&#61;array[1,1:2]
 dv> a&#61;array[1,1]
 dv> b&#61;array[2,2]
 dv> c&#61;array[2,1]
 
 # tes data
 dv> tes_emiss &#61; read($DV_EX+"/TES_emissivity_example.vic")
 
 # lab data
 dv> asu &#61; read($DV_EX+"/ASU_minlib.hdf")
 dv> spectrum &#61; asu.data[1]
 dv> lab_array &#61; asu.data[1:5]
 dv> spec_struc &#61; read($DV_EX+"/lab_structure.hdf")
 
 # make simulated aster data from lab - wavenumber x-axis
 dv> asterdata &#61; i2i(asu.data, from&#61;'lab1', to&#61;'aster', cm&#61;1)
 # make simulated tims data from lab - wavenumber x-axis
 dv> timsdata &#61; i2i(asu.data, from&#61;'lab1', to&#61;'tims', cm&#61;1)
 # make simulated aster data from lab - wavelength x-axis
 dv> asterdata &#61; i2i(asu.data, from&#61;'lab1', to&#61;'aster')
 
 # make some simulated themis data - wavelength x-axis
 dv> emthemis_cube &#61; i2i(asu.data[1:3], from &#61; 'lab1', to&#61;'themis')
 
 # convert lab1 data to aster, tims, and themis with wavenumber x-axis
 dv> a1 &#61; i2i(asu.data[1], from&#61;'lab1', to&#61;'aster', cm&#61;1)
 dv> t1 &#61; i2i(asu.data[1], from&#61;'lab1', to&#61;'tims', cm&#61;1)
 dv> th1 &#61; i2i(asu.data[1], from&#61;'lab1', to&#61;'themis', cm&#61;1)

'''Each image corresponds to the listed command'''

Simplest possible case: 1 vector, no colors, no labels - like xplot
 dv> pplot(spectrum)
 
 [[Image:pplot2_tutorial1.png|500px]]

2. Simple case: 1 vector with user-input x-axis, no colors, no labels
 dv> pplot(spectrum, xaxis&#61;xlab1)
 
 [[Image:pplot2_tutorial2.png|500px]]

3. Next simplest case: 1 array, user-input x-axis, no colors, no labels - will cycle colors automatically
 dv> pplot(lab_array, xaxis&#61;xlab1)
 
 [[Image:pplot2_tutorial3.png|500px]]

4. Array of spectra, 1 color (for all), 1 label, x and y axis labels, key in specfic location  Note: Order doesn't matter for colors and labels)
 dv> pplot(lab_array, 'Lab Spectra', 5, xaxis&#61;xlab1, xlabel &#61; 'Wavenumber', ylabel &#61; 'Emissivity', key&#61; '1400, .91')
 
 [[Image:pplot2_tutorial4.png|500px]]

5. Array of spectra, 1 color (for all), 1 label, standardized lab plot
 dv> pplot(lab_array, 5, 'Lab Spectra', lab1&#61;1)
 ...............................................................................
 
 [[Image:pplot2_tutorial5.png|500px]]

6. 1 array, 1 label, no colors specified - will auto cycle through colors and add index to label.  Key in bottom right.  Standard Mini-TES plot
 dv> pplot(array, 'Mini-TES Array', key &#61; 'bottom right', mtes&#61;1)
 
 [[Image:pplot2_tutorial6.png|500px]]

7. 2 different data elements with different x axes. x, y limits,  manually make lower x-axis wavelength but linear in cm-1. The x-axis is in wavenumber, but want a plot that is in wavelength on the lower axis. line weight 5, font_size 26, x-axis, y-axis, and plot labels
 dv> pplot({asudata[1:2], asterdata[1:2]}, xaxis&#61;{xlab1, xastercm}, {'Lab', 'ASTER'}, key &#61; 'outside', \
 wave&#61;1, x1&#61;500, x2&#61;1500,  y1&#61;.7, y2&#61;1.05, lw &#61; 5, font_size&#61;26, xlabel&#61;'Wavenumber', \
 ylabel&#61;'Emissivity', plot_title &#61; 'Lab and ASTER Example')
 
 [[Image:pplot2_tutorial7.png|500px]]

8. 2 data elements (1 array & 1 vector), 2 colors, 2 labels - single color and single label for each data element
 dv> pplot({array, avg(array, xy)}, {1,2}, {'Spectra', 'Ave. Spectrum'}, xaxis&#61;xmtes, x1&#61;500, \
 x2&#61;1800, y1&#61;.7, y2&#61;1.02, cm&#61;1)
 
 [[Image:pplot2_tutorial8.png|500px]]

9. ASU spectral structure example
 dv> pplot(spec_struc)
 
 [[Image:pplot2_tutorial9.png|500px]]

10. Resample example - lab, mtes, tes, and themis spectra
 dv> pplot({asu.data[1], a1, t1, th1},{1,2,3,4}, {'Lab','ASTER', 'TIMS', 'THEMIS'}, key&#61;'1450,.96', \
 xaxis&#61;{xlab1, xacm, xtcm, xthcm}, x1&#61;500,  x2&#61;1800, y1&#61;.88, y2&#61;1.01)
 doing broadband
 .........................................................
 doing broadband
 .........................................................
 doing broadband
 .........................................................
 
 [[Image:pplot2_tutorial10.png|500px]]

11. write to ps and pdf files. TES default plot, wavenumber on bottom x axis, ps, pdf
 dv> pplot({tes_emiss[1:3,1]}, {'TES Array'}, key &#61; '900, .9', tes&#61;1, emiss&#61;1, cm&#61;1, \
 ps&#61;$TMPDIR+"/junk.ps", pdf&#61;1)
 Outputting postscript to file: ./junk.ps
 Wavelength on the top
 converting to pdf file ./junk.ps.pdf
 Using postscript line colors from file: /u/cedwards/.Xdefaults
 
 [[Image:pplot2_tutorial11.png|500px]]

'''Now I'll demonstrate the more esoteric functions and features'''

3 vectors, 3 labels, no colors - will auto cycle through colors
 dv> pplot({a,b,c}, {'V1', 'V2', 'V3'}, mtes&#61;1)
 
 [[Image:pplot2_tutorial12.png|500px]]

2 data elements, 2 colors, 2 labels in structures - single color and single label for each data element
 dv> plot({array, xx}, {1,2}, {'a', 'b'}, mtes&#61;1)
 
 [[Image:pplot2_tutorial13.png|500px]]

Just 1 data element with 4 vectors in x, 4 labels, 4 colors
 dv> pplot(array[1:4, 1], {1,3,5,6}, {'a', 'b', 'b', 'd'}, mtes&#61;1)
 
 [[Image:pplot2_tutorial14.png|500px]]

3 data elements, colors WERE specified, one color for each and 1 label for 'a', 'xx', 'yy'
 dv> pplot({array,xx, yy}, {1,3,7}, {'a','xx', 'yy'}, mtes&#61;1)
 
 [[Image:pplot2_tutorial15.png|500px]]

THEMIS default plot, wavelength on bottom x axis
 dv> pplot(emthemis_cube, 'THEMIS', key &#61; '1100, .96', themis&#61;1, emiss&#61;1)
 
 [[Image:pplot2_tutorial16.png|500px]]
|
*[[pplot]]
*[[make_band]]
*[[i2i]]
*[[read]]
|
*[[pplot]]
*[[plot]]
*[[xplot]]
*[[triplot]]
}}
[[category:Procedures]]