如何用MODIS数据在TIMESAT中提取物候参数

2023-11-06

感觉研究植被物候的帖子好少呀,也找不到关于TIMESAT的操作教程。这里便记载我学习植被物候信息提取的学习笔记。下面是我综合了很多零零散散的信息学习到的东西,但是还是有很多看不明白的地方,比如一些参数的设置上面,希望有知道的可以让我学习一下。

一、TimeSat下载及matlab配置

准备:

1.TimeSat安装包,我使用TIMESAT3.3版本

2.MATLAB软件,我电脑是锐龙处理器,就安装的2020b版本

大家可以直接去官网(Welcome to the TIMESAT pages!)注册下载。注册时候需要验证码,验证码需要一点别的才能刷出来。

3.timesat3.3安装

  打开matlab2020b(b更稳定。)

  点击设置路径

  单击“添加并包含子文件夹”

找到timesat压缩包解压的文件夹路径

单击timesat-matlab文件夹

单击选择文件夹

安装成功如下:

  命令行窗口输入 TIMESAT  (大写) ,enter便可以调出TIMESAT窗口。

 一、TimeSat使用

  1.文件准备。

TIMESAT只能使用.dat和.img。这里使用ENVI5.3进行将tif数据批量格式转换。

2. 组成文件序列。

将所有需要进行时序处理的文件放在一个文件夹中,并建一个txt文件记录文件数和每个图像文件的绝对路径。如下所示,第一行写文件数目。需要注意的是,TIMESAT只能识别三个周期及以上的数据,所以要准备3年及以上的数据,或者你在开始和结束年进行复制一年的数据,周期要完整,中间不能缺少,缺少的需要进行复制补全。

3. 运行TIMESAT程序。

主菜单页面如下图所示。程序分为三个部分,分别为数据准备、数据处理、数据后处理,每个过程会生成不同的文件。

4. 显示并检查图像。运行TSM_imageview.

 No of rows in image代表图像的行数,columns代表列数

5. 处理时间序列数据。运行TSM_GUI。

这一步会生成一个初步的设置文件*.set。

TSM_GUI

选择不同的过滤方法Fitting method。

Common settings 部分设置数据值的范围

CLass-specific setting部分对提取方法进行设置,比如动态阈值还是固定值及其参数。

右下角的表格展示了当前像元(正上方标注)在这么多时相上表现出来的生长季节物候情况。

设置完之后点击上排setting保存成一个文件供下一步处理。

用Output可以输出某像素的时间序列到txt文件中

在导入文件列表时,进行以下操作。

TSM_GUI导入图像

TSM_GUI导入图像

注意:①timesat处理数据的范围要求缺失值或恒定值不超过20%,因此选取运算范围时应尽量选取有变化的区域。

②timesat有默认的数据范围,应注意元数据的最大最小值未超过在基础设置中现实的范围。若有个别极大或极小的异常值,可调整数据范围以去除异常值。

•  注意这里可以在Show mages里面选择一些小范围的像元来做实验,过于大的范围可能导致这一步无法加载数据

6. 调整设置文件。运行TSM_setting。

上面的设置是针对小范围数据的,在这里更改Rows/Columns to process将数据的处理范围设为全部图幅或需要的图幅。保存一个新的设置文件。还可以输入土地利用类型文件:不同的植被类型理论上应该使用不同的阈值,可以通过查阅相关文献拟定合适的阈值

在这里插入图片描述

job name 相当于项目名称

image mode 选择1 image files

Data series files 中选择到上面做好的txt list文件。

根据实际情况选择image file type, 按照我这种方法写下来是32位的

填上行列号

要处理的行列号,一般与上面一样

处的年份NO.of years

这里软件会自动读取txt里面第一行的数据,根据年份计算每一年是多少期数据,这里一定要确保核对,比如是16天一个数的话,一年就是23期。

7. 处理时间序列图像。运行TSF Process。

