Le portail francophone de la géomatique

est ce que tu as l'extension spatial analyst ?
si oui, tu peux faire ce genre de calcul avec

c'est je crois la manière la plus simple ou sinon tu passes par du code

Merci bien mais en faite j'avais vraiment besoin d'extraire le point bas d'une coupe pour reconstruire le fil d'eau de ma rivière. mais avec un code j'extrait un rectangle de mon raster et une zone tampon de la forme de ce rectangle, je fais une stat zonal et je compare le résultat de ce minimum avec mon raster rectangle en entrée, et là où ils sont égaux  j'extrais le pixel. Aprés je transforme en point et je relis les points. ça me semble bien tiré par les cheveux mais ça marche.
Dans de nombreuses applications j'ai envie de comparer un pixel avec son voisin, je ne comprends pas pourquoi c'est pas possible, avec la calculatrice raster par exemple?

Bonjour

Je suis pas sur d'avoir compris exactement ce que tu veux mais tu peux essayer ca dans la calculatrice raster:

(zonalmin(([raster] * 0),[raster],DATA) == [raster]) * [raster]

Explication :
- ([raster] * 0) cree un raster uniforme correspondant a ton raster avec la valeur 0 partout
- zonalmin(([raster] * 0),[raster],DATA) cree un raster qui contient pour chaque pixel de la zone la valeur minimum de cette zone. Les zones etant definies comme les pixels de meme valeur dans le premier raster en entree donc ici une seule zone puisque raster uniforme avec valeur 0 partout (voir point precedent) -> resultat un raster uniforme avec pour chaque pixel la valeur mini de ton raster.
- zonalmin(([raster] * 0),[raster],DATA) == [raster], ici je compare le raster que je viens de creer avec mon raster d'origine et la ou ils sont egaux j'ai ma valeur minimum, ceci cree un nouveau raster ayant pour valeur 1 aux pixels correspondant a la valeur minimum et 0 ailleurs.
- (zonalmin(([raster] * 0),[raster],DATA) == [raster]) * [raster], pour finir je multiplie le resultat par mon raster de depart pour recuperer mes pixels mini avec ma valeur mini. Resultat un raster avec 2 classes : 0 et ta valeur mini pour tes pixels mini.

Il ne reste plus ensuite qu'a faire un reclassify en remplacant ta classe 0 par NoData et puis une conversion en points et le tour est joue.

Attention ca ne marche que si le mini n'est pas egal a 0 et si tes valeurs sont des nombres entiers

A+

Rodolphe

Dernière modification par freesnowpetrel (Tue 14 October 2008 13:44)

C'est un peu ce que je fais sauf que je ne maitrise pas la calculatrice raster et j'ai fais un model builder. C'est dur de se faire comprendre...désolé par avance. En faite j'ai résolu mon problème c'est juste que je pensais qu'il devait avoir plus simple (un outil minimum tout fait) mais apparement non. Aprés j'élargissais ma question en demandant si on peut comparer deux pixels voisins. Et cela parce que à chaque fois que je cherche à faire qlqch, mon premier reflexe c'est de penser : si ce pixel a une valeur plus grande que son voisin alors ... si il a une valeur plus petite alors... Vous voyez? Une condition if then ça existe?