之前在運用哨兵2號數據做影像分析,發現網上很多操作的步驟比較零散,於是整理了一下整套的流程。
1 數據下載
https://scihub.copernicus.eu/數據下載網站
步驟1 框選需要下載影像的范圍、時間以及衛星類型,網站會自動選出這一時期覆蓋該區域的影響
步驟2:點擊其中的影像即可下載,Offline狀態的影響需要等其上線后才能下載
下載后的哨兵2號數據產品級別為Level-1C數據。要轉換成L2A級產品需要用戶自己進行處理生產。(L2A級數據主要包含經過輻射定標和大氣校正的大氣底層反射率數據(Bottom-of-Atmosphere corrected reflectance))
2 哨兵數據預處理(L1C級產品處理為L2A級)
L1C級數據轉L2A級數據主要是運用Sen2Cor插件,下載鏈接為:
http://step.esa.int/main/third-party-plugins-2/sen2cor/sen2cor_v2-8/
並且該插件必須安裝在用戶文件夾中
步驟3:下載好Sen2Cor插件后,雙擊L2A_Process.bat文件
步驟4:Win+R打開cmd,輸入插件的絕對路徑
步驟5:在cmd中輸入L2A_Process.bat --help,若出現以下內容,則為配置成功
步驟6:
檢驗配置環境好了以后,現在輸入L2A_Process "原始數據路徑"(就是L1C級產品數據的路徑),注意此處的原始數據路徑不能有空格和漢字,否則就會發生錯誤。
處理過程,Progress打到100%即為處理完成
最終得到的2A級產品
3 哨兵數據格式轉換,jp2轉TIFF
這一部分主要運用SNAP軟件完成。
SNAP下載地址:
http://step.esa.int/main/download/
步驟1:點擊file-open product,選擇需要轉格式輸出的數據.safe文件,選擇MTD開頭的.xml文件。
步驟2:右鍵點擊文件,選擇Open RGB window可以查看彩色圖像
步驟3:首先進行resampling。點擊raster- geometric operations- resampling,選擇I/0 parameters,確定保存格式和路徑;選擇resampling parameters,調節參數,點擊run。這一步主要是為了將所有波段的分辨率統一,因為哨兵2號的波段空間分辨率有10m、20m和60m,不便於后續分析。
設置重采樣的分辨率等
完成后SNAP左側會出現操作完成后的結果。選中它再點擊file - export- ENVI即可完成。轉出后的格式可以用ENVI打開。
轉換過程,要一定時間,而且轉出來的文件會非常大,差不多20個G。
但是轉出來的格式並不是TIFF,而是img光盤映像文件。因此還需要在ENVI中轉成TIFF格式數據。
4 ENVI處理哨兵數據(轉成tif格式)
步驟1:將上一部分中B1-B12.img文件導入至ENVI軟件
由於所有的波段是分開存儲的,這在后續分析中並不方便,因此需要合成為一個圖層導出為TIFF。
步驟2:
在toolbox中打開Layer Stacking,點擊Import File,將12個波段的圖層全部選中,點擊Reorder Files依據波段對圖片進行排序,選擇合適位置保存,點擊左下角OK。
生成時會有相當長一段時間的卡頓,因為生成了tif之后還會有一個ENP文件(金字塔影像)。生成了tif格式文件后,需要用特殊的方式打開,File、Open As、Optical Sensors、ADS40。
生成后的TIFF文件可以用ArcGIS等軟件打開進行后續分析。
作者:TuTu院長 https://www.bilibili.com/read/cv7764804/ 出處:bilibili