批量栅格投影（arcpy）

深藏阁楼爱情的钟 2022-09-17 04:23 244阅读 0赞

**栅格投影的意义：将分辨率单位为度的变为米  
定义投影的意义：将某个栅格的坐标系变为想要的坐标系**

批量栅格投影的代码如下：

# -*- coding: UTF-8 -*-
    #需要自行更改的如下：
    #inws：输入路径
    #outws：输出路径
    #Coordinate_System：目标投影坐标系
    import glob
    import os
    import arcpy
    
    # 输入路径  应该注意，中文路径，会导致读不出文件
    inws = r"------"
    
    # 参考文件路径    使用栅格数据集（从其导入方形像元大小）的路径
    Coordinate_System="PROJCS['UTM_Zone_46N',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['false_easting',500000.0],PARAMETER['false_northing',0.0],PARAMETER['central_meridian',93.0],PARAMETER['scale_factor',0.9996],PARAMETER['latitude_of_origin',0.0],UNIT['Meter',1.0]]"
    # 输出路径
    outws = r"----"
    
    # 利用glob包，将inws下的所有tif文件读存放到rasters中
    rasters = glob.glob(os.path.join(inws, "*.tif"))
    
    # 循环rasters中的所有影像，进行按掩模提取操作
    for raster in rasters:
        print str(raster)
        nameT = os.path.basename(raster).split(".")[0] + ".tif"	
        outname = os.path.join(outws, nameT)  # 合并输出文件名+输出路径
    
        arcpy.ProjectRaster_management(raster, outname, Coordinate_System, "NEAREST","#","#","#","#")
    
    print("OK!")

希望以上对大家有所帮助！

长路漫漫…  
唯有坚持…

发表评论取消回复

表情：

评论列表（有 0 条评论，244人围观）

还没有评论，来说两句吧...

相关阅读

相关 pyQt 栅格

!/usr/bin/python3 -- coding: utf-8 -- """ ZetCode PyQt5 tutorial

缺乏、安全感/ 2023年06月15日 04:59/ 0 赞/ 53 阅读

相关 arcpy栅格数据重命名

-- coding: utf-8 -- import arcpy wks = r'C:\Users\HP\Desktop\HB_Q

怼烎@/ 2022年12月11日 02:26/ 0 赞/ 141 阅读

相关 arcpy批量重命名

-- coding: utf-8 -- import arcpy wks = r'C:\Users\HP\Desktop\新建文件

￡神魔★判官ぃ/ 2022年12月06日 01:16/ 0 赞/ 264 阅读

相关批量栅格投影（arcpy）

栅格投影的意义：将分辨率单位为度的变为米定义投影的意义：将某个栅格的坐标系变为想要的坐标系批量栅格投影的代码如下： -- coding: UTF-8 --

深藏阁楼爱情的钟/ 2022年09月17日 04:23/ 0 赞/ 245 阅读

相关栅格数据做相关性分析前的预处理（批量定义投影、栅格投影、重采样）

最近在做两个栅格数据的相关性分析，这个最主要的是两个栅格数的行列号是一致的（也就是说坐标系、分辨率一致），那么标准的栅格数据的分辨率是以米的形式存在的，但是有一些栅格数据是以度

小鱼儿/ 2022年09月17日 04:22/ 0 赞/ 239 阅读

相关批量定义投影(arcpy)

今天突然发现一些影像拖到ArcGIS中和shp文件不重合，实际上是投影不一样的问题，这些问题很常见。那么如何将大量的影像进行定义投影呢？那就用arcpy!哈哈哈！代码如下

Dear 丶/ 2022年09月16日 13:27/ 0 赞/ 200 阅读

相关【Arcpy】批量合并重叠拓扑

需求图形存在大量重叠问题，需要对重叠做合并处理，有特定的对比规则（即判断应合并到哪个要素），人工处理耗时长，对因合并要素产生的图形变化而造成的误差不敏感。思路

「爱情、让人受尽委屈。」/ 2022年09月02日 13:55/ 0 赞/ 300 阅读

相关 PostGIS批量导入栅格数据

如果我们有一批以文件存储的影像数据如何利用PostGIS批量的导出到PostgreSQL数据库中进行管理呢？（单个数据的导入参见我上篇博文：[PostGIS导入导出栅格数

绝地灬酷狼/ 2022年07月12日 06:20/ 0 赞/ 684 阅读

相关栅格数据投影转换

栅格数据投影转换作者：阿振邮箱：tanzhenyugis@163.com 博客：[https://blog.csdn.net/theonegis/article/

本是古典何须时尚/ 2022年05月23日 04:35/ 0 赞/ 633 阅读

相关栅格数据裁剪

文章目录使用GDAL命令使用Python代码版权声明：本文为博主原创文章，转载请注明原文出处！写作时间：2019-03-22

缺乏、安全感/ 2022年03月01日 16:14/ 0 赞/ 343 阅读