ArcGIS中利用Arcpy批量对土地利用类型数据投影转换

最新推荐文章于 2024-04-05 09:47:30 发布

元宝kk

最新推荐文章于 2024-04-05 09:47:30 发布

阅读量761

点赞数

分类专栏： ArcGIS 土地利用数据 python 文章标签： arcgis python 经验分享

本文链接：https://blog.youkuaiyun.com/Yuanbaokk/article/details/129756151

版权

ArcGIS 同时被 3 个专栏收录

24 篇文章

订阅专栏

python

6 篇文章

订阅专栏

土地利用数据

4 篇文章

订阅专栏

该代码示例展示了如何利用Python的Arcpy模块对多个TIFF栅格文件进行批量投影转换。首先定义输入和输出文件夹，然后通过ListFiles获取所有TIFF文件，接着对每个文件进行坐标系定义和投影转换操作，使用NEAREST重采样方法。代码特别适用于跨地理坐标系的转换，如从GCS_WGS_1894到CGCS2000_3_Degree_GK_Zone_39。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

用模型构建器批量处理比较慢，因此可以利用Arcpy进行处理

# -*- coding: utf-8 -*-
import arcpy

# 多个shp文件所在文件夹
z_rasters = "H:/wd30/ww/sci/global30/z2001-2020/"
# 工作环境
arcpy.env.workspace = z_rasters
# 提取文件夹中的tif栅格文件
rasters = arcpy.ListFiles("*.tif")
# 输出raster所在路径
output_raster = "H:/wd30/ww/sci/global30/project/"
# 循环转换投影
for raster in rasters:
    # 获取名称
    catchment_name = raster[:-4]
    out_raster_name = output_raster + catchment_name + ".tif"
    # 投影坐标系WKID代码
    out_coor_system = 32639
    # 这个“30 30”是重采样的大小
    arcpy.ProjectRaster_management (raster, out_raster_name, out_coor_system, "NEAREST", "30 30")
    print "finished"

当数据需要进行跨地理坐标系进行投影转换时，例如我想从GCS_WGS_1894投影到CGCS2000_3_Degree_GK_Zone_39的坐标系，那么要首先进行一个地理坐标系的转换，这里可以直接进行定义坐标系，定义地理坐标系以后再进行投影转换。

# -*- coding: utf-8 -*-
import arcpy

# 多个shp文件所在文件夹
z_rasters = "H:/wd30/ww/sci/global30/z2001-2020/"
# 工作环境
arcpy.env.workspace = z_rasters
# 提取文件夹中的tif栅格文件
rasters = arcpy.ListFiles("*.tif")
# 输出raster所在路径
output_raster = "H:/wd30/ww/sci/global30/project/"
# 循环转换投影
for raster in rasters:
    # 获取名称
    catchment_name = raster[:-4]
    out_raster_name = output_raster + catchment_name + ".tif"
    # 这里的4490和4527是坐标系的WKID码
    arcpy.DefineProjection_management(raster, 4490)
    arcpy.ProjectRaster_management (raster, out_raster_name, 4527, "NEAREST", "30 30")
    print "finished"