MATLAB图像处理与应用

最新推荐文章于 2023-11-07 19:19:30 发布
最新推荐文章于 2023-11-07 19:19:30 发布
阅读量2k 收藏 33
点赞数
分类专栏： MATLAB学习笔记
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/qq_42672745/article/details/109955963
版权
本文介绍了如何使用MATLAB进行遥感图像的读写操作，包括头文件处理、数据类型判断和数据读取。重点讲解了直方图均衡化、指数变换、对数变换和线性拉伸等图像增强技术，并展示了小波变换和内置融合函数在遥感图像融合中的应用实例。
摘要由CSDN通过智能技术生成

一、遥感图像的读写

相关函数：multibandread与multibandwrite
1.遥感图像头文件读取

fid=fopen('OLI2017_sub.hdr','r');
info=fread(fid,'char=>char');%字符型或字符型
info=info';
display(info);%%fprintf(info);
fclose(fid);

2.读取列数、行数、波段数、数据类型

  • 列数
ac=strfind(info,'samples = ');		%查找字符串返回位置的套路
bc=length('samples = ');
cc=strfind(info,'lines ');
samples=[];
for i=ac+bc:cc-1
    samples=[samples,info(i)];
end
samples=str2num(samples);
  • 行数
ar=strfind(info,'lines   = ');
br=length('lines   = ');
cr=strfind(info,'bands ');
lines=[];
for i=ar+br:cr-1
    lines=[lines,info(i)];
end
lines=str2num(lines);
  • 波段数
ab=strfind(info,'bands   = ');
bb=length('bands   = ');
cb=strfind(info,'data type');%换汤不换药，总之就是查找搜寻
bands=[];
for i=ab+bb:cb-1
    bands=[bands,info(i)];
end
bands=str2num(bands);
  • 数据类型同理，其中的data type
%case=1 对应ENVI中的数据类型为Byte,对应MATLAB中的数据类型为unit8
%case=2 对应ENVI中的数据类型为Integer,对应MATLAB中的数据类型为unit16
%case=12 对应ENVI中的数据类型为Unsigned Int,对应MATLAB中的数据类型为unit16
%case=3 对应ENVI中的数据类型为Long Integer,对应MATLAB中的数据类型为unit32
%case=13 对应ENVI中的数据类型为Unsigned Long,对应MATLAB中的数据类型为unit32
%case=4 对应ENVI中的数据类型为Floating Point,对应MATLAB中的数据类型为float32
%case=5 对应ENVI中的数据类型为Double Precision,对应MATLAB中的数据类型为double
%otherwise 除以上几种常见的数据类型之外的数据类型视为无效的数据类型
  • 数据格式同上

3.遥感图像数据读取与显示

