"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) 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
"lon_max": 180., "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)
"lon_bin": 1., #HPC grid, probably would only want to change the bin sizes "hpcx_min": -1228.8, "hpcx_max": 1228.8, "hpcx_bin": 1.6, "hpcy_min": -1228.8, "hpcy_max": 1228.8, "hpcy_bin": 1.6 } #wave_maps = wave2d.simulate(params) 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) 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["epi_lon"], epi_lat=params["epi_lat"], xbin=1, ybin=0.05) ]
"lon_max": 180., "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)