Matlab可以处理数据

Question

我一直在处理855 4000x4000矩阵的处理2个列表。这是855个矩阵的列表，另一个是坐标列表（另外855 4000x4000矩阵）。在一个周期内执行此操作非常重要，以便不拥有成千上万的无用变量。对于每个文件，它都会剪切（读取不需要数据的NAN）协调数据，然后剪切与坐标相关的数据。然后，将所有值聚集到一个矩阵中。该代码是：

for x = 1:length(list_with_par)
    cd 'D:\Coord'
    par_lon = ncread(list_with_coordinates(x,:), 'longitude');
    par_lon(par_lon>=15) = nan;
    par_lon(par_lon<=-18) = nan;
    
    par_lat = ncread(list_with_coordinates(x,:), 'latitude');
    par_lat(par_lon>=84) = nan;
    par_lat(par_lon<=76) = nan;
    
    cd 'D:\Par'
    par = ncread(list_with_par(x,:), 'PAR');
    
    for i = 1:size(ncread(list_with_par(x,:),'PAR'),1) %size(,1)
        for z = 1:size(ncread(list_with_par(x,:),'PAR'),2) %size(,2)
            if isnan(par_lon(i,z))
                par(i,z) = nan;
            end
            if isnan(par_lat(i,z))
                par(i,z) = nan;
            end
        end
    end
    if size(par,2) < size(PAR_main,2)
        left_cells = size(PAR_main,2) - size(par,2);
        temp_cell = NaN(4865,left_cells);
        C2 = cat(2,par,temp_cell);
    end
    if size(par,2) == size(PAR_main,2)
        C2 = par(:,:,1);
    end
    PAR_main(:,:,x) = C2(:,:,1);
end

但是突然在处理4-5小时后会出现错误。

Error using netcdflib
The NetCDF library encountered an error during execution of 'open' function - 'HDF error (NC_EHDFERR)'.

Error in netcdf.open (line 67)
        [varargout{:}] = netcdflib ( 'open', filename, varargin{1} );

Error in internal.matlab.imagesci.nc/openToRead (line 1278)
            this.ncRootid = netcdf.open(this.Filename,'NOWRITE');

Error in internal.matlab.imagesci.nc (line 121)
                    this.openToRead();

Error in ncread (line 61)
ncObj   = internal.matlab.imagesci.nc(ncFile);

可能是一个问题？

Answer 1

我并不真正熟悉ncread（和关联的功能），但是有两件事跳出我的效率很低：

1. 在您的'i'和'z'中，是否有理由再次阅读数据以确定其大小，而不仅仅是使用您已经保存的“ PAR”变量？

   for i = 1：size（par，1）

   对于z = 1：size（par，2）

2. 为此，除非我缺少针对此功能的特定内容，否则您应该能够跳过循环'i'和'z'完全矢量化计算：

   par（isnan（par_lon））= nan;

   par（isnan（par_lat））= nan;

这无疑会大大减慢您的代码。除此之外，很难说，但是我绝对可以看到有数百万个无关紧要的文件读取可能会引起临时文件地址等的某些问题）或内存泄漏。