g.mapset mapset=zone location=local -c
v.patch -e input=donne1@zone,donne@zone output=donne_joint
v.surf.idw input="donne_joint@zone" output="interpolation_join" layer=1 column="N_SV_H" npoints=200 power=2
r.mapcalculator amap=zone_proba@zone bmap=zone_proba1@zone formula=A*B outfile=zone_croi
r.to.vect input=zone_croi@zone output=zone_croi feature=area
v.to.rast input="zone_croi@zone" layer=1 type="area" output="mask_croi" use="val" value=1
r.mapcalc "inter_croi= mask_croi * interpolation_join"
r.report map=inter_croi@zone units=c output=D:\protei\data\detection\jeu_donne\SET\stat_interpol_croi