其實DTW演算法更多應用于手勢識別、語音信號處理等領域，但是在近年來，研究發現該演算法在遙感時序資料處理方面，也具有一定的優勢，例如不同的作物，雖然其NDVI時序曲線的變化趨勢相同，但是不同地方播種時間會有所差異，而DTW演算法剛好可以彌補這個差異，以達到更好的分類于提取效果，
采用2018-2020年MODIS 植被指數合成產品MOD13A2，該資料時間間隔為16d，每年23期，三年共69期，空間解析度為1KM，（來源：https://ladsweb.modaps.eosdis.nasa.gov/ ，選擇資料產品，選擇時間和區域，就能夠下載HDF檔案）
MODIS資料產品預處理：
MODIS資料下載的格式均為HDF檔案，不便處理，因此需要通過其他工具統一拼接投影，并轉換成較易于處理和使用的TIF或DAT格式，可以采用NASA官網提供的MODIS Reprojection Tool（MRT）對MODIS資料進行批量鑲嵌、投影處理，并轉為TIF格式，
HDF資料拼接、轉換與投影：
首先打開modis tool工具，并匯入任意一期的資料，進行引數的預設定，選擇需要提取的指數型別（本次為ndvi）,設定輸出格式和投影資訊等（本次選擇了UTM投影，并定義坐標系為WGS84，選擇49帶），之后保存設定檔案，點擊“Save Parameter File”，保存prm格式的工程檔案（后面批量處理時要用到），最后點擊“RUN”，
這樣就得到了一期拼接投影后的tif影像：
HDF資料批處理：
實際操作程序中往往資料量比較大，這時候可以通過批處理命令批量呼叫MRT工具處理資料，第一種方方法是通過寫dat檔案，呼叫保存的設定批量完成處理，不過當時操作的時候，在投影階段出了一點問題，第二種：通過直接呼叫cmd命令，相比而言呼叫cmd命令代碼量較少，比較方便，本次通過呼叫cmd命令完成資料批量處理，
按下Win+R，打開運行界面，輸入cmd，打開命令提示符，按順序輸入以下代碼：
cd/d D:\MRT\MRT_download_Win\bin //MRT工具檔案夾
java -jar MRTBatch.jar -d C:\Users\shao\Desktop\modisdata2018\DATA -p C:\Users\shao\Desktop\modisdata2018\mod18.prm -o C:\Users\shao\Desktop\modisdata2018\modis_mosaic //選擇資料檔案夾、工程設定檔案夾、結果存放檔案夾
D:\MRT\MRT_download_Win\bin\MRTBatch.bat //呼叫MRTBatch工具
這里注意，其中第一行位置，是你安裝MRT的位置，第二行 -d后面，輸入存放原始資料的檔案夾；-p后面，輸入存放上一步生成的prm工程檔案位置；-o后面，是輸出的檔案夾，最后一行輸入MRTBatch.bat即可，運行后，就得到了所有NDVI影像：
到這里，MODIS資料預處理環節就告一段落了，下一步是對資料進行濾波除噪、例外值處理，以及裁剪目標區域影像，后續更新！
轉載請註明出處，本文鏈接：https://www.uj5u.com/qita/300317.html
標籤：其他