def test_EOF(self): x = np.random.random((self.D.nt, 20, 30)) self.D.data = np.ma.array(x, mask=x != x) self.D.cell_area = np.ones_like(self.D.data[0, :, :]) E = EOF(self.D) r = E.reconstruct_data() c = E.get_correlation_matrix() E.get_eof_data_correlation() #~ E.plot_channnel_correlations(100000) #slow!! E.plot_eof_coefficients(None, all=True) E._calc_anomalies() E.plot_EOF(None, all=True)
def test_EOF(self): x = np.random.random((self.D.nt, 20, 30)) self.D.data = np.ma.array(x, mask=x != x) self.D.cell_area = np.ones_like(self.D.data[0,:,:]) E = EOF(self.D) r = E.reconstruct_data() c = E.get_correlation_matrix() E.get_eof_data_correlation() #~ E.plot_channnel_correlations(100000) #slow!! E.plot_eof_coefficients(None, all=True) E._calc_anomalies() E.plot_EOF(None, all=True)
""" This file is part of pyCMBS. (c) 2012-2014 For COPYING and LICENSE details, please refer to the file COPYRIGHT.md """ """ EOF analysis """ from pycmbs.examples import download import matplotlib.pyplot as plt from pycmbs.diagnostic import EOF plt.close('all') air = download.get_sample_file(name='air') air.label = 'air temperature' # calculate climatological mean clim = air.get_climatology(return_object=True) # calculate EOF based on climatology because of performance issues for this example. E = EOF(clim) E.plot_EOF([0,1], use_basemap=True) # map_plot argument can be used here #~ E.plot_EOF(0,show_coef=False, use_basemap=True) # map_plot argument can be used here #~ E.plot_EOF(1,show_coef=False, use_basemap=True) plt.show()
# -*- coding: utf-8 -*- """ This file is part of pyCMBS. (c) 2012-2014 For COPYING and LICENSE details, please refer to the file COPYRIGHT.md """ """ EOF analysis """ from pycmbs.examples import download import matplotlib.pyplot as plt from pycmbs.diagnostic import EOF plt.close('all') air = download.get_sample_file(name='air') air.label = 'air temperature' # calculate climatological mean clim = air.get_climatology(return_object=True) # calculate EOF based on climatology because of performance issues for this example. E = EOF(clim) E.plot_EOF([0, 1], use_basemap=True) # map_plot argument can be used here #~ E.plot_EOF(0,show_coef=False, use_basemap=True) # map_plot argument can be used here #~ E.plot_EOF(1,show_coef=False, use_basemap=True) plt.show()
# -*- coding: utf-8 -*- """ This file is part of pyCMBS. (c) 2012-2014 For COPYING and LICENSE details, please refer to the file COPYRIGHT.md """ """ EOF analysis """ from pycmbs.examples import download import matplotlib.pyplot as plt from pycmbs.diagnostic import EOF air = download.get_sample_file(name='air') air.label = 'air temperature' # calculate climatological mean clim = air.get_climatology(return_object=True) # calculate EOF based on climatology because of performance issues for this example. E = EOF(clim) E.plot_EOF(0,show_coef=True, use_basemap=True) # map_plot argument can be used here E.plot_EOF(1,show_coef=True, use_basemap=True) plt.show()