how to merge two NetCDF4 files into one

Question

Sreeraj R 2021-10-11

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1560991-how-to-merge-two-netcdf4-files-into-one

评论： Sreeraj R 2021-10-11

Hi everyone, i am currently working on Sentinel 5P Tropomi data. the data is available in the format of NetCDF6. i am trying to merge two images into one. i need the lattitude, longitude and NO2 column data in the new merged file. is it possible to merge two files and have the merged file have all the varaiables which the two files have. if so please let me know how i can merge these two files in Matlab.

2 个评论
显示无隐藏无

KSSV 2021-10-11

Yes happily it can be merged.....shared your files so that we can have a look.

Sreeraj R 2021-10-11

编辑：KSSV 2021-10-11

yes sir....sorry there is a correction the data format is not NetCDF6 but NetCDF4 (not NetCDf4_classic) . i have attched the google drive link for the data

https://drive.google.com/drive/folders/1dbJs6Rb6evI3jvyT6BTSHrGF4O9ryqmC?usp=sharing

请先登录，再进行评论。

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

Answer 1

KSSV 2021-10-11

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/1560991-how-to-merge-two-netcdf4-files-into-one#answer_805846

在 MATLAB Online 中打开

file1 = 'tropomi1.nc' ;
file2 = 'tropomi2.nc' ;
x1 = ncread(file1,'/PRODUCT/longitude');
y1 = ncread(file1,'/PRODUCT/latitude');
z1 = ncread(file1,'/PRODUCT/nitrogendioxide_tropospheric_column') ;
x2 = ncread(file2,'/PRODUCT/longitude');
y2 = ncread(file2,'/PRODUCT/latitude');
z2 = ncread(file2,'/PRODUCT/nitrogendioxide_tropospheric_column') ;
% Make my gird 
xx0 = min(min(x1(:)),min(x2(:)));
xx1 = max(max(x1(:)),max(x2(:)));
yy0 = min(min(y1(:)),min(y2(:)));
yy1 = max(max(y1(:)),max(y2(:)));
% DEfine your resolutions needed 
dx = 0.5 ; dy = 0.5 ;
[X,Y] = meshgrid(xx0:dx:xx1,yy0:dy:yy1) ;
F1 = scatteredInterpolant(x1(:),y1(:),z1(:),'nearest','none') ;
Z1 = F1(X,Y) ; 
F2 = scatteredInterpolant(x2(:),y2(:),z2(:),'nearest','none') ;
Z2 = F2(X,Y) ; 
Z = (Z1+Z2)/2 ;
pcolor(X,Y,Z) 
shading interp 