python gdal 矢量_在Python中用GDAL实现矢量对栅格的切割实例

概述：

本文讲述如何在Python中用GDAL实现根据输入矢量边界对栅格数据的裁剪。

效果：

裁剪前

矢量边界

裁剪后

实现代码：

# -*- coding: utf-8 -*-

"""

@author lzugis

@date 2017-06-02

@brief 利用shp裁剪影像

"""

from osgeo import gdal, gdalnumeric, ogr

from PIL import Image, ImageDraw

import os

import operator

gdal.UseExceptions()

# This function will convert the rasterized clipper shapefile

# to a mask for use within GDAL.

def imageToArray(i):

"""

Converts a Python Imaging Library array to a

gdalnumeric image.

"""

a=gdalnumeric.fromstring(i.tobytes(),'b')

a.shape=i.im.size[1], i.im.size[0]

return a

def arrayToImage(a):

"""

Converts a gdalnumeric array to a

Python Imaging Library Image.

"""

i=Image.frombytes('L',(a.shape[1],a.shape[0]),

(a.astype('b')).tobytes())

return i

def world2Pixel(geoMatrix, x, y):

"""

Uses a gdal geomatrix (gdal.GetGeoTransform()) to calculate

the pixel location of a geospatial