matlab使用单像素成像方式重建三维物体不同深度图像的代码

以下是使用单像素成像方式重建三维物体不同深度图像的MATLAB代码示例：

% 定义物体的三维坐标和深度图像
[x,y,z] = meshgrid(-0.5:0.01:0.5,-0.5:0.01:0.5,-0.5:0.01:0.5);
depth1 = sqrt(x.^2+y.^2+z.^2);
depth2 = sqrt((x-0.1).^2+y.^2+z.^2);

% 定义相机参数
f = 1; % 焦距
pixel_size = 0.01; % 像素大小
image_size = 64; % 图像大小
u0 = image_size/2; % 像素中心横坐标
v0 = image_size/2; % 像素中心纵坐标

% 计算物体在图像上的投影
u1 = f*x./depth1/pixel_size+u0;
v1 = f*y./depth1/pixel_size+v0;
u2 = f*(x-0.1)./depth2/pixel_size+u0;
v2 = f*y./depth2/pixel_size+v0;

% 绘制图像
image1 = zeros(image_size,image_size);
image2 = zeros(image_size,image_size);
for i = 1:numel(u1)
    if round(u1(i))>0 && round(u1(i))<=image_size && round(v1(i))>0 && round(v1(i))<=image_size
        image1(round(v1(i)),round(u1(i))) = 1;
    end
    if round(u2(i))>0 && round(u2(i))<=image_size && round(v2(i))>0 && round(v2(i))<=image_size
        image2(round(v2(i)),round(u2(i))) = 1;
    end
end

% 显示结果
figure;
subplot(1,2,1);
imshow(image1);
title('深度为0的图像');
subplot(1,2,2);
imshow(image2);
title('深度为0.1的图像');

其中，`[x,y,z] = meshgrid(-0.5:0.01:0.5,-0.5:0.01:0.5,-0.5:0.01:0.5);` 定义了物体的三维坐标，`depth1` 和 `depth2` 分别是不同深度的图像。`f`、`pixel_size`、`image_size`、`u0` 和 `v0` 是相机参数，`u1`、`v1`、`u2` 和 `v2` 分别是物体在图像上的投影。最后通过循环遍历所有的投影点，将图像矩阵中对应的像素设置为1，最终得到两幅深度不同的图像。