Beispiel #1
0
    partRads=chan31[:max_rows,:max_cols]
    partMask=maskout[:max_rows,:max_cols]
    partLand=landout[:max_rows,:max_cols]
    partChan1rad=chan1[:max_rows,:max_cols]
    partChan1ref=chan1ref[:max_rows,:max_cols]
    
    numlatbins=800
    numlonbins=600

    bin_lats=fastbin(south,north,numlatbins,-999,-888)
    bin_lons=fastbin(west,east,numlonbins,-999,-888)

    lon_centers=bin_lons.get_centers()
    lat_centers=bin_lats.get_centers()

    new_hist=fh.pyhist(partLats,partLons,bin_lats,bin_lons)
    chan31_grid=new_hist.get_mean(partRads)
    mask_grid=new_hist.get_mean(partMask)
    land_grid=new_hist.get_mean(partLand)
    chan1rad_grid=new_hist.get_mean(partChan1rad)
    chan1ref_grid=new_hist.get_mean(partChan1ref)

    
    fig1=plt.figure(1)
    fig1.clf()
    axis1=fig1.add_subplot(111)
    im=axis1.pcolormesh(lon_centers,lat_centers,chan31_grid,cmap=cmap,\
                      norm=the_norm)
    cb=plt.colorbar(im,extend='both')
    the_label=cb.ax.set_ylabel('radiances ($W\,m^{-2}\,{\mu m}^{-1}\,sr^{-1}$)',rotation=270)
    axis1.set_title('MODIS channel 31 radiances (lat/lon binned)')
Beispiel #2
0
    north = lat_grid[0,0]
    south = lat_grid[2399,0]
    east = lon_grid[2399,2399]
    west = lon_grid[0,0]

    #lat/lon bins
    bin_lats=fastbin(south,north,numlatbins,-999,-888)
    bin_lons=fastbin(west,east,numlonbins,-999,-888)

    #lat/lon centers
    lon_centers=bin_lons.get_centers()
    lat_centers=bin_lats.get_centers()
    lon_centers_lst.append(lon_centers)
    lat_centers_lst.append(lat_centers)

    new_hist=fh.pyhist(lat_grid,lon_grid,bin_lats,bin_lons)
    lat_lon_counts=new_hist.get_hist2d()

    dirname='plots'
    make_dir(dirname)
    granule_info='count'

    #Figure 1: Lat/lon bin count
    fig1 = plt.figure(figcount)
    fig1.clf()
    cmap=cm.RdBu_r
    cmap.set_over('y')
    cmap.set_under('k')
    vmin= 0
    vmax= 20
    the_norm=Normalize(vmin=vmin,vmax=vmax,clip=False)
Beispiel #3
0
    #
    lon_hit = np.logical_and(lon_centers > -120, lon_centers < -115)
    lon_indices = np.where(lon_hit)[0]
    lat_hit = np.logical_and(lat_centers > 20, lat_centers < 25)
    lat_indices = np.where(lat_hit)[0]
    sub_lons = lon_centers[lon_indices]
    sub_lats = lat_centers[lat_indices]
    sub_chan1ref = chan1ref_grid[lat_indices[0] : lat_indices[-1], lon_indices[0] : lon_indices[-1]]
    sub_chan31 = chan31_grid[lat_indices[0] : lat_indices[-1], lon_indices[0] : lon_indices[-1]]

    bin_chan1ref = fastbin(0.05, 0.6, 50.0, -999, -888)
    bin_chan31 = fastbin(2.0, 18.0, 50.0, -999, -888)

    chan1_centers = bin_chan1ref.get_centers()
    chan31_centers = bin_chan31.get_centers()
    the_hist = fh.pyhist(sub_chan31, sub_chan1ref, bin_chan31, bin_chan1ref)

    counts = the_hist.get_hist2d()

    cmap = cm.RdBu_r
    cmap.set_over("y")
    cmap.set_under("k")
    vmin = 10.0
    vmax = 400.0
    the_norm = Normalize(vmin=vmin, vmax=vmax, clip=False)
    counts = counts.astype(np.float32)
    ## hit= (counts == 0)
    ## counts[hit] = 1.e-3
    ## log_counts=np.log10(counts)

    fig1 = plt.figure(1)
Beispiel #4
0
import numpy as np
import fasthist as fh

the_length = 15
test = np.ones([the_length], dtype=np.float32)
new_hist = fh.pyhist(test)
print new_hist.get_data()
print new_hist.indata
Beispiel #5
0
depolMPL = crosspolvals/copolvals

depolvals = depolMPL/(depolMPL+1)

copol_mean = np.mean(copolvals)
copol_std = np.std(copolvals)

copol_min = copol_mean-copol_std
copol_max = copol_mean+copol_std

bin_copol=fastbin(0.,0.002,100.,-999,-888)
bin_depol=fastbin(0.,2.0,100.,-999,-888)

copol_centers=bin_copol.get_centers()
depol_centers=bin_depol.get_centers()
the_hist=fh.pyhist(depolvals,copolvals,bin_depol,bin_copol)

counts=the_hist.get_hist2d()

cmap=cm.RdBu_r
cmap.set_over('y')
cmap.set_under('k')
counts=counts.astype(np.float32)
vmin= 0.
vmax= 4
the_norm=Normalize(vmin=vmin,vmax=vmax,clip=False)
hit= (counts <= 0)
counts[hit] = 1.e-3
log_counts=np.log10(counts)

olddir = os.getcwd()