Improve mask calculation doing it iteratively to remove effect of stars
Remove unnecessary print of fill value
Read images as float64 to avoid problems or precission when calculating statistics of images