矢量转换成栅格

将矢量图形转换成栅格图形是地理信息系统（GIS）中常见的任务，通常称为矢量转栅格。这个过程可以通过许多GIS软件完成，例如ArcGIS、QGIS等。下面是使用Python和GDAL库完成矢量转栅格的简单示例代码：

from osgeo import gdal, ogr
import numpy as np

# 打开矢量文件
vector_file = ogr.Open('vector_file.shp')

# 获取矢量图层
layer = vector_file.GetLayer()

# 获取矢量图层的空间参考信息
spatial_ref = layer.GetSpatialRef()

# 创建栅格图层
pixel_size = 0.1  # 栅格像元大小
target_ds = gdal.GetDriverByName('GTiff').Create('raster_file.tif', 100, 100, 1, gdal.GDT_Byte)
target_ds.SetGeoTransform((0, pixel_size, 0, 0, 0, -pixel_size))
target_ds.SetProjection(spatial_ref.ExportToWkt())

# 将矢量转换成栅格
gdal.RasterizeLayer(target_ds, [1], layer, burn_values=[1])

# 获取栅格数据
band = target_ds.GetRasterBand(1)
raster_data = band.ReadAsArray()

# 显示结果
print(raster_data)

这个程序首先打开一个矢量文件，然后获取其图层和空间参考信息。接着，使用GDAL库创建一个栅格图层，并设置其像元大小和空间参考信息。最后，使用`gdal.RasterizeLayer()`函数将矢量图层转换成栅格，并将结果保存到栅格图层中。最终，程序读取栅格数据并显示结果。

请注意，这只是一个简单的示例程序，实际的矢量转栅格过程可能需要更复杂的处理，例如在转换过程中应用插值算法来处理栅格的像元值。