def test_xsect_radanalmean(afile,detector,level1,level2):
"""given filename pattern or product object, detector number (0..3), level(0..7) and species? (0..5)
plot out a contour map on an image
"""
import map_tools as mt
prod = open_product(afile)
D = cloud_mask_images(prod,detector)
im = D.GCcsRadAnalMean[:,:,level1,level2]
m,x,y = mt.mapshow(D.longitude, D.latitude, im,resolution='c' )
vars(D).update( { 'x_transform': x, 'y_transform': y, 'basemap': m, 'image': im } )
return D
def test_xsect_sounder(afile,detector_number,wavenumber_index):
"""show sounder wavenumber cross section
afile: a glob pattern, filename, or eugene Product object
detector_number: 0..3
wavenumber_index: 0..8699
returns dictionary of an assortment of useful things
"""
prod = open_product(afile)
import map_tools as mt
lon,lat,im = sounder_image_at_wavenumber_index( prod, detector_number, wavenumber_index )
m,x,y = mt.mapshow(lon,lat,im,resolution='c')
return { 'latitude': lat, 'longitude': lon, 'x_transform': x, 'y_transform': y, 'basemap': m, 'image': im }
def test_xsect_cloudmask(afile,detector_number,level_number,entity='GCcsRadAnalWgt'):
"""show cloudmask level cross section
afile: a glob pattern, filename, or eugene Product object
detector_number: 0..3
level_number: 0..6
returns dictionary of an assortment of useful things
"""
import map_tools as mt
prod = open_product(afile)
lon,lat = image_location(prod,detector_number)
D = cloud_mask_images(prod,detector_number)
im = getattr(D,entity)[:,:,level_number]
m,x,y = mt.mapshow(D.longitude,D.latitude,im,resolution='c')
vars(D).update( { 'x_transform': x, 'y_transform': y, 'basemap': m, 'image': im } )
return D
def test_xsect_radanalmean(afile, detector, level1, level2):
"""given filename pattern or product object, detector number (0..3), level(0..7) and species? (0..5)
plot out a contour map on an image
"""
import map_tools as mt
prod = open_product(afile)
D = cloud_mask_images(prod, detector)
im = D.GCcsRadAnalMean[:, :, level1, level2]
m, x, y = mt.mapshow(D.longitude, D.latitude, im, resolution='c')
vars(D).update({
'x_transform': x,
'y_transform': y,
'basemap': m,
'image': im
})
return D
def test_xsect_sounder(afile, detector_number, wavenumber_index):
"""show sounder wavenumber cross section
afile: a glob pattern, filename, or eugene Product object
detector_number: 0..3
wavenumber_index: 0..8699
returns dictionary of an assortment of useful things
"""
prod = open_product(afile)
import map_tools as mt
lon, lat, im = sounder_image_at_wavenumber_index(prod, detector_number,
wavenumber_index)
m, x, y = mt.mapshow(lon, lat, im, resolution='c')
return {
'latitude': lat,
'longitude': lon,
'x_transform': x,
'y_transform': y,
'basemap': m,
'image': im
}
def test_xsect_cloudmask(afile,
detector_number,
level_number,
entity='GCcsRadAnalWgt'):
"""show cloudmask level cross section
afile: a glob pattern, filename, or eugene Product object
detector_number: 0..3
level_number: 0..6
returns dictionary of an assortment of useful things
"""
import map_tools as mt
prod = open_product(afile)
lon, lat = image_location(prod, detector_number)
D = cloud_mask_images(prod, detector_number)
im = getattr(D, entity)[:, :, level_number]
m, x, y = mt.mapshow(D.longitude, D.latitude, im, resolution='c')
vars(D).update({
'x_transform': x,
'y_transform': y,
'basemap': m,
'image': im
})
return D
