Function sharp_sepia 

fn sharp_sepia(
    field: &[f64],
    mask: &[u8],
    grid: &Grid,
    threshold: f64,
    radius: f64,
) -> (Vec<f64>, Vec<u8>)

Tian Liu’s SHARP.m as SEPIA runs it (steps 0-5 of the module docs).