ArcGIS应用(十九)Arcgis 统计分析计算多波段图像最大值、最小值、平均值等

1.计算目的

计算多波段图像中,每个象元对应位置的所有波段值的最大值,最小值平均值等参数信息。

使用场景,比如计算时间序列遥感指数的最大值。假如有三个时间的NDVI指数,我们需要计算三个时间的最大值图像,就需要这样计算每个象元位置三个时相的最大值。

2.实现方法

2.1 使用工具

使用Arcgis统计分析工具计算得到,打开工具箱:

打开Spatial Analyst Tools中的zonal工具:

如下图:

 使用ZonalStatistics工具进行统计分析.点击该工具,界面如下:

 参数说明(原版英文说明):

Input raster or feature zone data

Dataset that defines the zones. 这个参数用于定义统计区域,我们通常可以选择行政区域进行统计,这里我们选择需要统计的原图像,即统计原图像范围内极值等信息。

The zones can be defined by an integer raster or a feature layer.

Zone field

Field that holds the values that define each zone.

It can be an integer or a string field of the zone dataset.

Input value raster

Raster that contains the values on which to calculate a statistic. 

Output raster

The output zonal statistics raster.

Statistics type (optional)

Statistic type to be calculated.

  • MEAN — Calculates the average of all cells in the value raster that belong to the same zone as the output cell.
  • MAJORITY — Determines the value that occurs most often of all cells in the value raster that belong to the same zone as the output cell.
  • MAXIMUM — Determines the largest value of all cells in the value raster that belong to the same zone as the output cell.
  • MEDIAN — Determines the median value of all cells in the value raster that belong to the same zone as the output cell.
  • MINIMUM — Determines the smallest value of all cells in the value raster that belong to the same zone as the output cell.
  • MINORITY — Determines the value that occurs least often of all cells in the value raster that belong to the same zone as the output cell.
  • RANGE — Calculates the difference between the largest and smallest value of all cells in the value raster that belong to the same zone as the output cell.
  • STD — Calculates the standard deviation of all cells in the value raster that belong to the same zone as the output cell.
  • SUM — Calculates the total value of all cells in the value raster that belong to the same zone as the output cell.
  • VARIETY — Calculates the number of unique values for all cells in the value raster that belong to the same zone as the output cell.

Ignore NoData in calculations (optional)

Denotes whether NoData values in the Value input will influence the results of the zone that they fall within.

  • Checked—Within any particular zone, only cells that have a value in the input Value raster will be used in determining the output value for that zone. NoData cells in the Value raster will be ignored in the statistic calculation.
  • Unchecked—Within any particular zone, if any NoData cells exist in the Value raster it is deemed that there is insufficient information to perform statistical calculations for all the cells in that zone; therefore, the entire zone will receive the NoData value on the output raster.

2.2 参数设置与计算

打开图像数据, 在工具面板相关选项选择输入参数信息:

Zone Field这里默认Value会报错,目前还没找到具体原因.

 报错:

ERROR 010123: Could not get field name.
Failed to execute (ZonalStatistics).
Failed at Sun May 29 22:50:41 2022 (Elapsed Time: 0.06 seconds)

 将改选项改为count重新计算,可以得到计算结果(原图和结果图),结果图为单波段图像:

 

 3.结果分析

使用Identity工具查看象元属性发现,所计算结果只计算得到了前三个波段的最大值,而不是所有波段的最大值。这是这个工具的局限所在,可能和软件版本有关,以前版本针对raster数据计算也是只能计算前三个波段。后面将进一步分析原因。

实习内容: 练习1:使用缺省参数创建一个表面 练习2:数据检查 练习3:制作臭氧浓度图 练习4:模型比较 练习5:制作超出某一临界值的臭氧概率图 练习6:生成最终成果图 预备知识: 利用地统计分析模块,你可以根据一个点要素层中已测定 采样点、栅格层或者利用多边形质心,轻而易举地生成一个连续表面。这些采样点的值可以是海拔高度、地下水位的深度或者污染值浓度等。当与ArcMap一起使用时,地统计分析模块提供了一整套创建表面的工具,这些表面能够用来可视化、分析及理解各种空间现象。 美国环保局负责对加利佛尼亚州的大气臭氧浓度进行监测。臭氧浓度值是通过遍布全州的监测站来测定的。已经知道所有监测站的臭氧值浓度,但是我们还想知道加利佛尼亚州其他任一地方的臭氧浓度值,但是考虑到费用以及实用性问题,我们不可能在任何地方都建立监测站。地统计分析模块提供了许多工具,通过检测所有采样点之间的关系,生成一个关于臭氧浓度值、预测标差(不确定性)以及超出临界值的概率的连续表面,从而使对其他点的浓度值进行最佳预测成为可能。 在地统计分析模块的安装盘上,附有本次联系所需要的数据。这些数据包括: 数据集 描述 Ca_outline 加州轮廓图 Ca_ozone_pts 臭氧采样点数据(单位:ppm) Ca_cities 加州主要城市位置图 Ca_hilshade 加州山体阴影图 臭氧数据集(Ca_ozone_pts)表示的是1996年中每八个小时时段内的臭氧平均浓度的最大值,以ppm为单位。根据这些臭氧采样点的测量值,你可以生成两个连续表面(或地图),以现有的采样点数据为基础来预测加洲任意地方的臭氧浓度值。创建第一张图时你只需要简单地使用缺省选项即可,你会看到根据采样点数据生成表面是多么容易的事情。生成第二张图时,你可以较多的考虑采样点之间存在的空间关系,这时,你将要用到ESDA(空间数据探索分析)工具来检查你的数据。你还会学习到一些地统计选项,利用这些选项,你可以创建诸如剔除趋势并且模拟空间自相关的表面。利用ESDA工具,通过对各种地统计参数的操作,你可以创建出更加精彩的表面。 很多时候,我们关心的并不是那些有损于健康的风险的实际值,而是它是否超出某一有毒水平,如果确实超出,必须采取行动。你创建的第三个表面将用来评估臭氧浓度超过临界值的概率。 对于本次练习,临界值规定如下:如果最大值超过0.12ppm,那么这个地方就应该被严密监测。利用这个标准,你可以利用地统计分析模块来预测臭氧浓度超出临界值的概率大小。本练习被分为各个单独的任务,你可以根据自己的学习进度来逐步熟悉地统计分析模块个各项功能。 练习1带你进入地统计分析模块,通过创建臭氧浓度表面的过程,你会发现使用缺省参数创建表面是一件很容易的事情。 练习2教你创建表面之前如何对数据进行检查。数据检查的目的是为了找出数据中那些离群值并且发现数据中存在的趋势。 练习3你将创建第二个表面,这个表面更多地考虑了练习2中数据分析发现的空间关系,并且对练习1中生成的表面进行了改进。该练习还将向你介绍地统计学中的一些基本概念。 练习4教你如何对练习1和练习3中创建的表面进行比较,并判断哪个表面对未知值的预测更好。 练习5指导你创建臭氧浓度超出临界值的概率图,从而生成第三个表面。 练习6告诉你怎样利用ArcMap的功能将你在练习3和练习5中创建的表面放在一起做最终的显示。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

空中旋转篮球

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值