任务:将一个文件夹中的所有栅格边界批量转为矢量格式
说明:仅将栅格的边界转为矢量的shape格式,栅格内部不转,nodata值或0值也算做有效值,即得到的是矩形矢量,该操作在大量数据处理、深度学习样本处理中有较多应用。
代码:
代码文件下载:https://download.csdn.net/download/X_Cosmic/12838015
import gdal
import ogr
import osr
import fnmatch
import os
ogr.RegisterAll()
img_path="D:/img/" #影像所在的文件夹
shape_path="D:/shp" #shape输出位置
if not os.path.exists(shape_path):
os.mkdir(shape_path)
img_list = fnmatch.filter(os.listdir(img_path), '*.tif')
for img in img_list:
p_img=img_path+'/'+img
outfilename = shape_path+'/'+img[:-4]+".shp"
dataset = gdal.Open(p_img)
oDriver = ogr.GetDriverByName('ESRI Shapefile')
oDS = oDriver.CreateDataSource(outfilename)
srs = osr.SpatialReference(wkt=dataset.GetProjection())
geocd = dataset.GetGeoTransform()
oLayer = oDS.CreateLayer("polygon", srs, ogr.wkbPolygon)
oDefn = oLayer.GetLayerDefn()
row = dataset.RasterXSize
line = dataset.RasterYSize
geoxmin = geocd[0]
geoymin = geocd[3]
geoxmax = geocd[0] + (row) * geocd[1] + (line) * geocd[2]
geoymax = geocd[3] + (row) * geocd[4] + (line) * geocd[5]
ring = ogr.Geometry(ogr.wkbLinearRing)
ring.AddPoint(geoxmin, geoymin)
ring.AddPoint(geoxmax, geoymin)
ring.AddPoint(geoxmax, geoymax)
ring.AddPoint(geoxmin, geoymax)
ring.CloseRings()
poly = ogr.Geometry(ogr.wkbPolygon)
poly.AddGeometry(ring)
outfeat = ogr.Feature(oDefn)
outfeat.SetGeometry(poly)
oLayer.CreateFeature(outfeat)
outfeat = None
oDS.Destroy()