用python写一段矢量转栅格的代码

# 以下是一段用python写的将矢量转栅格的代码：from osgeo import gdal, ogr# 定义输入和输出矢量图层路径
in_vector = "input.shp"
out_raster = "output.tif"# 打开输入矢量图层
in_ds = ogr.Open(in_vector)
in_layer = in_ds.GetLayer()# 创建栅格数据集
out_driver = gdal.GetDriverByName("GTiff")
out_ds = out_driver.Create(out_raster, xsize, ysize, 1, gdal.GDT_UInt16)# 设置投影
out_ds.SetProjection(in_layer.GetSpatialRef().ExportToWkt)# 设置坐标系
out_ds.SetGeoTransform(in_layer.GetSpatialRef().GetGeoTransform())# rasterize
gdal.RasterizeLayer(out_ds, [1], in_layer, burn_values=[1])# 关闭数据集
out_ds = None
in_ds = None

相关问题

Python矢量转栅格代码

以下是一个Python矢量转栅格的代码示例：

import geopandas as gpd
import rasterio
from rasterio.features import rasterize

# 读取矢量数据
vector_data = gpd.read_file('path/to/vector_data.shp')

# 定义栅格化参数
width = 1000
height = 1000
transform = rasterio.transform.from_origin(vector_data.total_bounds[0], vector_data.total_bounds[1], width, height)

# 创建空白栅格
raster = rasterio.open('path/to/raster.tif', 'w', driver='GTiff', width=width, height=height, count=1, dtype=rasterio.uint8, crs=vector_data.crs, transform=transform)

# 栅格化矢量数据
burn_value = 255
rasterized = rasterize([(geom, burn_value) for geom in vector_data.geometry], out_shape=(height, width), transform=raster.transform, fill=0, all_touched=True)

# 将栅格数据写入栅格文件
raster.write(rasterized, 1)
raster.close()

这段代码使用了`geopandas`和`rasterio`库来进行矢量转栅格操作。首先，通过`gpd.read_file`函数读取矢量数据文件。然后，定义栅格化参数，包括栅格的宽度、高度和变换矩阵。接下来，使用`rasterio.open`函数创建一个空白的栅格文件，并指定其属性，如文件格式、宽度、高度、波段数、数据类型、坐标系和变换矩阵。然后，使用`rasterize`函数将矢量数据栅格化，并将结果写入栅格文件中。

利用python写一段代码将arcgis中shp矢量图转换为分辨率为100的栅格图

谢谢您的提问！以下是示例代码：

import arcpy

# 设置环境变量
arcpy.env.workspace = r"C:\data"

# 定义输出栅格图的参数
outRaster = r"C:\output\myraster.tif"
outputResolution = 100

# 定义要转换的矢量图的路径和文件名
inShapefile = r"C:\data\myshape.shp"

# 进行矢量转栅格操作
arcpy.FeatureToRaster_conversion(inShapefile, "myattribute", outRaster, outputResolution)

其中，`arcpy` 是 `ArcGIS` 的 `Python` 模块，`env.workspace` 设置工作空间，`FeatureToRaster_conversion` 将矢量图转换为栅格图，`inShapefile` 是输入的矢量文件路径，`myattribute` 是需要转换的属性字段，`outRaster` 是输出的栅格图，`outputResolution` 是输出栅格图的分辨率。需要根据实际情况修改路径和参数。希望对您有所帮助！