How to regrid 2 surface data to match each other

Question

Sarvesh 2024-8-14

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2145169-how-to-regrid-2-surface-data-to-match-each-other

评论： Sarvesh 2024-8-14

regrid_surface_data.mat

Hi,

Soo I have been going on for a bit with this issue. I have a set of wind data for the Australian region that I want to process. This data includes longitude, latitude, and wind speeds in both the x-direction and y-direction. Specifically, the x-direction data is represented by ulon, ulat, and uwind, while the y-direction data is represented by vlat, vlon, and vwind.

I want to calculate the resultant wind speed using the Pythagorean theorem. To do this, the dimensions of the u and v data need to be aligned, which is not the case currently.

How can I write a script to regrid the data so that the dimensions of the u variables match those of the v variables. And for a single lon/lat pair there should be a uwind and vwind.

The results should be stored as rlat, rlon, and rwind, where rlat and rlon are the regridded latitude and longitude, and rwind is the resultant wind speed.

thanks in advance.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Ronit 2024-8-14

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/2145169-how-to-regrid-2-surface-data-to-match-each-other#answer_1498764

在 MATLAB Online 中打开

regrid_surface_data.mat

Hi Sarvesh,

I understand you want to process a set of wind data, which includes latitude, longitude, and wind speeds in both the x-direction (u) and y-direction (v). To calculate the resultant wind speed using the Pythagorean theorem, the dimensions of the 'u' and 'v' data need to be aligned. Below is a MATLAB script that regrids the data so that the dimensions of the 'u' variables match those of the 'v' variables.

data = load("regrid_surface_data.mat");
% Create a common grid that spans the range of both ulat/vlat and ulon/vlon
commonLat = linspace(max(min(data.ulat), min(data.vlat)), min(max(data.ulat), max(data.vlat)), min(length(data.ulat), length(data.vlat)));
commonLon = linspace(max(min(data.ulon), min(data.vlon)), min(max(data.ulon), max(data.vlon)), min(length(data.ulon), length(data.vlon)));
% Create a common meshgrid
[CommonLon, CommonLat] = meshgrid(commonLon, commonLat);
% Store the results of common latitude and longitude
rlat = CommonLat(:,1);
rlon = CommonLon(1,:);

Now interpolate both 'uwind' and 'vwind' data onto the common grid using 'interp2' and then using the Pythagorean theorem to combine the 'u' and 'v' wind components, calculate the resultant wind speed.

% Interpolate uwind and vwind to the common grid
UwindInterp = interp2(data.ulon, data.ulat, data.uwind, rlon, rlat, 'linear');
VwindInterp = interp2(data.vlon, data.vlat, data.vwind, rlon, rlat, 'linear');
% Calculate the resultant wind speed
rwind = sqrt(UwindInterp.^2 + VwindInterp.^2);