don't use log, very unclear
loop over channels to find if there are in one specific regions
use caching with joblib, really really fast\!