How to create nc file having lat, long, time, and 2D matrix

13 次查看（过去 30 天）

显示更早的评论

Jitesh Dadich 2019-4-4

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/454392-how-to-create-nc-file-having-lat-long-time-and-2d-matrix

回答： ANKUR KUMAR 2019-5-9

采纳的回答： KSSV

Helo everyone,

I am trying to create a .nc file which will have lat, long, time and 2D matrix simultaneously.

I have loaded lat, long, and time parameters , but how will I put the corresponding 2D matrix for each lat long and time while creating .nc file.

The answer should come in a matrix which has time in 1st column, lat and long in 1st & 2nd rows, and parameter values corresponding to time and lat_long.

Thank you in advance

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

请先登录，再进行评论。

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

采纳的回答

KSSV 2019-4-4

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/454392-how-to-create-nc-file-having-lat-long-time-and-2d-matrix#answer_369056

Refer this demo example:

% nc filename to be written
file = 'myfile.nc' ;
%% Write lon and lat variables
% Get data
lon = 1:10 ;
lat = 1:10 ;
nx = length(lon) ;
nccreate(file,'lon','Dimensions',{'lon',1,nx},'DeflateLevel',7) ;
ny = length(lat) ;
nccreate(file,'lat','Dimensions',{'lat',1,ny},'DeflateLevel',7) ;
nccreate(file,'time','Dimensions',{'time',1,Inf},'DeflateLevel',7) ;
nccreate(file,'z','Dimensions',{'lon','lat','time'},'DeflateLevel',7) ;
for i = 1:10
    ncwrite(file,'time',i,i)   % write time 
    data = rand(10) ; 
    ncwrite(file,'z',data,[1,1,i]) ;   % write 3D data 
end

14 个评论
显示 12更早的评论隐藏 12更早的评论

KSSV 2019-4-5

Attach your xlsx file.

Jitesh Dadich 2019-4-5

HW.csv

Here I have attached a small file of my excel file.

Time is varying from 3rd row to 52 row in 1 st column and first 2 row has lon and lat.

请先登录，再进行评论。

更多回答（1 个）

ANKUR KUMAR 2019-5-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/454392-how-to-create-nc-file-having-lat-long-time-and-2d-matrix#answer_374337

I hope this helps.

[data]=xlsread('test.csv');
lat=data(2,2:end);
lon=data(1,2:end);
vardata=data(3:end,2:end);
vardata(:,size(vardata,2)+1)=nan;   %dummy column just to put nan over the points where data is not available
longitude=unique(lon);
latitude=unique(lat);
[lon_m,lat_m]=meshgrid(longitude,latitude);
index=arrayfun(@(x,y) find(lon==x & lat==y),lon_m,lat_m,'uni',0);
index(cellfun('isempty',index))={-9999.99};
datavar=str2double(string(index)); datavar(datavar==-9999.99)=786;
out=arrayfun(@(x) reshape(vardata(x,datavar),[48 51]),1:50,'uni',0);
outputdata=cat(3,out{:});
nccreatewrite('test.nc','lat',{'lat','c'},latitude')
nccreatewrite('test.nc','lon',{'lon','c'},longitude')
nccreatewrite('test.nc','TC',{'lat','lon','days'},outputdata)