arcpy使用RasterToNumPyArray和NumPyArrayToRaster处理栅格像元

最新推荐文章于 2025-02-09 14:08:37 发布

我是一只JU

最新推荐文章于 2025-02-09 14:08:37 发布

阅读量4.2k

点赞数 9

分类专栏： gis 文章标签： arcgis python

本文链接：https://blog.csdn.net/weixin_46668477/article/details/105701675

版权

本文介绍如何利用arcpy的RasterToNumPyArray和NumPyArrayToRaster函数进行栅格数据处理。通过将栅格转化为NumPy数组，可以方便地修改像元值或执行复杂运算，然后将处理后的数组重新转化为栅格。举例说明，代码示例展示了如何去除栅格中值为5和21的像元，实现自定义的数据处理。

摘要生成于 C知道，由 DeepSeek-R1 满血版支持，前往体验 >

使用arcpy提供的函数RasterToNumPyArray将栅格数据转为NumPy数组，得到像元值之后就可以更改各个像元值或者对整个栅格数据执行累积运算等相关功能操作，操作完成之后再使用函数NumpyArrayToRaster将NumPy数组转换为栅格数据并储存，完成整个流程。
比如下面这张图：

在这里插入图片描述
我们要将值为5和21的像元去除，参考以下代码：

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

tif_path = r'F:\All_Forest.tif\test.tif'  # 栅格文件路径

最低0.47元/天解锁文章

确定要放弃本次机会？

福利倒计时

: :

立减 ¥

普通VIP年卡可用

立即使用

我是一只JU

关注关注

9
点赞
踩
32

收藏

觉得还不错? 一键收藏
6
评论
分享

复制链接

分享到 QQ

分享到新浪微博

扫一扫
举报

举报

专栏目录

python numpy矩阵存储数据报 MemoryError

fu6543210的博客

06-07

5334

python调用np存储数据。突然间啥提示也没有，就MemoryError 就停了。网上有人说 python中的numpy 所创建的矩阵大小是有限的，无法创建上万行或者是列的矩阵，如下图错误想想我的代码，确实类似，我一次性存储4万图片，可能有点多了，于是改小一点，果然就不报这个错了。参考别人的建议方法：在用Python处理大数据时，本来16G的内存，内存还没使用四分之一...

arcpy的使用

WhatUwannadO的博客

06-01

727

使用矢量裁剪栅格，并统计栅格的和（排除noda与小于0的值）

6 条评论您还未登录，请先登录后发表或查看评论

关于栅格转数组的计算

John的博客

08-10

790

arcpy提供了一个函数用于栅格转数组，但是如果影像分辨率过高、数据量过大，经常出现报错而没法进行，因此写了一个分块读取的函数并合并到一起，供参考。这个函数是针对单波段栅格数据的。 """Created on Thu May 13 10:23:32 2021 @author: John_Huang @version:python 2.7""" import arcpy from arcpy import env from arcpy.sa import * import numpy def ra

arcpy批量重分类，分区统计

最新发布

2301_80740374的博客

02-09

691

output_raster = os.path.join(output_path, f'风险等级_{year}.tif')pop_raster = os.path.join(pop_path, f'学术严选_china_{year}.tif')print(f"\n处理第 {index}/{total_files} 个文件 ({progress:.1f}%)")print(f"\n处理第 {index}/{total_files} 个年份: {year}")

ArcPy常用函数梳理

吃好喝好，好好学习，好好工作

01-27

3561

1 地理处理环境及其设置 1) 获取和设置地理处理环境地理处理环境设置以env类的属性的方式访问。这些属性可用于检索或设置当前值。环境可作为环境类中的读/写属性进行访问，方法为arcpy.env.<环境名称>。示例1为设置地理处理环境。示例1：设置地理处理环境值 import arcpy # 设置地理处理环境值 acpy.env.worksapce = "c:/St_Johns/data.gdb" 2) 重新设置环境由于地理处理环境对工具操作和

【ArcPy】栅格与Numpy互相转换

Z_W_H_的博客

03-15

341

【代码】【ArcPy】栅格与Numpy互相转换。

10 栅格的使用（3）

weixin_43416590的博客

05-29

1333

10.9 使用 arcpy.sa 和 arcpy.ia 模块的类 arcpy.sa 和 arcpy.ia 模块包含几个主要用于定义栅格工具参数的类。通常，这些类用作工具参数的快捷方式，否则需要更复杂的字符串值。考虑重分类工具的示例。使用此工具，栅格像元会根据重映射表（也称为“重分类表”）赋予新值。重映射表显示每个旧值如何映射到新值。图中的工具对话框显示了土地利用栅格作为适宜性模型的一部分被重新分类为几个新值的示例。在此示例中，输入值由描述土地利用类型的字符串组成，但这些值也可以是数字。重分类工具的语

如何使用Python与ArcGIS库进行二进制数据的处理和编码操作？

12-18

