Convert sveral netCDF files to ASCII or CSV
12 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
Majid Mohamod
el 1 de Jun. de 2017
Respondida: javeed hamad
el 27 de Mzo. de 2022
HI all, I have 2000 files of netCDF (.nc). These files are for precipitation of my study area. I have the script which included below. I want to make this script to call all netCDF and output them as ASCII. Some netCDF files have attached too.. I really appreciate your help and time.
Thanks, Majid
10 comentarios
per isakson
el 5 de Jun. de 2017
Editada: per isakson
el 5 de Jun. de 2017
- What makes you think that this FEX-submission can read and convert your files?
- "ArcSWAT"   is that The Soil and Water Assessment Tool (SWAT)?
- What makes you think that the FEX-submission will write a csv-file that is appropriate for ArcSWAT ?
- NetCDFreader requires that the extension is .nc. Yours is .nc4.
Respuesta aceptada
KSSV
el 5 de Jun. de 2017
Dear Majid
Reading nc files is not a deal in MATLAB. It is pretty straight forward and easy. If you would have tried to read the documentation, by this time you could have solved your problem. No files is required to read the ncfiles. Just you need to know only few functions: ncdisp, ncreaad and ncreadatt.
In the below code, I am reading the precipitation data of the six ncfiles which you gave in the attachment. I am reading the begin, end date and time also. The output is taken into cells.
files = dir('*.nc4') ;
nfiles = length(files) ;
P = cell(nfiles,1) ; % precipitation of all files
date = cell(nfiles,2) ;
time = cell(nfiles,2) ;
for i = 1:nfiles
%%Read dat
date{i,1} = ncreadatt(files(i).name,'/','HDF5_GLOBAL.BeginDate') ;
date{i,2} = ncreadatt(files(i).name,'/','HDF5_GLOBAL.EndDate') ;
%%Read time
time{i,1} = ncreadatt(files(i).name,'/','HDF5_GLOBAL.BeginTime') ;
time{i,2} = ncreadatt(files(i).name,'/','HDF5_GLOBAL.EndTime') ;
%%REad precipitation
P{i} = ncread(files(i).name,'precipitation') ;
end
Now you can write the above data into what ever format you want.
5 comentarios
Más respuestas (2)
shobhit pipil
el 17 de Dic. de 2018
Editada: shobhit pipil
el 17 de Dic. de 2018
data = [];
Files=dir('*.nc');
for k=1:length(Files)
FileNames=Files(k).name
pr=ncread(FileNames,'rainfall_amount');%name of netCDF file; 'pr' is the precip variable in the netCDF file
long=ncread(FileNames,'x'); %'longitude' is the longitude variable in the netCDF file
lat=ncread(FileNames,'y'); %'latitude' is the latitude variable in the netCDF file
for j= 716:746(lat); %1:length(lat);
for i= 328:368(long); %1:length(long);
v=pr(i,j,:); %read through all lat,longs in netCDF
outfile=sprintf('%d_%d_PCP.txt',lat(j),long(i));%name of outputfile; format is 'LATITUDE_LONGITUDE_PCP.txt'
fid=fopen(outfile,'w+');
data2=(v);
data2(isnan(data2)) = -99.0 ; %cleanup in case the data is a NaN
data = {data;data2};
end
end
end
%fprintf(fid, '18900101\n');%#write the climate data starting date to header of outfile. Change '20200101' to the starting date of your dataset. For example, Jan. 1, 1950 is 19500101
fprintf(fid,'%5.1f\n',data); %formatting the data for SWAT specifications
fclose(fid);
disp([outfile 'created'])%display created file
fclose all
1 comentario
shobhit pipil
el 17 de Dic. de 2018
Would anybody please reply to me. I am trying to extract the NetCDF pricipitation data into SWAT format. The current script is not able to write data from NetCDF to txt file from each iteration. Please help.
javeed hamad
el 27 de Mzo. de 2022
HI every one!
I want to extract GHI (Global horizontal irradiation) data from NC file which include almost 3500 files for my study area, can someone please guide which script should use to extract multiple files simultaneously.
Thanks a lot, Javeed
0 comentarios
Ver también
Categorías
Más información sobre NetCDF en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!