• ArcGis Python脚本——批量对影像、要素类定义投影


    插个广告,制作ArcGIS的Tool工具学习下面的教程就对了:
    零基础学习Python制作ArcGIS自定义工具观看链接
    《零基础学习Python制作ArcGIS自定义工具》课程简介

    这一段是批量定义要素类(FeatureClasses)投影的ArcPy代码:

    把要处理的要素类塞进一个文件夹(工作空间workspace),然后将代码开头的路径换成这个“文件夹”的路径,处理完后再做其他处理(比如批量塞回原来文件夹,vba,python啥的代码、小工具处理即可)。

    define coordinate system那里把坐标系信息改为自己需要定义的坐标系信息。不知道哪儿获取?先定义好一个Shpfile的投影信息,然后把它的.prj文件用记事本打开。

    import arcpy
    # set environment workspace
    arcpy.env.workspace=r"C:UsersAdministratorDesktopArcCode新建文件夹"
    # get files
    clses = arcpy.ListFeatureClasses()  
    # define coordinate system
    Coordinate_System = "PROJCS['CGCS2000_3_Degree_GK_CM_114E',GEOGCS['GCS_China_Geodetic_Coordinate_System_2000',DATUM['D_China_2000',SPHEROID['CGCS2000',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Gauss_Kruger'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',114.0],PARAMETER['Scale_Factor',1.0],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]"
    # code source: https://www.cnblogs.com/yzhyingcool/
    # QQ:975601416
    for cls in clses: print str(cls) arcpy.DefineProjection_management(cls, Coordinate_System) print("OK")

    这是一段批量定义影像投影的Arcpy代码:

    影像图片格式:代码中定义为tif格式,其他格式自行更改。

    基本设置如上,不复述。

    import arcpy
    # set environment workspace
    arcpy.env.workspace=r"C:UsersAdministratorDesktop新建文件夹"
    # get files
    rasters = arcpy.ListRasters("*", "tif")   
    # define coordinate system
    Coordinate_System = "PROJCS['CGCS2000_3_Degree_GK_CM_114E',GEOGCS['GCS_China_Geodetic_Coordinate_System_2000',DATUM['D_China_2000',SPHEROID['CGCS2000',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Gauss_Kruger'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',114.0],PARAMETER['Scale_Factor',1.0],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]"
    #code source: https://www.cnblogs.com/yzhyingcool/ QQ:975601416
    for ras in rasters:
        print str(ras)
        arcpy.DefineProjection_management(ras, Coordinate_System)
    print("OK")
  • 相关阅读:
    with ,Row_Number,DateDiff,DateAdd用法学习
    jmeter 读取mysql数据库
    fidder 自动保存请求内容
    redis 常用方法整理
    解决:EXCEL复制粘贴,精度丢失
    MYSQL 创建常见问题
    MYSQL 存储过程、函数、临时表、游标
    MYSQL 测试常用语句使用技巧
    3-6
    selenium3 下载、配置
  • 原文地址:https://www.cnblogs.com/yzhyingcool/p/10274153.html
Copyright © 2020-2023  润新知