在向量資料的裁剪中，有兩種方式，一種是切割(clip)，一種是擦除(erase),形式如下，原始資料為a，裁剪資料為b，正向操作得到a∩b，反向操作得到a-b，如下圖
向量資料切割和刪除
在柵格資料的空間分析中，現有的工具只有類似於上述clip這種操作，即arcgis中的extract by mask操作，但是實際應用總會遇到需要對柵格資料做類似於erase的操作，即保留向量邊界以外的柵格資料。例如下面的中柵格資料為landsat8 遼東灣地區的影像，向量資料為陸地邊界，想要得到的結果是切割出海面上的資料。
目前有兩個思路。 思路1：先做正常的extract切割，向量範圍外的柵格資料設為無效值，然後讀取切割後的和原始的影像，只保留位置上切割後為無效值，原始為有效值的柵格，其餘設定為無效值。
python**實現
import os
import gdal
import numpy
GDAL 處理柵格資料
gdal 是柵格資料處理模組 ogr 是向量資料處理模組 gdal資料驅動，與ogr資料驅動類似，需要先建立某一型別的資料驅動，再建立響應的柵格資料集。一次性註冊所有的資料驅動，但是只能讀不能寫 gdal.allregister 單獨註冊某一型別的資料驅動，這樣的話可以讀也可以寫，可以新建資料集 d...
gdal 建立GeoTiff柵格資料
1 定義轉換引數 2private readonly double d transform 3 定義座標系統 4private readonly string d project geogcs wgs 84 datum wgs 1984 spheroid wgs 84 6378137,298.257...
ArcSDE柵格資料的儲存結構
乙個元資料表 sde raster columns 乙個業務表 business table 提供使用者檢視，具有柵格列 raster column，與向量資料業務表中的空間列對等 四個支撐表 對使用者遮蔽 1 sde blk 分塊表 定義柵格分塊 tile 2 sde bnd 波段表 定義波段 3...