Esempi in Python per CRD.heliographic

Esempio n. 1

File: zaw_pol_flux.py Progetto: amunozj/PolarFlux

def read_file(date):
	# Use date to find filename

	mgnt = CRD(filename)
	mgnt.heliographic(mgnt.im_raw.data)
	mgnt.magnetic_flux(mgnt.im_raw.data)
	mgnt.magnetic_flux(mgnt.im_raw.data, raw_field=True)

Esempio n. 2

def read_file(date):
    # Use date to find filename

    mgnt = CRD(filename)
    mgnt.heliographic(mgnt.im_raw.data)
    mgnt.magnetic_flux(mgnt.im_raw.data)
    mgnt.magnetic_flux(mgnt.im_raw.data, raw_field=True)

Esempio n. 3

File: zaw_util.py Progetto: ZachWerginz/PolarFlux

def CRD_read(date, instr):
    try:
        filename = search_file(date, instr)
    except IOError:
        return -1

    print(filename)
    
    try:
        mgnt = CRD(filename)
    except:
        return -1
    mgnt.heliographic()    
    mgnt.magnetic_flux()
    mgnt.magnetic_flux(raw_field=True)
    mgnt.date = mgnt.im_raw.date
    mgnt.md = date2md(date, instr)

    return mgnt

Esempio n. 4

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

Esempio n. 5

File: test.py Progetto: amunozj/PolarFlux

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