More advanced example

From CrewWiki

Revision as of 14:01, 18 February 2014 (view source) Hamann (Talk | contribs) (→dependency on viewing zenith angle) ← Older edit		Revision as of 14:09, 18 February 2014 (view source) Hamann (Talk | contribs) Newer edit →
Line 80:		Line 80:
	==Reading two datasets and have a look at the difference==		==Reading two datasets and have a look at the difference==
-	<code>
-	<nowiki>;</nowiki> define object <br>
-	o=obj_new('cws_read',xrange=[1500,2800],yrange=[2800,3400],scale=0.3)<br>
-
-	<nowiki>;</nowiki> define date and time<br>
-	o->set_date,2008,6,13,12,0<br>
-	<nowiki>;</nowiki> read coast vector<br>
-	coast=o->coast_line(/vector)<br>
-	<nowiki>;</nowiki> set cloud property<br>
-	o->set_product,'ctp'<br>
-	<nowiki>;</nowiki> set CREW algorithm<br>
-	o->set_group,'AWG'<br>
-	<nowiki>;</nowiki> read and display data<br>
-	ctp1=o->get_data()<br>
-	window,1<br>
-	view2d,ctp1,/cool,/colo,no_data_idx=where(ctp1 le 0),coast_vec=coast,title='AWG ctp', min=200<br>
-
-	<nowiki>;</nowiki> read and display cloud top pressure for another algorithm<br>
-	o->set_group,'CMS'<br>
-	ctp2=o->get_data()<br>
-	window,2<br>
-	view2d,ctp2,/cool,/colo,no_data_idx=where(ctp2 le 0),coast_vec=coast,title='CMS ctp',min=200<br>
-
-	<nowiki>;</nowiki> display the difference of the two datasets<br>
-	window,3<br>
-	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<br>
-
-	obj_destroy, o
-	</code>
-
-	[[File:difference_of_two_datasets.png|200px|difference of AWG and CMS cloud top pressure]]

---

Revision as of 14:09, 18 February 2014

Go back to CREW software package

Contents 1 Display a cloud property of one algorithm on the MSG disk 2 Display a cloud property of all algorithms on the MSG disk 3 1d histogram 4 dependency on latitude 5 dependency on viewing zenith angle 6 Simple statistics as a table 7 Thumbnails of all correlation plots 8 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


Display a cloud property of all algorithms on the MSG disk

o=obj_new('cws_stats')
o->set_product,'ctp'
o->set_date,2008,6,13,12,00
o->make_image, /no_track
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,'cth'
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


Simple statistics as a table

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


Thumbnails of all correlation plots

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


Reading two datasets and have a look at the difference