Binary Maps/Removing Transient Features

[DavidCnly]

Next step in segmentation process.

Current plan:

• create binary map for time = t
• create binary map for time = t - delta_t
• dilate both maps by 10 pixels
• subtract binary maps (Map_t - Map_delta_t)
• remove transient features from original map
• dilate final resultant map by 10 pixels

[DavidCnly]

Binary function will require a dilation function 'dilate_region' or similar.
Changes needed to differential rotation function (reproject) before beginning.
Will write tests for existing functions before starting this.

[samaloney]

Sounds good. I thought you eroded the maps first and then dilated?

I'd suggests functions something like

create_binary_map(amap: Map, threshold: Quantity)

dilate_binary_map(binary: Map, radius: Quantity[u.arcsec])

remove_transient(current: Map, previous: Map)