怎样用arcgis中的python进行批量投影

发布网友发布时间：2022-04-24 09:24

共2个回答

懂视网时间：2022-04-18 23:11

出图是项目里常见的任务，有的项目甚至会要上百张图片，所以批量出土工具很有必要。arcpy.mapping就是ArcGIS里的出图模块，能快速完成一个出图工具。

arcpy.mapping模块里常用的类有MapDocument、DataFrame、Layer、DataDrivenPages和TextElement。

MapDocument类是地图文档（.mxd文件）对应的类。初始化参数是一个字符串，一般是.mxd文件的路径：

 mxd=arcpy.mapping.MapDocument(r"F:GeoDataChinaAreaChinaVector.mxd")

DataFrame类用于操作地图内的Data Frame（即下图的Layers），能够控制地图的范围、比例尺等。用arcpy.mapping.ListDataFrames(map_document, {wildcard})函数获取。

df= arcpy.mapping.ListDataFrames(mxd)[0]

Layer类用于操作具体的图层。能够控制图斑的样式、可见性等。可以用.lyr文件的路径初始化，也可以通过arcpy.mapping.ListLayers(map_document_or_layer, {wildcard}, {data_frame})函数获取。

lyr1=arcpy.mapping.Layer(r" F:GeoDataChinaAreaProvince.lyr")

df.addLayer(lyr1)

lyr2=arcpy.mapping.ListLayer(mxd,"",df)[0]

DataDrivenPages类需要配合ArcMap中的Data Driven Pages工具使用。用于一个矢量文件内的全部或部分图斑每个出一张图的情况。

TextElement类用于操作地图上的文字，比如图名、页数。通过arcpy.mapping.ListLayoutElements (map_document, {element_type}, {wildcard})函数获取。

txtElm=arcpy.mapping.ListLayoutElements(mxd,"TEXT_ELEMENT")[0]

常见的出图模式有两种：一个矢量文件里每个图斑出一张图，一个文件夹下每个矢量文件出一张图。

每个图斑出一张图：

这种情况有Data Driven Pages工具配合最好。打开ArcMap的Customize->Toolbars->Data Driven Pages，设置好图层、名称字段、排序字段、显示范围和比例尺，保存地图。

# coding:utf-8

import arcpy

 

mxd=arcpy.mapping.MapDocument(r"F:GeoDataChinaAreaChinaVector.mxd")

for pageNum in range(1,mxd.dataDrivenPages.pageCount):

 mxd.dataDrivenPages.currentPageID=pageNum

 mapName=mxd.dataDrivenPages.pageRow.getValue(mxd.dataDrivenPages.pageNameField.name)

 print mapName

 arcpy.mapping.ExportToPNG(mxd,r"F:GeoDataChinaAreaProvince"+mapName+".png")

print 'ok'

一个文件夹下的每个矢量文件出一张图：

# coding:utf-8

import arcpy

import os

 

def GetShpfiles(shpdir):

 shpfiles=[]

 allfiles=os.listdir(shpdir)

 for file in allfiles:

  if os.path.isfile(file):

  if file.endswith('.shp'):

   shpfiles.append(file)

  else:

  shpfiles.extend(GetShpfiles(file))

 return shpfiles

 

allshps=GetShpfiles(r"F:GeoDataChinaAreaProvince")

mxd=arcpy.mapping.MapDocument(r"F:GeoDataChinaAreaChinaVector.mxd")

lyr=arcpy.mapping.ListLayer(mxd)[0]

for shp in allshps:

 paths=os.path.split(shp)

 print paths[1]

 lyr.replaceDataSource(paths[0],"SHAPEFILE_WORKSPACE",paths[1])

 arcpy.mapping.ExportToPNG(mxd,r"F:GeoDataChinaAreaProvince"+paths[1]+".png")

print 'ok'

热心网友时间：2022-04-18 20:19

接到的需求是对文件夹下包括子目录所有影像进行投影变换，在网上找了找，大家的代码好像和需求不太一致，自己写了一个小工具。

思路是遍历文件夹所有.tif文件。

1.用os.walk()函数来获取路径和文件名

2.用 if 语句*读取后缀为.tif的文件

3.将文件路径和文件名组合起来

4.调用arcpy.ProjectRaster_management语句投影变换并输出

由于刚刚接触python，难点主要在于os.walk()目录的读取这个函数的了解上。

以下是代码

#coding = utf-8
import os
import arcpy
from arcpy import *

#定义一个函数，第一个参数是一个集合，用来接收所有.tif的路径。第二个函数是指定的文件夹。
def huoqu(jieshou,lujing):
    for dirpath,dirnames,filenames in os.walk(lujing):
        for filename in filenames:
            if  os.path.splitext(filename)[1] == '.tif':
                file1 = dirpath+'\\'+filename
                jieshou.append(file1)

#这里变量的设置，GetParameterAsText函数来获取对应在python工具箱中的位置
env.workspace = GetParameterAsText(0)
work = env.workspace
ref = GetParameterAsText(1)
tifz = []

#调用自定义的函数
huoqu(tifz,work);

#对所有的影像进行循环投影变换，并且将输出路径设置为指定的文件夹+原本的名字
for file in tifz:
    out = GetParameterAsText(2)+"/"+ os.path.basename(file)
    ProjectRaster_management(file,out,ref)

全部栏目

怎样用arcgis中的python进行批量投影