Resampling non-uniformly sampled 2D Surface

19 ビュー (過去 30 日間)
Ranjan Sonalkar
Ranjan Sonalkar 2020 年 12 月 24 日
コメント済み: Star Strider 2020 年 12 月 29 日
I have a 2D surface that is created from orthographic transformation of uniformly sampled surface in latitude/longitude coordinates. So, the transformed surface is non-uniformly sampled in orthographic coordinates (x, y). I would like to resample it to a specific uniform resolution in the orthographic axes.
Is there a function that can do this?
  4 件のコメント
Ranjan Sonalkar
Ranjan Sonalkar 2020 年 12 月 29 日
It looks like all I need is griddata, since I already have the desired uniformly spaced vectors, x, y at which I want the surface to be sampled at. x,y, are the nonuniformly sampled vectors and Z is the corresponding matrix of surface values. So, it appears that the following three lines of code will do the job.
[X, Y] = meshgrid(x, y);
[Xq,Yq] = meshgrid(x_i, y_i); % Create Interpolation Matrices With Uniform Grid Spacing
ZqNatural = griddata(X, Y, Z, Xq, Yq, 'natural'); % Interpolate Matrices To New Regular Grid - natural
Star Strider
Star Strider 2020 年 12 月 29 日
That will work!
There was a bit of a discrepancy between the title of your post: Resampling non-uniformly sampled 2D Surface and ‘... orthographic transformation of uniformly sampled surface ...’ so I went with the non-uniform situation in my Answer, because it will cover both possibilities.

サインインしてコメントする。

採用された回答

Star Strider
Star Strider 2020 年 12 月 24 日
Not a single function, however there are a group of functions that used in concert can do it.
Try this example:
x = sort(rand(1, 10))*10; % Latitude Vector
y = sort(rand(1, 12))*10; % Longitude Vector
[X,Y] = meshgrid(x, y); % Create Matrices
Z = sin(X*pi/2) .* cos(Y*pi/2); % Create Matrices
figure
surf(x, y, Z)
xlabel('x')
ylabel('y')
zlabel('Z')
title('Surface With Non-Uniform Grid Spacing')
xv = linspace(min(x), max(x), 40); % Create New Uniformly-Spaced Vector From Original ‘x’
yv = linspace(min(y), max(y), 50); % Create New Uniformly-Spaced Vector Vector From Original ‘y’
[Xq,Yq] = ndgrid(xv, yv); % Create Interpolation Matrices With Uniform Grid Spacing
Zq = griddata(X(:),Y(:),Z(:), Xq, Yq, 'natural'); % Interpolate Matrices To New Regular Grid
figure
surf(Xq, Yq, Zq)
xlabel('Xq')
ylabel('Yq')
zlabel('Zq')
title('Interpolated Surface With Uniform Grid Spacing')
.

その他の回答 (0 件)

カテゴリ

Help Center および File ExchangeSurface and Mesh Plots についてさらに検索

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by