tif数据与某省shp矢量文件有一致的坐标系

ronvicki2023-11-13 22:28

可以使用GDAL库中的Python绑定来读取tif数据并指定其坐标系为与某省.shp文件一致的坐标系。具体步骤如下：

1.导入必要的库

python 复制代码

from osgeo import gdal, osr

2.读取省.shp文件并获取其坐标系

python 复制代码

# 读取某省.shp文件
province_ds = ogr.Open('province.shp')
province_layer = province_ds.GetLayer()

python 复制代码

# 获取某省.shp文件的坐标系
province_srs = province_layer.GetSpatialRef()

3.读取tif数据并获取其坐标系

python 复制代码

# 读取tif数据
tif_ds = gdal.Open('data.tif')

python 复制代码

# 获取tif数据的坐标系
tif_srs = osr.SpatialReference()
tif_srs.ImportFromWkt(tif_ds.GetProjection())

4.将tif数据的坐标系转换为与省.shp文件一致的坐标系

python 复制代码

# 创建坐标系转换对象
transform = osr.CoordinateTransformation(tif_srs, province_srs)

进行坐标系转换

python 复制代码

transformed_tif_ds = gdal.AutoCreateWarpedVRT(tif_ds, None, tif_srs.ExportToWkt())

5.使用转换后的tif数据进行后续操作

python 复制代码

# 读取转换后的tif数据
transformed_tif_data = transformed_tif_ds.ReadAsArray()

6.进行后续操作