Combining Maps
Maps are combined using a minimum intensity merge method. Maps are originally created for each instrument image individually then merged. There can be some points of overlap between different instrument maps and this is resolved by taking the data with the minimum intensity from each point of overlap. Coronal Hole data is then chosen based off which data points are used to create the original EUV maps. This method ensures that resulting maps are more continuous at seams. We also use a cutoff mu value to limit limb data distortion. In merging regions of overlap, we use data with a mu value greater than the cutoff value. In areas without overlap, any data available is used (mu cutoff of 0).
Combine Maps Function
The combine maps function can be found here.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 |
|
- 1.)
map_list[ii].data[map_list[ii].mu < mu_cutoff] = map_list[ii].no_data_val
- for all pixels with mu < mu_cutoff, set intensity to no_data_val
- 2.1)
np.logical_or(np.logical_and(overlap[:, :, ii], mu_array[:, :, ii] >= mu_cut_over), np.logical_and( data_array[:, :, ii] != map_list[0].no_data_val, mu_array[:, :, ii] >= mu_cutoff))
- to determine "good indices" based of Caplan et. al. 2016:
- check for overlap, in areas of overlap choose data where mu > mu_cut_over
- in areas of no overlap, choose data where mu > mu_cutoff
- to determine "good indices" based of Caplan et. al. 2016:
- 2.2)
mu_array[:, :, ii] > (max_mu - del_mu)
- determines "good indices" based off the value of del_mu
- 3.)
data_array[np.logical_not(good_index)] = float_info.max, data_array[data_array == map_list[0].no_data_val] = float_info.max
- make poor mu pixels unusable to merge, make no_data_vals unusable to merge
- 4.)
map_index = np.argmin(data_array, axis=2)
- find minimum intensity of remaining pixels
- 5.)
keep_data = data_array[row_index, col_index, map_index],
keep_data = data_array[row_index, col_index, map_index]- choose data to use for the EUV and CHD map
- 6.)
psi_d_types.PsiMap
- create new PsiMap object for both EUV and CHD combined maps