【GIS】arcpy对栅格的常见操作——拼接裁剪代数计算

HappyFinch

已于 2022-03-12 21:16:57 修改

阅读量1.2k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： GIS与python 文章标签： python

于 2022-03-12 21:12:12 首次发布

本文链接：https://blog.youkuaiyun.com/weixin_48337229/article/details/123449504

最近在做一个项目，主要任务是对栅格数据进行各种操作，将利用arcpy进行的各项操作记录于此帖中【持续更新】

下载栅格数据

主要用到MODIS数据，下载链接：MODIS数据下载

注意事项：

MODIS数据的官网下载需要科学上网。在我们国家的一些地理数据平台上下载MODIS数据不需要科学上网，但是上面的数据不是很全。
关于下载格式：① 网站默认的MODIS下载格式为hdf，这个格式可以在ArcGIS与ENVI中打开；② 在提供的链接下载地址中下载MODIS数据时，第五步 REVIEW&ORDER 页面中提供了下载的格式选择，可以勾选该选项，选择下载格式为GeoTiff，之后便可以下载tif格式的栅格卫星图了。该格式好处是arcpy处理tif格式的栅格文件更便捷。

栅格图像拼接

# -*- coding: utf-8 -*-
import arcpy

arcpy.CheckOutExtension("Spatial") 
 
#设置工作空间
arcpy.env.workspace = r'C:\Users\Hqy\Desktop\data'

#files列表里是需要拼接成一张图的所有栅格文件
files=[];
for ras in arcpy.ListRasters("*.tif"):
    files.append(ras)

base = files[0]  #设置一个基准的文件 （此处我选择了第一个栅格）

out_coor_system = arcpy.Describe(base).spatialReference 
dataType = arcpy.Describe(base).DataType  
piexl_type = arcpy.De