#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)