larrys_script module
Multi-resolution filtering of radio images technique described in Rudnick, 2002 https://iopscience.iop.org/article/10.1086/342499/pdf larry@umn.edu – please contact for assistance, as needed code below courtesy of Viral Parekh, SARAO , vparekh@ska.ac.za
technique creates a diffuse emission map, called “open”; for small scale features, filtered = original_map - open
pick a box size 3x the beam-size or 3x size of features you want to remove
open map has an offset zero level - determine it and correct
open map is at original resolution - units are the same in Jy/beam
open map will show sharp edges to diffuse regions, but is boxy; convolve for aesthetics