Matlab实现灰度图m*n分割

今天（5月7日）进行了第一次线上考试，开卷第一题就是编程题，出乎意料，想了一会儿先跳过做后面的，结果发现不能回退检查修改，所以开卷20分直接送走了…
考后再做这道题，发现了很多问题，但最终也成功编写出了代码，作此总结。

题目要求：编写一个函数，对输入的图像做2*2维分割，并返回图像的矩阵。
个人答案：用matlab实现，首先编写：

1 主文件read.m

img = imread("lena.jpg");
m = 2;
n = 2;

% RGB图转GRAY图
img = rgb2gray(img);
% 调用自己编写的imgSeg()函数做图像分割，将img分割成m*n的小图
[newimg, imgnum] = imgSeg(img, m, n);

for i = 1:imgnum
    % mat2gray做图像矩阵的归一化操作，否则保存的图为全白
    imwrite(mat2gray(newimg(:, :, i)), strcat('img', num2str(i), '.bmp')); 
end

2 分割函数文件imgSeg.m

function [newImg, imgNum] = imgSeg(img, segM, segN)
%UNTITLED Summary of this function goes here
%   Detailed explanation goes here
%   将输入的img分割成segM*segN维，保存到newImg并返回
% 获取图像尺寸
[w, h] = size(img);
% 定义newImg保存分割后的pixels
newImg = zeros(fix(w/segN), fix(h/segM), segM*segN);
imgNum = 1;
for i = 1:segN
    for j = 1:segM
        newImg(:, :, imgNum) = img(1+(i-1)*fix(w/segN):i*fix(w/segN), 1+(j-1)*fix(h/segM):j*fix(h/segM));
        % 调试点，查看当前分割的图像
        % imshow(img(1+(i-1)*fix(w/segN):i*fix(w/segN), 1+(j-1)*fix(h/segM):j*fix(h/segM)));
        imgNum = imgNum + 1;
    end
end

% imgNum在循环结束时多加了一次，故-1，使之=segM*segN
imgNum = imgNum - 1;

end

运行效果为：
输入：

输出：

3 问题总结

3.1 灰度图与彩色图

最开始我只使用了imread()，没有转灰度图，所以得到的img的RGB图的结构，对其操作就出现了错误。

3.2 高维矩阵降维

注意，在matlab中，若A是一个333的三维矩阵，则：
（1）A(1, :, : )是3D的；
（2）A(:, 1, : )是3D的；
（3）A(:, :, 1 )是2D的；
而用imwrite保存图像时，其矩阵应为2维，所以newImg如果采用了（1）或（2）方式需要做降维处理，例如squeeze(newimg(i, :, : ))，参考：matlab squeeze

3.3 图像矩阵归一化

没有用mat2gray()直接保存图像会发现图像全白，这是因为操作像素时将值变成了double型，而保存图像时double取值在0.0~1.0，大于1就白色，故需将其转换成uint8型，参考：matlab mat2gray()

今后要多多编程练习啊，不然这样期末考试真的会崩！！！

附工程文件下载地址：
下载地址1
下载地址2