Le portail francophone de la géomatique

Bonjour,

Il n'est pas possible dans Arcgis d'éditer un raster comme ça.
Une solution (pas top) qui me vient à l'idée, est de convertir ton MNT en fichier ASCII (Conversion Tools -> From Raster -> Raster to ASCII).
Tu obtiens un fichier texte avec entête, chaque valeur dans le corps du fichier correspond  à un pixel de ton raster.
Le problème est si ton raster est grand pour trouver les valeurs que tu veux changer. Enfin, si tu y arrives tu fais la démarche inverse et convertis ton ASCII en raster.

Pierre

Dernière modification par Pierre Vernier (Fri 14 August 2009 05:48)

Hello,

Pourquoi ne pas utiliser les opérateurs spatiaux ? On pourrait par exemple utiliser un masque pour mettre les valeurs à modifier = 1, puis de multiplier le MNT par un raster ayant des valeurs sur ces zones pour obtenir un MNT modifié (en veillant à la résolution, l'emprise et la superposition à chaque fois, qui se règlent dans les propriété de l'environnement dans la boite d'outil).

Pour créer les zones "masquées", on peut multiplier les zones à modifier par 100 (ou 1000 si besoin) et utiliser un outil Booléen (dans Math/logique/) pour mettre ces zones à 1 ou même juste une reclassification simple d'ailleurs.

Ensuite, il suffira de multiplier le MNT "masqué" par un MNT constitué de 1 pour les zones qui ne changent pas et de valeurs pour les zones modifiées


Exemple avec des valeurs de pixel :

A. MNT initial (coupe d'une butte que l'on va aplanir sur les pixels 3,4,5 avec un Z = 14) :
9-10-12-15-12-10-9

B. Raster "nouveau" avec alti modifiées créé à partir des points topo nouveaux :
Rien-Rien-14-14-14-Rien-Rien

C. MNT nécessaire pour créer le masque, obtenu par reclassification du précédent, Rien =1 ; Autres =100
1-1-100-100-100-1-1

D. Multiplication A*C :
(9*1)-(10*1)-(12*100)-(15*100)-(12*100)-(10*1)-(9*1)
Résultat : 9-10-1200-1500-1200-10-9

E. Reclassification du raster D : Zones > 1000 = 1 :
9-10-1-1-1-10-9

F. Reclassification du raster B : NoData = 1
1-1-14-14-14-1-1

G. Multiplication E*F :
(9*1)-(10*1)-(1*14)-(1*14)-(1*14)-(10*1)-(9*1)
Résultat : 9-10-14-14-14-10-9

Peut être compliqué, mais au moins ça devrait marcher si on respecte bien la résolution, l'emprise et la superposition des rasters.

Il y a sûrement d'autres pistes, auxquelles je n'ai pas pensé.

Tiens par exemple, peut être avec l'opérateur "CON" => conditional qui permettrait de raccourcir tout ça à partir de la phase D (Si MNT_D > 1000, appliquer RASTER_TOPO (on va de la phase D à G direct).

Robin.

Bonjour et merci pour toutes ces réponses !
Je n'ai pas tout testé, mais la solution de bayle_o (passer en vecteur - point ou polygone - pour la modification et rebasculer en raster ensuite) fonctionne très bien et très simplement !!
bonne journée
Paul