数学软件及应用.docVIP

数学软件及应用 ——图像的缩放 学院：理学院 专业：信息与计算科学 学号：050108010023 姓名：董玉波 日期：2008-12-24 Matlab图像缩放算法 算法的目的： 根据图像上的关键部分和不重要区域，对图片进行拉伸或压缩。本文实现的是图像的纵向和横向seam的删除。 1、特性和功能 该程序允许用户调整图像，从纵向或横向消除连续路径像素。纵向seam的定义是一个持续的路径运行的像素从顶部向底部的一个形象的情况下，而横向seam是一个持续不断的像素线横跨从左至右的形象。 该程序的主要功能包括：通过给定的列数运行程序，可以得到图像的横向压缩效果，适合于由普通尺寸的图片转向各种尺寸的图像压缩。 2、详细的算法实现 第一步计算删除或插入的seam涉及到计算原始图像的梯度图像。梯度图像是一种常见的图像，是用来在横向和纵向seam的计算，然后平均每三个梯度图像。在这个程序中用Sobel算子计算图像梯度，但其他的梯度也可使用。 （ MATLAB的功能： findEnergy ） 像素指数 (i-1,j-1) (i-1,j) (i-1,j+1) (i,j-1) (i,j) (i,j+1) (i+1,j-1) (i+1,j) (i+1,j+1) 一旦确定了图像梯度的计算方法，下一步就是计算energymap图像。 energymap图像需要分开计算垂直或横向的seam，而且也需要将seam删除后重新计算。纵向seam的计算方法如下：对梯度图像中的每个像素(i,j)，在energymap中(i，j)的值是现在值和前一行中min((i-1,j-1),(i-1,j),(i-1,j+1))的和。对i=1(初始行)，它们的值被设成梯度图像中相对应的值。（ MATLAB的功能：findSeamImg ） 完成了energymap的计算方法后，该方法寻找最佳seam是先找到最后一排的最小值(i,j)，然后程序倒退找到3相邻像素(i-1,j-1),(i-1,j),(i-1,j+1)中最低的一个，并且保存它的路径。这个过程反复进行，直至第一行完成。（ MATLAB的功能： findSeam ） 优化后，发现seam的路径，使像素的seam从梯度图像和原始的RGB图像中删除，其后的像素依次回退一个像素，形成了一个连续的形象。 对于插入seam的情况（增加图像大小） ，可以计算出沿某一方向的一个seam，插入相邻两个像素的平均值。如果所需的图像尺寸是要增加在一个特定的方向上增加像素，前N条被删除的seam的计算必须首先完成（ MATLAB的功能： removalMap ） ，然后平均像素插入每次连续seam，从而限制了最大增加图像尺寸。 3，实现算法的关键 （1）. 寻找图像上能量最小的缝(Seam)： Vertical seam: Horizontal seam: 最简单的能量定义： （图像x,y方向的梯度和） 利用动态规划寻找能量最小的缝，主要是计算累计能量： （2）. 行缝和列缝的处理顺序： 主要是引入 4、运行程序 （i）列缝处理 打开SandboxScript.m文件，按F5运行，等待大约6分钟后，将出现经过压缩的图像。 列缝处理的具体MATLAB函数和代码如下： （1）、主函数： close all clear all X=imread(pagodas_small.jpg); X=double(X)/255; [rows cols dim]=size(X); Y=rgb2hsv(X); E1=findEnergy(X); figure(1) imshow(E1) S1=findSeamImg(E1); figure(2) imshow(S1,[min(S1(:)) max(S1(:))]) SeamVector=findSeam(S1); %显示叠加着seam的图像。 SeamedImg=SeamPlot(E1,SeamVector); SeamCutImg=SeamCut(X,SeamVector); Figure(3) imshow(SeamCutImg) M=removalMap(X,cols); %cols表示删除的像素数。 MSeamedImg=SeamCut(X,M); Figure(4) %显示压缩后的图像。 imshow(MSeamedImg) （2）、findEnergy函数从给定的RGB或灰白图像中创建了一幅梯度图。假如RGB图像被输入，可以用Sobel算子找到那水平和垂直的梯度。 function Emean=findEnergy(x) [rows cols dim]=size(x); Grd=[ -1 -2 -1; 0 0 0; 1 2