matlab矩阵缩小和放大

Matlab里的imresize函数可以对图像放大和缩小，但这同时也会改变图像矩阵的大小，如果想要上图所示的结果，需要再进行一些处理，处理代码如下所示。

clc

close all

% 从当前目录下打开一张图片

[filename, filepath] = uigetfile({'*.jpg*.ppmjpeg *.*.bmp*.png'},'Choose Input Image')

if isequal(filename,0) || isequal(filepath,0)

disp('User pressed cancel')

return

else

fullfp = fullfile(filepath, filename)

end

image = imread(fullfp) %代表要处理的图像

mysize = size(image)

%把图像转换成灰度图

if numel(mysize) >2

image = rgb2gray(image)

end

r_e = mysize(1)

c_e = mysize(2)

subplot(2,3,1)imshow(image,[])title('Input Image')

temp1 = imresize(image,2)%表示把图像放大到销雀陆原来的两倍，但同时图像矩阵也是变成了原来的两倍

[r_t1,c_t1] = size(temp1)

s = temp1(round(r_t1/2)-floor(r_e/2) : round(r_t1/2)+ceil(r_e/2)-1, round(c_t1/岁芹2)-floor(c_e/2) : round(c_t1/2)+ceil(c_e/2)-1)

subplot(2,3,2)imshow(s,[])title('Magnification')

temp2 = imresize(image,0.5)%表示把图像缩小到原来的一半，但同时图像矩阵也变成了原来的一半

[r_t2,c_t2] = size(temp2)

temp3 = zeros(r_e,c_e)

temp3(round(r_e/2)-floor(r_t2/2) : round(r_e/2)+ceil(r_t2/2)-1, round(c_e/2)-floor(c_t2/2) : round(c_e/2)+ceil(c_t2/2)-1) = temp2

ss = temp3

subplot(2,3,3)imshow(ss,[])title('Minification')

登录后复制

ps：直接调用imresize函数而不进行处理的效果如下图所示：

看起来图像没怎么变化，但存储图像的矩阵已经变大或变小了，如下图所示：

image为原亏顷始图像矩阵，temp1为放大后的图像矩阵，temp2为缩小后的图像矩阵

function [l,d]=ldlt(a)

n=length(a(:,1))

l(1,1)=1d(1,1)=a(1,1)

l(2,1)=a(2,1)/d(1,1)

l(2,2)=1d(2,2)=a(2,2)-l(2,1)*d(1,1)*l(2,1)

for i=3:n,

l(i,1)=a(i,1)/d(1,1)

l(i,i)=1

for j=2:i-1,

l(i,j)=(a(i,j)-l(i,1:j-1)*d(1:j-1,1:j-1)*l(j,1:j-1)')/d(j,j)

d(i,i)=a(i,i)-l(i,1:i-1)*d(1:i-1,1:i-1)*l(i,1:i-1)'

end

end

I=imread('cameraman.tif')

I=im2double(I) %转换图像矩阵为双精度型。

T=dctmtx(8) %产生二维DCT变换矩阵

%计算二维DCT，矩档锋键阵T及其转置是DCT函数P1*X*P2的参数

B=blkproc(I,[8,8],'P1*x*P2',T,T') %二值掩模，用来压缩DCT系数，只留下DCT系数中左上角的10个

mask=[1 1 1 1 0 0 0 0

1 1 1 0 0 0 0 0

1 1 0 0 0 0 0 0

1 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0

0 0 0 0 0 0 0 0 ]

B2=blkproc(B,[8 8],'P1.*x',mask) %只保留DCT变换的10个系数

I2=blkproc(B2,[8 8],'P1*x*P2',T',T) %重构图像

subplot(1,2,1)

imshow(I)

title('原始图像')

subplot(1,2,2)

imshow(I2)

title('压缩图像')

我这个程序也是改了后的，不过我运行过了，图片压缩效果都能出来，就是压缩比没有变小，我也不知道基戚为什么，如果有人给决绝了麻烦通知一下行巧，我邮箱sunnyxiaoxu@126.com

---