Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Code Block
titleodb ls
lat@hdr         lon@hdr         obsvalue@body
84.559998       -44.100006      100180.000000
84.349998       -46.989990      100140.000000
83.610001       -89.299988      100380.000000
83.419998       -71.970001      100140.000000
83.300003       -69.040009      100090.000000
82.480003       -93.190002      100350.000000
82.449997       -170.309998     101510.000000
82.260002       -128.949997     99180.000000
81.449997       -144.850006     100180.000000
81.419998       -144.619995     100180.000000
81.400002       -145.550003     100250.000000
80.709999       -110.500000     100180.000000
80.320000       -179.860001     102969.992188
80.019997       -151.399994     100440.000000

  

In this example, we will ask Magics to load this ODB file, and plot the position of each observation using the simple marker. WE have to inform Magics about the name of the columns to use to find the latitude, and longitude information.

Image Added

Code Block
languagepy
titleSimple symbol plotting
# importing Magics module
from Magics.macro import *

# Setting of the output file name
output = output(output_formats=['png'],
                output_name_first_page_number='off',
                output_name="odb_step1")

# Background Coastlines
background = mcoast(
    map_coastline_sea_shade_colour='white',
    map_coastline_land_shade_colour='cream',
    map_grid='on',
    map_coastline_land_shade='on',
    map_coastline_sea_shade='on',
    map_label='on',
    map_coastline_colour='tan',
    )
# Import odb data
odb = odb_geopoints(odb_filename='geo.odb',
                    odb_latitude_variable='lat@hdr',
                    odb_longitude_variable='lon@hdr',
                    )
# Define the symbol plotting
symbol = msymb(symbol_type='marker',
               symbol_colour='navy',
               symbol_marker_index=3,
               symbol_height=0.4,
            )
# Add a title 
lines = ['Using odb...', 
         'select lat@hdr, lon@hdr, obsvalue@body where (source='ispdv2.2') and (varno=110)',
         ]

title = mtext(
    text_lines=lines,
    text_justification='left',
    text_font_size=0.7,
    text_colour='charcoal',
)

#Create the plot
plot(output, background, odb, symbol, title)