ArcGIS中重采樣栅格像元對齊問題

我們通常要進行基于像元的運算，往往我們的研究中涉及到多源資料，是以就需要對資料進行地理配準、空間配準、重采樣等操作。那麼當不同來源，不同分辨率的資料重采樣為同一空間分辨率之後，各個像元不一一對應有偏移，該怎麼辦呢？

轉發自我的知乎文章

我們通常要進行基于像元的運算，往往我們的研究中涉及到多源資料，是以就需要對資料進行地理配準、空間配準、重采樣等操作。

一開始，我認為相同的地理橢球與投影坐标系下，不同來源，不同分辨率的資料重采樣為同一空間分辨率之後，各個像元會一一對應。然而有的時候卻會遇到重采樣之後的像元大小雖然一樣，但像元之間會錯位，并不完全比對。

具體的原因我不太清楚，我猜測可能是重采樣算法内部的格網生成機制、原始資料分辨率與目标分辨率的比對（比如1000米重采樣成500米，因為除得盡，直接分成四個就行，但是如果重采樣成300米則除不盡，重采樣後的格網就不一定劃成什麼樣）

那麼這個問題該怎麼解決呢？目前在ArcGIS中進行重采樣操作時（resample 或者project raster）需要設定Environment --> Processing Extent --> Snap Raster 為基準栅格資料，這樣重采樣後資料的像元就會和基準栅格資料的像元一一對應，沒有位置偏移。

我的這篇部落格介紹了gdal中此類問題的解決辦法。

希望小夥伴有其它關于ENVI、GDAL等工具中此類問題的解決辦法進行分享哦！

參考：

https://link.zhihu.com/?target=http%3A//zhihu.esrichina.com.cn/question/13814

https://link.zhihu.com/?target=http%3A//zhihu.esrichina.com.cn/question/18622