Interpolation mittels griddata

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

Partner:

Forum

Option

[Erweitert]

• Diese Seite per Mail weiterempfehlen

Gehe zu:

RoggerThat

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 18.04.2016, 19:48 Titel: Interpolation mittels griddata

Hallo!

Ich habe folgendes Problem: ich möchte unregelmäßig verteilte Messdaten auf ein regelmäßiges Gitter interpolieren! Da wollte ich griddata verwenden! Geht zwar prinzipiell auch, aber die Dimensionen der Ergebnismatrix waren...nun ja, unerwartet ^^ Deshalb gehe ich stark davon aus dass ich den Befehl nicht richtig benutze!

Anbei mal wie ich es versucht habe:

Code:

%Testmatrix
data=rand(896,360);

%ich definiere mir die Grenzen meines Gitters:
x_start=-224;
x_end=224;

y_start=-224;
y_end=224;

%definiere mir die Schrittgröße meines "neuen" grids in x und y Richtung:
s_x=0.1;
s_y=0.1;

% definiere ein Gitter für x und y in der gewünschten Größe :
[newX, newY] = meshgrid([x_start:s_x:x_end],[y_start:s_y:y_end])

% evaluiere Daten auf meinem neuen Gitter:
newDataGrid = griddata(oldX, oldY, data, newX, newY);

Funktion ohne Link?

Die Matrizen mit oldX und oldY haben übrigens die selben Dimensionen wie meine data-Matrix.

Ich hätte jetzt offen gestanden als Ergebnis eine Matrix erwartet, die in etwa 4480x 4480 Einträge hat. Bekommen habe ich aber eine mit 2238x40 Einträgen. Tja. Offen gestanden weiß ich nicht so recht warum...

Ich hoffe sehr jemand kann mir da weiterhelfen, stehe gerade ziemlich auf der Leitung!

Lg RoggerThat

Winkow

Moderator


	Beiträge: 3.842

	Anmeldedatum: 04.11.11

	Wohnort: Dresden

	Version: R2014a 2015a

Verfasst am: 18.04.2016, 20:42 Titel:

da dein beispiel nicht lauffähig ist kann ich das nicht reproduzieren. bei mir klappts mit testdaten

Code:

%Testmatrix
data=rand(896,360);
%ich definiere mir die Grenzen meines Gitters:
x_start=-224;
x_end=224;

y_start=-224;
y_end=224;
[oldX,oldY]=meshgrid(linspace(x_start,x_end,360),linspace(y_start,y_end,896));

%definiere mir die Schrittgröße meines "neuen" grids in x und y Richtung:
s_x=0.1;
s_y=0.1;

% definiere ein Gitter für x und y in der gewünschten Größe :
[newX, newY] = meshgrid([x_start:s_x:x_end],[y_start:s_y:y_end]) ;

% evaluiere Daten auf meinem neuen Gitter:
newDataGrid = griddata(oldX, oldY, data, newX, newY);

Funktion ohne Link?

_________________

richtig Fragen

Harald

Forum-Meister


	Beiträge: 24.495

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 18.04.2016, 20:44 Titel:

Hallo,

was sind oldX und oldY?

Die Matrizen werden extrem groß, daher habe ich es mal mit

Code:

s_x=1;
s_y=1;

Funktion ohne Link?

versucht.

Dann liefert

Code:

newDataGrid = griddata(1:360, 1:896, data, newX, newY);

Funktion ohne Link?

ein newDataGrid derselben Größe wie newX und newY.

Grüße,
Harald

RoggerThat

Gast


	Beiträge: ---

	Anmeldedatum: ---

	Wohnort: ---

	Version: ---

Verfasst am: 18.04.2016, 22:16 Titel:

Hallo,

erst einmal vielen dank für eure Antworten!

Ich glaube ich habe da generell einen Denkfehler:

Also, ich habe 896 mal 360 Messwerte. Jeder dieser Messwerte hat eine x-Koordinate und eine y-Koordinate. Also habe ich folglich auch 2 Matrizen mit 896 mal 360 x-Werten (=oldX) und 896 mal 360 y-Werten (=oldY) .

Bei euren Beispielen habt ihr aber nur 360 x-Werte und 896 y-Werte die ihr einlest. Das verstehe ich nicht.

Tut mir leid, ist wahrscheinlich total offensichtlich..

Lg RoggerThat

Harald

Forum-Meister


	Beiträge: 24.495

	Anmeldedatum: 26.03.09

	Wohnort: Nähe München

	Version: ab 2017b

Verfasst am: 19.04.2016, 06:15 Titel:

Hallo,

das bedeutet, dass die Daten jeweils als Kombination für x- und y-Werte interpretiert werden. Stelle doch bitte dein oldX und oldY zur Verfügung, damit man das Problem nachvollziehen kann.

Grüße,
Harald


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.