今天(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()
今后要多多编程练习啊,不然这样期末考试真的会崩!!!