示例#1
0
    "hpcx_max": 1228.8,
    "hpcx_bin": 2.4,
    "hpcy_min": -1228.8,
    "hpcy_max": 1228.8,
    "hpcy_bin": 2.4
}



wave_maps = wave2d.simulate(waveparams, verbose = True)
for i in range(len(wave_maps)):
    wave_maps[i] = sunpy.map.Map.superpixel(wave_maps[i], (4,4), method = 'sum')

    #To get simulated HG' maps (centered at wave epicenter):
wave_mapsp_raw = wave2d.simulate_raw(waveparams)
wave_mapsp_raw_noise = wave2d.add_noise(waveparams, wave_mapsp_raw)
   
#visualize(wave_maps)
   
import util
   
new_wave_maps = []
   
for wave in wave_maps:
    print("Unraveling map at "+str(wave.date))
    new_wave_maps += [util.map_hpc_to_hg_rotate(wave, epi_lon = waveparams.get('epi_lon'), epi_lat = waveparams.get('epi_lat'), lon_bin = 5, lat_bin = 0.2)]


input_maps = new_wave_maps

#wave_maps = wave2d.simulate(params, verbose = True)
示例#2
0
    "lon_bin": 5.,
    
    #HPC grid, probably would only want to change the bin sizes
    "hpcx_min": -1228.8,
    "hpcx_max": 1228.8,
    "hpcx_bin": 2.4,
    "hpcy_min": -1228.8,
    "hpcy_max": 1228.8,
    "hpcy_bin": 2.4
}

wave_maps = wave2d.simulate(params, verbose = True)
    
    #To get simulated HG' maps (centered at wave epicenter):
wave_maps_raw = wave2d.simulate_raw(params)
wave_maps_raw_noise = wave2d.add_noise(params, wave_maps_raw)
    
visualize(wave_maps)
    
import util
    
new_wave_maps = []
    
for wave in wave_maps:
    print("Unraveling map at "+str(wave.date))
    new_wave_maps += [util.map_hpc_to_hg_rotate(wave, epi_lon = params.get('epi_lon'), epi_lat = params.get('epi_lat'), xbin = 5, ybin = 0.2)]

input_maps = new_wave_maps

#wave_maps = wave2d.simulate(params, verbose = True)
#visualize(wave_maps)
示例#3
0
    "lon_bin": 5.,
   
    #HPC grid, probably would only want to change the bin sizes
    "hpcx_min": -1228.8,
    "hpcx_max": 1228.8,
    "hpcx_bin": 2.4,
    "hpcy_min": -1228.8,
    "hpcy_max": 1228.8,
    "hpcy_bin": 2.4
}

wave_maps = wave2d.simulate(waveparams, verbose = True)

    #To get simulated HG' maps (centered at wave epicenter):
wave_maps_raw = wave2d.simulate_raw(waveparams)
wave_maps_raw_noise = wave2d.add_noise(waveparams, wave_maps_raw)
   
#visualize(wave_maps)

import util
   
new_wave_maps = []
   
for wave in wave_maps:
    print("Unraveling map at "+str(wave.date))
    new_wave_maps += [util.map_hpc_to_hg_rotate(wave, epi_lon = waveparams.get('epi_lon'), epi_lat = waveparams.get('epi_lat'), lon_bin = 5, lat_bin = 0.2)]

input_maps = new_wave_maps

#wave_maps = wave2d.simulate(params, verbose = True)
#visualize(wave_maps)
示例#4
0
    "lon_bin": 5.,

    #HPC grid, probably would only want to change the bin sizes
    "hpcx_min": -1228.8,
    "hpcx_max": 1228.8,
    "hpcx_bin": 2.4,
    "hpcy_min": -1228.8,
    "hpcy_max": 1228.8,
    "hpcy_bin": 2.4
}

wave_maps = wave2d.simulate(params, verbose=True)

#To get simulated HG' maps (centered at wave epicenter):
wave_maps_raw = wave2d.simulate_raw(params)
wave_maps_raw_noise = wave2d.add_noise(params, wave_maps_raw)

visualize(wave_maps)

import util

new_wave_maps = []

for wave in wave_maps:
    print("Unraveling map at " + str(wave.date))
    new_wave_maps += [
        util.map_hpc_to_hg_rotate(wave,
                                  epi_lon=params.get('epi_lon'),
                                  epi_lat=params.get('epi_lat'),
                                  xbin=5,
                                  ybin=0.2)