数字图像处理:实验八 遥感图像增强
实验八 遥感图像增强
一、实验目的
- 熟悉MATLAB中图像增强的原理;
- 掌握遥感图像灰度变换,直方图调整等增强方法。
二、实验仪器及设备
计算机、Matlab图像处理软件、Landsat8遥感图像
三、实验内容及步骤
- 遥感图像头文件读取
图像对比度增强的方法可以分为两类,一类是直接对比增强方法,另一类是间接对比增强方法,直方图均衡化和线性拉伸是两种最常见的间接对比度增强方法,直方图均衡是通过使用累计函数对灰度值进行调整以实现对比度的增强。线性拉伸通过对比度拉伸对直方图进行调整,从而扩大前景和背景灰度的差别,以达到增强对比度的目的,这种方法可以利用线性或非线性的方法来实现,如指数变换,对数变换和线性拉伸等。
1.直方图均衡化
直方图均衡化处理的中心思想是将原始图像的灰度直方图从比较集中的某个灰度区间变成在全部灰度范围内的均匀分布,通过对图像进行非线性拉伸,重新分配图像像素值,使一定灰度范围内的像素数据大致相同,从而使得给定图像的直方图分布变成均匀分布。2.指数变换
指数变换可以将遥感图像的灰度值部分进行压缩,高灰度值部分进行扩展,以达到增强遥感影像高灰度部分的目的。3.对数变换
对数变换可以将遥感图像的低灰度值部分进行扩展,高灰度值部分进行压缩,以达到增强遥感影像低灰度部分的目的。4.线性拉伸
线性拉伸是为了突出感兴趣的目标或灰度区间相对抑制那些不感兴趣的灰度区域。
四、实验内容及步骤
- 直方图均衡化
src_img=data(:,:,6);
%转为8bit图像,以波段6为例
min_val=min(min(src_img));
max_val=max(max(src_img));
for i=1:linesfor j=1:samplessrc_img(i,j)=(src_img(i,j)-min_val)/(max_val-min_val)*255;end
end
src_img=uint8(src_img);
figure(1);
subplot(2,2,1),imshow(src_img),title('原图像');
subplot(2,2,2),imhist(src_img),title('原图像直方图');
matlab_eq=histeq(src_img);
subplot(2,2,3),imshow(matlab_eq),title('matlab直方图均衡化原图像');
subplot(2,2,4),imhist(matlab_eq),title('matlab直方图均衡化直方图');
图一 直方图均衡化
2. 指数变换
gamma=0.4;
c=1;
g2=data(:,:,6);
g2=c*(g2.^gamma);
%归一化
min_val=min(min(g2));
max_val=max(max(g2));
for i=1:linesfor j=1:samplesg2(i,j)=(g2(i,j)-min_val)/(max_val-min_val)*255;end
end
g2=uint8(g2);
figure(1);
subplot(2,2,1),imshow(src_img),title('原图像');
subplot(2,2,2),imhist(src_img),title('原图像直方图');
subplot(2,2,3),imshow(g2),title('增强图像');
subplot(2,2,4),imhist(g2),title('增强图像直方图');
图二 指数变换
3. 对数变换
c1=1.0;
g=double(src_img);
v1=10;
v2=100;
v3=200;
g11=c*log2(1+v1*g)/log2(v1+1);
g22=c*log2(1+v2*g)/log2(v2+1);
g33=c*log2(1+v3*g)/log2(v3+1);
%归一化
min_val1=min(min(g11));
max_val1=max(max(g11));
min_val2=min(min(g22));
max_val2=max(max(g22));
min_val3=min(min(g33));
max_val3=max(max(g33));
for i=1:linesfor j=1:samplesg11(i,j)=(g11(i,j)-min_val1)/(max_val1-min_val1)*255;g22(i,j)=(g22(i,j)-min_val2)/(max_val2-min_val2)*255;g33(i,j)=(g33(i,j)-min_val3)/(max_val3-min_val3)*255;end
end
g=uint8(g);
g11=uint8(g11);
g22=uint8(g22);
g33=uint8(g33);
figure(1);
subplot(2,2,1),imshow(g),title('原图像');
subplot(2,2,2),imshow(g11),title('增强图像v=10');
subplot(2,2,3),imshow(g22),title('增强图像v=100');
subplot(2,2,4),imshow(g33),title('增强图像v=2000');
图三 对数变换
4. 线性拉伸
g=double(src_img);
[m,n,w]=size(g);
figure(1);
imshow(src_img);title('原图像');
mid=mean(mean(g));
%横轴
fa=20;fb=120;
%纵轴
ga=100;gb=255;
[height,width]=size(g);
dst_img=uint8(zeros(height,width));
g=double(g);
%三段斜率
k1=ga/fa;
k2=(gb-ga)/(fb-fa);
k3=(255-ga)/(255-fa);
for i=1:heightfor j=1:widthif g(i,j)<=fadst_img(i,j)=k1*g(i,j);elseif fa<g(i,j) && g(i,j)<=fbdst_img(i,j)=k2*(g(i,j)-fa)+ga;elsedst_img(i,j)=k3*(g(i,j)-fb)+gb;endend
end
dst_img=uint8(dst_img);
j=dst_img;
figure(2);
imshow(j);title('线性拉伸图像');pixel_f=1:256;
pixel_g=zeros(1,256);
%三段斜率,小于1表示该段江北压缩
k1=ga/fa;
k2=(gb-ga)/(fb-fa);
k3=(255-ga)/(255-fa);
for i=1:256if i<= fapixel_g(i)=k1*i;elseif fa<i && i<=fbpixel_g(i)=k2*(i-fa)+ga;elsepixel_g(i)=k3*(i-fb)+gb;end
end
figure(3);
plot(pixel_f,pixel_g);
图四 线性拉伸
五、实验心得
1.学会了遥感影像的增强处理。
2. 通过合适的增强操作,可凸显出图像的一些细节信息,使得影像信息提取、图像解译变得方便。
数字图像处理:实验八 遥感图像增强相关推荐
- 数字图像处理实验八图像的傅里叶变换
自选一幅灰度图像,编写MATLAB程序,完成图像的傅里叶变换并显示图像傅里叶变换谱. 1.快速图像傅里叶变换函数 Y = fft2(X) 2.将图像频谱零频分量移动到图像频谱中心 Y = fftshi ...
- 数字图像处理实验三图像增强
一.实验目的 (1)了解图像增强的目的及意义,加深对图像增强的 感性认识,巩固所学的图像增强的理论知识和相 关算法. (2)熟练掌握直方图均衡化和直方图规定化的计算过 程. (3)熟练掌握空域滤波中常 ...
- 图像处理边缘增强matlab,数字图像处理实验 matlab 图像增强 边缘检测 图像操作.doc...
数字图像处理实验 matlab 图像增强 边缘检测 图像操作 实验1 点运算和直方图处理 实验目的 1. 掌握利用Matlab图像工具箱显示直方图的方法 2. 掌握运用点操作进行图像处理的基本原理. ...
- 数字图像处理实验(总计23个)汇总
以下这些实验中的代码全部是我自己编写调试通过的,到此,最后进行一下汇总. 数字图像处理实验(1):PROJECT 02-01, Image Printing Program Based on Half ...
- 数字图像处理matlab实验对图像复原,数字图像处理实验07图像的复原处理
数字图像处理实验 一.数字图像处理实验 实验七 图像的复原处理 一.实验目的 熟悉几种在实际应用中比较重要的图像复原技术,学会用MATLAB复原函数对退化图像进行复原处理. 二.实验内容 1.用点扩散 ...
- 数字图像处理实验四图像频域增强
一.实验目的 (1)了解图像增强的目的及意义,加深对图像增强的感性认识,巩固所学的图像增强的理论知识和相关算法. (2)熟练掌握低通.高通.带通.同态滤波器的使用方法,明确不同性质的滤波器对图像的影响 ...
- matlab数字图像实验报告,数字图像处理实验报告(matlab)
数字图像处理实验报告(matlab) 学院:自动化学院 班级:电081班 姓名:李林树 学号:40850099 2011年10月 实验一 直方图均衡化 一. 实验目的: 1. 熟悉图像数据在计算机中的 ...
- 【计算机视觉】数字图像处理(四)—— 图像增强
数字图像处理(四)-- 图像增强 图像增强的定义 图像增强方法 一.图像增强的点运算 (一)灰度变换 1. 线性变换 2. 分段线性变换 3. 非线性灰度变换 对数变换 指数变换 (二) 直方图修整法 ...
- 数字图像处理实验——Python语言实现
数字图像处理实验--Python语言实现 实验一:数字图像处理入门 实验二:直方图均衡 实验三:线性平滑和锐化--掩模法 实验四:非线性平滑--中值滤波 实验五:非线性锐化--梯度法 GitHub地址 ...
- 数字图像处理实验5图像复原
一.实验目的 (1)了解图像复原的目的及意义,加深对图像复原理论的认识. (2)掌握维纳滤波复原基本原理. (3)掌握约束最小二乘方复原方法. (4)掌握盲解卷积复原方法 二.实验内容 (1)维纳滤 ...
最新文章
- 第四章:2.串 -- 串的模式匹配算法(KMP)
- idea自动生成not null判断语句
- 达文西画中的数学密码
- 深拷贝的缺点_拷贝?还傻傻分不清深浅?
- 《恋上数据结构第1季》动态扩容数组原理及实现
- iOS开发初探篇——UITabBarController
- [Android]解决ClickableSpan中点击后ListView中item的长按冲突的问题
- 惯导系统模型及其仿真(五)
- Python乌龟吃鱼小游戏
- 实验吧ctf题库:这个看起来有点简单!
- 网易的猪场有多豪?网友:请你低调一点
- 不愿意跟同事分享的强大办公软件,超超超超级好用
- FT-TRN-BEG-C安装教程及问题解决
- input框监听输入法输入中文
- “爱心帮”APP走进厦门高校
- 静态HTML网页设计作品——生鲜超市网站设计(5页)HTML+CSS+JavaScript 学生DW网页设计作业成品 美食站
- Flash中使用Filereference上传文件的一些注意事项
- 使用easyExcel批量导入导出数据
- C语言中BCL的作用,恶性淋巴瘤基因检测结果中Bcl-2、Bcl-6和C-myc代表什么?什么是“双重打击”?...
- 崇尚个人英雄主义很难铸就优秀的网站设计公司