ENVI/IDL编程:相对辐射校正-直方图匹配
參考博文http://blog.sina.com.cn/s/blog_764b1e9d0102vqws.html?提供的envimosaicraster方法,編寫批處理代碼
效果:測試結果與該博文的完全相同
主要的方法:ENVIMosaicRaster
數據格式:ENVI( image與 image.hdr)
PRO ZY3histogrammatch
? ? ;relative normalization
? ? ;全色影像之間進行相對輻射校正 直接用直方圖匹配的方法 NAD FWD BWD 因為是一個時相的數據
? ? ;每個時期應該有一個標準,大家都配準到某一個尺度上去
? ? ;20190214 修改代碼直方圖匹配。
? ? COMPILE_OPT IDL2
? ? e =ENVI()
? ? ;Determine input scenes?
? ? DIRPATH='待校正影像路徑'
? ? RESPATH='校正后影像存放路徑'
? ? refname='參考影像路徑\名稱'
? ? filelist = FILE_SEARCH(DIRPATH,'*.hdr')
? ? FOR i=0,n_elements(filelist)-1 DO BEGIN
? ? ? filename=STRMID(filelist[i],0,STRLEN(filelist[i])-4)
? ? ? basename=FILE_BASENAME(filename)?
? ? ? resname=respath+basename+'_hm'
? ? ? ;Test whether the current image has been processed?
? ? ? IF FILE_TEST(resname) EQ 1 THEN CONTINUE
? ? ? ? adjustraster=e.OPENRASTER(filename, DATA_IGNORE_VALUE=0)
? ? ? ? referraster=e.OPENRASTER(refname,DATA_IGNORE_VALUE=0)
? ? ? ? ;Create the mosaic raster
? ? ? ? Scenes = [adjustraster, referraster]
? ? ? ? MosaicRaster=ENVIMosaicRaster(Scenes)
? ? ? ? MosaicRaster.COLOR_MATCHING_METHOD = 'histogram matching'
? ? ? ? MosaicRaster.COLOR_MATCHING_STATS = 'entire scene'
? ? ? ? MosaicRaster.COLOR_MATCHING_ACTIONS = ['adjust','reference']
? ? ? ? ;Save it as ENVI format
? ? ? ? MosaicRaster.Export, resname, 'ENVI'
? ? ? ? ENDIF
? ? ENDFOR
? ? ;Close files in memory
? ? fids = ENVI_GET_FILE_IDS()
? ? size = SIZE(fids) ?;dimensions,respective length,data type, num_of_elements total
? ? length = size[1]
? ? FOR i = 0L, length-1 DO BEGIN
? ? ? ENVI_FILE_MNG,id = fids[i],/remove
? ? ENDFOR
END
總結
以上是生活随笔為你收集整理的ENVI/IDL编程:相对辐射校正-直方图匹配的全部內容,希望文章能夠幫你解決所遇到的問題。
- 上一篇: [react] React中怎么操作虚拟
- 下一篇: Spring Boot开发MongoDB