IDL实现TM遥感影像直方图统计（中值、均值、方差、众数及峰度系数计算）

1 代码

近红外波段：

Pro statistic_tmfn=dialog_pickfile(title='Plesse choose image:')tm_img=read_image(fn)sz=size(tm_img)n_columns=sz[2] & n_raws=sz[3];window,0,xsize=n_columns,ysize=n_raws;tvscl,tm_img,/true;tm_ref=reform(tm_img,)tm_rgb=tm_img[0:2,*,*]tm_nir=tm_img[3,*,*]tm_r=tm_img[2,*,*]tm_g=tm_img[1,*,*]tm_b=tm_img[0,*,*]tm_sf=[tm_nir,tm_r,tm_g];szr=size(tm_r);n_c=szr[2] & n_r=szr[3];statasticMED=median(tm_nir)MEA=mean(tm_nir)VAR=variance(tm_nir)KUR=kurtosis(tm_nir)His=histogram(tm_nir,nbins=100,locations=locations,omin=omim,omax=omax)Sel=where(His EQ max(His))MODE=tm_nir[Sel[0]]MED_label='Median='+string(MED,format='(f5.2)')MEA_label='Mean='+string(MEA,format='(f5.2)')VAR_label='Variance='+string(VAR,format='(f7.2)')KUR_label='Kurtosis='+string(KUR,format='(f5.2)')MOD_label='Mode='+string(MODE,format='(f5.2)')p1=plot(locations,His,linestyle=0,color='red',thick=2,dimensions=[600,400],$title='Histogram Of NIR Band',window_title='Histogram',xrange=[0,200],yrange=[0,10000])t1=text(0.20,0.80,MED_label,target=p1,font_size=12)t2=text(0.20,0.75,MEA_label,target=p1,font_size=12)t3=text(0.20,0.70,VAR_label,target=p1,font_size=12)t4=text(0.20,0.65,KUR_label,target=p1,font_size=12)t5=text(0.20,0.60,MOD_label,target=p1,font_size=12);window,0,xsize=n_columns,ysize=n_raws;tvscl,tm_img,/true;imgshow_r=image(tm_r,dimension=[n_c,n_r],margin=0,title='TM Image of R',window_title='True R');imgshow_false=image(tm_sf,dimension=[n_columns,n_raws],margin=0,title='TM Image of Standared False',window_title='False Image');imgshow_true.save,'true.jpeg',resolution=600;imgshow_false.save,'false.jpeg',resolution=600
End

注：以上代码是TM近红外波段统计与直方图显示，其他波段类似只需换一下如下代码，如统计红波段改为：

 ;statasticMED=median(tm_r)MEA=mean(tm_r)VAR=variance(tm_r)KUR=kurtosis(tm_r)His=histogram(tm_r,nbins=100,locations=locations,omin=omim,omax=omax)Sel=where(His EQ max(His))MODE=tm_r[Sel[0]]

p1=plot(locations,His,linestyle=0,color='red',thick=2,dimensions=[600,400],$title='Histogram Of R Band',window_title='Histogram',xrange=[0,200],yrange=[0,12000])

2 结果：

近红外波段直方图统计：

简书分享

此时的你，

无与伦比！

IDL实现TM遥感影像直方图统计（中值、均值、方差、众数及峰度系数计算）相关推荐

计算遥感影像变异系数以及像元值变化趋势的小软件
计算遥感影像变异系数以及像元值变化趋势的小软件最近由于学习上的需要,需要计算很多遥感影像的相关系数,使用envi软件波段计算比较繁琐,就在网上找了不少计算相关指数的python代码,并且自己写了一些 ...
【Matlab图像去噪】中值+均值+Lee+Kuan图像滤波【含源码 1179期】
一.代码运行视频(哔哩哔哩) [Matlab图像去噪]中值+均值+Lee+Kuan图像滤波[含源码 1179期] 二.matlab版本及参考文献 1 matlab版本 2014a 2 参考文献 [1] ...
使用numpy计算平均值，中值，方差，标准差
使用numpy计算平均值,中值,方差,标准差 import numpy as np dataset = np.random.rand(10,12) # 输出前两行 dataset[0:2] # 平均值 ...
图像局部均值、中值、方差求取结合opencv
结合opencv求取图像的局部均值.中值.方差: //src 为待处理图像 //indexrows 为图像遍历的行数 //indexcols 为图像遍历的列数 //meanv 保存均值 //ker 窗 ...
TM遥感影像波段/通道bands
遥感影像波段band 一.遥感影像波段 1.原理 2.举例说明二.TM影像各波段简介 1.TM影像概述 2.各波段影像特征 3.波段组合 4.类型提取 5.光谱差异三.遥感图像--多波段数据存储的 ...
（三）Landat_5 TM 遥感影像辐射标定以及大气表观辐射率计算
前言一.辐射定标辐射定标:将记录的原始DN值转换为辐射值,也就是将图像的亮度灰度值转换为绝对辐射亮度. 如下便是通过imtool函数查看Landsat_5 TM band5的原始亮度灰度值,其整体 ...
图像的熵、灰度平均值、灰度中值、方差
//图像的熵========================================================================================= ...
（二）Landsat_5 TM 遥感影像波段合成真假彩色影像
前言回顾以下Landsat系列卫星的发展史,可以发现Landsat_5卫星已经退役,作为世界上寿命最长的地球观测卫星,其在影像方面对人类的发展具有很深刻的意义.Landsat_5 TM共有7个波段, ...
ArcMap|珞珈一号遥感影像在ArcMap中的处理
1.珞珈一号下载北京地区夜光影像下载地址:http://59.175.109.173:8888/index.html (原来压缩包可以批量解压) 成像时间14点多是因为"卫星数据的时间一般 ...
（五）Landat_5 TM 遥感影像计算NDVI、MNDWI、NDBI以及地表温度反演
流程概况图一.植被覆盖指数(NDVI) 所谓植被指数,就是利用多波段遥感图像的可见光波段以及近红外波段相组合,形成能够反映突出植被分布状况和强度的指数.植被指数的求解方法不唯一,定义的植被指数已经多 ...

IDL实现TM遥感影像直方图统计（中值、均值、方差、众数及峰度系数计算）

1 代码

2 结果：

IDL实现TM遥感影像直方图统计（中值、均值、方差、众数及峰度系数计算）相关推荐

最新文章

热门文章