Merging Multiple text files

7 visualizzazioni (ultimi 30 giorni)
minh lan
minh lan il 20 Mag 2019
Commentato: KSSV il 21 Mag 2019
I have multiple text files, each file is data for 1 day with 3 columns, 1 columns for longtitude, 1 columns for latitude and 1 columns for precip value.
It looks like that
108.650 13.950 0.000
108.450 13.383 0.000
105.833 22.150 0.000
106.217 21.300 -99.000
104.283 22.533 0.000
105.717 9.283 0.100
The value of longtitude and latitude is same for all file.
Now, I want to merge all file to 1 netcdf or text file with dimension is 455x455x365 (455 is the number of long and lat, 365 is precip value for each day of a year)
How can I do that?
  2 Commenti
Geoff Hayes
Geoff Hayes il 21 Mag 2019
minh - do you have 365 files (one for each day of the year)? Does each file have 455 rows? Please clarify.
minh lan
minh lan il 21 Mag 2019
Geoff Hayes,
Yes, I have 365 files and each file has 455 rows and 3 columns,

Accedi per commentare.

Risposta accettata

KSSV
KSSV il 21 Mag 2019
Modificato: KSSV il 21 Mag 2019
txtfiles = dir('*.txt') ;
ncfile = 'myfile.nc' ;
N = length(txtfiles) ;
nx = 455 ; ny = 455 ; % number of lat, lon
for i = 1:N
data = importdata(txtfiles(i).name) ;
lon = data(:,1) ;
lat = data(:,2) ;
p = data(:,3) ;
if i == 1
% Do interpolation
idx = boundary(lon,lat) ;
xi = linspace(min(lon),max(lon),nx) ;
yi = linspace(min(lat),max(lat),ny) ;
[X,Y] = meshgrid(xi,yi) ;
idx = inpolygon(X,Y,lon(idx),lat(idx)) ;
Z = griddata(lon,lat,p,X,Y) ;
Z(~idx) = NaN ;
% Longitude
nccreate(ncfile,'lon','Dimensions',{'lon',1,nx}) ;
ncwrite(ncfile,'lon',xi) ;
ncwriteatt(ncfile,'lon','long_name','longitude');
% Latitude
nccreate(ncfile,'lat','Dimensions',{'lat', 1, ny})
ncwrite(ncfile,'lat',yi)
ncwriteatt(ncfile,'lat','long_name','latitude');
% Precipitation
nccreate(ncfile,'time','Dimensions',{'time',1,Inf}) ;
nccreate(ncfile,'p','Dimensions',{'lon',nx,'lat',ny,'time',Inf})
ncwriteatt(ncfile,'p','long_name','precipitation');
ncwrite(ncfile,'time',i,i) ;
ncwrite(ncfile,'p',Z,[1,1,i])
end
Z = griddata(lon,lat,p,X,Y) ;
Z(~idx) = NaN ;
ncwrite(ncfile,'time',i,i) ;
ncwrite(ncfile,'p',Z,[1,1,i])
end
  11 Commenti
minh lan
minh lan il 21 Mag 2019
One more question,
If I want to have scatter data to compare the different map between the scattered data and the data after interpolation?
How can I change your code?
Thank you.
KSSV
KSSV il 21 Mag 2019
USe lon,lat,p for scattered data.....use X,Y,Z for grid data.

Accedi per commentare.

Più risposte (0)

Tag

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by