[GRASS-de] Fragen zu r.mapcalc

Thomas Zimmermann thomaszimmermann12 at web.de
Do Okt 9 10:21:53 CEST 2003


Hallo Grass-User!

Ich hätte da mal ein paar Fragen zu r.mapcalc, zu denen ich in keinem der Online-Beschreibungen Antworten fand.

1.: Wie kann man eine Rasterzelle durch ihre geogr. Länge bzw. Breite genau definieren, oder wie kann man diesen Wert in eine Formel mit einbauen?
       
zB: 
r.mapcalc>neuekarte=karte.r*("geographische Breite von")vegetation.r

"neuekarte=100+row()+col()" macht sowas ja für "y/x", wie aber bei Longitude/Latitude?
       
Konkret geht es u.A. darum, die geographische Entfernung zweier Punkte in units zu messen.


2.,und die wichtigste Frage: Kann man bei Nachbarschaftsbeziehungen die Pixelangaben durch Parameter ersetzen? 
Es soll beliebig weit in eine Richtung gesucht werden, um die genau nächste Rasterzelle zu finden, die eine bestimmte Bedingung erfüllt. 


(Hier: "0", und linker nachbar:"1" ZB.:

     r.mapcalc> neuekarte = eval (if((test.r>test.r[0,-1]),"a",null()),a=min(if((test.r<test.r[b,c]),abs(formel),null()),b=n*0.5,c=n*0.5,0<=n<=1000) #oder so ähnlich.

Ist dazu dieses "x()" und "y()" irgendwie geeignet? Hab sonst keine Anleitung gefunden. 

Gibt`s bei 1000 auch einen Wert für 'Unendlich'? Es soll eben beliebig weit bis zum Ende der Rasterkarte abgefahren und geprüft werden. Das Verhältnis von n=a+b stellt hierbei einfach den "Winkel" dar, der hier (0,5:0,5) 45 Grad betragen würde.


3.Kann man auch das Erstellungsdatum eines Bildes in die Rechnung integrieren? 

Gruß, und Danke für alle Antworten, auch wenn was oben so nicht geht!

Thomas Zimmermann


______________________________________________________________________________
Horoskop, Comics, VIPs, Wetter, Sport und Lotto im WEB.DE Screensaver1.2
Kostenlos downloaden: http://screensaver.web.de/?mc=021110