MATLAB中不规则极坐标网格转规则网格插值失败求助

一、先明确网格坐标类型 #

首先确认你的Xi、Yi是极坐标（r, θ）还是已转换的笛卡尔坐标（x, y）：

• 若为极坐标：先将所有网格点转笛卡尔坐标：

  x = Xi .* cos(Yi);
  y = Xi .* sin(Yi);
  

• 若已是笛卡尔坐标：直接进入散点插值环节。

二、排查散点插值失败的常见问题及修复 #

你尝试scatteredInterpolant、griddata失败，通常是以下原因导致，对应处理：

1. 存在重复/重合网格点
   不规则极坐标网格容易出现重合点（比如极径为0时，所有极角对应同一点），会触发插值函数报错。处理方式：

   % 把坐标和数据扁平化并合并
   all_points = [Xi(:), Yi(:), DATA(:)];
   % 保留唯一的(x,y)点，稳定排序避免打乱原始顺序
   [~, unique_idx] = unique(all_points(:,1:2), 'rows', 'stable');
   unique_data = all_points(unique_idx,:);
   
   x_unique = unique_data(:,1);
   y_unique = unique_data(:,2);
   data_unique = unique_data(:,3);
   

2. 规则网格超出原始数据凸包范围
   griddata默认对凸包外的点返回NaN，若你的规则网格大量点在原始数据覆盖范围外，会被误认为插值失败。处理方式：

   • 先可视化原始数据的凸包，确认规则网格范围：

     hull_idx = convhull(x_unique, y_unique);
     plot(x_unique(hull_idx), y_unique(hull_idx), 'r--'); hold on;
     plot(x_unique, y_unique, 'b.'); hold off;
     

   • 若必须对凸包外点插值，用scatteredInterpolant设置外插方法：

     interp_fun = scatteredInterpolant(x_unique, y_unique, data_unique, 'linear', 'nearest');
     DATA_interp = interp_fun(Xr, Yr);
     

三、极坐标网格的专属优化方案 #

如果原始数据是不规则极坐标网格（Xi为极径网格，Yi为极角网格），可以先在极坐标下插值到规则极网格，再转笛卡尔坐标，能减少边缘误差：

1. 生成规则极坐标网格：

   R_range = linspace(min(Xi(:)), max(Xi(:)), Nr);
   Theta_range = linspace(min(Yi(:)), max(Yi(:)), Mr);
   [Rr, Thetar] = meshgrid(R_range, Theta_range);
   

2. 极坐标下插值：

   DATA_polar = griddata(Xi(:), Yi(:), DATA(:), Rr, Thetar);
   

3. 转换为笛卡尔规则网格：

   Xr = Rr .* cos(Thetar);
   Yr = Rr .* sin(Thetar);
   DATA_interp = DATA_polar;
   

四、验证插值结果 #

用绘图确认插值合理性：

figure;
subplot(1,2,1);
pcolor(Xi, Yi, DATA); shading flat; axis equal; title('原始不规则网格');
subplot(1,2,2);
pcolor(Xr, Yr, DATA_interp); shading flat; axis equal; title('插值后规则网格');

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