python tif栅格数据投影转换

ronvicki

于 2023-11-13 20:32:41 发布

阅读量239

点赞数

分类专栏： Python 文章标签： python 开发语言

本文链接：https://blog.csdn.net/ronvicki/article/details/134385902

版权

Python 专栏收录该内容

26 篇文章 0 订阅

订阅专栏

本文介绍了在Python2.7环境下使用ArcGIS库进行地理信息系统中的操作，具体展示了如何获取tif文件列表，设置投影坐标系，以及使用ProjectRaster_management函数对tif文件进行投影转换的过程。

摘要由CSDN通过智能技术生成

python 2.7下运行

import arcpy
from arcpy import env

# 设置工作空间
env.workspace = r"H:\data"

# 获取所有tif文件的路径
rasters = arcpy.ListRasters("*", "TIF")

# 设置投影坐标系
out_coordinate_system = arcpy.SpatialReference("CGCS_2000")
# out_coordinate_system = arcpy.SpatialReference("WGS 1984 UTM Zone 11N")

# 循环处理每个tif文件
for raster in rasters:
    # 构建输出路径和文件名
    out_raster = arcpy.Describe(raster).baseName + "_projected.tif"
    out_path = r"H:\course\RSwork\1yearr\ndvi\0ori\5py\ndvi\newproj"
    out_raster_path = out_path + "\\" + out_raster

    # 投影转换
    arcpy.ProjectRaster_management(raster, out_raster_path, out_coordinate_system)