数字图像处理——月亮图像的识别与裁剪

已知一张带有月亮的图片(左)，截取只有月亮的区域(右)，效果如图：

预备知识：
find函数作用：
如图运行结果所示，a是一个二维矩阵，find函数作用是找到矩阵中所有不为零的点的坐标。坐标编号是从最左边一列开始，从上到下依次为1,2,3。。。换列时数字继续累加。

方法一：find法

img=imread('moon.png');%读取moon图像
W=im2bw(img);%把moon图像转化为二值图
[r,c]=find(W);%用find函数找到二值图中从上到下，从左至右所有像素点中值不为零的坐标
figure;%绘制一个窗口
subplot(121),imshow(img);%显示moon图像的原图
subplot(122),imshow(img(min(r):max(r),min(c):max(c),:));%从上到下，从左至右从前到后显示截取后的图像

方法二：投影法，利用sum函数

img=imread('moon.png');
W=im2bw(img);
c=find(sum(W));把二值图中的每个列相加，得到一个一维的矩阵，即得到一个原图中从上到下的一个投影。再用find函数找到这个一维数组中所有值不为0的值的坐标
r=find(sum(W'));通过图