matlab 多个波段,MatLab读取ENVI图像统计多波段图像信息
在ENVI统计遥感多波段图像中每个波段的均值、方差、最大值、最小值是比较容易办到的,但是如果要处理多批的数据就没有那么方便了,这里转载一个MatLab读取ENVI图像(img+hdr)的程序,并且计算了相关系数。
之前我在利用MatLab读取ENVI图像里分享了一个MatLab读取ENVI图像的函数,这里可以用上这个函数,具体的请查看上篇文章。
在新的.m文件里面实现批处理程序,代码如下:
clc;
clear;
filestring='E:\郭\实验数据2\*.img';%计算不同的路径中的图像,只需更改这里
%______以下,对于没有显式扩展名的情况,添加扩展名'.img'_________________%
subfile=filestring;
changefile=strrep(subfile,'*.img','');%找到图像所在文件夹
cfile=dir(changefile);%找到该文件夹下所有文件
for k=1:length(cfile)
if (size(strfind(cfile(k).name,'.'))==0)%判断文件名中有没有'.',如果没有意味着没有扩展名
copyfile(strcat(changefile,cfile(k).name),strcat(changefile,strcat(cfile(k).name,'.img')));
%上面这句,为没有'.img'的添加上
%delete(strcat(changefile,cfile(k).name));%删除多余数据,可选
end
end
%_________________________________________________________________%
filedir=strrep(filestring,'*.img','');
file=dir(filestring);
filenum=length(file);
%___以下代码计算每个图像在每个波段的均值,以及每个波段上所有图像均值的均值___%
fid=fopen(strcat(filedir,'均值.txt'),'wt');
fprintf(fid,'%s', '每个图像在每个波段的均值');
total=zeros(52,14);
for k=1:filenum
tempfile=strcat(filedir,file(k).name);
imgdata=read_ENVIimagefile(tempfile);
fprintf(fid,'\n%s',strcat('file',cfile(k).name,':'));
for i=1:14%事先已经知道共有14个波段,这点不好,需要对read_ENVIimagefile进行改造
meanvalue=[];
imgband=imgdata(:,:,i);
[x,y]=size(imgband);
reband=reshape(imgband,x*y,1);
countmean=mean(reband);
total(k,i)=countmean;
countmean=num2str(countmean);
meanvalue=[meanvalue countmean ' '];
meanvalue=num2str(meanvalue);
fprintf(fid,'%s',meanvalue);
end
end
fclose(fid);
fid1=fopen(strcat(filedir,'所有图像均值的均值.txt'),'wt');
fprintf(fid1,'所有图像均值的均值:');
fprintf(fid1,'%s',num2str(mean(total)));
fclose(fid1);
%___以下代码计算每个图像在每个波段的最大值,以及每个波段上所有图像的最大值的均值___%
fid=fopen(strcat(filedir,'最大值.txt'),'wt');
total=zeros(52,14);
for k=1:filenum
tempfile=strcat(filedir,file(k).name);
imgdata=read_ENVIimagefile(tempfile);
fprintf(fid,'\n%s',strcat('file',cfile(k).name,':'));
for i=1:14
maxvalue=[];
imgband=imgdata(:,:,i);
[x,y]=size(imgband);
reband=reshape(imgband,x*y,1);
countmax=max(reband);
total(k,i)=countmax;
countmax=num2str(countmax);
maxvalue=[maxvalue countmax ' '];
maxvalue=num2str(maxvalue);
fprintf(fid,'%s',maxvalue);
end
end
fclose(fid);
fid1=fopen(strcat(filedir,'所有图像最大值的均值.txt'),'wt');
fprintf(fid1,'所有图像最大值的均值:');
fprintf(fid1,'%s',num2str(mean(total)));
fclose(fid1);
%___以下代码计算每个图像在每个波段的最小值,以及每个波段上所有图像的最小值的均值___%
fid=fopen(strcat(filedir,'最小值.txt'),'wt');
total=zeros(52,14);
for k=1:filenum
tempfile=strcat(filedir,file(k).name);
imgdata=read_ENVIimagefile(tempfile);
fprintf(fid,'\n%s',strcat('file',cfile(k).name,':'));
for i=1:14
minvalue=[];
imgband=imgdata(:,:,i);
[x,y]=size(imgband);
reband=reshape(imgband,x*y,1);
countmin=min(reband);
total(k,i)=countmin;
countmin=num2str(countmin);
minvalue=[minvalue countmin ' '];
minvalue=num2str(minvalue);
fprintf(fid,'%s',minvalue);
end
end
fclose(fid);
fid1=fopen(strcat(filedir,'所有图像最小值的均值.txt'),'wt');
fprintf(fid1,'所有图像最小值的均值:');
fprintf(fid1,'%s',num2str(mean(total)));
fclose(fid1);
%___以下代码计算每个图像在每个波段的方差,以及每个波段上所有图像的方差的均值___%
fid=fopen(strcat(filedir,'方差.txt'),'wt');
total=zeros(52,14);
for k=1:filenum
tempfile=strcat(filedir,file(k).name);
imgdata=read_ENVIimagefile(tempfile);
fprintf(fid,'\n%s',strcat('file',cfile(k).name,':'));
for i=1:14
meanvalue=[];
imgband=imgdata(:,:,i);
[x,y]=size(imgband);
reband=reshape(imgband,x*y,1);
countvar=var(reband);
total(k,i)=countvar;
countvar=num2str(countvar);
varvalue=[meanvalue countvar ' '];
varvalue=num2str(varvalue);
fprintf(fid,'%s',varvalue);
end
end
fclose(fid);
fid1=fopen(strcat(filedir,'所有图像方差的均值.txt'),'wt');
fprintf(fid1,'所有图像方差的均值:');
fprintf(fid1,'%s',num2str(mean(total)));
fclose(fid1);
matlab 多个波段,MatLab读取ENVI图像统计多波段图像信息相关推荐
- matlab 读取ENVI中*.img格式图片代码
function data=read_ENVIimagefile_img(imgfilename) %本函数读取img格式,前提是img图像显式带有'.img'后缀名. if length(imgfi ...
- matlab数组保存envi,MATLAB 读取envi img格式
function data=read_ENVIimagefile_img(imgfilename) %本函数读取img格式,前提是img图像显式带有'.img'后缀名. if length(imgfi ...
- 【Matlab】Matlab读取dcm图像的函数以及CT值失真的问题处理
Matlab版本:2020a 一.dicomread函数 Matlab读取dcm图像的函数是dicomread,根据dicomread的帮助文档,该函数有四种参数输入方式: X = dicomread ...
- 【MATLAB】读取序列图像raw文件并求序列均值
读取序列图像raw文件并求序列均值 读取16位raw图像,图像大小为640*512,也可自行调整. matlab代码如下: %% 初始化 foldname='C:\Users\admin\Deskto ...
- 光谱分类算法 matlab,Matlab K-means聚类算法对多光谱遥感图像进行分类(一)
Matlab K-means聚类算法对多光谱遥感图像进行分类 作者: 白艺亭 测试了下matlab自带kmeans函数,作者编写函数,以及ENVI下的Kmeans方法,对比其效果,代码及结果图展示见下 ...
- Matlab图像处理创新实践-实验3【图像锐化】
Matlab图像处理创新实践-实验1[图像滤波基础(1)] Matlab图像处理创新实践-实验2[图像滤波基础(2)] Matlab图像处理创新实践-实验3[图像锐化] Matlab图像处理创新实践- ...
- Matlab图像处理创新实践-实验2【图像滤波基础(2)】
Matlab图像处理创新实践-实验1[图像滤波基础(1)] Matlab图像处理创新实践-实验2[图像滤波基础(2)] Matlab图像处理创新实践-实验3[图像锐化] Matlab图像处理创新实践- ...
- Matlab图像处理创新实践-实验1【图像滤波基础(1)】
Matlab图像处理创新实践-实验1[图像滤波基础(1)] Matlab图像处理创新实践-实验2[图像滤波基础(2)] Matlab图像处理创新实践-实验3[图像锐化] Matlab图像处理创新实践- ...
- JAVA imread_利用Matlab中的imread怎么读取图片
展开全部 方法/步骤 1.在matlab软件中,读取32313133353236313431303231363533e59b9ee7ad9431333337613137图像数据(载入)利用的是imre ...
最新文章
- 2018总结及2019计划
- 大数据针对高校_“大数据+移动支付” 安徽高校食堂智慧升级
- java.io.eof_java.io.IOException: java.io.EOFException: Unexpected end of input stream错误
- python历史以及基础知识
- Redis笔记4-持久化方案
- 专家:金融科技发展需行业自律和监管合力推动
- 怎样判断ajax请求,如何判断一个请求为ajax请求?
- ICEFaces Note(7)
- 统计学作业4 概念
- 学生成绩互帮互助分组C语言,C语言实现—学生成绩管理系统(Linux下运行)
- 语录(心灵鸡汤来一波)
- 轻量级web富文本框——wangEditor使用手册(3)——如何自定义配置菜单 demo
- 中国移动计算机岗专业知识,【中国移动】计算机技术岗面经
- c#.net配置mysql数据库连接池_c# 创建 mysql数据库连接池
- AcrGIS 做成本距离分析时提示ERROR 999999:无法启动配置 RasterCommander.ImageServer
- vs2019报错 声明与xxxx不兼容
- 前端工程化-我们需要做什么
- 不只卖电脑的联想,还有什么“新看点”?
- mysql on.000002_mysql | 同乐学堂
- 如何将网站发布上线,并配置自定义域名
热门文章
- 看看这几个版本Linux系统总有一款适合你
- http / 关于长连接和短链接的理解
- 数据结构与算法 / 霍夫曼树、霍夫曼编码和解码
- 乌班图linux配置yum仓库,Linux系统的仓库配置(yum、epel)与软件安装(wireshark,wget,nginx)...
- php实现附件上传下载,PHP实现文件上传与下载
- 计算机专业c语言题库,计算机专业C语言考试题.doc
- rj45接口引脚定义_简单修复数控主机网络接口通讯故障
- layui select 修改_layui修改select的值的方法
- server2016安装mysql_windows server2016安装MySQL5.7.19解压缩版教程详解
- echarts symbol 回调函数_凹函数和凸函数到底什么样?傻傻分不清楚