针对你需要从带453×453二维经纬矩阵的NetCDF文件中，提取纬度36°、经度9.15°位置全时段pr变量的需求，这里提供一个实用的MATLAB实现方案：

核心思路 #

因为你的经纬坐标是二维网格矩阵（而非一维的规则序列），每个网格点的经纬度是独立的，所以需要先找到与目标点（36°N, 9.15°E）空间距离最近的网格单元，再提取该单元对应的全时段降水数据。

完整代码实现 #

% 1. 定义目标坐标
target_lat = 36;
target_lon = 9.15;

% 2. 读取NetCDF文件（替换为你的文件路径）
ncfile = 'your_file.nc';

% 加载二维经纬矩阵和pr变量
lat_matrix = ncread(ncfile, 'lat'); % 453x453
lon_matrix = ncread(ncfile, 'lon'); % 453x453
pr_data = ncread(ncfile, 'pr');     % 假设pr的维度是 [x, y, time]，需根据实际调整

% 3. 计算每个网格点与目标点的球面距离（Haversine公式，精度更高）
lat_rad = deg2rad(lat_matrix);
lon_rad = deg2rad(lon_matrix);
target_lat_rad = deg2rad(target_lat);
target_lon_rad = deg2rad(target_lon);

dlat = lat_rad - target_lat_rad;
dlon = lon_rad - target_lon_rad;
a = sin(dlat/2).^2 + cos(target_lat_rad) .* cos(lat_rad) .* sin(dlon/2).^2;
c = 2 .* atan2(sqrt(a), sqrt(1-a));
distance = 6371 * c; % 地球半径6371km，结果单位为km

% 4. 找到距离最小的网格索引
[min_dist, idx] = min(distance(:));
[y_idx, x_idx] = ind2sub(size(lat_matrix), idx); % MATLAB矩阵索引是(row,col)对应(y,x)

% 5. 提取该网格点的全时段pr数据
% 注意：需根据pr实际维度顺序调整索引，比如维度为[time,x,y]则写成(:,x_idx,y_idx)
target_pr = pr_data(x_idx, y_idx, :); 

% 可选：验证提取的网格点实际经纬值
fprintf('最接近目标点的网格坐标：lat=%.4f, lon=%.4f，距离=%.4fkm\n', ...
    lat_matrix(y_idx,x_idx), lon_matrix(y_idx,x_idx), min_dist);

关键注意事项 #

• 维度顺序：一定要确认pr变量的维度顺序，可通过ncinfo(ncfile, 'pr')查看变量的维度信息，再对应调整代码中的索引位置。
• 简化计算：如果研究区域范围很小（几十公里内），可以用平面距离简化计算：distance = sqrt((lat_matrix - target_lat).^2 + (lon_matrix - target_lon).^2)，结果为度数差，足够用于筛选最近邻网格。
• 数据还原：如果pr变量带有scale_factor或add_offset属性，提取后记得还原原始值：target_pr = target_pr * scale_factor + add_offset;，同样可通过ncinfo(ncfile, 'pr')查看属性。

内容的提问来源于stack exchange，提问作者laykeh llamesi