Tach Leute!
Bin grad dabei ne Habitatmodellierung zu machen.
Über ne multiple logistische Regression habe ich nun meine Regressionskoeffizienten für verschiedenen Biotop-und Bodentypen sowie für Höhe und Hangneigung.
Hat jemand Erfahrung, wie diese nun elegant ins GIS gefüttert werden, um mir meine Vorkommenswahrscheinlichkeiten in einer Habitateignungskarte (aus den 4 vereinigten thematischen Karten) darzustellen?
Bin grad dabei ne Habitatmodellierung zu machen.
Über ne multiple logistische Regression habe ich nun meine Regressionskoeffizienten für verschiedenen Biotop-und Bodentypen sowie für Höhe und Hangneigung.
Hat jemand Erfahrung, wie diese nun elegant ins GIS gefüttert werden, um mir meine Vorkommenswahrscheinlichkeiten in einer Habitateignungskarte (aus den 4 vereinigten thematischen Karten) darzustellen?
- Anmelden oder Registieren, um Kommentare verfassen zu können
Gespeichert von MartinD am Fr., 06.03.2009 - 17:38
PermalinkVersuch ichs halt mal im nächsten Forum...
Tach Leute!
Bin grad dabei ne Habitatmodellierung zu machen.
Über ne multiple logistische Regression habe ich nun meine Regressionskoeffizienten für verschiedenen Biotop-und Bodentypen sowie für Höhe und Hangneigung.
Hat jemand Erfahrung, wie diese nun elegant ins GIS gefüttert werden, um mir meine Vorkommenswahrscheinlichkeiten in einer Habitateignungskarte (aus den 4 vereinigten thematischen Karten) darzustellen?
Gespeichert von Johannes Kamp am Di., 17.03.2009 - 15:54
Permalinksofern Du mit ArcView 3.x arbeitest, gibt es die Moeglichkeit, Regressionskoeffizienten multivariater Modelle auf Rasterdaten flaechig vorliegender Modellvariablen (=Habitatparameter) anzuwenden. Fuer 9.x gibt es sicher auch Moeglichkeiten, habe ich aber noch nie probiert.
Dazu kann man die Regressionsgleichung ueber den Map calculator (bei geladenem Spatial Analyst) direkt eingeben und die vier Variablen durch Klick auf das entsprechende shape im 'layer'-Menue auswaehlen, so sie als Rasterkarte vorliegen.
Das ist alles andere als komfortabel, da der Map calculator so seine Tuecken hat. Es gibt aber eine Schritt-fuer-Schritt-Anleitung, die ich Dir bei Bedarf mailen kann.
Beste Gruesse
Johannes
Gespeichert von MartinD am Do., 19.03.2009 - 15:12
Permalink...habe mein Problem -so glaube ich zumindest - so gelöst,
dass ich die verschiedenen thematischen Karten in ArcMap mithilfe von Select by Attributes und dem Field Calculator erstmal die jeweiligen Regressionskoeffiezienten zugewiesen habe,
die einzelnen Karten mit Analyst Tool - Overlay - Intersect zusammengelegt habe, und danach die einzelnen Regressionskoeffizienten wieder mit dem Field Calculator summiert habe...
...müsste doch auch so passen, oder??
Gespeichert von Johannes Kamp am Di., 07.04.2009 - 14:35
Permalinksorry fuer die spaete Antwort.
Der Knackpunkt liegt glaube ich darin, was Du unter 'Regressionskoeffizienten zugewiesen' verstehst.
Das Ziel ist ja, fuer jeden Punkt der Karte (i.d.R. jeden Rasterpunkt) eine Vorkommenswahrscheinlichkeit zwischen 0 und 1 zu erzeugen.
Nach Deiner mail gehe ich davor aus, dass Du ein multivariates logistisches Modell (i.e. GLM mit logit-link) gebastelt hast, dass vier Praediktorvariablen enthaelt, die Dir alle flaechendeckend als shapes vorliegen.
Daher muss die Zuweisung (=Modellgleichung) wie folgt aussehen:
f(y=1)= 1/(1+(e^([intercept]+[regkoeff1]*[predvar1]+...+[r egkoeff4]*[predvar4])))
Das mathematische 'e' ist im GIS meist 'exp'.
Wenn Du einfach nur die shapewerte mit dem Regressionskoeffizienten multiplizierst und das ganze dann verschneidest, fehlt dir der Rest der Formel der logistischen Regression (incl. intercept=Achsenabschnitt und saemtlicher Interaktionen), und Du duerftest kaum Werte zwischen 0 und 1 angezeigt bekommen, oder?
Beispiel: Du hast eine Karte mit Vegetationshoehen. Wenn Du diese (sagen wir 20-120 cm) einfach mit dem Regressionskoeff. multiplizierst (nehmen wir an, der sei beta=-0.25), bekommst Du eine Karte mit Werten zwischen -5 und -30 cm, aber noch keine Karte der Vorkommenswahrscheinlichkeit.
Beste Gruesse
Johannes
Gespeichert von MartinD am Di., 07.04.2009 - 17:10
Permalinkdass Du an mich denkst/mir antwortest!
Stimmt, hab mithilfe eines GLM mit link-logit die Regressionskoeffizienten ausspucken lassen!
Dann denn einzelnen themat. Karten diese Regressionskoeffizienten zugewiesen,
wobei ich mich (hoffentlich richtigerweise) dafür entschieden habe, den kategorialen Variablen (Bodentyp, Biotoptyp) pro Kategorie den jeweiligen RegKoeff zuzuweisen (ohne Multiplikation mit der jeweiligen Biotop-bzw.Bodentyp-Schlüsselnr...würde ja auch keinen Sinn machen) die numerischen Variablen (Höhe, Hangneigung, Expo) hingegen mit dem jeweiligen RegKoeff multipliziert habe.
Danach dann das "intersect-tool" bei dem ja nur alle die Punkte bleiben, die in allen Karten vorhanden sind (ohne die "missing values")???
In dem entstandenen Attribut-table in einem neuen Feld die RegKoeffs gemäss "Deiner" Formel dann verrechnet,
so komm ich dann auf Werte/Vorkommenswahrscheinlichkeiten zwischen 0 und 1...........die dann rein von "meiner Methodik" her hoffentlich korrekt sind??...oder
Gruss,
Martin
Gespeichert von Johannes Kamp am Mi., 08.04.2009 - 15:22
Permalinkdanke fuer die ergaenzenden Anmerkungen. Ich denke, das muesste so in Ordnung sein! Kann Dir jedenfalls bzgl. kategorialer Variablen und der weiteren Vorgehensweise folgen, und die Ergebniskarten scheinen ja auch sinnvoll zu sein.
Vielleicht ist es trotzdem ratsam, es nochmal einem Statistiker zu zeigen - aus der Ferne ist uebersieht man manchmal doch was.
Beste Gruesse,
Johannes