More advanced example
From CrewWiki
(Difference between revisions)
| Line 6: | Line 6: | ||
o->set_group,'AWG'<br> | o->set_group,'AWG'<br> | ||
ctp1=o->get_data()<br> | ctp1=o->get_data()<br> | ||
| + | window,1 | ||
view2d,ctp1,/cool,/colo,no_data_idx=where(ctp1 le 0),coast_vec=coast,title='EUM ctp', min=200<br> | view2d,ctp1,/cool,/colo,no_data_idx=where(ctp1 le 0),coast_vec=coast,title='EUM ctp', min=200<br> | ||
longitude=o->lon()<br> | longitude=o->lon()<br> | ||
| Line 13: | Line 14: | ||
o->set_group,'GSF'<br> | o->set_group,'GSF'<br> | ||
ctp2=o->get_data()<br> | ctp2=o->get_data()<br> | ||
| + | window,2 | ||
view2d,ctp2,/cool,/colo,no_data_idx=where(ctp2 le 0),coast_vec=coast,title='GSF ctp',min=200<br> | view2d,ctp2,/cool,/colo,no_data_idx=where(ctp2 le 0),coast_vec=coast,title='GSF ctp',min=200<br> | ||
| + | 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-GSF',min=-150,max=150<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-GSF',min=-150,max=150<br> | ||
obj_destroy, o | obj_destroy, o | ||
</code> | </code> | ||
Revision as of 09:51, 14 February 2014
o=obj_new('cws_read',xrange=[1500,2800],yrange=[2800,3400],scale=0.3)
ls = o->land_sea()
coast=o->coast_line(/vector)
o->set_product,'ctp'
o->set_group,'AWG'
ctp1=o->get_data()
window,1
view2d,ctp1,/cool,/colo,no_data_idx=where(ctp1 le 0),coast_vec=coast,title='EUM ctp', min=200
longitude=o->lon()
latitude =o->lat()
rgb = o->rgb() & view2d,rgb,/cool,/no_axes
o->set_date,2008,6,13,12,0
o->set_group,'GSF'
ctp2=o->get_data()
window,2
view2d,ctp2,/cool,/colo,no_data_idx=where(ctp2 le 0),coast_vec=coast,title='GSF ctp',min=200
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-GSF',min=-150,max=150
obj_destroy, o
