Octave调用griddata插值环面帽数据时大量返回NaN值问题

griddata返回大量NaN是两个问题共同导致的：

• 传入的采样点包含大量无效NaN值：你执行z(z<0) = NaN后，没有同步过滤对应位置的x、y坐标，相当于给插值函数传入了近一半配对值为NaN的无效点，这部分点会被griddata直接丢弃，实际参与插值的有效点数量远低于预期。
• 插值网格范围和有效点实际坐标范围不匹配：你设置的[-5,5]网格区间远小于当前环面参数下的有效点分布范围，绝大多数有效采样点都落在网格边界外，外插区域默认返回NaN。
• 额外注意：你当前写的环面参数方程默认对称轴为x轴，而非z轴，如果需要生成对称轴沿z轴的标准环面，需要调整x/y/z和参数Theta、Phi的对应关系。

1. 先过滤所有含NaN的无效采样点，只保留坐标和值均有效的点传入插值函数
2. 按有效点的实际坐标范围调整插值网格边界，不要硬设不符合数据分布的区间
3. 插值方法优先选自然邻域插值natural，比默认线性插值对非均匀分布点的适配性更好，空值更少

R = 2;
r = 7;

alpha_degree = 90;
alpha_rad = alpha_degree*pi/180;

theta = linspace(0,2*pi,100);
phi = linspace(0,pi,100);
[Phi,Theta] = meshgrid(phi,theta); 

% 当前参数方程生成的是对称轴沿x轴的环面
x = r.*sin(Theta);
y = (R+r.*cos(Theta)).*cos(Phi);
z = (R+r.*cos(Theta)).*sin(Phi);

% 裁剪z<0的部分
z(z<0) = NaN;

% 坐标旋转
x_rot = x.*cos(alpha_rad) - y.*sin(alpha_rad);
y_rot = y.*cos(alpha_rad) + x.*sin(alpha_rad);
z_rot = z;

% 过滤所有NaN无效点
x_flat = x_rot(:);
y_flat = y_rot(:);
z_flat = z_rot(:);
valid_mask = ~isnan(z_flat);
x_valid = x_flat(valid_mask);
y_valid = y_flat(valid_mask);
z_valid = z_flat(valid_mask);

% 按有效点范围生成规则网格，留1单位余量避免边缘截断
x_range = [min(x_valid)-1, max(x_valid)+1];
y_range = [min(y_valid)-1, max(y_valid)+1];
[xx,yy] = meshgrid(linspace(x_range(1),x_range(2),100), linspace(y_range(1),y_range(2),100));

% 自然邻域插值
zz = griddata(x_valid, y_valid, z_valid, xx, yy, 'natural');

% 可直接用surf查看插值结果
% surf(xx,yy,zz, 'EdgeColor','none')

如果需要旋转后再沿其他轴裁剪环面帽，把z(z<0) = NaN的裁剪步骤移到坐标旋转操作之后，按需要裁剪对应坐标轴的负值即可。

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