本文介绍基于Python中ArcPy模块，对大量长时间序列栅格遥感影像文件的每一个像元进行多时序平均值的求取。
在遥感应用中，我们经常需要对某一景遥感影像中的全部像元的像素值进行平均值求取——这一操作很好实现，基于ArcMap软件或者简单的Python代码就可以实现；但有时候，我们会需要结合同一地区、不同时相的多景遥感影像，求取每一个像元在全部时相中像素值的平均值——这一需求的实现较之前者就有些麻烦，本文对此加以介绍。
首先，我们来明确一下本文的具体需求。现有一个存储有大量.tif格式遥感影像的文件夹，其中每一个遥感影像的文件名中都包含有该图像的成像时间，如下图所示。且其中除了.tif格式的遥感影像文件外，还具有其它格式的文件。
我们希望，对于同一年成像的遥感影像进行逐像元平均值的求取。例如，上图中具有2001年第185天成像、第193天成像、第201天成像……等等遥感影像8幅，每一幅都是这一年不同时间在同一空间位置的成像；同时，还有2005年不同时间成像的遥感影像9幅。我们希望，首先将2001年成像的8幅遥感影像加以逐像元平均值的求取，即求取每一个像元在这8景图像中像素值的平均；随后再对2005年成像的9幅遥感影像加以逐像元平均值的求取，以此类推。
明确了需求后，我们就可以开始具体的操作。首先，本文所需用到的代码如下。
# -*- coding: utf-8 -*-
"""
Created on Sat Apr 16 10:48:37 2022
@author: fkxxgis
"""
import arcpy
from arcpy.sa import *
tif_file_path="E:/LST/Data/MODIS/05_Resample/"
average_file_path="E:/LST/Data/MODIS/06_Average/"
arcpy.env.workspace=tif_file_path
tif_file_name=arcpy.ListRasters("*","tif")
tif_file_year=tif_file_name[0][0:4]
one_year_tif_list=[]
sum_pic=0
for tif_file in tif_file_name:
if tif_file[0:4]==tif_file_year:
one_year_tif_list.append(tif_file)
tif_file_temp=tif_file
if tif_file==tif_file_name[len(tif_file_name)-1]:
pic_num=len(one_year_tif_list)
for tif_file_new in one_year_tif_list:
sum_pic=sum_pic+Raster(tif_file_new)
(sum_pic/pic_num).save(average_file_path+tif_file_year+"_Ave.tif")
else:
pic_num=len(one_year_tif_list)
for tif_file_new in one_year_tif_list:
sum_pic=sum_pic+Raster(tif_file_new)
(sum_pic/pic_num).save(average_file_path+tif_file_year+"_Ave.tif")
one_year_tif_list=[]
sum_pic=0
one_year_tif_list.append(tif_file)
tif_file_year=tif_file[0:4]
其中，tif_file_path是原有计算平均值前遥感图像的保存路径，average_file_path是我们新生成的求取平均值后遥感影像的保存路径，也就是结果保存路径。
在这里，和我们前期的博客Python ArcPy批量拼接长时间序列栅格图像类似，需要首先在资源管理器中，将tif_file_path路径下的各文件以“名称”排序的方式进行排序；随后，利用arcpy.ListRasters()函数，获取路径下原有的全部.tif格式的图像文件，并截取第一个文件的部分文件名，从而获取其成像时间的具体年份。
接下来，遍历tif_file_path路径下全部.tif格式图像文件。其中，我们通过一个简单的判断语句if tif_file[0:4]==tif_file_year:，来确定某一年的遥感影像是否已经读取完毕——如果已经读取完毕，例如假如2001年成像的8幅遥感影像都已经遍历过了，那么就对这8景遥感影像加以逐像元的平均值求取，并开始对下一个年份（即2005年）成像的遥感影像继续加以计算；如果还没有读取完毕，例如假如2001年成像的8幅遥感影像目前仅遍历到了第5幅，那么就不求平均值，继续往下遍历，直到遍历完2001年成像的8幅遥感影像。
这里相信大家也看到了为什么我们要在前期先将文件夹中的文件按照“名称”排序——是为了保证同一年成像的所有遥感影像都排列在一起，遍历时只要遇到一个新的年份，程序就知道上一个年份的所有图像都已经遍历完毕了，就可以将上一个年份的所有栅格图像加以平均值求取。
在这里，逐像元的平均值求取其实也非常简单——我们对每一个像元分别执行以下操作：首先将该像元在当前年份里所有遥感影像的像素值相加，随后除以这一年份的遥感影像的数量，得到的就是该像元在这一年中像素值的平均值。
最后，通过if tif_file==tif_file_name[len(tif_file_name)-1]:这个判断，来确认是否目前已经遍历到文件夹中的最后一个图像文件。如果是的话，就需要将当前成像年份的所有图像进行平均值的求取，并宣告代码完成运行。
在 IDLE (Python GUI) 中运行代码。代码运行完毕后，我们看一下结果文件夹。可以看到，其中的图像已经是按照成像时间，分别完成平均值求取后的结果了。
在最后，还需要说明一点——用以上代码来求取长时间序列遥感影像的像元平均值，对于任意一个像元，只要该像元在任意一个时相的图像中是无效值（即为NoData），那么该像元在最终求出的平均值结果图中，像素值也将会是无效值NoData。针对这一问题的解决，我们将在下一篇博客中介绍。
这里是Ruby新手。完成一些练习后碰壁了。练习:计算一系列成绩的字母等级创建一个方法get_grade来接受测试分数数组。数组中的每个分数应介于0和100之间，其中100是最大分数。计算平均分并将字母等级作为字符串返回，即“A”、“B”、“C”、“D”、“E”或“F”。我一直返回错误:avg.rb:1:syntaxerror,unexpectedtLBRACK,expecting')'defget_grade([100,90,80])^avg.rb:1:syntaxerror,unexpected')',expecting$end这是我目前所拥有的。我想坚持使用下面的方法或.join，
项目介绍随着我国经济迅速发展,人们对手机的需求越来越大,各种手机软件也都在被广泛应用,但是对于手机进行数据信息管理,对于手机的各种软件也是备受用户的喜爱小学生兴趣延时班预约小程序的设计与开发被用户普遍使用,为方便用户能够可以随时进行小学生兴趣延时班预约小程序的设计与开发的数据信息管理,特开发了小程序的设计与开发的管理系统。小学生兴趣延时班预约小程序的设计与开发的开发利用现有的成熟技术参考,以源代码为模板,分析功能调整与小学生兴趣延时班预约小程序的设计与开发的实际需求相结合,讨论了小学生兴趣延时班预约小程序的设计与开发的使用。开发环境开发说明：前端使用微信微信小程序开发工具：后端使用ssm：VU
我对如何计算通过{%assignvar=0%}赋值的变量加一完全感到困惑。这应该是最简单的任务。到目前为止，这是我尝试过的:{%assignamount=0%}{%forvariantinproduct.variants%}{%assignamount=amount+1%}{%endfor%}Amount:{{amount}}结果总是0。也许我忽略了一些明显的东西。也许有更好的方法。我想要存档的只是获取运行的迭代次数。 最佳答案 因为{{incrementamount}}将输出您的变量值并且不会影响{%assign%}定义的变量，我
给定一个nxmbool数组:[[true,true,false],[false,true,true],[false,true,true]]有什么简单的方法可以返回“该列中有多少个true？”结果应该是[1,3,2] 最佳答案 使用转置得到一个数组，其中每个子数组代表一列，然后将每一列映射到其中的true数:arr.transpose.map{|subarr|subarr.count(true)}这是一个带有inject的版本，应该在1.8.6上运行，没有任何依赖:arr.transpose.map{|subarr|subarr.in
给定两个大小相等的数组，如何找到不考虑位置的匹配元素的数量？例如:[0,0,5]和[0,5,5]将返回2的匹配项，因为有一个0和一个5共同；[1,0,0,3]和[0,0,1,4]将返回3的匹配项，因为0有两场，1有一场；[1,2,2,3]和[1,2,3,4]将返回3的匹配项。我尝试了很多想法，但它们都变得相当粗糙和令人费解。我猜想有一些不错的Ruby习惯用法，或者可能是一个正则表达式，可以很好地回答这个解决方案。 最佳答案 您可以使用count完成它:a.count{|e|index=b.index(e)andb.delete_at
Ruby中如何“一般地”计算以下格式(有根、无根)的JSON对象的数量？一般来说，我的意思是元素可能不同(例如“标题”被称为其他东西)。没有根:{[{"title":"Post1","body":"Hello!"},{"title":"Post2","body":"Goodbye!"}]}根包裹:{"posts":[{"title":"Post1","body":"Hello!"},{"title":"Post2","body":"Goodbye!"}]} 最佳答案 首先，withoutroot代码不是有效的json格式。它将没有包
目标我正在尝试计算自给定日期以来周的距离，而无需跳过任何步骤。我更喜欢用普通的Ruby来做，但ActiveSupport无疑是一个可以接受的选择。我的代码我写了以下内容，这似乎可行，但对我来说似乎还有很长的路要走。require'date'DAYS_IN_WEEK=7.0defweeks_sincedate_stringdate=Date.parsedate_stringdays=Date.today-dateweeks=days/DAYS_IN_WEEKweeks.round2endweeks_since'2015-06-15'#=>32.57ActiveSupport的#weeks
我想使用Ruby检查数千对文件中的每对文件是否包含相同的信息。有人能指出我正确的方向吗？ 最佳答案 require'fileutils'FileUtils.compare_file('file1','file2')当且仅当文件file1和file2相同时返回true。 关于ruby-如何批量检查文件内容是否相同，我们在StackOverflow上找到一个类似的问题： https://stackoverflow.com/questions/33769865/
技术选型1，前端小程序原生MINA框架cssJavaScriptWxml2，管理后台云开发Cms内容管理系统web网页3，数据后台小程序云开发云函数云开发数据库（基于MongoDB）云存储4，人脸识别算法基于百度智能云实现人脸识别一，用户端效果图预览老规矩我们先来看效果图，如果效果图符合你的需求，就继续往下看，如果不符合你的需求，可以跳过。1-1，登录注册页可以看到登录页有注册入口，注册页如下我们的注册，需要管理员审核，审核通过后才可以正常登录使用小程序1-2，个人中心页登录成功以后，我们会进入个人中心页我们在个人中心页可以注册人脸，因为我们做人脸识别签到，需要先注册人脸才可以进行人脸比对，进
如何计算两个字符串之间的字符交集？例如(假设我们有一个名为String.intersection的方法):"abc".intersection("ab")=2"hello".intersection("hallo")=4好的，男孩女孩们，感谢你们的大量反馈。更多示例:"aaa".intersection("a")=1"foo".intersection("bar")=0"abc".intersection("bc")=2"abc".intersection("ac")=2"abba".intersection("aa")=2一些补充说明:维基百科定义intersection如下:Int