鴿了很久的up主今天又更新啦，這次給大家?guī)淼氖怯嘘P(guān)于如何用ENVI實(shí)現(xiàn)modis影像文件批量的幾何校正處理問題。

引言

幾何校正是指遙感成像過程中，受多種因素的綜合影響，原始圖像上地物的幾何位置、形狀、大小、尺寸、方位等特征與其對(duì)應(yīng)的地面地物的特征往往是不一致的，這種不一致就是幾何變形，也稱幾何畸變。

幾何校正是遙感中的專業(yè)名詞。一般是指通過一系列的數(shù)學(xué)模型來改正和消除遙感影像成像時(shí)因攝影材料變形、物鏡畸變、大氣折光、地球曲率、地球自轉(zhuǎn)、地形起伏等因素導(dǎo)致的原始圖像上各地物的幾何位置、形狀、尺寸、方位等特征與在參照系統(tǒng)中的表達(dá)要求不一致時(shí)產(chǎn)生的變形。

----------引自百度百科

舉個(gè)例子

?

在有關(guān)于遙感圖像的處理當(dāng)中，大家往往需要對(duì)某些遙感圖像進(jìn)行幾何校正。像我本人做有關(guān)于modis遙感影像的科創(chuàng)中，遇到的首要問題就是要對(duì)遙感圖像進(jìn)行幾何校正。在ENVI中，實(shí)現(xiàn)modis單幅圖像的幾何校正其實(shí)是一個(gè)相對(duì)簡(jiǎn)單的問題，但是如果想實(shí)現(xiàn)modis的批量ENVI幾何校正是較為困難的。因此我以MYD02HKM文件為大家進(jìn)行一個(gè)說明與介紹。（也許對(duì)真正的大佬而言批量的幾何校正也許是一個(gè)簡(jiǎn)簡(jiǎn)單單的問題，但對(duì)于像up主這樣的新手來說，這確實(shí)是一個(gè)很大的困難，因此，這篇專欄更多的在于分享自己的科創(chuàng)經(jīng)歷，學(xué)習(xí)成果之類的，希望各位大佬不喜勿噴。）

?

準(zhǔn)備工作

準(zhǔn)備工作方面僅需要大家的ENVI安裝MCTK。（自己網(wǎng)上找一下，有問題評(píng)論區(qū)里留言，與此同時(shí)MCTK是有用戶手冊(cè)的，雖然是全英文的，但是如果有耐心去看的話，在文末有使用IDL編程實(shí)現(xiàn)調(diào)用MCTK對(duì)單幅圖像進(jìn)行處理，本文也是基于這一方法實(shí)現(xiàn)的批量幾何校正。）

代碼

代碼簡(jiǎn)介及解析

我們?cè)谶@一部分里主要是為大家解釋一下代碼，這樣子大家遇到Bug或者是要進(jìn)行修改會(huì)方便一些。

1.?

這兩句代碼沒有太多意義，第一個(gè)的filesearch指向的是輸入文件夾（即存放數(shù)據(jù)的文件夾），以此類推第二個(gè)的filesearchx指向的是輸出文件夾（即輸出幾何校正后文件的文件夾）。
警告：輸入和輸出文件夾的路徑不宜過深，不然可能出現(xiàn)運(yùn)行無結(jié)果的情況。

?

2.?

這句代碼的用處是讀取出filesearch這一變量指向的文件夾中的所有hdf文件的文件名。

?

3.?

為什么要單獨(dú)把這一行代碼拎出來講呢，這句代碼的用處是讀取filesearchx這一變量指向的文件夾中所有文件名帶有a的文件名。這邊是為了實(shí)現(xiàn)每一份文件校正后直接存入單獨(dú)的文件夾。（如圖）

4.?

這句代碼是記錄該文件夾中的文件數(shù)，用于確定循環(huán)次數(shù)。

5.?

這一部分無太多需注意的，可以自己修改

6.?

這是指輸入MCTK的參數(shù)，下面的以此類推，具體可以在MCTK的說明書里找。

?

尾聲

本次的分享就到此為止了，大家如果有什么更好的代碼也希望能夠分享出來，大家互相學(xué)習(xí)促進(jìn)呀。希望大家多多支持我呀，有什么建議也給up提一下哈！

?