GIS数据共享:官方网站

GIS

当前位置:首页 > 3S > GIS

ArcGIS 二次开发系列之Python版示例代码——矢量栅格数据批量投影转化

ArcGIS 二次开发系列之Python版示例代码——矢量栅格数据批量投影转化

说明:本示例代码使用最简单的python语法,简要介绍如何使用Python做ArcGIS二次开发。适合初学者阅读。

# -*- coding: utf-8 -*-
 
import arcpy
import os
 
input_folder = "H:\\Original\\qhnu\\wangxch20170318\\MAP"
fileExtent = [".shp",".tif",".tiff",".img"]
 
for files in os.walk(input_folder ):
    for filename in files[2]:
        for Extent in fileExtent:
            if filename.endswith(Extent):
                inputData = files[0] + "\\" + filename
                oupPutData = files[0].replace("wangxch20170318","modif") + "\\r" + filename
                print inputData,oupPutData
                if os.path.exists(files[0].replace("wangxch20170318","modif")) == False:
                    os.mkdir(files[0].replace("wangxch20170318","modif"))
                if Extent == ".shp":
                    arcpy.Project_management(inputData, oupPutData, "PROJCS['WGS_1984_Albers',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Albers'],PARAMETER['False_Easting',4000000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',96.5],PARAMETER['Standard_Parallel_1',34.0],PARAMETER['Standard_Parallel_2',37.0],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]", "", "PROJCS['WGS_1984_Albers',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Albers'],PARAMETER['False_Easting',4000000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',105.0],PARAMETER['Standard_Parallel_1',25.0],PARAMETER['Standard_Parallel_2',47.0],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]", "NO_PRESERVE_SHAPE", "", "NO_VERTICAL")
                else:
                    arcpy.ProjectRaster_management(inputData, oupPutData, "PROJCS['WGS_1984_Albers',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Albers'],PARAMETER['False_Easting',4000000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',96.5],PARAMETER['Standard_Parallel_1',34.0],PARAMETER['Standard_Parallel_2',37.0],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]", "NEAREST", "1000 1000", "", "", "PROJCS['Albers_Conical_Equal_Area',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Albers'],PARAMETER['false_easting',4000000.0],PARAMETER['false_northing',0.0],PARAMETER['central_meridian',105.0],PARAMETER['standard_parallel_1',25.0],PARAMETER['standard_parallel_2',47.0],PARAMETER['latitude_of_origin',0.0],UNIT['Meter',1.0]]")

在上面的代码中:


    input_folder定义了输入文件的存放路径,读者可按自己实际情况进行修改。


    fileExtent对input_folder文件夹中的数据进行筛选。


    下面的三个for循环是从input_folder文件夹中查找指定的文件,如果读者对这段代码不熟悉,需要学习一下python基础知识。 


    inputData定义了从input_folder文件夹找到的文件的完整路径,oupPutData定义了做投影转换后文件的保存位置(因为我的输入文件夹中有“wangxch20170318”文件夹,所以使用python的replace方法将输入路径中的“angxch20170318”文件夹替换为“modify”文件夹,读者可按实际情况进行改写)。


    print语句下面的if语句是判断输出文件夹是否存在,如果不存在则创建(这里有两点说明:1. 因为前面讲输入文件的“wangxch20170318”文件夹修改为“modify”文件夹,在python中只是对字符串进行编辑,实际上这个文件夹不存在,所以需要新建文件夹。2. 如果读者新建的目录较深,os的mkdir会报错,可以使用os的makedirs方法)。


    ArcMap对矢量和栅格的投影转换方式在不同的toolbox,所以此处我们也需要对矢量和栅格使用不同的函数(具体函数的参数读者请参考ArcMap自带的help,找到对应的工具说明)。

扫码查看

上一篇:已经是第一篇

下一篇:ArcMap各个版本下载

相关内容

文章评论

表情

共 0 条评论,查看全部
  • 这篇文章还没有收到评论,赶紧来抢沙发吧~

热门标签