%读取图像
data=multibandread('OLI2017_sub.dat',[lines,samples,bands],'double',0,'bil','ieee-le');%有问题啊。。。把[lines,samples,bands]参数值调小点就可以了
% % 错误使用 multibandread (line 118)
% % 文件太小，无法包含指定的数据。请检查大小、偏移和精度参数
data=double(data);
%数值转化为0~255的整型用于显示
data_unit8=data;
for k=1:bands
	min_val=min(data(:,:,k));
	max_val=max(data(:,:,k));
	for i=1:lines
		for j=1:samples
		data_unit8(i,j,k)=unit8((data_unit8(i,j,k)-min_val)/(max_val-min_val)*255;
		end
	end
end
%单波段遥感图像显示
im1=data_unit8(:,:,k);
figure;imshow(im1);
%真彩色显示
im3=data_unit8(:,:,6:8);
figure;imshow(im3);
%假彩色显示
i3m=data_unit8(:,:,5:7);
figure;inshow(i3m);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%另一种方法，用ENVI将影像转化为.tif的格式，然后
I=imread('sub_01.tif');
im1=I(:,:,3);
imshow(im1);

另：参考博客

二、遥感图像增强

  1. 直方图均衡化
  2. 指数变换
  3. 对数变换
  4. 线性拉伸
  • 直方图均衡化
I=imread('sub_01.tif');
im1=I(:,:,1);
[r, c, b]=size(im1);
%转化为8bit图像，以1波段为例
min_val=min(min(im1));
max_val=max(max(im1));
for i=1:r
    for j=1:c
        im1(i,j)=255*(im1(i,j)-min_val)/(max_val-min_val);
    end
end
im1=uint8(im1);
figure(1);
subplot(2,2,1),imshow(im1),title('原图');
subplot(2,2,2),imhist(im1),title('原图直方图');
eq=histeq(im1);
subplot(2,2,3),imshow(eq),title('直方图均衡化');
subplot(2,2,4),imhist(eq),title('直方图均衡化后的');
  • 指数变换
I=imread('sub_01.tif');
im1=I(:,:,1);
im0=im1;
[r, c, b]=size(im1);
Gamma=0.4;
C=1;
im1=C*(double(im1).^Gamma);     %double类型运算
%归一化
min_val=min(min(im1));
max_val=max(max(im1));
for i=1:r
    for j=1:c
        im1(i,j)=255*(im1(i,j)-min_val)/(max_val-min_val);
    end
end
im1=uint8(im1);
figure(1);
subplot(2,2,1),imshow(im0),title('原图');
subplot(2,2,2),imhist(im0),title('原图直方图');
subplot(2,2,3),imshow(im1),xlabel('\gamma=0.4');
subplot(2,2,4),imhist(im1),title('直方图均衡化后的');
  • 对数变换
I=imread('sub_01.tif');
im1=double(I(:,:,1));
[r, c, b]=size(im1);
C=1.0;
v1=10;v2=100;v3=200;
f1=C*log2(1+v1*im1)/log2(v1+1);     %double类型运算
f2=C*log2(1+v2*im1)/log2(v2+1);
f3=C*log2(1+v3*im1)/log2(v3+1);
%归一化
min_val1=min(min(f1));
max_val1=max(max(f1));
min_val2=min(min(f2));
max_val2=max(max(f2));
min_val3=min(min(f3));
max_val3=max(max(f3));
for i=1:r
    for j=1:c
        f1(i,j)=255*(f1(i,j)-min_val1)/(max_val1-min_val1);
        f2(i,j)=255*(f2(i,j)-min_val2)/(max_val2-min_val2);
        f3(i,j)=255*(f3(i,j)-min_val3)/(max_val3-min_val3);
    end
end
f1=uint8(f1);
f2=uint8(f2);
f3=uint8(f3);
figure(1);
subplot(2,2,1),imshow(I(:,:,1)),title('原图');
subplot(2,2,2),imshow(f1),title('v=10');
subplot(2,2,3),imshow(f2),xlabel('v=100');
subplot(2,2,4),imshow(f3),title('v=200');
  • 线性拉伸
I=imread('sub_01.tif');
im1=I(:,:,1);
[r, c, b]=size(im1);
%遥感图像读入的是uint16的需要转化为uint8然后再显示
imshow(im1);title('原图像');
mid=mean(mean(im1));
%横轴
fa=20;fb=120;
%纵轴
ga=100;gb=255;
dst_img=uint8(zeros(r,c));
im1=double(im1);
%三段斜率
k1=ga/fa;
k2=(gb-ga)/(fb-fa);
k3=(255-gb)/(255-fb);
for i=1:r
    for j=1:c
        if im1(i,j)<=fa
            dst_img(i,j)=k1*g(i,j);
        elseif fa<im1(i,j)&&im1(i,j)<=fb
            dst_img(i,j)=k2*(g(i,j)-fa)+ga;
        else
            dst_img(i,j)=k3*(im1(i,j)-fb)+gb;
        end
    end
end
dis_img=uint8(dst_img);
J=dst_img;
figure(2);imshow(J);title('线性变换图像');

三、遥感图像融合

%%1.小波变换融合
X1=data(:,:,6);
min_val=min(min(X1));
max_val=max(max(X1));
for i=1:lines
    for j=1:samples
        X1(i,j)=255*(X1(i,j)-min_val)/(max_val-min_val);
    end
end
X1=double(X1)/256;
figure;imshow(X1),title('高分辨率');
axis square;
X2=data(:,:,10);
min_val=min(min(X2));
max_val=max(max(X2));
for i=1:lines
    for j=1:samples
        X2(i,j)=255*(X2(i,j)-min_val)/(max_val-min_val);
    end
end
X2=double(X2)/256;
figure;imshow(X2),title('低分辨率');
axis square;
[c1,s1]=wavedec2(X1,2,'syms4');%将x1进行2维，使用'syms'进行变换
sizec1=size(c1);
for i=1:sizec1(2)
    c1(i)=1.2*c1(i);    %将分解后的值扩大1.2倍
end
[c2,s2]=wavedec2(X2,2,'syms4');%将x1进行2维，使用'syms'进行变换
c=(c1+c2)/2;
s1=(s1+s2)/2;
xx=waverec2(c,s,'sym4');    %进行重构
figure;imshow(xx),title('融合后的遥感图像');
axis square;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%2.利用自带融合函数实现融合
X1=data(:,:,6);
min_val=min(min(X1));
max_val=max(max(X1));
for i=1:lines
    for j=1:samples
        X1(i,j)=255*(X1(i,j)-min_val)/(max_val-min_val);
    end
end
X2=data(:,:,10);
min_val=min(min(X2));
max_val=max(max(X2));
for i=1:lines
    for j=1:samples
        X2(i,j)=255*(X2(i,j)-min_val)/(max_val-min_val);
    end
end
XFUS=wfusimg(X1,X2,'sym4',5,'max','max');
min_val=min(min(XFUS));
max_val=max(max(XFUS));
for i=1:lines
	for j=1:samples
		XFUS(i,j)=(XFUS(i,j)-min_val)/(max_val-min_val)*255;
	end
end
%显示融合后图像
X1=uint8(X1);
X2=uint8(X2);
XFUS=uint8(XFUS);
subplot(1,3,1),imshow(X1),axis square,title('高分辨率');
subplot(1,3,2),imshow(X2),axis square,title('低分辨率');
subplot(1,3,3),imshow(XFUS),axis square,title('融合后');

四、图像处理常用函数

  1. 读写图像文件

imread 、imwrite 、imfinfo 、multibandread 、multibandwrite

  1. 图像的显示

imshow 、colorbar 、figure

  1. 图像的变换

fft2二维傅里叶变换
ifft2二维傅里叶反变换
fft2可以进行卷积，然后用conv2校验

  1. 模拟噪声生成函数和预定义滤波器

imnoise 、fspecial

  1. 图像的增强

imhist直方图
histeq直方图均衡化
log对数变换
filter2基于卷积的图像滤波函数
conv2线性滤波
medfilt2中值滤波
锐化：Sobel算子、拉普拉斯算子

  1. 图像的分析

P=impixel(I)%交互式获取图像像素值
P=impixel(I,c,r)%指定点坐标像素值，c,r为行坐标和列坐标

im2uint8,im2double和uint8,double

在这里插入图片描述

PRSer_Carrot
关注
  • 0
    点赞
  • 踩
  • 33
    收藏
    觉得还不错? 一键收藏
  • 8
    评论
专栏目录
Matlab 图像处理
gotomic的专栏
10-20 2万+
 MATLAB知识var stattitle=MATLAB知识; yaaan17:40 PM fnChangeFloat(left,oImg)实验一 Matlab语言、数字图象基本操作一、实验目的1、复习MATLAB语言的基本用法；2、掌握MATLAB语言中图象数据与信息
matlab图像处理
dianmao0917的专栏
12-04 1936
Matlab中的图像数据类型转换  MATLAB中读入图像的数据类型是uint8，而在矩阵中使用的数据类型是double因此 I2=im2double(I1) :把图像数组I1转换成double精度类型；如果不转换，在对uint8进行加减时会产生溢出，可能提示的错误为：Function '*' is not defined for values of class 'uint8'。
【Matlab】读取".hdr"文件
zhangyonghui007的博客
05-08 1万+
写在最前 本文主要介绍了 Matlab 中函数 multibandread 的用法。 该函数可以用于读取 “.hdr” 文件。 一、函数：multibandread() % 该函数语法规则为：X = multibandread(filename, size, precision, offset, interleave, byteorder) % filename： 文件名 % size：图像尺寸和...
【MATLAB Image Processing Toolbox 入门教程三】快速入门之“在多光谱图像中寻找植被”
weixin_43024526的博客
10-30 1772
【MATLAB Image Processing Toolbox 入门教程三】本篇摘要一、从多光谱图像文件导入彩色红外通道二、构建近红外光谱散射图三、计算植被系数并显示其定位四、综合实例部分总结 本篇摘要 本篇演示了如何使用MATLAB数组对图像数据进行处理，并将数据绘制出来。图像数据差异可用于区分图像的不同表面特征，这些特征在不同光谱通道之间具有不同的反射率。示例部分采用了三维图像阵列，三个平面来自电磁频谱的不同信号，包含可见红光和近红外（NIR）通道，通过查找可见红光和NIR通道之间的差异，确定了示例
数字图像处理：实验七 遥感图像读取
chengzilhc的博客
05-03 4301
实验七 遥感图像读取 一、实验目的 熟悉MATLAB中遥感图像读写的格式； 掌握多波段遥感图像的读取与存储； 掌握多波段遥感图像中单波段、不同波段的组合显示。 二、实验仪器及设备 计算机、Matlab图像处理软件、Landsat8遥感图像 三、实验基础 1.实验遥感图像简介 2013年2月11日Landsat8卫星成功发射，设计寿命为五年，星上携带两个主要设备，陆地成像仪和热红外传感器。陆地成像仪共11个波段，其中可见光，近红外，短波红外波段空间分辨率为30米，热红外波段空间分辨率为100米，全色波段分
matlab 读写ENVI标准数据
最新发布
心网千劫
11-07 1592
本博客主要讲解如何读、生成ENVI标准格式的数据。主要分为四部分：读取ENVI头文件、读取ENVI数据、写入ENVI头文件、生成ENVI标准数据，最后附加讲解了本人写的生成hdr文本文件代码。
在做文件上传时.使用了MultipartFile时,文件太小导致错误的解决方案.
weixin_30517001的博客
08-16 514
今天出现了这么一个问题...当用户上传的文件太小时.导致文件上传失败，然后一顿乱找原因...发现 --- 当文件上传时在tomcat的临时路径里不存在.( 上传文件tomcat会存放在一个他自己的临时路径中,当代码结束时tomcat会自动销毁 )，所以我认为文件上传的前端出现了问题 因为( 虽然在后边拿到了文件名等数据，但是既然不存在于tomcat的临时路径，那么就没上传成功... )最终把问题锁...
MATLAB图像处理与应用案例
11-05
本资源“MATLAB图像处理与应用案例”提供了一系列的源代码，旨在帮助学习者深入理解和掌握MATLAB在图像处理中的应用。 一、图像基本操作 MATLAB提供了丰富的图像读取、显示和保存函数，如`imread`用于读取图像，`...
Matlab图像处理与应用
10-24
Matlab图像处理与应用 黄长虹 数字图像处理
MATLAB图像处理实例详解_MATLAB图像处理实例详解_matlab图像处理_
10-02
本文将深入探讨“MATLAB图像处理实例详解”这一主题，旨在帮助读者理解并掌握MATLAB在图像处理中的应用。 首先，我们要了解MATLAB的基础知识。MATLAB全称是“矩阵实验室”（Matrix Laboratory），它支持数值计算、...
matlab读取ENVI .img文件
03-01
高光谱图像中，matlab利用multibandread函数读取ENVI .img文件，并以.mat文件存储数据，用于我们后面的图像处理。
Matlab数据导入导出方法
01-18
Matlab数据导入导出方法,告诉你怎么从matlab中导入导出数据。
matlab读取并显示遥感dat图像
10-17
主要使用multibandread函数读取dat文件，显示单波段图像，多波段可以按波段显示，详细方法前面有解释，multibandread函数的参数可以根据自己的hdr文件输入
《MATLAB图像处理——能力提高与应用案例》程序.zip
05-18
《MATLAB图像处理——能力提高与应用案例》程序是一个压缩包，其中包含了多个与数字图像处理相关的MATLAB代码示例。这个资源对于那些希望通过实践来提升MATLAB在图像处理方面技能的学习者来说是非常宝贵的。MATLAB是...
img转换matlab,matlab怎么导入img格式数据
weixin_34627088的博客
03-17 478
哎，靠天靠地还是靠自己。。。。。。。。matlab读入img多波段图像数据用multibandread函数，语法为：X = multibandread(filename, size, precision, offset, interleave, byteorder)X = multibandread(...,subset1,subset2,subset3)ArgumentDescriptionfi...
MATLAB高光谱图像处理基础
热门推荐
TxyITxs的博客
06-25 1万+
介绍高光谱图像的基本知识，便通过MATLAB对高光谱图像进行基本的处理。 文章目录前言一、高光谱图像二、MATLAB高光谱图像处理1.加载.MAT文件数据2.图像的显示3.图像维度变换总结 前言 高光谱图像是一个立方体结构，维度为M x N x B，M为水平方向像素数量，N为垂直方向像素数量，B为每个像素对应的波段数量。 实验数据集： http://www.ehu.eus/ccwintco/index.php?title=Hyperspectral_Remote_Sensing_Scenes#Pavi
fcm模糊聚类matlab实例_MATLAB高光谱影像处理
weixin_39968946的博客
11-28 317
高光谱影像由于含有近百个波段，用matlab自带的图像读写函数imread和imwrite往往不能直接操作，以下针对高光谱影像读取、保存图像数组以及影像分块予以说明：影像读取：利用multibandread函数，此函数用于读取多波段二进制影像文件。至于何为二进制影像文件，简单地说ENVI主菜单file—save file as—envi standard得到的就是二进制影像文件，有时甚至会看到后缀...
使用MATLAB自动生成遥感影像数据的HDR文件
Zephyr Hou 的博客
07-17 3445
使用MATLAB自动生成遥感影像数据的HDR文件前言1. 头文件信息2. matlab 程序 欢迎学习交流： zengf.hou@bit.edu.cn https://www.jianshu.com/p/cb6d9c634380 声明：算法仅供参考，需要根据自身需求进行适当更改，不要拿来主义！ 前言 遥感影像数据的常用格式有很多，其中使用ENVI读取数据时，我们经常会发现很多数据都会带有 .HDR 后缀的文件，而且这个文件名称与影像数据同名，如果删除之后，影像就无法读取。其实这是遥感影像数据中的hdr头文件
[转载]matlab高光谱影像读取、保存、分块
weixin_30911451的博客
12-21 1713
原文地址：matlab高光谱影像读取、保存、分块作者：wind_静水流深_cloud 高光谱影像由于含有近百个波段，用matlab自带的图像读写函数imread和imwrite往往不能直接操作，以下针对高光谱影像读取、保存图像数组以及影像分块予以说明： 影像读取： 利用multibandread函数，此函数用于读取多波段二进制影像文件。至于何为二进制影像文件，简单地说ENVI主菜单fi...
评论 8
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值