Binary Maps/Removing Transient Features

TCDSolar / SMARTpy

SolarMonitor Active Region Tracking

0 stars 2 forks source link

Binary Maps/Removing Transient Features #4

Open DavidCnly opened 2 months ago

DavidCnly commented 2 months ago

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 commented 2 months ago

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 commented 2 months ago

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)