Surface_rec has a feature of located reconstruction (see
mask.hbf at Data description). With a boolean intersection of all of the previous indicators, we can create such a file and reconstruct only on this pixels.
This may not be the optimal strategy but it may be interesting to compare it with the initial one.
For the sake of keeping a maximum of informations on the surface, it is interesting to create not only the intersection mask we just talked about but also an identification map that will tell us for each pixel how many indicators were valid or not and which combinaison of them.