如何用Matlab将.mat转.nc？解决元素数量不匹配错误

问题描述 #

需将.mat数据转换为NetCDF格式，编写MATLAB脚本时写入数据步骤失败，报错：

Error using matlab.internal.imagesci.netcdflib
The number of input elements does not match the variable size.

待转换的.mat为1×1结构体，包含字段：

• pld：4492500×1 double
• release_date：4492500×1 double
• x：4492500×72 double
• y：4492500×72 double

原代码如下：

% Define the path to the folder containing .mat files
folder_path = 'D:/.../output_test';

% Get a list of .mat files in the folder
mat_files = dir(fullfile(folder_path, '*.mat'));


% Loop through each .mat file
for i = 1:numel(mat_files)
    mat_file = mat_files(i);
    mat_data = load(fullfile(folder_path, mat_file.name));
    
    % Extract data from the loaded struct
    pld = mat_data.pld;
    release_date = mat_data.release_date;
    x = mat_data.x;
    y = mat_data.y;
    
    % Create a NetCDF file
    nc_file = strrep(mat_file.name, '.mat', '.nc');
    ncid = netcdf.create(nc_file, 'NETCDF4');
    
    % Define dimensions
    dim_x = netcdf.defDim(ncid, 'x', size(x, 2));  
    dim_y = netcdf.defDim(ncid, 'y', size(x, 2));  
    
    % Define variables
    var_pld = netcdf.defVar(ncid, 'pld', 'double', [dim_x, dim_y]);
    var_release_date = netcdf.defVar(ncid, 'release_date', 'double', [dim_x, dim_y]);
    var_x = netcdf.defVar(ncid, 'x', 'double', [dim_x, dim_y]);
    var_y = netcdf.defVar(ncid, 'y', 'double', [dim_x, dim_y]);
    
    % Complete the definition of the NetCDF file
    netcdf.endDef(ncid);
    
    % Write data to variables
    netcdf.putVar(ncid, var_pld, pld);
    netcdf.putVar(ncid, var_release_date, release_date);
    netcdf.putVar(ncid, var_x, x);
    netcdf.putVar(ncid, var_y, y);
    
    % Close the NetCDF file
    netcdf.close(ncid);
end

错误原因 #

核心问题是维度定义与数据实际结构完全不匹配：

1. 原代码定义的x和y维度长度均为72（取size(x,2)），并给所有变量指定[dim_x, dim_y]的二维结构（即72×72），但实际数据：
   • pld和release_date是4492500×1的一维数组，元素总数4492500，远大于72×72=5184，直接写入必然触发元素数量不匹配错误。
   • x和y是4492500×72的二维数组，元素总数为4492500×72，同样和72×72的变量维度不匹配。
2. 从数据结构看，应该定义两个核心维度：一个代表样本/轨迹数（4492500），另一个代表时间步/序列长度（72）。

修正后的代码 #

% Define the path to the folder containing .mat files
folder_path = 'D:/.../output_test';

% Get a list of .mat files in the folder
mat_files = dir(fullfile(folder_path, '*.mat'));

% Loop through each .mat file
for i = 1:numel(mat_files)
    mat_file = mat_files(i);
    mat_data = load(fullfile(folder_path, mat_file.name));
    
    % Extract data from the loaded struct
    pld = mat_data.pld;
    release_date = mat_data.release_date;
    x = mat_data.x;
    y = mat_data.y;
    
    % Get actual data dimensions
    num_samples = size(pld, 1);    % 4492500
    num_timesteps = size(x, 2);    % 72
    
    % Create a NetCDF file (save to the same folder as original .mat)
    nc_file = strrep(mat_file.name, '.mat', '.nc');
    ncid = netcdf.create(fullfile(folder_path, nc_file), 'NETCDF4');
    
    % Define dimensions
    dim_samples = netcdf.defDim(ncid, 'samples', num_samples);
    dim_timesteps = netcdf.defDim(ncid, 'timesteps', num_timesteps);
    
    % Define variables with matching dimensions
    var_pld = netcdf.defVar(ncid, 'pld', 'double', dim_samples);
    var_release_date = netcdf.defVar(ncid, 'release_date', 'double', dim_samples);
    var_x = netcdf.defVar(ncid, 'x', 'double', [dim_samples, dim_timesteps]);
    var_y = netcdf.defVar(ncid, 'y', 'double', [dim_samples, dim_timesteps]);
    
    % Complete definition phase
    netcdf.endDef(ncid);
    
    % Write data to variables
    netcdf.putVar(ncid, var_pld, pld);
    netcdf.putVar(ncid, var_release_date, release_date);
    netcdf.putVar(ncid, var_x, x);
    netcdf.putVar(ncid, var_y, y);
    
    % Close the NetCDF file
    netcdf.close(ncid);
    
    fprintf('Converted %s to %s successfully.\n', mat_file.name, nc_file);
end

补充说明 #

1. 维度命名改为samples和timesteps，避免与变量x/y混淆，更贴合数据含义。
2. 修正了NetCDF文件的保存路径，确保生成的.nc文件与原.mat文件在同一文件夹。
3. 增加转换完成提示，方便批量处理时查看进度。

内容的提问来源于stack exchange，提问作者Charlotte Sève