matlab练习程序（随机直线采样）

我只是感觉好玩，写了这样一段程序。

原理就是先随机生成两个点，然后根据这两个点画直线，最后在直线上的像素保留，没在直线上的像素丢弃就行了。

最后生成了一幅含有很多空洞的图像。

当然，对含有空洞的图像是可以用修复算法修复的。

我也尝试修复了一下，用的算法我过去也写过，可以看这里。

这一次就不贴修复代码了，那段程序中的输入图像img、mask和这里的输出图像img、mask是一模一样的。

原图：

采样后：

修复后：

matlab代码如下：

main.m:

clear all; close all;clc;

img=double(imread('lena.jpg'));
imshow(img,[]);

[h,w]=size(img);

n=;
pix=rand(n,);
pix(:,:)=round(pix(:,:)*(w-))+;
pix(:,:)=round(pix(:,:)*(h-))+;
mask=zeros(h,w);

for i=:n
    x1=pix(i,);x2=pix(i,);
    y1=pix(i,);y2=pix(i,);
    mask=drawline(mask,x1,y1,x2,y2);
end

ind=find(mask~=);
img(ind)=;
figure;imshow(img,[])

drawline.m:

function img=drawline(img,x1,y1,x2,y2)   %因为图像坐标和数学函数坐标y轴朝向相反，所以这里所有y变量取相反数
    [h w]=size(img);
    y1=-y1;
    y2=-y2;

    if x1~=x2
        k=(y2-y1)/(x2-x1);
        b=y1-k*x1;
        for i=:w
            yy=-round(i*k+b);
            xx=i;
            if yy>= && yy<=h && xx>= && xx<=w
                img(yy,xx)=;
            end
        end
    end

    if y1~=y2
        k=(x2-x1)/(y2-y1);
        b=x1-k*y1;
        for i=-h:
            yy=-i;
            xx=round(i*k+b);
            if yy>= && yy<=h && xx>= && xx<=w
                img(yy,xx)=;
            end
        end
    end

end