Download source and data
ODB Radiance Map Example
#Metview Macro # **************************** LICENSE START *********************************** # # Copyright 2019 ECMWF. This software is distributed under the terms # of the Apache License version 2.0. In applying this license, ECMWF does not # waive the privileges and immunities granted to it by virtue of its status as # an Intergovernmental Organization or submit itself to any jurisdiction. # # ***************************** LICENSE END ************************************ # # ------------------------------------------------------------------ # Tags: ODB, Map # Title: ODB satellite radiance map # Description: Demonstrates how to plot AMSU-A radiances from ODB # ------------------------------------------------------------------ # define channel channel=5 #5,6 # read ODB amsua_odb = read("amsua.odb") # filter channel and define plot data vis = odb_visualiser( odb_where : "vertco_reference_1@body = " & channel, odb_data : amsua_odb ) # find out value range f = odb_filter(odb_query : "select obsvalue@body as val where vertco_reference_1@body = " & channel, odb_data : amsua_odb) val = values(f,"val") v1 = minvalue(val) v2 = maxvalue(val) # symbol plotting (adjusted to the actual min and max) symbol_fixed = msymb( legend : "on", symbol_type : "marker", symbol_table_mode : "advanced", symbol_advanced_table_selection_type : "interval", symbol_advanced_table_min_value : v1, symbol_advanced_table_max_value : v2, symbol_advanced_table_interval : 5, symbol_advanced_table_max_level_colour : "red", symbol_advanced_table_min_level_colour : "blue", symbol_advanced_table_colour_direction : "clockwise", symbol_advanced_table_marker_list : 3, symbol_advanced_table_height_list : 0.15 ) # coastlines coast = mcoast( map_coastline_colour : "RGB(0.5,0.5,0.5)", map_coastline_resolution : "low", map_coastline_land_shade : "on", map_coastline_land_shade_colour : "RGB(0.21,0.21,0.21)", map_coastline_sea_shade : "on", map_coastline_sea_shade_colour : "RGB(0.67,0.69,0.69)", map_grid_colour : "RGB(0.38,0.37,0.37)" ) # legend with histogram legend = mlegend( legend_display_type : "histogram" ) # title title = mtext( text_line_count : 2, text_line_1 : "Date: 2010-12-22 Time: 00", text_line_2 : "Sat: METOP-A Sensor: AMSU-A Channel: " & channel & " Param: Tb [K]" ) # define the output plot file setoutput(pdf_output(output_name : 'odb_radiance')) # plot plot(coast, vis, symbol_fixed, legend, title)
ODB Radiance Map Example
# **************************** LICENSE START *********************************** # # Copyright 2019 ECMWF. This software is distributed under the terms # of the Apache License version 2.0. In applying this license, ECMWF does not # waive the privileges and immunities granted to it by virtue of its status as # an Intergovernmental Organization or submit itself to any jurisdiction. # # ***************************** LICENSE END ************************************ # import metview as mv # define channel channel = 5 #5,6 # read ODB amsua_odb = mv.read("amsua.odb") # filter channel and define plot data vis = mv.odb_visualiser( odb_where = "vertco_reference_1@body = " + str(channel), odb_data = amsua_odb ) # find out value range f = mv.odb_filter(odb_query = "select obsvalue@body as val where vertco_reference_1@body = " + str(channel), odb_data = amsua_odb) val = mv.values(f, "val") v1 = mv.minvalue(val) v2 = mv.maxvalue(val) # symbol plotting (adjusted to the actual min and max) symbol_fixed = mv.msymb( legend = "on", symbol_type = "marker", symbol_table_mode = "advanced", symbol_advanced_table_selection_type = "interval", symbol_advanced_table_min_value = v1, symbol_advanced_table_max_value = v2, symbol_advanced_table_interval = 5, symbol_advanced_table_max_level_colour = "red", symbol_advanced_table_min_level_colour = "blue", symbol_advanced_table_colour_direction = "clockwise", symbol_advanced_table_marker_list = 3, symbol_advanced_table_height_list = 0.15 ) # coastlines coast = mv.mcoast( map_coastline_colour = "RGB(0.5,0.5,0.5)", map_coastline_resolution = "low", map_coastline_land_shade = "on", map_coastline_land_shade_colour = "RGB(0.21,0.21,0.21)", map_coastline_sea_shade = "on", map_coastline_sea_shade_colour = "RGB(0.67,0.69,0.69)", map_grid_colour = "RGB(0.38,0.37,0.37)" ) # legend with histogram legend = mv.mlegend( legend_display_type = "histogram" ) # title title = mv.mtext( text_line_count = 2, text_line_1 = "Date: 2010-12-22 Time: 00", text_line_2 = "Sat: METOP-A Sensor: AMSU-A Channel: " + str(channel) + " Param: Tb [K]" ) # define the output plot file mv.setoutput(mv.pdf_output(output_name = 'odb_radiance')) # plot mv.plot(coast, vis, symbol_fixed, legend, title)