from matplotlib import use
use('agg')
import matplotlib.pyplot as plt
import numpy as np
from os import path
from temperature import TemperatureMap as tmap


#data_dir = path.expanduser('~/CoronaTemps/')
maps_dir = path.expanduser('~/CoronaTemps/')
data_dir = '/imaps/sspfs/archive/sdo/aia/fulldisk/data/'
thismap1 = tmap('2011-02-15', data_dir=data_dir, maps_dir=maps_dir)
thismap3 = tmap('2011-02-15', data_dir=data_dir, maps_dir=maps_dir, n_params=3, verbose=True)

thismap1.save()
print thismap3.shape
print np.nanmin(thismap3.data[..., 0]), np.nanmean(thismap3.data[..., 0]), np.nanmax(thismap3.data[..., 0])

fig = plt.figure(figsize=(14, 14))
thismap1.plot()
plt.colorbar()
plt.savefig(path.expanduser('~/CoronaTemps/1paramtemps'))
plt.close()

fig = plt.figure(figsize=(14, 14))
thismap3.plot()
plt.colorbar()
plt.savefig(path.expanduser('~/CoronaTemps/fullgausstemps'))
plt.close()

from matplotlib import use
use('agg')
import matplotlib.pyplot as plt
import numpy as np
from os import path
from temperature import TemperatureMap as tmap


#data_dir = path.expanduser('~/CoronaTemps/')
maps_dir = path.expanduser('~/CoronaTemps/')
data_dir = '/imaps/sspfs/archive/sdo/aia/fulldisk/data/'
#thismap1 = tmap('2011-02-15', data_dir=data_dir, maps_dir=maps_dir)
thismap3 = tmap('2011-02-15', data_dir=data_dir, maps_dir=maps_dir, n_params=3, verbose=True)

#thismap1.save()
print thismap3.shape
print np.nanmin(thismap3.data[..., 0]), np.nanmean(thismap3.data[..., 0]), np.nanmax(thismap3.data[..., 0])

"""fig = plt.figure(figsize=(14, 14))
thismap1.plot()
plt.colorbar()
plt.savefig(path.expanduser('~/CoronaTemps/1paramtemps'))
plt.close()"""

fig = plt.figure(figsize=(14, 14))
thismap3.plot()
plt.colorbar()
plt.savefig(path.expanduser('~/CoronaTemps/fullgausstemps'))
plt.close()

Comment out 1-param stuff