This adds the functionality to exclude pixels flagged in band mask. To ensure all neighboring pixels are removed that may be saturated or anomalous, the band mask is dilated with iterations = 10 in spicy.ndimage.binary_dilation(). This worked ok for the two test cases but is not optimized for all cases.
Test cases: emit20230326t112533_o08508_s003 and emit20230326t112544_o08508_s003
Notes:
1) ghg_process.py called local_surface_control_simple which crashed for me. The corresponding line without _simple is uncommented here. Not sure how to resolve this.
2) The setting for iterations = 10 could be made an input parameter, but not sure if this is desired.
This adds the functionality to exclude pixels flagged in band mask. To ensure all neighboring pixels are removed that may be saturated or anomalous, the band mask is dilated with iterations = 10 in spicy.ndimage.binary_dilation(). This worked ok for the two test cases but is not optimized for all cases.
Test cases: emit20230326t112533_o08508_s003 and emit20230326t112544_o08508_s003
Notes: 1) ghg_process.py called local_surface_control_simple which crashed for me. The corresponding line without _simple is uncommented here. Not sure how to resolve this.
2) The setting for iterations = 10 could be made an input parameter, but not sure if this is desired.