Split the mask computation and the area computation to avoid double code
and to keep the threshold in the stats.cxx only
Properly interpolate the scalar mask, the kend field was missing. Now
stats look smooth