咨询：笛卡尔转极坐标图像空白区域填充方案（用于OpenGL纹理）

问题原因 #

空白区域是因为原代码中warp函数的正向映射（笛卡尔→极坐标）无法覆盖极坐标空间的所有像素位置，部分区域没有对应原图像的像素，导致出现空白。

解决方案：反向映射（极坐标→笛卡尔） #

这种方式主动为每个极坐标输出像素寻找原图像的对应位置，通过插值确保所有区域都有有效像素，彻底避免空白。

方法1：使用interp2手动实现

% 加载图像
img = imread('my_image.png');
[h,w,ch] = size(img);

% 定义极坐标输出尺寸（可与原图像一致，或按需调整）
out_h = h;
out_w = w;

% 生成极坐标网格：角度0~2π，半径0到原图像对角线的一半
max_rho = sqrt((w/2)^2 + (h/2)^2);
[theta_out, rho_out] = meshgrid(linspace(0, 2*pi, out_w), linspace(0, max_rho, out_h));

% 转换为笛卡尔坐标，以原图像中心为原点
X_out = rho_out .* cos(theta_out) + w/2;
Y_out = rho_out .* sin(theta_out) + h/2;

% 双线性插值生成极坐标图像，边界空白用0填充（可改为邻近像素值）
polar_img = zeros(out_h, out_w, ch, class(img));
for c = 1:ch
    polar_img(:,:,c) = interp2(double(img(:,:,c)), X_out, Y_out, 'linear', 0);
end
polar_img = cast(polar_img, class(img));

% 查看结果
figure;
subplot(121); imshow(img); title('原图像');
subplot(122); imshow(polar_img); title('极坐标图像（无空白）');

方法2：使用imwarp简化实现

img = imread('my_image.png');
[h,w] = size(img);
center = [w/2, h/2]; % 图像中心坐标
max_rho = sqrt((w/2)^2 + (h/2)^2);

% 定义极坐标到笛卡尔坐标的逆变换函数
inv_transform = @(rho_theta) ...
    [rho_theta(:,2).*cos(rho_theta(:,1)), rho_theta(:,2).*sin(rho_theta(:,1))] + center;

% 生成极坐标输出网格
[theta_grid, rho_grid] = meshgrid(linspace(0, 2*pi, w), linspace(0, max_rho, h));
output_points = [theta_grid(:), rho_grid(:)];

% 执行变换，双线性插值+边界填充
polar_img = imwarp(img, inv_transform, ...
    'OutputPoints', output_points, ...
    'Interpolation', 'linear', ...
    'FillValues', 0);
polar_img = reshape(polar_img, h, w, []);

% 显示结果
imshow(polar_img); title('极坐标图像（无空白）');

备选方案：事后填充空白 #

如果坚持使用原warp方法，可通过后处理填充空白区域：

% 运行原代码生成warp图像并捕获
figure;
warp(theta, rho, Z, img); view(2); axis square;
frame = getframe(gcf);
polar_img = frame.cdata;
close(gcf);

% 识别空白区域（假设空白为低亮度值，可根据实际调整阈值）
gray_img = rgb2gray(polar_img);
mask = gray_img < 10;

% 填充空白孔洞
filled_img = imfill(polar_img, mask, 'holes');

% 查看填充结果
imshow(filled_img);

注意事项 #

• 优先选择反向映射方法，生成的图像像素无缺失，更适合作为OpenGL纹理使用。
• 可调整插值方式（如'nearest'或'cubic'）和填充值，匹配你的图像需求。

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