More advanced example
From CrewWiki
(→dependency on latitude) |
(→dependency on viewing zenith angle) |
||
| Line 55: | Line 55: | ||
</code><br> | </code><br> | ||
[[File:viewing_zenith_angle.jpg|200px|cloud top pressure per viewing zenith angle]]<br> | [[File:viewing_zenith_angle.jpg|200px|cloud top pressure per viewing zenith angle]]<br> | ||
| + | <br> | ||
| + | |||
| + | ==Simple statistics as a table== | ||
| + | <code> | ||
| + | o=obj_new('cws_stats') <br> | ||
| + | o->set_product,'ctp'<br> | ||
| + | o->set_date,2008,6,13,12,00<br> | ||
| + | o->make_tables<br> | ||
| + | obj_destroy, o<br> | ||
| + | </code><br> | ||
| + | [[File:ctp_table.jpg|200px|cloud top pressure statistics]]<br> | ||
<br> | <br> | ||
Revision as of 14:01, 18 February 2014
Go back to CREW software package
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
; 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
