I understand that you wish to unify the coordinate systems in the Arctic from two different datasets and calculate the difference between the velocity vectors in both datasets.
The following steps can be followed to achieve the required results:
Step-1: Interpolate the NSIDC data for the coordinates in the ECCO grid system using "interp2" function.
load('NSIDC_vs_ECCO.mat');
% Find the range of latitude and longitude
% Initialize interpolated NSIDC data
VelocityNSIDC_interpolated = zeros(90, 90, 31);
latVec = linspace(min(min(latitude)), max(max(latitude)), 361);
lonVec = linspace(min(min(longitude)), max(max(longitude)), 361);
[newLon, newLat] = meshgrid(lonVec, latVec);
for day = 1:31
% Interpolate each day's data
VelocityNSIDC_interpolated(:,:,day) = interp2(newLon, newLat, seaice_NSIDC(:,:,day), Y, X,'linear');
end
F = 5;
Step-2: Calculate the difference as required.
for day = 1:31
figure(1), clf; clc;
t = tiledlayout(1,2); t.TileSpacing = 'tight'; t.Padding = 'compact';
% NSIDC
nexttile
pcolor(X,Y,VelocityNSIDC_interpolated(:,:,day));
shading interp; hold on;
title('NSIDC Polar Pathfinder','FontSize',F);
% ECCO
nexttile
pcolor(X,Y,seaice_ECCO(:,:,day)); colorbar; clim([0 0.55]);
shading interp; hold on;
title('ECCO sea ice velocity','FontSize',F);
% sgtitle([datestr(NSIDC_datetime(day), 'yyyy-mm-dd')],'FontSize',F+5);
diff(:,:,day) = VelocityNSIDC_interpolated(:,:,day) - seaice_ECCO(:,:,day);
end
Step-3: Plot the difference data.
for day = 1:31
figure(1), clc;
t = tiledlayout(1,2); t.TileSpacing = 'tight'; t.Padding = 'compact';
% ECCO
nexttile
pcolor(X,Y,diff(:,:,day)); colorbar; clim([0 0.55]);
shading interp; hold on;
title('ECCO sea ice velocity','FontSize',F);
% sgtitle([datestr(NSIDC_datetime(day), 'yyyy-mm-dd')],'FontSize',F+5);
end
This process will provide the required result.
For more information related to the "interp2" function, please refer the following documentation link: https://www.mathworks.com/help/matlab/ref/interp2.html
Hope this helps!
Regards
Arun
