Python利用arcpy模块实现栅格的创建与拼接

与天争命

本文介绍基于Python语言

1. arcpy

复制代码

模块，实现栅格影像图层建立与多幅遥感影像数据批量拼接（Mosaic）的操作。
首先，相关操作所需具体代码如下：

1. import os
   
2. import arcpy
   
3. 
   
4. file_path="G:/Postgraduate/LAI_Glass_RTlab/A2018161_Dif/DRT/"
   
5. out_file_path="G:/Postgraduate/LAI_Glass_RTlab/A2018161_Dif/DRT/"
   
6. out_file_name="Global.tif"
   
7. 
   
8. file_name_list=os.listdir(file_path)
   
9. 
   
10. tif_file_path=file_path+file_name_list[0]
    
11. cell_size_x=arcpy.GetRasterProperties_management(tif_file_path,"CELLSIZEX")
    
12. cell_size=cell_size_x.getOutput(0)
    
13. value_type=arcpy.GetRasterProperties_management(tif_file_path,"VALUETYPE")
    
14. describe=arcpy.Describe(tif_file_path)
    
15. spatial_reference=describe.spatialReference
    
16. 
    
17. arcpy.CreateRasterDataset_management(out_file_path,out_file_name,cell_size,"16_BIT_SIGNED",
    
18.                                      spatial_reference,"1")
    
19. 
    
20. out_file=out_file_path+out_file_name
    
21. for file in file_name_list:
    
22.     file_path_name=file_path+file
    
23.     print(file_path_name)
    
24.     arcpy.Mosaic_management([file_path_name],out_file)

复制代码

其中，

1. file_path

复制代码

为存放有多景初始遥感影像的路径格式为

1. .tif

复制代码

栅格文件（如果不是

1. .tif

复制代码

格式，例如是

1. .hdf

复制代码

等文件，需首先进行文件格式的转换）；

1. out_file_path

复制代码

为拼接后所得结果栅格图层的存放路径；

1. out_file_name

复制代码

为拼接后所得结果栅格图层的文件名称，其可选格式有很多，如下图所示。

在这里，我们默认所得拼接结果图层为一个（也就是

1. file_path

复制代码

文件夹中全部的待处理遥感影像最终全拼接在一起）；如果大家需要使得拼接结果图层是多幅（也就是

1. file_path

复制代码

文件夹中待处理遥感影像依据区域、时间等分为很多不同的部分，每一部分拼接在一起），可以参考Python GDAL读取栅格数据并基于质量评估波段QA对指定数据加以筛选掩膜，利用其中的循环方式实现需求。
随后，通过

1. os.listdir()

复制代码

函数获取

1. file_path

复制代码

路径下的栅格文件，并存储于

1. file_name_list

复制代码

列表中。
接下来需要创建一个新的栅格图层。之所以要进行这一步骤，是因为本文后期选择用

1. arcpy.Mosaic_management()

复制代码

函数进行栅格的批量拼接，因此需要首先创建一个新的、空的栅格图层作为拼接的基准。如果大家的需求不是批量拼接栅格数据，而是单纯想利用

1. arcpy

复制代码

进行新栅格的创建，那就只看这一部分的代码即可。
在这里，我们选择用

1. file_path

复制代码

路径下的第一个栅格数据（下称“第一栅格”）作为新栅格图层中各项属性（例如像素边长、像素数据格式等）的依据。首先，

1. arcpy.GetRasterProperties_management()

复制代码

函数获取第一栅格的像素

1. x

复制代码

边边长；因为一般栅格数据中像素都是正方形，因此我们就通过

1. cell_size=cell_size_x.getOutput(0)

复制代码

将第一栅格的像素

1. x

复制代码

边边长作为新栅格图层像素

1. x

复制代码

边与

1. y

复制代码

边二者的边长。再利用

1. arcpy.GetRasterProperties_management()

复制代码

函数获取第一栅格的数据格式；最后利用中间变量

1. describe

复制代码

获取第一栅格的空间参考信息。
完成以上步骤后，将已获取的第一栅格的各类信息通过函数

1. arcpy.CreateRasterDataset_management()

复制代码

带入新栅格中。在这里需要注意：尽可能在将要拼接时选择新栅格为

1. "16_BIT_SIGNED"

复制代码

及以下的数据格式（具体数据格式类别如下图），且将

1. file_path

复制代码

路径下待拼接的栅格数据的数据格式也全部修改为这一格式；否则可能会由于数据量大而导致拼接过程极慢。我之前就是由于选用了

1. 32 bit float

复制代码

格式的栅格数据进行拼接，导致全球范围的MODIS一个植被产品数据拼接花了将近一天的时间。如果大家的栅格像素数据包含小数，可以通过乘上一个缩放系数的方式进行数据整数化。

代码最后的一个

1. for

复制代码

循环，就是遍历

1. file_name_list

复制代码

中的各个栅格数据，并通过

1. arcpy.Mosaic_management()

复制代码

函数加以拼接即可。
以上，便完成了本次批量拼接的操作。这里还有一点需要注意：由于

1. arcpy

复制代码

模块的限制，如果大家的Python版本是

1. 3.0

复制代码

及以上，往往不能直接运行上述代码，最好是在ArcMap的Python运行框或其对应IDLE（如下图所示）中运行。

至此，大功告成。
到此这篇关于Python利用arcpy模块实现栅格的创建与拼接的文章就介绍到这了,更多相关Python arcpy栅格内容请搜索脚本之家以前的文章或继续浏览下面的相关文章希望大家以后多多支持脚本之家！

来源:https://www.jb51.net/article/276370.htm
免责声明：由于采集信息均来自互联网，如果侵犯了您的权益，请联系我们【E-Mail:cb@itdo.tech】 我们会及时删除侵权内容，谢谢合作！