我以' .hdf'的形式提供数据。各种变量的文件,我必须使用它来制作网格数据库。 我的数据具有相同数量的纬度,经度和数据值,如下所示
~~~~~~~
lat lon data1 data2 .... data125
x x x x ....... x
x x x x ........ x
~~~~~~~~
我现在已经复制了部分代码
这里,data1,2 ... 125用于不同高度的数据值。 我正在使用' gridfit'用于使网格超出lat和lon值,并将数据(data1)放在网格中的适当位置。
还有其他替代网格吗? griddata不适合我的应用程序。如何检查网格值是对还是错?
有人可以帮忙吗?
使用的部分代码:
[Attr_name,lat,lon,D1,D2] = data_read_cloudsat_ppt_col(dir(fullfile(files,'/*.hdf'))); %the data_read_cloudsat_ppt_col is a function to read the .hdf files
%% grid making from swaths
for hh=1:size(D1,2)
x=lat;
y=lon;
z=D1(:,hh);
cellsize=0.1;
minx=min(x);
maxx=max(x);
miny=min(y);
maxy=max(y);
xi=(minx:cellsize:maxx);
yi=(miny:cellsize:maxy);
[X,Y]=meshgrid(xi,yi);
[m,n]=size(X);
ZI{hh} = gridfit(y',x',z,yi,xi);
ZZZ_water=cat(3,ZI{:}); %concatenates the 2D files in 3D
clear ZI
end