arcpy对栅格数据进行坐标转换——arcpy在处理路径中踩得坑

虽然在之前说的在FME中,处理数组我经常用pythonCaller,但是我真不会python,这也并不是我的学习重心,算是买一赠一的赠品,只不过这赠品有的使用还挺好用。
最近尝试着用arcpy写了一个小脚本,一路踩了各种坑,简直是坑的我死去活来,记录一下,做个备忘。
(里面的代码做了删减,运行不起来的,需要添加你自己的目标坐标系定义,源坐标系因为我的数据已经有了,也没有进行设置。)

这算是我写的第一个比较完整的包含了路径处理的小脚本了。
想要实现数据的批量处理,其实arcpy的难点并不多。主要就两个方法:
arcpy.ListRasters()和arcpy.ProjectRaster_management()
关键在于是中文格式编码的处理,和路径的拼接。

1.引入文件路径
在处理中文的时候首先需要声明
#--coding=utf-8--
其次,由于路径中存在中文,对于引入的中文路径,需要做decode
这步是极其关键的一步,我就是处理好中文路径,坑了自己好长时间。

#-*-coding=utf-8-*-
import arcpy,os,os.path
from arcpy import env
input_path = r'I:\源坐文件'.decode("utf-8")
output_path = r'I:\转换完成后'.decode("utf-8")
Coordinate_System = '这里写你自己的目标坐标系'

2.用os.walk()方法,遍历文件夹,获取所有的路径及栅格文件,并以此为基础进行坐标转换

#遍历文件夹,获取所有的栅格数据
for parent, dirnames, filenames in os.walk(input_path):#三个参数:1.父目录 2.所有文件夹名字(不含路径)3.所有文件名字
    for dirname in  dirnames:#输出文件夹信息
        env.workspace = os.path.join(parent,dirname)
        rasters = arcpy.ListRasters()
        for raster in rasters:
            output_dir = os.path.join(output_path,dirname)
            out_raster = os.path.join(output_dir,raster)
            #print "output_dir is" + output_dir
            print "output_dir is" + output_dir
            isExists=os.path.exists(output_dir)
            # 判断结果
            if not isExists:
                # 如果不存在则创建目录
                os.makedirs(output_dir)
            print 'start to project '+output_dir+"\\"+raster               
            try:
            	arcpy.ProjectRaster_management(raster,out_raster,Coordinate_System,"NEAREST","#","#","#","#")
            except Exception as e:
                print e
            else:
                print output_dir+" has projected "+raster

没了,没什么好写的,我只是记录一下自己踩的坑。

  • 3
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 4
    评论
### 回答1: arcpy是一个强大的Python库,可以用于处理GIS数据。如果要遍历栅格数据并将数据复制到同名文件夹内,可以使用以下代码: ```python import arcpy import os # 设置工作目录 arcpy.env.workspace = r"path/to/workspace" # 获取所有的栅格数据 raster_list = arcpy.ListRasters() # 遍历所有栅格数据 for raster in raster_list: # 获取栅格数据路径和名称 raster_path = os.path.dirname(raster) raster_name = os.path.basename(raster) # 创建同名文件夹 output_folder = os.path.join(raster_path, raster_name) if not os.path.exists(output_folder): os.makedirs(output_folder) # 复制栅格数据到同名文件夹内 arcpy.CopyRaster_management(raster, os.path.join(output_folder, raster_name)) ``` 这段代码首先设置工作目录,然后获取所有的栅格数据。接下来,遍历所有的栅格数据,获取栅格数据路径和名称,并创建同名文件夹。最后,使用`arcpy.CopyRaster_management`将栅格数据复制到同名文件夹内。 ### 回答2: arcpy是一个用于处理地理空间数据的Python模块。要遍历栅格数据并将数据复制到同名文件夹内,可以按照以下步骤进行操作: 1. 首先,需要导入arcpy模块: import arcpy 2. 然后,使用arcpy.ListRasters函数获取当前工作空间的所有栅格数据: rasters = arcpy.ListRasters() 3. 创建一个循环遍历所有栅格数据: for raster in rasters: # 获取当前栅格数据的名称,并去掉文件扩展名 raster_name = arcpy.Describe(raster).baseName # 创建同名文件夹 folder_path = arcpy.env.workspace + "\\" + raster_name if not arcpy.Exists(folder_path): arcpy.CreateFolder_management(arcpy.env.workspace, raster_name) # 复制栅格数据到同名文件夹内 arcpy.CopyRaster_management(raster, folder_path + "\\" + raster_name) 以上代码arcpy.Describe函数用于获取栅格数据的属性信息,arcpy.env.workspace表示当前工作空间路径arcpy.Exists函数用于判断文件夹是否存在,arcpy.CreateFolder_management函数用于创建文件夹,arcpy.CopyRaster_management函数用于复制栅格数据到指定路径。 这样,就可以通过arcpy遍历栅格数据并将数据复制到同名文件夹内了。 ### 回答3: arcpy是ArcGIS软件的一个Python模块,它可以用于地理空间数据的处理和分析。如果想要遍历栅格数据并将数据复制到同名文件夹内,可以使用arcpy模块来完成这个任务。 首先,需要导入arcpy模块,然后使用ListRasters函数来获取栅格数据的列表。这个函数可以接收一个输入文件夹路径,返回一个包含所有栅格数据的列表。 接下来,遍历栅格数据列表。对于每一个栅格数据,首先获取它的文件名(包括扩展名),然后使用os模块的splitext函数将文件名和扩展名分开。再次使用os模块的join函数将栅格数据的文件名与输入文件夹路径拼接起来,得到同名文件夹的路径。 然后,使用arcpy的CopyRaster函数来复制栅格数据。这个函数需要提供栅格数据的输入路径、输出路径和可选的数据格式。 最后,重复上述步骤直到遍历完所有栅格数据。 下面是一个示例代码: ```python import arcpy import os input_folder = "输入文件夹路径" # 获取栅格数据列表 raster_list = arcpy.ListRasters("*", "ALL", input_folder) # 遍历栅格数据列表 for raster in raster_list: # 获取栅格数据的文件名和扩展名 file_name, file_extension = os.path.splitext(raster) # 构建同名文件夹的路径 output_folder = os.path.join(input_folder, file_name) # 如果同名文件夹不存在,则创建 if not os.path.isdir(output_folder): os.mkdir(output_folder) # 复制栅格数据到同名文件夹 arcpy.CopyRaster_management(os.path.join(input_folder, raster), os.path.join(output_folder, raster), "", "", "", "NONE", "NONE", "") ``` 以上代码会遍历输入文件夹的所有栅格数据,并将每一个栅格数据复制到同名的文件夹内。注意要替换"输入文件夹路径"为实际的输入文件夹路径

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 4
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值