选择上一步中的*.set文件,会开始按行处理。处理完后会显示如下。在目录里会生成两个文件,*.tts/*tpa文件。

TSM_progress

9. 生成拟合图像。运行TSF_fit2img。

打开原始的拟合数据*.tts,并按照以下对话输入。

TSF_fit2img

>>设置缺失像元的值

0

>>设置输出的文件名(无扩展名)

Whatever

>>设置输出的数据格式

1=8位,2=16位,3=32位

3

>>设置需要生成图像的时间点(若输入-1即全部输出)

54

10. 在ENVI/ARCGIS中打开拟合图像。

处理之后的图像是无格式的,使用原图像的头文件,修改为拟合图像的文件名,若图像大小有改变,记得在头文件中修改为正确的行列数。然后在ENVI中,OPENAS->binary->选中物候期数据->header info->input header infofrom->open->原始文件头文件->OK。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

如何用MODIS数据在TIMESAT中提取物候参数 的相关文章

  • 微信公众号网页开发之拍照、上传本地图片

    微信网页开发 JS SDK说明文档 https developers weixin qq com doc offiaccount OA Web Apps JS SDK html 0 绑定域名 登录微信公众平台进入 公众号设置 gt 功能设置
  • 前端blob下载文件

    在一个项目中 需要下载excel pdf word zip等数据模板 而后端返回的是一个文件流 前端就需要使用blob来下载 1 首先要在接口中添加 responseType blob 模板下载 export function downLo

随机推荐

  • wazuh常用内容、防御sql注入

    目录 安装wazuh 常用内容 检测sql注入 主动响应 安装wazuh 本地测试的话建议用ova文件 直接导入虚拟机就能用了 官网 Virtual Machine OVA Installation alternatives 常用内容 目录
  • (字典树)acwing835. Trie字符串统计 算法基础班第二讲

    题目 维护一个字符串集合 支持两种操作 I x 向集合中插入一个字符串 x Q x 询问一个字符串在集合中出现了多少次 共有 N 个操作 输入的字符串总长度不超过 105 字符串仅包含小写英文字母 输入格式 第一行包含整数 N 表示操作数
  • 【目标检测】20、ATSS: bridging the gap between anchor-based and anchor-free detection via ATSS

    文章目录 一 背景和动机 二 方法 2 1 分析不同 2 2 方法 2 3 嵌入 FCOS 三 效果 四 代码 本文贡献点 指出了 anchor based 和 anchor free 方法性能差别的决定性因素 如何定义正样本和负样本 提出
  • WPF编程,Live Charts使用说明(4)——主题

    可以为应用程序中的任何图表设置默认样式 必须根据需要选择颜色 大小和笔触 色彩 这组颜色将是系列中的默认颜色 当系列数大于主题中的颜色数时 这些颜色将重复 材质 基于Google的材质设计 基于MoedernUi设计的Metro 蓝色蓝调
  • 打印沙漏 C语言

    题目 本题要求你写个程序把给定的符号打印成沙漏的形状 例如给定17个 要求按下列格式打印 所谓 沙漏形状 是指每行输出奇数个符号 各行符号中心对齐 相邻两行符号数差2 符号数先从大到小顺序递减到1 再从小到大顺序递增 首尾符号数相等 给定任
  • android Intent常用标识

    Intent常用标识 FLAG ACTIVITY BROUGHT TO FRONT 这个标志一般不是由程序代码设置的 如在launchMode中设置singleTask模式时系统帮你设定 FLAG ACTIVITY CLEAR TOP 如果
  • 简单有限状态机代码实现

    学习链接 Unity有限状态机编写 木子微冷 博客园 实现思路 1 创建状态基类BasState 所有具体状态类继承该类 基类有3个方法 进入状态 状态中 离开状态 2 创建管理类来管理状态机 StateMgr 并实现状态机的各个方法 状态
  • MariaDB数据库

    MariaDB数据库是MySQL数据库的一个分支 仍是一个开放的数据库 数据和表定义文件是兼容的 所有客户端 API 和协议都是兼容的 文件名 二进制文件和路径在MySQL和MariaDB上是相同的 端口和套接字是相同的 所有MySQL连接
  • 接口测试问题以及经验合集

    一 接口测试常见问题 前景提要 由于本人测试小白 可能所遇问题都较为基础 测试小白可以参考 Q1 postman会报connect ECONNREFUSED jemeter会报org apache http conn HttpHostCon
  • 14:00面试,14:06就出来了,问的问题有点变态。。。

    从小厂出来 没想到在另一家公司又寄了 到这家公司开始上班 加班是每天必不可少的 看在钱给的比较多的份上 就不太计较了 没想到5月一纸通知 所有人不准加班 加班费不仅没有了 薪资还要降40 这下搞的饭都吃不起了 还在有个朋友内推我去了一家互联
  • async_retrying

    from async retrying import retry import aiohttp import asyncio retry attempts 6 async def fetch print 1 async with aioht
  • 南京特殊教育师范:校外如何登陆Webplus站群发布文章

    因为网络安全需要 Webplus站群2021年完成三级等保工作 为了满足三级等保整改要求 网站群后台禁止外网访问 如果在家需要发布文章 可以使用学校的VPN进入后台 发布文章 1 访问学校的WebVPN https vpn njts edu
  • javascript小说阅读器分页算法的实现

    好久好久不写代码了 也好久没更新博客了 这次就和大家分享一个电子书阅读器分页的算法吧 像一些主流的阅读器 如QQ阅读 iReader等 都实现了txt文档分页显示的功能 打开一个txt文档可以快速把文档分割成若干页 每页文字正好铺满屏幕 点
  • ctfshow web入门 命令执行34-38

    1 web34 和web33差不多 过滤了冒号 第一个为官方解法 payload c include GET 1 gt 1 php filter read convert base64 encode resource flag php c
  • 立一个小小的flag

    距离上次写博客已经很久了 下个月也是我入职一周年了 很久没有沉下心好好写一些文章 内心还有有点不安的 准备接下来花点时间好好写点东西 把工作学习中的事情记录下来 希望对你我都有帮助 17年还是大三的时候 每天都会坚持学习写博客 沉淀自己学到
  • 基础实验篇

    导读 uORB是PX4 Pixhawk系统中非常重要且关键的模块之一 是用于无人机模块间通信的协议机制 本篇将详细介绍uORB并详细拆解uORB消息读写与自定义实验 一 基础实验篇 uORB消息读写与自定义实验 一 ect Request
  • ctfshow--web入门(web101--web115&web123&web125-web133)

    ctfshow web入门 web101 web115 web123 web125 web133 web101
  • 什么是CSRF

    目录 1 CSRF是什么 2 CSRF漏洞检测 3 防御CSRF攻击 3 1 验证 HTTP Referer 字段 3 2 在请求地址中添加 token 并验证 3 3 在 HTTP 头中自定义属性并验证 1 CSRF是什么 CSRF Cr
  • 结合高阶函数聊聊useMemo和useCallback

    使用useMemo可以实现useCallback useCallback fn deps 其实相当于 useMemo gt fn deps 所以说 使用useMemo完全可以实现useCallback useCallback functio
  • 如何用MODIS数据在TIMESAT中提取物候参数

    感觉研究植被物候的帖子好少呀 也找不到关于TIMESAT的操作教程 这里便记载我学习植被物候信息提取的学习笔记 下面是我综合了很多零零散散的信息学习到的东西 但是还是有很多看不明白的地方 比如一些参数的设置上面 希望有知道的可以让我学习一下