Ministack limitations

[falkamelung]

The ministack feature as currently implemented has several limitations. First it creates annual ministacks starting in January which causes problems for reference images in early or late in the year as shown in the network below. Second, it has no options on how to connect the ministacks which may result in connections with long spatial or temporal baseline which risks additional unwrap errors to occur.

It may be desirable to connect ministacks with only one connection. A version of find_short_baselines.py that does this is here. The required code changes are listed below. This creates this network (early period is the same network as above, but different reference month)

diff find_short_baselines.py  /Users/famelung/code/rsmas_insar/tools/MiaplPy/src/miaplpy/find_short_baselines.py
191,195d190
< def find_short_tbaseline_pair(date_list, first_index):
<     # starting with second entry (i>0)
<     j = first_index - 1
<     pair = (date_list[first_index ], date_list[j])
<     return pair
197d191
<
205,206c199,200
<     l_ind = np.zeros(indices_first.shape, dtype=int)
<     l_ind[0:-1] = np.array(f_ind[1::]).astype(int)
---
>     l_ind = np.zeros(indices_first.shape, dtype=np.int)
>     l_ind[0:-1] = np.array(f_ind[1::]).astype(np.int)
225,227c219
<             #pairs.append(find_short_pbaseline_pair(bperp, date_list, ministack_size, l_ind[i]))
<             pairs.append(find_short_tbaseline_pair(date_list,f_ind[i]))
<             #pass
---
>             pairs.append(find_short_pbaseline_pair(bperp, date_list, ministack_size, l_ind[i]))
230,232c222,223
<         #pairs.append((date_list[ref_inds[i]], date_list[ref_inds[i+1]]))
<         #pairs.append((date_list[l_ind[i]-1], date_list[l_ind[i]]))
<         pass
---
>         pairs.append((date_list[ref_inds[i]], date_list[ref_inds[i+1]]))
>         pairs.append((date_list[l_ind[i]-1], date_list[l_ind[i]]))