前言
最近辦公室的師兄在處理橫斷山區水土耦合的時候,遇到乙個問題,需要對柵格資料進行擦除,也就是反向提取。如果是向量資料,arcgis中提供直接的工具可以進行,但是柵格資料並沒有直接的工具。如果將柵格資料轉換為向量資料,則會遇到座標系、轉換方法和像元大小確定的問題,還會造成精度損失。網上搜尋的方法不甚明確,且較為複雜,現在提供一種直接簡便的方法。
第一步
首先確定需要被裁剪的柵格和擦除柵格圖層,我們旨在獲取擦除圖層空白部分所對應的被擦除圖層部分。
第二步
開啟柵格計算器,路徑為空間分析工具——地圖代數——柵格計算器**如下:
con(isnull(『擦除圖層』),『被擦除柵格圖層』,0)
**含義為當擦除圖層為空值時,賦值為被擦出圖層,其餘部分賦值為0。
第三步柵格計算的結果如下圖,可以發現除去我們柵格本身的值外,還有0值,其實這部分就是擦除柵格部分被賦值為0,所以我們需要刪除這部分柵格像元,才能得到最終的我們需要的柵格影象。
第四步
我們根據提取工具來選擇我們需要的柵格影象,工具路徑如下:空間分析工具——提取分析——按屬性提取,截圖如下:
可以發現,去除0值部分柵格,就得到了擦除反向裁剪的柵格影象。如圖:
注意事項
操作中,應該統一座標系和像元大小,並且去除擦除柵格和被擦除柵格的背景像元,令其為nodata。
使用gdal實現柵格資料的擦除(反裁剪)操作
在向量資料的裁剪中,有兩種方式,一種是切割 clip 一種是擦除 erase 形式如下,原始資料為a,裁剪資料為b,正向操作得到a b,反向操作得到a b,如下圖 向量資料切割和刪除 在柵格資料的空間分析中,現有的工具只有類似於上述clip這種操作,即arcgis中的extract by mask操...
ArcGIS批量裁剪柵格資料
arcpy批量計算柵格資料平均值 arcgis批量裁剪柵格資料 arcpy批量柵格重取樣 arcpy批量裁剪柵格資料 idl多程序批處理遙感資料 arcpy批量拼接柵格資料 做遙感資料處理,經常會用到,長時間序列的某些影像。這些資料的特點是 相同區域,不同時間,要使用同乙個shp來進行柵格裁剪。關於...
IDL實現向量 shp 裁剪柵格TASK(一
隨著envi idl版本的更新,idl對向量和柵格資料的處理也變得越來越簡單化。其提供了很多方便的介面,使得使用者呼叫和學習練習便捷成為了可能。最近接觸idl,發現好多網上的 都是延後的,新的介面 理解和編寫起來都比較方便,尤其是在做大量資料研究和應用時,使用批處理的方式顯得尤其重要。新的介面還在摸...