数字图像处理技术与应用练习题
说明：不允许直接使用MATLAB （或者OPENCV 等等软件）所带的图像图像
函数，重点考察大家是否理解了各种处理算法，算法可用伪代码描述。
算法应较详细。
1、 设一幅大小为M ×N 的灰度图像I 中，灰度为g 的像素数为h(g)， 0255g ≤≤。
请写出对图像I 进行直方图均衡化，得到图像J 的计算方法。
clc;
clear;
I= imread('a.jpg');
I= rgb2gray(I) ; %将图像转换为灰度图像
J= histeq( I) ; %对I 进行直方图均衡化
subplot( 121) ,imshow(I) ,title('原始图像') ;
subplot (122), imshow(J), title('直方图均衡化后的图像');
figure( 2) ;
subplot( 121) ,imhist(I, 64), title( '原始的直方图');
subplot( 122) ,
imhist(J,64) ,
title(' 均衡化后的直方图');
2、 设一幅大小为M ×N 的灰度图像I 中，现要变成（放大或缩小）为 P ×Q 的
图像J ，请写出J 的生成算法（可以使用近邻插值）。
I=imread('f.jpg');%读入图像
%图像属性
% Filename: 'f.jpg'
% FileModDate: '24-Aug-2008 16:50:30'
% FileSize: 20372
% Format: 'jpg'
% FormatVersion: ''
% Width: 480
% Height: 640
% BitDepth: 8
% ColorType: 'grayscale'
% FormatSignature: ''
% NumberOfSamples: 1
% CodingMethod: 'Huffman'
% CodingProcess: 'Sequential'
% Comment: {}
[rows,cols]=size(I);
K1 = str2double(inputdlg('请输入行缩放倍数', 'INPUT scale factor', 1, {'0.6'}));%行默认变为原来的0.6倍
K2 = str2double(inputdlg('请输入列缩放倍数', 'INPUT scale factor', 1, {'0.4'}));%列默认变为原来的0.4倍
width = K1 * rows;
height = K2 * cols;
im2 = uint8(zeros(width,height)); %定义输出图像矩阵
widthScale = rows/width;
heightScale = cols/height;
for x = 6:width - 6 %为防止矩阵溢出而选择的参数6
for y = 6:height - 6
oldX = x * widthScale; %oldX，oldY为原坐标，x，y为新坐标
oldY = y * heightScale;
if (oldX/double(uint16(oldX)) == 1.0) & (oldY/double(uint16(oldY)) == 1.0) im2(x,y) = I(int16(oldX),int16(oldY));
else
a = double(round(oldX));
b = double(round(oldY)); %若不是整数四舍五入后把临近值赋过去
im2(x,y) = I(a,b);
end
end
end
imshow(I); %输出原图像
figure;
imshow(im2); %输出缩放后图像
3、设一幅大小为M×N的灰度图像I中，现要将其逆时针旋转A度，得到图
像J，请写出J的生成算法（可以使用近邻插值）。
clear;%此题是用最近邻域法实现图像旋转
im1=imread('b.jpg');
[m,n,p]=size(im1);
% 将图像旋转30度
a=0.5; %a=sin30=0.5
b=0.866; %b=cos30=0.866
row=n*a+m*b;
col=n*b+m*a;
for i=1:row %先把图象填充成全黑for j=1:col
im2(i,j,:)=uint8(0);
end
end
for i=1:m %把原图象像素点旋转后变为新图象点
for j=1:n
xx=round(abs((i-m/2)*b-(j-n/2)*a+row/2));
yy=round(abs((i-m/2)*a+(j-n/2)*b+col/2));
for k=1:3
im2(xx,yy,k)=im1(i,j,k);
end
end
end
temp1=uint8(0);
temp2=uint8(0);
temp3=uint8(0);
for i=1:row %把画面上的空点按照最近邻插值法填充
temp1=uint8(0);
temp2=uint8(0);
temp3=uint8(0);
for j=1:col %找到最右的图象边界点
if (im2(i,j,:)==uint8(0))
else
kk=j;
end
end
for j=1:kk
if (im2(i,j,:)==uint8(0))
im2(i,j,1)=temp1;
im2(i,j,2)=temp2;
im2(i,j,3)=temp3;
else
temp1=im2(i,j,1);
temp2=im2(i,j,2);
temp3=im2(i,j,3);
end
end
end
imshow(im1);
figure;
imwrite(im1,'5.jpg'); %保存原图像
imshow(im2);
imwrite(im2,'6.jpg');%保存旋转后图像
4、请写出生成（2N+1）×（2N+1）大小的高斯模板H（方差为sigma）的方
法。
5、请写出生成（2N+1）×（2N+1）大小的高斯一阶导数模板HX(水平方向的
梯度)、HY（垂直方向的梯度）（高斯的方差为sigma）的方法。
6、请写出使用大小为（2N+1）×（2N+1）模板H对图像I进行滤波，生成图
像J的方法。
7、请写出使用大小为3×3的模板对图像I进行中值滤波，生成图像J的方法。clear;
A=imread('num22','bmp');
subplot(1,2,1);
B=rgb2gray(A);
subimage(B);
title('处理前的图');
C=B;
xsize=size(B);
for k=2:(xsize(1)-1)
for j=2:(xsize(2)-1)
t=B(k-1:k+1,j-1:j+1);
C(k,j)=median(t(1:9));
end
end
subplot(1,2,2);
subimage(C);
title('处理后的图');
8、请写出求Otsu阈值（即最大类间距准则）的计算方法。
试证明采用最大类间距准则计算出的阈值与采用最小类内距准则计算出的阈值相同。
9、设有一幅二值图像（元素取值为0或1），请生成该图像的标记图像。（即
第一个连通区域中的每一个白色像素的值都置为1，第二个连通区域中的每一个白色像素的值都置为2，依此类推。区域编号可不考虑顺序）
clear all;
close all;
clc;
img=imread('liantong.bmp');
imgn=img>128;
s=uint8(1-imgn);
%{
s=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0; %这个矩阵是维基百科中的矩阵
1 1 1 1 1 1 1 1 0 0 1 1 1 1 0 0;
0 0 1 1 1 1 0 0 0 1 1 1 1 0 0 0;
0 1 1 1 1 0 0 0 1 1 1 0 0 1 1 0;
1 1 1 0 0 1 1 0 0 0 1 1 1 0 0 0;
0 1 1 0 0 0 0 0 1 1 0 0 0 1 1 0;
0 0 0 0 0 1 1 1 1 0 0 1 1 1 1 0;
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
%}
imshow(mat2gray(s));
[m n]=size(s);
tmp=zeros(m,n); %标记图像
label=1;
queue_head=1; %队列头
queue_tail=1; %队列尾
neighbour=[-1 -1;-1 0;-1 1;0 -1;0 1;1 -1;1 0;1 1]; %和当前像素坐标相加得到八个邻域坐标
for i=2:m-1
for j=2:n-1
if s(i,j)==1 && tmp(i,j) ==0
tmp(i,j)=label;
q{queue_tail}=[i j]; %用元组模拟队列，当前坐标入列
queue_tail=queue_tail+1;
while queue_head~=queue_tail
pix=q{queue_head};
for k=1:8 %8邻域搜索
pix1=pix+neighbour(k,:);
if pix1(1)>=2 && pix1(1)<=m-1 && pix1(2) >=2 &&pix1(2)<=n-1
if s(pix1(1),pix1(2)) == 1 && tmp(pix1(1),pix1(2)) ==0 %如果当前像素邻域像素为1并且标记图像的这个邻域像素没有被标记，那么标记
tmp(pix1(1),pix1(2))=label;
q{queue_tail}=[pix1(1) pix1(2)];
queue_tail=queue_tail+1;
end
end
end
queue_head=queue_head+1;
end
clear q; %清空队列，为新的标记做准备
label=label+1;
queue_head=1;
queue_tail=1;
end
end
end
10、设一幅二值图像中，只有一个白色区域，试给出求该区域外围轮廓线的方法
（要求按顺时针的顺序给出各点的坐标，即行/列号）。
轮廓提取的算法很多，本题对于二值图像可以用一个简单的算法进行处理，首先遍历图像中的每一个像素点，查看该点像素值是否为白色像素点；若是再查看该点周围8邻域像素点的灰度值之和为2040；若是，则该点是内部点，输出图像相应位置置为黑色。若不是，则输出该点的坐标信息。
11、设有一幅二值图像，采用3×3的结构元（每个元素均为1）对其进行腐蚀
操作，试写出得到结果图像的方法。
%A是处理的矩阵B是结构体，此题为3x3的元素均为1的结构体，用来对A 进行扫描的
%实现的是简单的二值图像的处理灰度图像有待研究
%实现图像的腐蚀已完成未测试~
function ans =mimerode(A,B)
%初始化
[ha,wa] =size(A);
[hb,wb] =size(B);
origin =[0,0];
ans =zeros(ha,wa);
for i = 1:hb
for j = 1:wb
if (B(i,j)== 1)
origin = [i,j];
break;
end
end
end
%对矩阵A进行边扩张以0补充
new_A =zeros(ha+hb*2,wa+wb*2);
for i = 1:ha
for j = 1:wa
new_A(i+hb,j+wb) = A(i,j);
end
end
%对矩阵A进行遍历
for i = 1:ha
for j = 1:wa
if A(i,j) == 0
continue;
end
%把B的origin点移到当前处理点
%寻找B的影响区域块
begin_i = i + hb - origin(1);
begin_j = j + wb - origin(2);
%对矩阵new_A中和重合的B中的元素对比
for i1 =begin_i:begin_i+hb-1
for j1 =begin_j:begin_j+wb-1
flag = true;
if (B(i1-begin_i+1,j1-begin_j+1)== 1 & new_A(i1,j1) == 0)
flag = false
break;
end
end
end
if flag
ans(i,j) = 1;
else
ans(i,j) = 0;
end
end
end
12、写出二维离散傅立叶变换、反变换的计算公式，给出基于频域滤波的基本步
骤。
二维离散傅里叶变换公式：假设以正方形网格采样得到的图像用f(x,y)来表示，其大小为N行M列，则f(x,y)的二维离散傅里叶变换可以表示为：Fu,v=x=0M-1y=0N-1fx,ye-j2π(uxM+vyN), u=0,1,…M-1, v=0,1,…N-1
其反变换为：
fx,y= 1MNu=0M-1v=0N-1Fu,ve-j2π(uxM+vyN), x=0,1,…M-1,
y=0,1,…N -1
基于频域滤波的基本步骤：
（1） 用(-1)x+y 乘以图像来进行中心变换
（2） 由（1）计算图像的DFT （离散傅里叶变换），即F(u,v)
（3） 用滤波器函数H(u,v)乘以F(u,v)
（4） 计算（3）中结果的反DFT
（5） 得到（4）中结果的实部
（6） 用(-1)x+y 乘以（5）中的结果
13. 请写出C 均值（K-means ）聚类分割的基本步骤。
14. 请写出Canny 算子检测边缘的详细步骤。
算法过程
1高斯函数
2221G(x,y)exp ||,,22y m
x m x m y m x y n m n σ++--+-=-=∑∑表示 高斯滤波器窗口大小
2偏导数:使用微分算子求出偏导数
x y 111111G ,G 111122-????=?=?????---????
B arctan y
x B B ==梯度大小
3非极大值抑制：
沿幅角方向检测模值的极大值点，即边缘点，遍历8个方向图像像素，把每个像素偏导值与相邻像素的模值比较，取其MAX 值为边缘点，置像素灰度值为0.
3双阈值检测：
由于单阈值处理时，合适的阈值选择较 困难，常常需要采用反复试验，因此采用双阈值检测算法。
对经过非极大值抑制后的图像作用两个阈值th1，th2，th1=0.4th2,两个阈值作用后得到两个图像1、2，较大阈值检测出的图像2去除了大部分噪声，但是也损失了有用的边缘信息。
较小阈值检测得到的图像1则保留着较多的边缘信息，以此为基础，补充图像2中的丢失的信息，连接图像边缘。
4链接边缘的具体步骤如下：
? 对图像2进行扫描，当遇到一个非零灰度的像素p(x,y)时，跟踪以p(x,y)
为开始点的轮廓线，直到轮廓线的终点q(x,y)。
? 考察图像1中与图像2中q(x,y)点位置对应的点s(x,y)的8邻近区域。如
果在s(x,y)点的8邻近区域中有非零像素s(x,y)存在，则将其包括到图像2
中，作为r(x,y)点。从r(x,y)开始，重复第一步，直到我们在图像1和图像
2中都无法继续为止。
? 当完成对包含p(x,y)的轮廓线的连结之后，将这条轮廓线标记为已经访问。
回到第一步，寻找下一条轮廓线。重复第一步、第二步、第三步，直到图
像2中找不到新轮廓线为止。
? 至此，完成canny 算子的边缘检测。
15、试给出一种图像增强方法，使得图像中暗区域中的内容的对比度加大，以便
更好的分辨暗区域中的内容。
16、试简述中值滤波的特性和适用场合。
中值滤波是一种非线性滤波，由于它在实际运算过程中并不需要图像的统计特性，所以比较方便。中值滤波首先是被应用在一维信号处理技术中，后来被二维图像信号处理技术所应用。在一定的条件下，可以克服线性滤波器所带来的图像细节模糊，而且对滤除脉冲干扰及图像扫描噪声最为有效。中值滤波的目的是保护图像边缘的同时去除噪声。
17、试写出孔洞填充的算法。对二值图像中所有被白色区域包围（封闭）的黑色
像素即为孔洞。
18、设有两个白色区域，被一条细小的白线所连接，试设计一种算法，消除两个
区域之间的细线，使两个区域分开。
练习题 1、图像灰度量化用6比特编码时，量化等级为（ B ） A 32个 B 64个 C128个 D 256个 2.下面说法正确的是：（ B ） A、基于像素的图像增强方法是一种线性灰度变换； B、基于像素的图像增强方法是基于空间域的图像增强方法的一种； C、基于频域的图像增强方法由于常用到傅里叶变换和傅里叶反变换，所以总比基于图像 域的方法计算复杂较高； D、基于频域的图像增强方法比基于空域的图像增强方法的增强效果好。 3、采用幂次变换进行灰度变换时，当幂次取大于1时，该变换是针对如下哪一类图像进行增强。（ B ） A 图像整体偏暗 B 图像整体偏亮 C图像细节淹没在暗背景中D图像同时存在过亮和过暗背景 4、采用模板［-1 1］T主要检测（ A ）方向的边缘。 A.水平 B.45 C.垂直 5、下列算法中属于图象锐化处理的是（ C ） A.低通滤波 B.加权平均法 C.高通滤波 D. 中值滤波 6、维纳滤波器通常用于（C ） A、去噪 B、减小图像动态范围 C、复原图像 D、平滑图像 7、彩色图像增强时，（ C ）处理可以采用RGB彩色模型。 A. 直方图均衡化 B. 同态滤波 C. 加权均值滤波 D. 中值滤波 8、（ B ）滤波器在对图像复原过程中需要计算噪声功率谱和图像功率谱。 A. 逆滤波 B. 维纳滤波 C. 约束最小二乘滤波 D. 同态滤波 9、高通滤波后的图像通常较暗，为改善这种情况，将高通滤波器的转移函数加上一常数量以便引入一些低频分量。这样的滤波器叫B ）。 A. 巴特沃斯高通滤波器 B. 高频提升滤波器 C. 高频加强滤波器 D. 理想高通滤波器
《遥感数字图像处理》习题与答案 第一部分 1.什么是图像？并说明遥感图像与遥感数字图像的区别。 答：图像（image）是对客观对象的一种相似性的描述或写真。图像包含了这个客观对象的信息。是人们最主要的信息源。 按图像的明暗程度和空间坐标的连续性划分，图像可分为模拟图像和数字图像。模拟图像（又称光学图像）是指空间坐标和明暗程度都连续变化的、计算机无法直接处理的图像，它属于可见图像。数字图像是指被计算机储存，处理和使用的图像，是一种空间坐标和灰度都不连续的、用离散数字表示的图像，它属于不可见图像。 2.怎样获取遥感图像？ 答：遥感图像的获取是通过遥感平台搭载的传感器成像来获取的。根据传感器基本构造和成像原理不同。大致可分为摄影成像、扫描成像和雷达成像三类。 m= 3.说明遥感模拟图像数字化的过程。灰度等级一般都取2m（m是正整数），说明8时的灰度情况。 答：遥感模拟图像数字化包括采样和量化两个过程。 ①采样：将空间上连续的图像变换成离散点的操作称为采样。空间采样可以将模拟图像具有的连续灰度（或色彩）信息转换成为每行有N个像元、每列有M个像元的数字图像。 ②量化：遥感模拟图像经离散采样后，可得到有M×N个像元点组合表示的图像，但其灰度（或色彩）仍是连续的，不能用计算机处理。应进一步离散、归并到各个区间，分别用有限个整数来表示，称为量化。 m=时，则得256个灰度级。若一幅遥感数字图像的量化灰度级数g=256级，则灰当8 度级别有256个。用0—255的整数表示。这里0表示黑，255表示白，其他值居中渐变。由于8bit就能表示灰度图像像元的灰度值，因此称8bit量化。彩色图像可采用24bit量化，分别给红，绿，蓝三原色8bit，每个颜色层面数据为0—255级。 4.什么是遥感数字图像处理？它包括那些容？ 答：利用计算机对遥感数字图像进行一系列的操作，以求达到预期结果的技术，称作遥感数字图像处理。 其容有： ①图像转换。包括模数（A/D）转换和数模（D/A）转换。图像转换的另一种含义是为使图像处理问题简化或有利于图像特征提取等目的而实施的图像变换工作，如二维傅里叶变换、沃尔什-哈达玛变换、哈尔变换、离散余弦变换和小波变换等。 ②数字图像校正。主要包括辐射校正和几何校正两种。 ③数字图像增强。采用一系列技术改善图像的视觉效果，提高图像的清晰度、对比度，突出所需信息的工作称为图像增强。图像增强处理不是以图像保真度为原则，而是设法有选择地突出便于人或机器分析某些感兴趣的信息，抑制一些无用的信息，以提高图像的使用价值。 ④多源信息复合（融合）。 ⑤遥感数字图像计算机解译处理。 5.说明遥感数字图像处理与其它学科之间的关系。 答：应具备的基础理论知识有：数学、地学、信息论、计算机、GIS、现代物理学。 6.说明全数字摄影测量系统的任务和主要功能。目前，比较著名的全数字摄影测量系统有哪些？
数字图像处理 实验指导书 学院：通信与电子工程学院 专业：电子信息工程 班级： 学号： 姓名： XX理工大学
实验一 MATLAB数字图像处理初步 一、实验目的与要求 1．熟悉及掌握在MATLAB中能够处理哪些格式图像。 2．熟练掌握在MATLAB中如何读取图像。 3．掌握如何利用MATLAB来获取图像的大小、颜色、高度、宽度等等相关信息。 4．掌握如何在MATLAB中按照指定要求存储一幅图像的方法。 5．图像间如何转化。 二、实验原理及知识点 1、数字图像的表示和类别 一幅图像可以被定义为一个二维函数f(x,y),其中x和y是空间(平面)坐标，f 在任何坐标处(x,y)处的振幅称为图像在该点的亮度。灰度是用来表示黑白图像亮度的一个术语，而彩色图像是由单个二维图像组合形成的。例如，在RGB彩色系统中，一幅彩色图像是由三幅独立的分量图像(红、绿、蓝)组成的。因此，许多为黑白图像处理开发的技术适用于彩色图像处理，方法是分别处理三副独立的分量图像即可。 图像关于x和y坐标以及振幅连续。要将这样的一幅图像转化为数字形式，就要求数字化坐标和振幅。将坐标值数字化成为取样；将振幅数字化成为量化。采样和量化的过程如图1所示。因此，当f的x、y分量和振幅都是有限且离散的量时，称该图像为数字图像。 作为MATLAB基本数据类型的数值数组本身十分适于表达图像，矩阵的元素和图像的像素之间有着十分自然的对应关系。 图1 图像的采样和量化 根据图像数据矩阵解释方法的不同，MA TLAB把其处理为4类： 亮度图像(Intensity images) 二值图像(Binary images) 索引图像(Indexed images) RGB图像(RGB images)
名词：*数字图像，数字图像处理，图像采样，线性拉伸，高通滤波，低通滤波，中值滤波，特征空间，图像分析，图像分割 问答题：1、设一幅图像有如图所示直方图，对该图像进行直方图均衡化，写出均衡化过程，并画出均衡化后的直方图。若在原图像一行上连续8个像素的灰度值分别为：0、1、2、3、4、5、6、7，则均衡后，他们的灰度值为多少？ 如图为一幅16级灰度的图像。请写出均值滤波和中值滤波的3x3滤波器；说明这两种滤波器各自的特点；并写出两种滤波器对下图的滤波结果（只处理灰色区域，不处理边界）。（15分） 设一幅灰度图像，其目标和背景的像素点灰度呈正态分布，灰度直方图如图所示。其中：、分 别为目标点的灰度分布密度函数、均值；、分别为背景点的灰度分布密度函数、均值。并设目标点和背景点的方差均为，目标点个数和图像总像点数的比为1:2。T是根据最小误差准则确定的最佳阈值。（15分） 试证明：
1．根据所学过的图像处理和分析方法，设计一套算法流程来实现汽车牌照的定位和数字的识别（给出设计思想即可）。 1、如图所示，A和B的图形完全一样，其背景与目标的灰度值分别标注于图中， 请问哪一个目标人眼感觉更亮一些？为什么？（10分） 选择题： 图像灰度方差说明了图像哪一个属性。（B ） A 平均灰度 B 图像对比度 C 图像整体亮度D图像细节 下列算法中属于图象锐化处理的是：( C ) A.低通滤波 B.加权平均法 C.高通滤波 D. 中值滤波 高通滤波后的图像通常较暗，为改善这种情况，将高通滤波器的转移函数加上一常数量以便引入一些低频分量。这样的滤波器叫 B 。 A. 巴特沃斯高通滤波器 B. 高频提升滤波器 C. 高频加强滤波器 D. 理想高通滤波器 ( )7.下列图象边缘检测算子中抗噪性能最好的是： a.梯度算子 b.Prewitt算子 c.Roberts算子 d. Laplacian算子
第一章． 遥感概念 遥感（Remote Sensing，简称RS），就是“遥远的感知”，遥感技术是利用一定的技术设备和系统，远距离获取目标物的电磁波信息，并根据电磁波的特征进行分析和应用的技术。 遥感技术的原理 地物在不断地吸收、发射（辐射）和反射电磁波，并且不同物体的电磁波特性不同。 遥感就是根据这个原理，利用一定的技术设备和装置，来探测地表物体对电磁波的反射和地物发射的电磁波，从而提取这些物体的信息，完成远距离识别物体。 图像 人对视觉感知的物质再现。图像可以由光学设备获取，如照相机、镜子、望远镜、显微镜等；也可以人为创作，如手工绘画。图像可以记录、保存在纸质媒介、胶片等等对光信号敏感的介质上。随着数字采集技术和信号处理理论的发展，越来越多的图像以数字形式存储。因而，有些情况下“图像”一词实际上是指数字图像。 物理图像：图像是人对视觉感知的物质再现 数字图像：图像以数字形式存储。 图像处理 运用光学、电子光学、数字处理方法，对图像进行复原、校正、增强、统计分析、分类和识别等的加工技术过程。 光学图像处理 应用光学器件或暗室技术对光学图像或模拟图像(胶片或图片)进行加工的方法技术 数字图像处理 是通过计算机对图像进行去除噪声、增强、复原、分割、提取特征等处理的方法和技术。图像处理能做什么？（简答） 是通过计算机对图像进行去除噪声、增强、复原、分割、提取特征等处理的方法和技术。数字图像处理主要目的：提高图像的视感质量，提取图像中所包含的某些特征或特殊信息，进行图像的重建，更好地进行图像分析，图像数据的变换、编码和压缩，更好图像的存储和传输。数字图像处理在很多领域都有应用。 遥感图像处理(processing of remote sensing image data )是对遥感图像进行辐射校正和几何纠正、图像整饰、投影变换、镶嵌、特征提取、分类以及各种专题处理的方法。常用的遥感图像处理方法有光学的和数字的两种。
《数字图像处理》实验指导书 编写: 罗建军 海南大学三亚学院 10月
目录 一、概述 ....................................................................... 错误!未定义书签。 二、建立程序框架 ....................................................... 错误!未定义书签。 三、建立图像类 ........................................................... 错误!未定义书签。 四、定义图像文档实现图像读/写.............................. 错误!未定义书签。 五、实现图像显示 ....................................................... 错误!未定义书签。 六、建立图像处理类................................................... 错误!未定义书签。 七、实现颜色处理功能............................................... 错误!未定义书签。 (一) 亮度处理................................................................. 错误!未定义书签。 (二) 对比度处理............................................................. 错误!未定义书签。 (三) 色阶处理................................................................. 错误!未定义书签。 (四) 伽马变换................................................................. 错误!未定义书签。 (五) 饱和度处理............................................................. 错误!未定义书签。 (六) 色调处理................................................................. 错误!未定义书签。 八、实现几何变换功能............................................... 错误!未定义书签。 (一) 图像缩放................................................................. 错误!未定义书签。 (二) 旋转......................................................................... 错误!未定义书签。 (三) 水平镜像................................................................. 错误!未定义书签。 (四) 垂直镜像................................................................. 错误!未定义书签。 (五) 右转90度................................................................. 错误!未定义书签。 (六) 左转90度................................................................. 错误!未定义书签。 (七) 旋转180度............................................................... 错误!未定义书签。 九、实现平滑锐化功能............................................... 错误!未定义书签。 十、图像处理扩展编程............................................... 错误!未定义书签。
第一章引言 一.填空题 1. 数字图像是用一个数字阵列来表示的图像。数字阵列中的每个数字，表示数字图像的一个最小单位，称为__________。1. 像素 2. 数字图像处理可以理解为两个方面的操作：一是从图像到图像的处理，如图像增强等；二是____________________，如图像测量等。2. 从图像到非图像的一种表示 3. 数字图像处理可以理解为两个方面的操作：一是__________________，如图像增强等；二是从图像到非图像的一种表示，如图像测量等。3. 从图像到图像的处理 5. 数字图像处理包含很多方面的研究内容。其中，________________的目的是根据二维平面图像数据构造出三维物体的图像。5. 图像重建 二.简答题 1. 数字图像处理的主要研究内容包含很多方面，请列出并简述其中的4种。 ①图像数字化：将一幅图像以数字的形式表示。主要包括采样和量化两个过程。 ②图像增强：将一幅图像中的有用信息进行增强，同时对其无用信息进行抑制，提高图像的可观察性。 ③图像的几何变换：改变图像的大小或形状。 ④图像变换：通过数学映射的方法，将空域的图像信息转换到频域、时频域等空间上进行分析。 ⑤图像识别与理解：通过对图像中各种不同的物体特征进行定量化描述后，将其所期望获得的目标物进行提取，并且对所提取的目标物进行一定的定量分析。 2. 什么是图像识别与理解？ 图像识别与理解是指通过对图像中各种不同的物体特征进行定量化描述后，将其所期望获得的目标物进行提取，并且对所提取的目标物进行一定的定量分析。比如要从一幅照片上确定是否包含某个犯罪分子的人脸信息，就需要先将照片上的人脸检测出来，进而将检测出来的人脸区域进行分析，确定其是否是该犯罪分子。 3. 简述数字图像处理的至少3种主要研究内容。 ①图像数字化：将一幅图像以数字的形式表示。主要包括采样和量化两个过程。 ②图像增强：将一幅图像中的有用信息进行增强，同时对其无用信息进行抑制，提高图 像的可观察性。 ③图像的几何变换：改变图像的大小或形状。 ④图像变换：通过数学映射的方法，将空域的图像信息转换到频域、时频域等空间上进 行分析。 5. 简述图像几何变换与图像变换的区别。 ①图像的几何变换：改变图像的大小或形状。比如图像的平移、旋转、放大、缩小等， 这些方法在图像配准中使用较多。 ②图像变换：通过数学映射的方法，将空域的图像信息转换到频域、时频域等空间上进
遥感 第一章 1遥感数字图像；遥感数字图像的分类方式和对应类别。 （1）定义：遥感数字图像是数字形式的遥感图像。不同的地物能够反射或辐射不同波长的电磁波，利用这种特性，遥感系统可以产生不同的遥感数字图像。 （2）可见图像和不可见图像 单波段和多波段，超波段 数字图像和模拟图像 2遥感图像的成像方式（三大种：摄影、扫描、雷达）。 （1）摄影，扫描属于被动遥感 雷达属于主动遥感 （2）摄影：根据芦化银物质在关照条件下回发生分解这一机制，将卤化银物质均匀涂在片基上，制成感光胶片 扫描：扫描类遥感传感器逐点逐行地以时序方式获取的二维图像 雷达：由发射机向侧面发射一束窄波段，地物反射的脉冲，由无线接收后被接收机接收 3遥感图像的数字化（模数转换）过程——两大过程：采样、量化，名词解释。 采样：将空间上连续的图像变换成离散点的操作称为采样，即：图像空间位置的数字化。采样是空间离散。 量化：遥感模拟图像经离散采样后，可得到由Ｍ×Ｎ个像素点组合表示的图像，但其灰度（或彩色）仍是连续的，还不能用计算机处理。它们还要进一步离散并归并到各个区间，分别用有限个整数来表示，这称之为量化，即：图像灰度的数字化。量化属于亮度属性离散。 遥感图像数字化过程两个特点：亮度和空 4遥感数字图像的存储空间大小的计算。 图像的灰度级有：2,64,128,256 存储一幅大小为M*N，灰度量化位数G的图像，所需要的存储空间（图像数据量）为M*N*G(bit) 1B=8bit 1KB=1024B 1MB=1024KB 1GB=1024MB TM空间分辨：1,2,3,4,5,7为30米，6为120米 5遥感数字图像的分辨率（时间、空间、光谱、辐射分辨率）； （1）时间分辨率：指对同一地点进行遥感采样的时间间隔即采样的时间频率，也称重访周期空间分辨率：指图像像素所代表的相应地面范围的大小，空间分辨率愈高，像素所代表的范围愈小 光谱分辨率：光谱分辨率是指成像的波段范围，分得愈细，波段愈多，光谱分辨率愈高 辐射分辨率：是传感器区分反射或发射的电磁波辐射强度差异的能力。高辐射分辨率可以区分信号强度的微小差异。 （2）常见传感器和空间分辨率书17-18页 6遥感数字图像的数据（数据级别、数据存储格式、元数据定义） （1）数据级别： 0级产品：未经过任何校正的原始图像数据 1级产品：经过了初步辐射校正的图像校正 2级产品：经过了系统级的几何校正，即根据卫星的轨道和姿态等参数以及地面系统中的有关参数对原始数据进行几何校正。产品的几何精度由上述参数和处理模型决定。 3级产品：经过几何精校正，即利用地面控制点对图像进行了校正，使之具有了更精确的地理坐标信息。产品的几何精度要求在亚像素量级上。 不同点：不同级别的产品使用条件不同，但是他们都是数据的集合，是信息量的汇总。一般来说，都是由元数据和图像基本数据两部分数据汇总的结果。
数字图像处理程序
数字图像处理实验 图像处理实验（一）直方图 灰度变换是图像增强的一种重要手段，使图像对比度扩展，图像更加清晰，特 征更加明显。 灰度级的直方图给出了一幅图像概貌的描述，通过修改灰度直方图来得到图像 增强。 1、灰度直方图 （1）计算出一幅灰度图像的直方图 clear close all I=imread('004.bmp'); imhist(I) title('实验一（1）直方图'); （2）对灰度图像进行简单的灰度线形变换， figure subplot(2,2,1) imshow(I); title('试验2-灰度线性变换'); subplot(2,2,2) histeq(I); （3）看其直方图的对应变化和图像对比度的变化。 原图像 f(m,n) 的灰度范围 [a,b] 线形变换为图像 g(m,n),灰度范围[a’,b’]公式：g(m,n)=a’+(b’-a’)* f(m,n) /(b-a) figure subplot(2,2,1) imshow(I) J=imadjust(I,[0.3,0.7],[0,1],1); title(' 实验一（3）用g(m,n)=a’+(b’-a’)* f(m,n) /(b-a)进行变换 '); subplot(2,2,2) imshow(J) subplot(2,2,3) imshow(I) J=imadjust(I,[0.5 0.8],[0,1],1); subplot(2,2,4) imshow(J) (4) 图像二值化（选取一个域值，(5) 将图像变为黑白图像） figure subplot(2,2,1)
第六章图像的锐化处理 一.填空题 1. 在图像的锐化处理中，通过一阶微分算子和二阶微分算子都可以进行细节的增强与检测。垂直方向的微分算子属于________________。（填“一阶微分算子”或“二阶微分算子”） 2. 在图像的锐化处理中，通过一阶微分算子和二阶微分算子都可以进行细节的增强与检测。Roberts交叉微分算子属于________________。（填“一阶微分算子”或“二阶微分算子”） 3. 在图像的锐化处理中，通过一阶微分算子和二阶微分算子都可以进行细节的增强与检测。Sobel 微分算子属于________________。（填“一阶微分算子”或“二阶微分算子”） 4. 在图像的锐化处理中，通过一阶微分算子和二阶微分算子都可以进行细节的增强与检测。Priwitt微分算子属于________________。（填“一阶微分算子”或“二阶微分算子”） 5. 在图像的锐化处理中，通过一阶微分算子和二阶微分算子都可以进行细节的增强与检测。Laplacian微分算子属于________________。（填“一阶微分算子”或“二阶微分算子”） 6. 在图像的锐化处理中，通过一阶微分算子和二阶微分算子都可以进行细节的增强与检测。Wallis 微分算子属于________________。（填“一阶微分算子”或“二阶微分算子”） 7. 在图像的锐化处理中，通过一阶微分算子和二阶微分算子都可以进行细节的增强与检测。水平方向的微分算子属于________________。（填“一阶微分算子”或“二阶微分算子”） 8. 图像微分______________了边缘和其他突变的信息。（填“增强”或“削弱”） 9. 图像微分______________了灰度变化缓慢的信息。（填“增强”或“削弱”） 10. 图像微分算子______________用在边缘检测中。（填“能”或“不能”） 四.简答题 1. 图像中的细节特征大致有哪些？一般细节反映在图像中的什么地方？ 2. 一阶微分算子与二阶微分算子在提取图像的细节信息时，有什么异同？ 3. 简述水平方向的微分算子的作用模板和处理过程。 4. 简述垂直方向的微分算子的作用模板和处理过程。 5. 已知Laplacian微分算子的作用模板为：，请写出两种变形的Laplacian算子。解答: 1. 图像的细节是指画面中的灰度变化情况，包含了图像的孤立点、细线、画面突变等。孤 立点大都是图像的噪声点，画面突变一般体现在目标物的边缘灰度部分。 2. 一阶微分算子获得的边界是比较粗略的边界，反映的边界信息较少，但是所反映的边界 比较清晰；二阶微分算子获得的边界是比较细致的边界。反映的边界信息包括了许多的细节 信息，但是所反映的边界不是太清晰。 五.应用题 1. 已知Roberts算子的作用模板为：，Sobel算子的作用模板为： 。 设图像为：
█一、叙述常见的数字图像处理开发工具有哪些？各有什么特点？（10分） 答．目前图像处理系统开发的主流工具为Visual C++（面向对象可视化集成工具）和MATLAB的图像处理工具箱（Image Processing Tool box）。两种开发工具各有所长且有相互间的软件接口。 Microsoft公司的VC++是一种具有高度综合性能的面向对象可视化集成工具，用它开发出来的Win 32程序有着运行速度快、可移植能力强等优点。VC++所提供的Microsoft 基础类库MFC对大部分与用户设计有关的Win 32应用程序接口API进行了封装，提高了代码的可重用性，大大缩短了应用程序开发周期，降低了开发成本。由于图像格式多且复杂，为了减轻程序员将主要精力放在特定问题的图像处理算法上，VC++ 6.0提供的动态链接库ImageLoad.dll支持BMP、JPG、TIF等常用6种格式的读写功能。 Microsoft公司的VC++是一种具有高度综合性能的面向对象可视化集成工具，开发出来的Win 32程序有着运行速度快、可移植能力强等优点。VC++对WIN32应用程序接口API 进行了封装，提高了代码的重用性，缩短了应用程序开发周期，降低了开发成本。VC++6.0提供的动态链接库imageload.dll支持BMP、JPG、TIF等六种格式的读写功能。MATLAB的图像处理工具箱MATLAB是由MathWorks公司推出的用于数值计算的有力工具，是一种第四代计算机语言，它具有相当强大的矩阵运算和操作功能，力求使人们摆脱繁杂的程序代码。MATLAB图像处理工具箱提供了丰富的图像处理函数，灵活运用这些函数可以完成大部分图像处理工作，从而大大节省编写低层算法代码的时间，避免程序设计中的重复劳动。MATLAB图像处理工具箱涵盖了在工程实践中经常遇到的图像处理手段和算法，如图形句柄、图像的表示、图像变换、二维滤波器、图像增强、四叉树分解域边缘检测、二值图像处理、小波分析、分形几何、图形用户界面等。但是，MATLAB也存在不足之处限制了其在图像处理软件中实际应用。首先，强大的功能只能在安装有MATLAB系统的机器上使用图像处理工具箱中的函数或自编的m文件来实现。其次，MATLAB使用行解释方式执行代码，执行速度很慢。第三，MATLAB擅长矩阵运算，但对于循环处理和图形界面的处理不及C++等语言。为此，通应用程序接口API和编译器与其他高级语言（如C、C++、Java等）混合编程将会发挥各种程序设计语言之长协同完成图像处理任务。API支持MATLAB与外部数据与程序的交互。编译器产生独立于MATLAB环境的程序，从而使其他语言的应用程序使用MATLAB。 MATLAB的图像处理工具箱MATLAB是由MathWorks 公司推出的用于数值计算的有力工具，是一种第四代计算机语言，它具有相当强大的矩阵运算和操作功能。MATLAB 图像处理工具箱提供了丰富的图像处理函数，灵活运用这些函数可以完成大部分图像处理工作。MATLAB图像处理工具箱涵盖了在工程实践中经常遇到的图像处理手段和算法。它的缺点有：1、强大的功能只能在安装有MATLAB 系统的机器上使用图像处理工具箱中的函数或自编的m文件来实现。2、MATLAB使用行解释方式执行代码，执行速度很慢。3、MATLAB擅长矩阵运算，但对于循环处理和图形界面的处理不及C++等语言。█二、叙述常见的数字图像应用软件有哪些？各有什么特 点？（10分） 答：图像应用软件是可直接供用户使用的商品化软件。用户 从使用功能出发，只要了解软件的操作方法就可以完成图像 处理的任务。对大部分用户来说，商品化的图像应用软件无 需用户进行编程，操作方便，功能齐全，已经能满足一般需 求，因而得到广泛应用。常用图像处理应用软件有以下几种： 1）PHOTOSHOP：当今世界上一流的图像设计与制作工具。 PHOTOSHOP已成为出版界中图像处理的专业标准。高版 本的PHOTOSHOP支持多达20多种图像格式和TWAIN接 口，接受一般扫描仪、数码相机等图像输入设备采集的图像。 PHOTOSHOP支持多图层的工作方式，只是PHOTOSHOP 的最大特色。使用图层功能可以很方便地编辑和修改图像， 使平面设计充满创意。利用PHOTOSHOP还可以方便地对 图像进行各种平面处理、绘制简单的几何图形、对文字进行 艺术加工、进行图像格式和颜色模式的转换、改变图像的尺 寸和分辨率、制作网页图像等。 2）CorelDRAW：一种基于矢量绘图、功能强大的图形图像 制作与设计软件。矢量式图像以几何、色彩参数描述图像， 其内容以线条和色块为主，数据量较小。CorelDraw是当今 流行的图像处理软件中为数不多的特点明显、功能强大的基 于矢量绘图的软件包。利用它，可以方便地制作精美的名片、 贺卡、书签、图书封面、广告、宣传画等作品。 3）ACDSee：快速、高性能的看图程序，是目前最享盛名 的图片浏览器。它能广泛应用于图片的获取、管理、浏览和 优化，支持BMP、GIF、JPG、TGA、TIF等超过50种常见 的图形文件格式，图片打开速度极快，可以直接查看动画 GIF，处理如Mpeg之类常用的视频文件，还可以为每一个 目录建立一个相册。ACDSee可以从数码相机和扫描仪高效 获取图片，并进行便捷的查找、组织和预览。ACDSee可以 轻松处理数码影像。 █三、傅里叶变换、加窗傅里叶变换和小波变换的时间－频 率特性有什么不同？（10分） 解：傅里叶变换使得时间信号变成了频域信号，加窗傅里叶 变换使得时间信号变成了时频信号，但是窗口是固定的，小 波变换同样变成了视频信号，但是时频的窗口是变化的。 █四、用JPEG标准，对于576行×720列的CCIR601建议分辨 率的彩色图像，其亮度分量可分割成多少个子块，而两个 色差分量可分别分割成多少子块？（10分） 解：对于576行×720列的CCIR601建议分辨率的彩色图像， JPEG将其亮度分量分割成（576/8）×（720/8）=6480块。 两个色差分量都可分割成两组：（576/8）×（360/8）=3240 块。 █五、二维傅里叶变换的分离性有什么实际意义？（10分） 解：该性质表明，一个二维傅里叶变换可由连续两次 一维傅里叶变换来实现。实现的方法如下图所示： █六、有了离散傅里叶及其快速算法FFT，为什么还要提出 离散余弦算法DCT及其快速算法？为什么许多视频国际标 准将DCT作为帧内编码的基本压缩算法？（10分） 答：在所有的变换编码方案中，离散K-L变换是最佳变换， 理论价值较高，常常作为对其他变换特性进行评价的标准。 但此变换没有快速算法，在工程应用中受到限制。在次最 佳变换算法中，DFT和DCT都是常用的变换编码方法，它们 分别有快速算法：FFT和FCT。这两种方法相比较，DFT涉及 到复数运算，而DCT是实数变换具有十分吸引人的一些特 点：它是一种实数变换，计算量较小，其变换矩阵的基向 量很好地描述了人类视觉的相关性，且对于大多数图像来 说，该变换的压缩性能很接近离散K-L变换，而且其变换矩 阵与图像内容无关，另外由于它构造对称的数据序列，避 免了在图像边界处的跳跃及所引起的Gibbs效应，并且也有 快速算法，因而得到广泛的应用。作为准最佳变换，它已 成为一些静态图像、视频压缩国际标准（或建议）中的基 本处理模块。 █七、扫描仪的光学分辨率是600×1200线，一个具有5000 个感光单元的CCD器件，用于A4幅面扫描仪，A4幅面的纸张 宽度是8.3英寸，该扫描仪的光学分辨率是多少dpi？（10 分） 解：（1）600×1200线，其中前一个数字代表扫描仪的横向 分辨率，后一数字则代表纵向分辨率。 （2）dpi是指单位面积内像素的多少，也就是扫描精度， 目前国际上都是计算一英寸面积内像素的多少。光学分辨 率是扫描仪的光学部件在每平方英寸面积内所能捕捉到的 实际的光点数,是指扫描仪CCD 的物理分辨率，也是扫描仪 的真实分辨率，它的数值是由CCD的像素点除以扫描仪水平 最大可扫尺寸得到的数值。 每一个感光单元对应一个像素。由于CCD感光单元个 数为5000， 5000/8.3＝602 (dpi) █八、直方图均衡，若一个64×64的离散图像，灰度分成8 层，其灰度rk的值和分布情况如下：请绘制该图像的直方 图，并求经过直方图均衡后的图像的直方图。 （20分） 解： █九、一图像大小为640×480，256色。用软件工具SEA （version 1.3）将其分别转成24位色BMP，24位色JPEG， GIF（只能转成256色）压缩格式，24位色TIFF压缩格式， 24位色TGA压缩格式，得到的文件大小分别为：921,654字 节；17,707字节；177,152字节；923,044字节；768,136字 节。分别计算每种压缩图像的压缩比。（10分） 解：不计算较小的文件头和彩色查找表（LTU）的数据量， 原始图像的数据量为： 640×480×1 byte=307,200 byte。 经转换后各种格式的压缩比如下： 24位色BMP格式： 307,200/921,654=0.333（增加了冗余度） 24位色JPEG格式： 307,200/17,707=17.35 GIF压缩格式： 307,200/177,152=1.73 24位色TIFF压缩格式： 307,200/923,044=0.333（增加了 冗余度） 24位色TGA压缩格式： 307,200/768,136=0.400（增加了冗 余度） █十、用JPEG标准，对于576行×720列的CCIR601建议分辨 率的彩色图像，其亮度分量可分割成多少个子块，而两个 色差分量可分别分割成多少子块？（10分） 解：对于576行×720列的CCIR601建议分辨率的彩色图像， JPEG将其亮度分量分割成（576/8）×（720/8）=6480块。 两个色差分量都可分割成两组：（576/8）×（360/8）=3240 块。 █十一、对下面的图像采用基于区域灰度差进行区域增长， 给出灰度差值 T①＝1；T②＝2；T③＝3三种情况下的分割图像。（10分） 解： 十二、用4连通或8连通准则，判断如下图像中的目标。 (10 分) 1 1 1 0 0 0 0 0 解： 1 1 1 0 1 1 0 0 1 1 0 0 1 1 0 0 1 1 0 0 0 1 1 0 0 1 0 0 0 0 1 0 1 0 1 0 0 0 1 0 1 1 1 0 0 1 1 0 1 1 1 0 0 0 0 0 X 1) 2) 3) 四邻域：L4 = 1 1 1 0 0 0 0 0 1 1 1 0 3 3 0 0 1 1 0 0 3 3 0 0 1 1 0 0 0 3 3 0 0 1 0 0 0 0 3 0 2 0 2 0 0 0 3 0 2 2 2 0 0 3 3 0 2 2 2 0 0 0 0 0 八邻域：L8 = 1 1 1 0 0 0 0 0 1 1 1 0 2 2 0 0 1 1 0 0 2 2 0 0 1 1 0 0 0 2 2 0 0 1 0 0 0 0 2 0 1 0 1 0 0 0 2 0 1 1 1 0 0 2 2 0 1 1 1 0 0 0 0 0
遥感数字图像处理教程 第一章概论 1.1图像和遥感数字图像 1.1.1图像和数字图像 本书定义图像为通过镜头等设备得到的视觉形象 根据人眼的视觉可视性可将图像分为可视图像和不可视图像。可视图像有图片、照片、素描和油画等，以及用透镜、光栅和全息技术产生的各种可见光图像。不可见图像包括不可见光成像和不可测量值 按图像的明暗程度和空间坐标的连续性，可将图像分为数字图像和模拟图像。数字图像是指用计算机存储和处理的图像，是一种空间坐标和灰度不连续、以离散数字原理表达的图像。在计算机，数字图像表现为二维阵列，属于不可见图像。模拟图像指空间坐标和明暗程度连续变化的、计算机无法直接处理的图像，属于可见图像。 利用计算机技术，可以实现模拟图像和数字图像之间相互转换。把模拟图像转化为数字图像成为模/数转换，记作A/D转换； 数字图像最基本的单位是像素。像素是A/D转换中国的取样点，是计算机图像处理的最小单位；每个像素具有特定的空间位置和属性特征。 1.1.2遥感数字图像 遥感数字图像时数字形式的遥感图像。不同的地物能够反射或辐射不同长波的电磁波，利用这种特性，遥感系统可以产生不同的遥感数字图像。 遥感数字图像中的像素成为亮度值。亮度值的高低由遥感传感器所探测到的地物电磁波的辐射强度决定。由于地物反射或辐射电磁波的性质不同受大气的影响不同，相同地点不同图像的亮度值可能不同。 图像的每个像素对应三维世界中的一个实体、实体的一部分或多个实体。在太阳照射下，一些电磁波被这个实体反射，一些被吸收。反射部分电磁波到达传感器被记录下来，成为特定像素点的值。 1.2压感数字图像处理 1.2.1遥感数字图像处理概述 遥感数字图像处理是利用计算机图像处理系统对遥感图像中的像素进行系列操作的过程。遥感数字图像处理主要包括三个方面 1.图像增强，使用多种方法，如：灰度拉伸、平滑、瑞华、彩色合成、主成分变换K-T变换、代数运算、图像融合等压抑、去除噪声、增强整体图像或突出图像中的特定地物的信息，是图像更容易理解、解释和判读、 图像增强着重强调特定图像特征，在特征提取、图像分析和视觉信息的显示很有用。 2.图像校正：图像校正也成图像回复、图像复原，主要是对传感器或环境造成的退化图像进行模糊消除、噪声滤除、几何失真或非线性校正。 信息提取：根据地物光谱特征和几何特征，确定不同地物信息的提取规则。 1.2.2 遥感数字图像处理系统 数字图像处理需要借助数字图像处理系统来完成。一个完整的遥感数字图像处理系统包括硬件系统和软件系统两大部分。 1.硬件系统 包括计算机、数字化设备、大容量存储、显示器和输出设备以及操作台 1）计算机 是图像处理核心，大的存和高的CPU速度有助于加快处理的进度。 2）数字化设备
数字图像处理实验报告
实验一数字图像处理编程基础 一、实验目的 1. 了解MA TLAB图像处理工具箱； 2. 掌握MA TLAB的基本应用方法； 3. 掌握MA TLAB图像存储/图像数据类型/图像类型； 4. 掌握图像文件的读/写/信息查询； 5. 掌握图像显示--显示多幅图像、4种图像类型的显示方法； 6. 编程实现图像类型间的转换。 二、实验内容 1. 实现对图像文件的读/写/信息查询，图像显示--显示多幅图像、4种图像类型的显示方法、图像类型间的转换。 2. 运行图像处理程序，并保存处理结果图像。 三、源代码 I=imread('cameraman.tif') imshow(I); subplot(221), title('图像1'); imwrite('cameraman.tif') M=imread('pout.tif') imview(M) subplot(222), imshow(M); title('图像2'); imread('pout.bmp') N=imread('eight.tif') imview(N) subplot(223), imshow(N); title('图像3'); V=imread('circuit.tif') imview(V) subplot(224), imshow(V); title('图像4');
N=imread('C:\Users\Administrator\Desktop\1.jpg') imshow(N); I=rgb2gary(GRB) [X.map]=gary2ind(N,2) RGB=ind2 rgb(X,map) [X.map]=gary2ind(I,2) I=ind2 gary(X,map) I=imread('C:\Users\dell\Desktop\111.jpg'); subplot(231),imshow(I); title('原图'); M=rgb2gray(I); subplot(232),imshow(M); [X,map]=gray2ind(M,100); subplot(233),imshow(X); RGB=ind2rgb(X,map); subplot(234),imshow(X); [X,map]=rbg2ind(I); subplot(235),imshow(X); 四、实验效果
数字图像处理 一、填空题 1、数字图像的格式有很多种，除GIF格式外，还有jpg 格式、tif 格式。 2、图像数据中存在的有时间冗余、空间冗余、结构冗余、信息熵冗余、知识 冗余、视觉冗余。 3、在时域上采样相当于在频域上进行___延拓。 4、二维傅里叶变换的性质___分离性、线性、周期性与共轨对称性、__位 移性、尺度变换、旋转性、平均值、卷积。（不考） 5、图像中每个基本单元叫做图像元素；在早期用picture表示图像时就称为 像素。 6、在图象处理中认为线性平滑空间滤波器的模板越大,则对噪声的压制越 好 ;但使图像边缘和细节信息损失越多; 反之, 则对噪声的压制不好 ,但对图像的细节等信息保持好。模板越平,则对噪声的压制越好 ,但对图像细节的保持越差;反之,则对噪声的压制不好,但对图像细节和边缘保持较好。 7、哈达玛变换矩阵包括___+1 和___—1 两种矩阵元素。（不要） 8、对数变换的数学表达式是t = Clog ( 1 + | s | ) 。 9、傅里叶快速算法利用了核函数的___周期性和__对称性。（不要） 10、直方图均衡化的优点是能自动地增强整个图像的对比度。（不要） 二、选择题 ( d )1.一幅灰度级均匀分布的图象，其灰度范围在[0，255]，则该图象的信息量为： a. 0 .255 c ( c )2.采用模板［-1 1］主要检测____方向的边缘。 a.水平 b.45 c.垂直 ( c )3. 下列算法中属于图象平滑处理的是： a.梯度锐化 b.直方图均衡 c. 中值滤波增强 ( b )4.图象与灰度直方图间的对应关系是： a.一一对应 b.多对一 c.一对多 d.都不对 ( a )5.对一幅图像采样后，512*512的数字图像与256*256的数字图像相比较具有的细节。 a.较多 b.较少 c.相同 d.都不对 ( b )6.下列算法中属于点处理的是： a.梯度锐化 b.二值化 c.傅立叶变换 d.中值滤波 ( d )7.二值图象中分支点的连接数为： .1 c ( a )8.对一幅100100像元的图象，若每像元用８bit表示其灰度值，经霍夫曼编码后压缩图象的数据量为40000bit，则图象的压缩比为： :1 :1 c.4:1 :2 ( d )9.下列算法中属于局部处理的是： a.灰度线性变换 b.二值化 c.傅立叶变换 d.中值滤波 ( b )10.下列图象边缘检测算子中抗噪性能最好的是： a.梯度算子算子算子d. Laplacian算子