More advanced example

From CrewWiki

Revision as of 15:23, 17 February 2014 (view source) Hamann (Talk | contribs) ← Older edit		Revision as of 15:59, 17 February 2014 (view source) Hamann (Talk | contribs) (→1d histogram) Newer edit →
Line 19:		Line 19:
	o->set_date,2008,6,13,12,00<br>		o->set_date,2008,6,13,12,00<br>
	o->hist1d<br>		o->hist1d<br>
		+	obj_destroy, o<br>
		+	</code><br>
		+	[[File:single_image.png|200px|cloud top pressure histogram]]<br>
		+	<br>
		+
		+	==dependency on latitude==
		+	<code>
		+	o=obj_new('cws_stats') <br>
		+	o->set_product,'ctp'<br>
		+	o->set_date,2008,6,13,12,00<br>
		+	o->dep<br>
		+	obj_destroy, o<br>
		+	</code><br>
		+	[[File:single_image.png|200px|cloud top pressure histogram]]<br>
		+	<br>
		+
		+	==dependency on viewing zenith angle==
		+	<code>
		+	o=obj_new('cws_stats') <br>
		+	o->set_product,'ctp'<br>
		+	o->set_date,2008,6,13,12,00<br>
		+	o->dep, /vza<br>
	obj_destroy, o<br>		obj_destroy, o<br>
	</code><br>		</code><br>

---

Revision as of 15:59, 17 February 2014

Go back to CREW software package

Contents 1 Display a cloud property of one algorithm on the MSG disk 2 1d histogram 3 dependency on latitude 4 dependency on viewing zenith angle 5 Reading two datasets and have a look at the difference

Display a cloud property of one algorithm on the MSG disk

o=obj_new('cws_stats')
o->set_product,'cth'
o->set_group,'CMS'
o->set_date,2008,6,13,12,15
o->make_image, /single
obj_destroy, o


1d histogram

o=obj_new('cws_stats')
o->set_product,'ctp'
o->set_date,2008,6,13,12,00
o->hist1d
obj_destroy, o


dependency on latitude

o=obj_new('cws_stats')
o->set_product,'ctp'
o->set_date,2008,6,13,12,00
o->dep
obj_destroy, o


dependency on viewing zenith angle

o=obj_new('cws_stats')
o->set_product,'ctp'
o->set_date,2008,6,13,12,00
o->dep, /vza
obj_destroy, o


Reading two datasets and have a look at the difference

; define object
o=obj_new('cws_read',xrange=[1500,2800],yrange=[2800,3400],scale=0.3)


; define date and time
o->set_date,2008,6,13,12,0
; read coast vector
coast=o->coast_line(/vector)
; set cloud property
o->set_product,'ctp'
; set CREW algorithm
o->set_group,'AWG'
; read and display data
ctp1=o->get_data()
window,1
view2d,ctp1,/cool,/colo,no_data_idx=where(ctp1 le 0),coast_vec=coast,title='AWG ctp', min=200


; read and display cloud top pressure for another algorithm
o->set_group,'CMS'
ctp2=o->get_data()
window,2
view2d,ctp2,/cool,/colo,no_data_idx=where(ctp2 le 0),coast_vec=coast,title='CMS ctp',min=200


; display the difference of the two datasets
window,3
view2d,ctp1-ctp2,/cool,/colo,no_data_idx=where((ctp1 le 0) or (ctp2 le 0) or (abs(ctp1-ctp2) gt 150)),coast_vec=coast,title='difference ctp AWG-CMS',min=-150,max=150


obj_destroy, o