Commit 90a2e24f authored by Arve Kylling's avatar Arve Kylling
Browse files

First version

parent bbb7da6b
This diff is collapsed.
Known problems:
import IASITools as IASITools
from glob import glob
import os
def ReadAndPlotSingleIASIFile(fn, Projection='lcc', PlotFile=None):
Scene = IASITools.Scene()
Scene.ReadBufr(fn, test=False, verbose=False)#True)#
Scene.IfCrossDateLineCorrectLongitude()
Scene.CalculateFootprint()
Scene.ConvertRadiance2BT()
# Scene.CalculateSEVIRIChannels()
# Type is plot type
# dBT is brightness temperature difference between two IASI channels
Scene.PlotLatLonVals(Type='dBT', #'BT', # 'SEVIRIdBT',
Projection=Projection,
# llcrnrlon=181,
# llcrnrlat=50,
# urcrnrlon=120,
# urcrnrlat=70,
verbose=False,
Title=False,
test=True, PlotFile=PlotFile)#,
return
#############################################################################################
if __name__ == "__main__":
path = '/home/aky/NILU/xnilu_wrk/sat_data/events/2011-02_SaharanDust/eumetcast/EPS-11/'
fs = path+"iasi_*_eps_o.l1_bufr"
# All file names
fns = sorted(glob(fs))
ii=0
for fn in fns:
fnbits = fn.split('/')
fnpng = fnbits[len(fnbits)-1]
fnpng = 'Plots/'+fnpng.replace('.l1_bufr','.png')
if os.path.exists(fnpng):
continue
lonmin, lonmax, latmin, latmax = IASITools.GetLatLonRegion(fn)
print lonmin, lonmax, latmin, latmax, fnpng
if latmin < -75:
Projection='spstere'
elif latmax > 75:
Projection='npstere'
else:
Projection='lcc'
ReadAndPlotSingleIASIFile(fn, Projection=Projection, PlotFile=fnpng)
# exit()
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment