g.mapset mapset=zone location=local -c
r.recode input="interpolation@zone" output="inter_int" rules="D:\protei\data\detection\jeu_donne\SET\recode_inter.txt"
r.reclass input="inter_int@zone" output="reclass" rules="D:\protei\data\detection\jeu_donne\SET\reclass.txt"
r.to.vect input=reclass@zone output=mask feature=area
r.mapcalculator amap=reclass@zone bmap=inter_int@zone formula=A*B outfile=cone
r.mapcalculator amap=zone_impossible@zone bmap=cone@zone formula=A*B outfile=cone2
r.mapcalculator amap=zone_etude_mask@zone bmap=cone2@zone formula=A*B outfile=cone3
r.report map=cone3@zone units=c output=D:\protei\data\detection\jeu_donne\SET\stat_cone