在Python中，我们可以利用ArcGIS库（arcpy和geopandas等模块）来处理地理空间的二进制数据，比如栅格数据。以下是一个基本步骤： 1. 安装必要的库：首先确保你已经安装了`arcgis`, `pandas`, 和 `numpy` 这些库，...

计算多个栅格图层的sum

05-19

3. 在sum运算中，对于每个像元，将所有栅格图层的值相加。 4. 将结果保存到输出栅格图层中。下面是一个使用Python和numpy库计算多个栅格图层的sum的示例代码： ```python import numpy as np import arcpy # ...

基于宁波DEM高程和宁波人口栅格数据利用八邻域种子蔓延洪水淹没算法的python代码如何计算宁波各区县淹没面积和受灾人口

06-13

# 遍历淹没DEM和人口栅格数据的每个像元 for y in range(pop.shape[0]): for x in range(pop.shape[1]): # 如果该像元在该区县内，并且淹没了，则将该像元的受灾人口数量加入popAffected数组 if mask[y, x] != ...

在 ArcGIS 中使用 NumPy.docx

12-07

ArcPy 提供了两个函数，`RasterToNumPyArray` 和 `NumPyArrayToRaster`，它们可以将 ArcGIS 栅格转换为 NumPy 数组，并且也能将 NumPy 数组转换回栅格格式。这种转换特别适用于执行那些可以应用于 NumPy 数组上的...

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

luobengge的博客

08-22

1427

虽然在之前说的在FME中，处理数组我经常用pythonCaller，但是我真不会python，这也并不是我的学习重心，算是买一赠一的赠品，只不过这赠品有的使用还挺好用。最近尝试着用arcpy写了一个小脚本，一路踩了各种坑，简直是坑的我死去活来，记录一下，做个备忘。（里面的代码做了删减，运行不起来的，需要添加你自己的目标坐标系定义，源坐标系因为我的数据已经有了，也没有进行设置。）这算是我写的第一个比较完整的包含了路径处理的小脚本了。想要实现数据的批量处理，其实arcpy的难点并不多。主要就两个方法：

Arcpy进行地图代数栅格运算

I have a dream, here my dreams come true!

10-07

1万+

Python利用Arcpy进行栅格运算ArcGIS栅格计算器介绍Arcpy对栅格文件的读取、创建和写入Arcpy进行栅格的数学分析和逻辑运算归纳总结环境;ArcGIS、Python 2.7 ArcGIS栅格计算器介绍 ArcGIS中的工具箱有栅格计算器工具，它提供的时地图代数功能，可以实现很多复杂的栅格代数运算处理，但栅格计算器工具专门用于 ArcGIS for Desktop 应用程序（仅作为 GP 工具对话框）或模型构建器。它不适用于脚本的编写，而且也不能用于 ArcPy Spatial Analys

Arcpy神器之NumPyArrayToFeatureClass

Anywalker

12-06

1695

Arcpy神器之NumPyArrayToFeatureClass 之前在“ofo数据获取&amp;坐标转换&amp;生成shp数据”一文中涉及到用Arcpy创建点要素数据，最近再次使用时发现怎么那么麻烦于是有另辟新路，翻阅Arcpy接口文档，终于找到了一个神器NumPyArrayToFeatureClass。我们先看看原来的写法 #创建shp点数据 def CreateFeaturc...

关于arcpy.da.FeatureClassToNumPyArray的一些应用

白氏可乐python学习归纳区

03-14

2450

import arcpy import numpy as np fc=ur"E:\test.shp" array=arcpy.da.FeatureClassToNumPyArray(fc,("Name_CHN","Price")) print array[array["Name_CHN"]==u"北京市"]["Price"].ptp() print array[array["Name_CHN"]

[ArcPy] 遥感影像去黑边-第六届全国大学生GIS技能大赛试题

geodoer

12-21

5535

去除遥感影像黑边： 1. 将边界进行矢量化，然后用淹膜提取 2. 边界像元值为0，将0值重分类为NoData 3. 法一法二联合做就可以解决他们两个的缺点 4. 使用ArcPy仅将边缘像元为0的地方赋予NoData，里面的地方不处理 5. 使用边界跟踪算法得出边界，存储成shp格式，再进行淹膜

[空间分析] 低分辨率转高分辨率-DEM与栅格图像配准-空间内插

geodoer

10-10

8636

【问题背景】两幅不同的分辨率的DEM（DEM1，DEM2），将分辨率统一【思路】低分辨率的DEM1 --&gt;转成点要素 --&gt; 插值（按照DEM2的范围、DEM2的分辨率）–&gt; 得到产品raster 【原理】一幅栅格图=矩阵+栅格大小（分辨率）+栅格的范围（左下角的经纬度） --&gt; 栅格大小、栅格范围一定，栅格的框框位置也一定，也对的上【判断分辨率是否转成一样】肉眼看...

Python学习：numpy库数据量太大出现Memory Error问题的解决方法汇总