#Magnetic flux testing
# fluxpix = x.magnetic_flux(238, 238)
# flux = x.magnetic_flux(x.im_raw.data)
# print(flux[750,750])
#print(area.nansum())
#KPVT
# kpvt.heliographic(kpvt.im_raw.data)
# kpvt.los_corr(kpvt.im_raw.data)
#kpvt.eoa(kpvt.im_raw.data)
# kpvtflux = kpvt.magnetic_flux(kpvt.im_raw.data)
#print(np.nanmax(kpvtarea))
#print(kpvt.area.nansum())
#SPMG
# spmg.heliographic(spmg.im_raw.data)
# spmg.los_corr(spmg.im_raw.data)
# spmgarea = spmg.eoa(spmg.im_raw.data)
# spmg.magnetic_flux(spmg.im_raw.data)
# #print(np.nanmax(spmgarea))
# print(spmgarea.nansum())
#HMI
hmi.heliographic(hmi.im_raw.data)
#hmi.los_corr(hmi.im_raw.data)
hmiarea = hmi.eoa(hmi.im_raw.data)
#hmi.magnetic_flux(hmi.im_raw.data)
print(np.nansum(hmi.area))
stop = timeit.default_timer()
print("Time = ", stop - start)
from zaw_coord import CRD
import numpy as np
import matplotlib.mlab as mlab
import matplotlib.pyplot as plt
x = CRD('MDI\\fd_M_96m_01d.1222.0005.fits')
area = x.eoa(x.im_raw.data).value.flatten()
area = area[~np.isnan(area)]
loc = np.arange(0, 5e17, 5e17/999)
n, bins, patches = plt.hist(area, loc, log=True)
print(n.max())
plt.axis([0,5e17, 0, 60000])
l = plt.plot(bins, 'r--', linewidth=1)
plt.show()
from zaw_coord import CRD
import numpy as np
import matplotlib.mlab as mlab
import matplotlib.pyplot as plt
x = CRD('MDI\\fd_M_96m_01d.1222.0005.fits')
area = x.eoa(x.im_raw.data).value.flatten()
area = area[~np.isnan(area)]
loc = np.arange(0, 5e17, 5e17 / 999)
n, bins, patches = plt.hist(area, loc, log=True)
print(n.max())
plt.axis([0, 5e17, 0, 60000])
l = plt.plot(bins, 'r--', linewidth=1)
plt.show()
# fluxpix = x.magnetic_flux(238, 238)
# flux = x.magnetic_flux(x.im_raw.data)
# print(flux[750,750])
#print(area.nansum())
#KPVT
# kpvt.heliographic(kpvt.im_raw.data)
# kpvt.los_corr(kpvt.im_raw.data)
#kpvt.eoa(kpvt.im_raw.data)
# kpvtflux = kpvt.magnetic_flux(kpvt.im_raw.data)
#print(np.nanmax(kpvtarea))
#print(kpvt.area.nansum())
#SPMG
# spmg.heliographic(spmg.im_raw.data)
# spmg.los_corr(spmg.im_raw.data)
# spmgarea = spmg.eoa(spmg.im_raw.data)
# spmg.magnetic_flux(spmg.im_raw.data)
# #print(np.nanmax(spmgarea))
# print(spmgarea.nansum())
#HMI
hmi.heliographic(hmi.im_raw.data)
#hmi.los_corr(hmi.im_raw.data)
hmiarea = hmi.eoa(hmi.im_raw.data)
#hmi.magnetic_flux(hmi.im_raw.data)
print(np.nansum(hmi.area))
stop = timeit.default_timer()
print ("Time = ", stop - start)
