利用ENVI的ROI统计影像信息

首先在ENVI当中建立ROI，然后保存成XML格式

利用下面代码的

raster = e.OpenRaster(file)和

raster.ExportRois,infilepath+'incidence_angle_sample.csv' , rois, 'CSV'

就可以实现按照ROI统计的影像信息的CSV

打开生成的CSV后得到

RO1 #1 类别 npts（像素个数）: 320

RO1 #2 类别 npts（像素个数）: 320

B1、B2、B3为影像的波段1、2和3的像素值

最后用READ_CSV打开该文件

hh=roicsv.FIELD7; 读取B1
hv=roicsv.FIELD8;读取B2
angle=roicsv.FIELD9;

pro data_select_by_roiCOMPILE_OPT IDL2; Start the applicatione = ENVI(/HEADLESS)ys=2016;ye=2021;ms=1;me=4;for yr=ys,ye do beginfor mon=ms,me do begininfilepath='H:\Lead\data\sample\'+STRING(yr, FORMAT='(I04)')+'\'+STRING(mon, FORMAT='(I01)')+'\'cd,infilepathtifdata = file_search('S1*.tif',count = num);infilepath,if num eq 0 then beginprint,infilepath,' has no tif file.'continueendif; Open an input rasterfile = FILEPATH(tifdata[0], ROOT_DIR=infilepath);raster = e.OpenRaster(file); Open an ROI file
;      roiname=infilepath+'incidence_angle_sample.xml'file = FILEPATH('incidence_angle_sample.xml', ROOT_DIR=infilepath);, SUBDIRECTORY = ['data']rois = e.OpenRoi(file)if file_test(infilepath+'incidence_angle_sample.csv') eq 0 then beginraster.ExportRois,infilepath+'incidence_angle_sample.csv' , rois, 'CSV'print,'文件新建'endifroicsv = READ_CSV( infilepath+'incidence_angle_sample.csv' ,N_TABLE_HEADER=13,TYPES='Float')filenum=N_elements(roicsv.FIELD1);x=roicsv.FIELD1;y=roicsv.FIELD2;hh=roicsv.FIELD7;hv=roicsv.FIELD8;angle=roicsv.FIELD9;flag=0for i=0,filenum-1 do begin;搜集数据if i eq 0 or ((i mod 16) eq 0) then begintemp1=[hh[i]]temp2=[hv[i]]temp3=[angle[i]]endif else begintemp1=[temp1,hh[i]]temp2=[temp2,hv[i]]temp3=[temp3,angle[i]]endelse;判断是否需要提交结果if i ne 0 and  (N_elements(temp1) eq 16) then begin ;判断是否与前一个不为同ROI  abs(x[i-1]-x[i]) gt 5 and abs(y[i-1]-y[i])if flag eq 0 then begin ;判断是否第一次建立结果数组result1=mean(temp1)result2=mean(temp2)result3=mean(temp3)flag=1endif else beginresult1=[result1,mean(temp1)]result2=[result2,mean(temp2)]result3=[result3,mean(temp3)]          endelse      endif   endforif yr gt 2017 or (yr eq 2017 and mon ge 3) then beginfinal_ice=[[result1[0:19]],[result2[0:19]],[result3[0:19]]]final_ice=Transpose(final_ice)final_ow=[[result1[20:-1]],[result2[20:-1]],[result3[20:-1]]]final_ow=Transpose(final_ow)endif else beginfinal_ow=[[result1[0:19]],[result2[0:19]],[result3[0:19]]]final_ow=Transpose(final_ow)final_ice=[[result1[20:-1]],[result2[20:-1]],[result3[20:-1]]]final_ice=Transpose(final_ice)endelseWRITE_CSV,infilepath+'HH_HV_incidence_angle_statistics_ow.csv',final_owWRITE_CSV,infilepath+'HH_HV_incidence_angle_statistics_ice.csv',final_iceprint,'hello',infilepath+'HH_HV_incidence_angle_statistics.csv'endforendforend

利用ENVI的ROI统计影像信息相关推荐

利用电子计算机处理数字化的影像信息,影像诊断习题题库适合复习考试使用
<医学影像诊断学>试题库第一篇总论一.单选题(每题仅一个最佳答案) 1.X线在人体内的透过率从大到小,其正确排列为 ( E ) A.气体.液体及软组织.脂肪.骨 B.骨.脂肪.液体及 ...
利用电子计算机处理数字化的影像信息,医学影像学总资料
一.名词解释 1. 医学影像学:指包括常规X线检查.超声检查.核素显像诊断.CT和MRI诊断在内的医学影像. 2. 自然对比:人体组织自然存在的密度差别称自然对比. 3. 人工对比:对于缺乏自然对比 ...
利用ENVI裁剪sentinel 1A影像并生成ArcGIS可直接检索影像值的数据
对sentinel 1A雷达影像进行处理后会生成vv和vh两种数据,在ENVI中打开后利用cursor工具可查看相应影像值,但该数据格式不支持在ArcGIS中打开并查看. 处理好的雷达影像还需裁剪,使 ...
利用ENVI和ARCGIS做相关性分析的两种方法
一.ENVI里如何对两个波段进行相关性分析 1.首先获取要统计的两个波段的信息,这里用的ndvi和地表反照率(aldeo),两个波段值经过归一化后进行计算 2.首先利用layer stacking工具 ...
利用ENVI对LANDSAT系列数据进行预处理步骤
数据来源:地理空间数据云首先,利用envi打开 File>Open 下载好的Landsat数据中的MTL.txt文件进行辐射定标: 在工具栏Toolbox中搜索 Radiometric Co ...
基于微信小程序云开（统计学生信息并导出excel）2.0版
前言随着移动端的不断发展,人们大部分的办公及生活应用都开始趋向于移动端.然而在2017年"微信之父"张小龙带领团队,开发了一款叫做微信小程序的东西,它的出现打破了人们认识移动端的 ...
利用SNMP获取H3C交换机端口信息
利用SNMP获取H3C交换机端口信息 1.什么是Snmp SNMP是英文"Simple Network Management Protocol"的缩写,中文意思是"简单网 ...
利用envi对Landsat8图像进行NDVI处理
利用envi对Landsat8图像进行NDVI提取提取NDVI 选择一定范围的NDVI 提取NDVI 版本:envi4.7 第一步:File -> Open image File 第二步: B ...
DepthInSpace：多帧影像信息在单目结构光深度估计中的应用（ICCV2021）
点击上方"3D视觉工坊",选择"星标" 干货第一时间送达作者| flow,3D视觉开发者社区签约作者,毕业于中国科学院大学,人工智能领域优质创作者. 编辑| ...

利用ENVI的ROI统计影像信息

利用ENVI的ROI统计影像信息相关推荐

最新文章

热门文章