Satellitenswath auf Weltkarte darstellen

WICHTIG: Der Betrieb von goMatlab.de wird privat finanziert fortgesetzt. - Mehr Infos...

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

kilkenny

Forum-Anfänger


	Beiträge: 13

	Anmeldedatum: 19.01.11

	Wohnort: ---

	Version: ---

Verfasst am: 09.11.2011, 13:58 Titel: Satellitenswath auf Weltkarte darstellen

Hallo zusammen,

ich sitze jetzt schon seit Stunden an meinem Problem und komme einfach nicht weiter:
Meine Aufgabe besteht darin aus Satellitendaten die gemessenen Windgeschwindigkeiten in eine Weltkarte darzustellen.
Das ganze lässt sich auch wunderbar auslesen, sodass ich am Ende jeweils eine

Code:

Funktion ohne Link?

36*860 Matrix

Code:

Funktion ohne Link?

mit allen Werten bekomme.

Das wirkliche problem tauch jetzt auf, wenn ich die Daten plotten möchte:
Da es sich bei den Daten nur um einen Swatch handelt, das heißt der Satellit nur ein bestimmtes Gebiet überfliegt (u.a. auch Land) stimmt die Platzierung der Werte in der Matrix natürlich nicht mit der tats"achlichen Breite und L"ange überein sondern muss aus der Longitude- und Latitudematrix vor dem plotten erst ausgelesen werden.
Hier ich komme einfach nicht weiter!!!
Ans"atze mit contourm scheitern im Ansatz...
Wäre echt super, wenn ihr mir helfen könntet!

Anbei mein Code bis jetzt:

Code:

wind_file = h5read('S1L2B2011196_09574_09575.h5','/science_data/Wind_speed');
lon_file = h5read('S1L2B2011196_09574_09575.h5','/science_data/Longitude');
lat_file = h5read('S1L2B2011196_09574_09575.h5','/science_data/Latitude');

%initialisiere swath_karte
for mm=1:36 %Breite
for kk = 1:860 %Laenge
swath_karte(mm,kk,[1:3]) = 0;
end
end

for ii =1:36;
for jj = 1:860;

swath_karte(ii,jj,1)= floor(lon_file(ii,jj)*0.01); %Laenge
swath_karte(ii,jj,2)= floor(lat_file(ii,jj)*0.01); %Breite
swath_karte(ii,jj,3)= (wind_file(1,ii,jj)*0.01); %Windspeed

end
end
axesm('mapprojection','mercato', 'maplonlimit', [181 180], 'frame', 'on');
load coast
plotm(lat, long);
framem; gridm; mlabel; plabel;
% plot data on map
lat = linspace(-90,90,180);
long = linspace(0, 360,180);
for ii =1:36;
for jj = 1:86;
contourfm(swath_karte(ii,jj,1),swath_karte(ii,jj,2),swath_karte(ii,jj,3),'r');

end
end
framem; gridm; mlabel; plabel;
tightmap

Funktion ohne Link?

Als Fehler bekomme ich:

Code:

Warning: Contour lines not rendered for uniform Z.
> In mapdisp/private/parseContourInputs>checkZ at 73
In mapdisp/private/parseContourInputs at 47
In contourm at 98
In contourfm at 37
In h5_read_test at 27
Operation terminated by user during save__listener__
(line 16)

In contourm>addColorPropertyListeners (line 422)
axListener = addlistener(ax,'CLim','PostSet',
@setColormapsCallback);

In contourm (line 134)
addColorPropertyListeners(h)

In contourfm (line 37)
contourm(varargin{:},'Fill','on','DefaultLineColor','black');

Funktion ohne Link?

Viele Grüße und Danke schon einmal für eure Hilfe!!!!! Idea


Einstellungen und Berechtigungen
Beiträge der letzten Zeit anzeigen:	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

Impressum | Nutzungsbedingungen | Datenschutz | FAQ |

RSS

Hosted by:

Copyright © 2007 - 2025 goMatlab.de | Dies ist keine offizielle Website der Firma The Mathworks

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.