from lunarspec import *

lst1=Glob('DATA/h*')

for it1 in lst1:
	it2 = str(it1)
	it=it2[5:]
	loadspec(it)

Figure('./raw-h8xxxx.png',
	'''
	sgfieldout fields=wavelength,reflectance,uncerntainty
	| sgsort key=group,id
	| sgplotps
	''')


Figure('./filter-h8xxxx.png', 
	'''
	sgfieldout fields=wavelength,reflectance
	| sgsort key=group,id
	| sggraphps 
	  x=wavelength y=reflectance reflectance.color=id
	''')