Verfasst am: 26.06.2013, 13:56
Titel: Farbliche Darstellung für X,Y-Punkt mit Wert
Hallo,
ich habe Probleme eine gewünschte Grafik zu erzeugen.
Gegeben:
- ein Haufen Postleitzahlen (ca. 4000) in Deutschland
- Latitude und Longitude der Zentroide
- Farblich darzustellender Wert für jede PLZ
Aktuell bekomme ich zwar eine Grafik erstellt, allerdings habe ich Zweifel an der Richtigkeit. Zudem würde ich am liebsten das ganze noch georeferenziert haben. Sodass man größere Städte/Flüsse einzeichnen könnte. (Worlmaps?!)
Matlab schlägt mir auch triscatteredinterp statt griddata vor, allerdings bekomm ich das gar nicht zum Laufen.
Fig2 = figure;
hold on
axis xy
axis([013040750])% without lat and long imagesc(Z)% without lat and long
%axis([minlat maxlat minlong maxlong])
%imagesc(lat,long,Z) colorbar('location','southoutside') hold off
Verfasst am: 26.06.2013, 15:38
Titel: Re: Farbliche Darstellung für X,Y-Punkt mit Wert
Hallo qwerob,
Du hast erklärt, was gegeben ist, aber nicht, was gesucht ist.
Ich weiß nicht, was eine "Zentroide" und was ein "Farblich darzustellender Wert für jede PLZ" ist.
Was hast Du denn für Zweifel an der Richtigkeit? Da wir ja nich wissen, was Du genau erreichen möchtest, können wir dieses wichtige Detail nicht erraten.
Was hast Du denn mit "triscatteredinterp" genau versucht und was bedeutet "bekomm ich das gar nicht zum Laufen" genau? Fehlermeldungen oder unerwartete Ergebnisse?
Gruß, Jan
qwerob
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am: 27.06.2013, 11:55
Titel:
Naja die Werte stellen meinet wegen die Höhe über Null dar, sodass du danach nen Höhenprofil von Deutschland erhälst. Zentroide heißt schlicht der Flächenschwerpunkt des PLZ-Bezirks. Sind ja leider keine schön, regelmäßig angeordneten Flächen. Genau das stellt auch das größe Problem dar.
Imagesc will eine regelmäßige Matrix zum plotten und die erzeuge (interpoliere) ich gerade mit "griddata".
Letzteres habe ich lediglich versucht durch den von Matlab vorgeschlagenen Befehl zu ersetzen. Das funktioniert leider nicht.
(Die genaue Fehlermeldung kann ich posten, wenn ich wieder am Arbeitsplatz bin)
Bin ich vll auf dem Holzweg und es geht einfacher? zB ohne imagesc und mit worldmaps oder so?
wenn du eine Karte plotten möchtest, (was du ja möchtest) und das wohlmöglich noch georeferenziert haben möchtest, musst du dir auch immer über die Wahl der Projektion im Klaren sein.
Es gibt für dein Problem übrigens schon perfekt zusammengeschnittene Lösungen in der m_map toolbox (frei, zum Beispiel hier: http://www.eos.ubc.ca/~rich/map.html )
Wenn du die Daten korrekt vorliegen hast, solltest du (nach dem korrekten "installieren", d.h. die Toolbox zum Suchpfad hinzufügen) der Toolbox nur noch 2 Zeilen Code benötigen, um das graphisch darzustellen. Man muss ja das Rad nicht neu erfinden.
Klingt auf jeden Fall schonmal besser. Danke dir!!
Nach (sehr mühsamer) Installation der Toolbox verlangt aber auch pcolor eine Matrix und nicht einen Wert für ein X,Y-Paar.
Eieieiei, das war mir nicht klar. Eigentlich hätte es mir aber auffallen können, da pcolor ja Matrizen plotted. Tut mir Leid. Aber es ist ja nicht ausgeschlossen, dass du diese Toolbox trotzdem später für die Graphiken nutzen möchtest (-:.
Ich bin mir aber trotzdem sehr sicher, dass du Dein Problem damit lösen kannst, da ich eigentlich das gleiche schon mal gemacht habe. Und tatsächlich haben wir da auch erst interpoliert.
Code:
x = lat;
y = long;
z = data;
F = TriScatteredInterp(x,y,z);
% x,y Richtungen und Gittergröße (hier musst du vielleicht was anderes Einstellen, damit es schön wird)
tx = linspace(minlat-0.1,maxlat+0.1,100);
ty = linspace(minlong-0.1,maxlong)+0.1,100);
% Gritter erzeugen [mx,my] = meshgrid(tx,ty);
mz = F(mx,my);
Versuch mal bitte, ob du damit schon mal grunsätzlich so einen Plot bekommst, wie du ihn haben möchtest. Durch den Contourplot erhältst du dann Höhenlinien mit den interpolierten Werten zwischen den gegebenen Messwerten.
Viele Grüße,
Nras.
Edit: Übrigens kannst du natürlich auch andere Projektionen auswählen. Ich glaube, dass auch oft 'Mercator' benutzt wird, gerade wenn man nicht die ganze Erde darstellen will, sondern zum Beispiel nur Deutschland.
qwerob
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am: 28.06.2013, 11:22
Titel:
Super! Keine Fehlermeldung.
Jetzt weist mich matlab nur darauf hin, dass ich auch triscattered nicht verwenden soll, weil es in nem "future release" durch "scatteredInterpolant" ersetzt wird. Allerdings erzeugt das bei mir nur ein wildes Mandala.
Da triscattered ja aktuell noch funktioniert, belass ichs einfach ma dabei.
Was mich wundert ist, wenn ich F im variables editor ansehe, ist x ein 4000x2 double und y 4000x1.
lat und long sind aber beide schlicht 4000x1.
Nun zu der eigentlichen Krux: Wie bekomme ich Städte und Flüsse eingezeichnet? Köln und der Rhein wären mir wichtig.
für meinen Fall brauchte ich keine Flüsse/Städe oder Grenzen. Diese sind da aber sicher auch irgendwo drin zu finden. Eine Google-Suche nach "m_map river" führt einen auf das hier:
Zitat:
New in release 1.4e are
GSHHS Now comes with the WDB rivers and borders database, and links to them have been added via m_gshhs.m
Das heißt, dass da auf jeden Fall eine Datenbank zu existiert, die man nutzen kann. Da musst du dich aber nun wohl alleine durchwühlen. Du hast da auch sicher mehr Motivation zu als ich (-:.
Nras
qwerob
Gast
Beiträge: ---
Anmeldedatum: ---
Wohnort: ---
Version: ---
Verfasst am: 09.07.2013, 15:28
Titel:
Soooo.. habe endlich eine Lösung gefunden.
Vielleicht gibt es ja noch Leute, die Ähnliches versuchen und denen ich hiermit weiterhelfen kann:
Code:
lat = lat*180/pi; % Umrechnung in Grad
long = long*180/pi;
Du kannst Beiträge in dieses Forum schreiben. Du kannst auf Beiträge in diesem Forum antworten. Du kannst deine Beiträge in diesem Forum nicht bearbeiten. Du kannst deine Beiträge in diesem Forum nicht löschen. Du kannst an Umfragen in diesem Forum nicht mitmachen. Du kannst Dateien in diesem Forum posten Du kannst Dateien in diesem Forum herunterladen
MATLAB, Simulink, Stateflow, Handle Graphics, Real-Time Workshop, SimBiology, SimHydraulics, SimEvents, and xPC TargetBox are registered trademarks and The MathWorks, the L-shaped membrane logo, and Embedded MATLAB are trademarks of The MathWorks, Inc.
RSS