Repository: climate Updated Branches: refs/heads/master 1c85945fa -> 9deb819f7
CLIMATE-861 - An example to plot aerosol optical depth climatology from NASA MISR instrument - example/draw_climatology_map_MISR_AOD.py has been added. Project: http://git-wip-us.apache.org/repos/asf/climate/repo Commit: http://git-wip-us.apache.org/repos/asf/climate/commit/992ab591 Tree: http://git-wip-us.apache.org/repos/asf/climate/tree/992ab591 Diff: http://git-wip-us.apache.org/repos/asf/climate/diff/992ab591 Branch: refs/heads/master Commit: 992ab591b8be01061340697750bf66d508adc301 Parents: 818784e Author: huikyole <huiky...@argo.jpl.nasa.gov> Authored: Wed Aug 24 16:28:54 2016 -0700 Committer: huikyole <huiky...@argo.jpl.nasa.gov> Committed: Wed Aug 24 16:28:54 2016 -0700 ---------------------------------------------------------------------- examples/draw_climatology_map_MISR_AOD.py | 34 ++++++++++++++++++++++++++ 1 file changed, 34 insertions(+) ---------------------------------------------------------------------- http://git-wip-us.apache.org/repos/asf/climate/blob/992ab591/examples/draw_climatology_map_MISR_AOD.py ---------------------------------------------------------------------- diff --git a/examples/draw_climatology_map_MISR_AOD.py b/examples/draw_climatology_map_MISR_AOD.py new file mode 100644 index 0000000..1af834d --- /dev/null +++ b/examples/draw_climatology_map_MISR_AOD.py @@ -0,0 +1,34 @@ +import ocw.dataset as ds +import ocw.data_source.local as local +import ocw.dataset_processor as dsp +import ocw.plotter as plotter + +import numpy as np +import numpy.ma as ma + + +''' data source: https://dx.doi.org/10.6084/m9.figshare.3753321.v1''' +dataset = local.load_file('/home/huikyole/climate/examples/AOD_monthly_2000-MAR_2016-FEB_from_MISR_L3_JOINT.nc', + 'nonabsorbing_ave') +''' Subset the data for East Asia''' +Bounds = ds.Bounds(lat_min=20, lat_max=57.7, lon_min=90, lon_max=150) +dataset = dsp.subset(dataset, Bounds) + +'''The original dataset includes nonabsorbing AOD values between March 2000 and February 2015. +dsp.temporal_subset will extract data in September-October-November.''' +dataset_SON = dsp.temporal_subset(dataset, month_start=9, month_end=11, average_each_year=True) + +ny, nx = dataset_SON.values.shape[1:] + +# multi-year mean aod +clim_aod = ma.zeros([3, ny, nx]) + +clim_aod[0,:] = ma.mean(dataset_SON.values, axis=0) # 16-year mean +clim_aod[1,:] = ma.mean(dataset_SON.values[-5:,:], axis=0) # the last 5-year mean +clim_aod[2,:] = dataset_SON.values[-1,:] # the last year's value + +# plot clim_aod (3 subplots) +plotter.draw_contour_map(clim_aod, dataset_SON.lats, dataset_SON.lons, + fname='nonabsorbing_AOD_clim_East_Asia_Sep-Nov', + gridshape=[1,3],subtitles=['2000-2015: 16 years','2011-2015: 5 years', '2015: 1 year'], + clevs=np.arange(21)*0.02)
