基于地球发动机的ndvi-et模型

openet-ndvi的Python项目详细描述


======
openet-ndvi-et模型
=====


基于python api的h引擎实现了一个简单的模型,用于计算蒸散量(et)作为归一化差分植被指数(ndvi)的线性函数。在这个模型中,参考et(etf)的分数被计算为:

。数学:

etf=m*ndvi+b

,其中m和b的默认值分别为1.25和0.0。实际ET计算如下:

…数学:

et=etf*etr


input collections
=etf*etr


ndvi-et模型目前用于以下地球引擎图像集合:


landsat-toa
*landsat/lc08/c01/t1-toa或landsat/lc08/c01/t1-rt-toa
*landsat/le07/c01/t1-toa或是地球卫星/LE07/C01/C01/T1/U/TOA
*陆地卫星/LT05/LT05/C01/C01/T1/TOA
*陆地卫星/LT04/LT04/C01/T1/T1/TOA
>陆地卫星sr
*陆地卫星/LC08/C01/T1/U sr
*陆地卫星/LE07/C01/T1/T1/U sr
*陆地卫星/LT05/C01/T1/T1/U sr
*陆地卫星/LT04/C01/C01/T1/T1/U sr
第2号哨兵到A
*第2号卫星/br/>*第2号哨兵到A
*第2号卫星/2号哨兵到A

模型结构
=====


与ndvi-et模型交互的主要方式是通过"collection"和"image"类。

==


collection类应用于生成et的图像集合(以及其他模型"变量")。这些集合可以仅用于图像"立交桥"日期,也可以插入到每日、每月或每年的时间步。可以为多个输入集合类型构建集合,例如合并Landsat 8和Sentinel 2。

r/>地球引擎集合ID的列表(请参见"输入集合"。
直接传递给collection.filterdate()调用的开始日期
ISO格式的开始日期字符串(即yyyy-mm-dd)。
结束日期
直接传递给.filterdate()调用的ISO格式的结束日期字符串。结束日期必须是独占的(即,数据将更新到此日期,但不包括该日期)。
传递给集合的geometry
e e.geometry()调用。
将包括具有与几何体相交的示意图的所有图像。

---——————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————————集合将被筛选为云覆盖率小于此值的图像。
可选,默认值为70%。
筛选器参数
输入集合的自定义筛选器参数。
此参数尚未完全实现。
Model参数
传递到映像类初始化的参数。
此参数尚未完全实现。

overpass method
----


variables
要计算/返回的变量列表。


interpolate method
----

要计算/返回的变量列表。
t_interval
用于插值和聚合值的时间间隔。
选项:"每日"、"每月"、"每年"、"自定义"
可选,默认值为"自定义"。
interp_method
插值方法。
选项:"线性"
可选,默认值为"线性"。
插值计算中要包括的开始日期之前和结束日期之后的额外天数。
可选,默认值为"线性"。它是32。


集合示例代码块::python

65,38.7399),
etr_source='idaho_epscor/gridmet',
etr_band='etr')\
.立交桥(变量=['et','etr','etf'])


月结=ndvi_et.集合(
集合=['landsat/lc08/c01/t1_toa',
开始日期='2017-06-01',
结束日期='2017-09-01',
几何体=ee.几何体.点(-121.5265,38.7399),
etr_source='idaho_epscor/gridmet',
etr_band='etr')\
.插值(变量=['et','etr','etf']t_interval='monthly')

image
=

图像类应用于处理单个图像、具有自定义过滤功能的图像集合,或将自定义参数应用于集合中的每个图像。


通常,使用下面列出的集合/传感器特定的帮助器方法之一初始化ndvi-et图像(见下文)。这些方法将波段重命名为通用命名方案,应用基本的云掩蔽,并且可以通过将其中一个辅助方法映射到图像集合上来构建。

请参阅"image mapping<;examples/image_mapping.ipynb>;``示例笔记本"以了解更多详细信息。

图像
-----------------------------------


要将陆地卫星集合1的类实例化为一个图像,请使用图像.from_landsat_c1_to a()方法。


输入的陆地卫星图像必须具有以下波段和属性:

==========================一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个又一个的5 b1、b2、b3、b4、b5、b7、b6、bqa
陆地卫星7 b1、b2、b3、b4、b5、b7、b6_VCID_1, BQA
LANDSAT_8 B2, B3, B4, B5, B6, B7, B10, BQA
================= =============================================

================= =============================================
Property Description
================= ===================================
系统:索引-陆地卫星场景ID
-必须采用地球引擎格式(例如LC08_)
系统:时间开始图像日期1970年以来以毫秒为单位的时间
航天器ID-用于确定哪种陆地卫星类型(用于波段重命名)
-必须是:landsat U 4、landsat U 5、landsat U 7或landsat U 8
要实例化该图像,则必须是:landsat U 4、landsat U 5、landsat U 7或landsat U 8
E类对于Landsat集合1 SR图像,请使用该图像。from_Landsat_c1_sr()方法。


输入的Landsat图像必须具有以下波段和属性:

=====================================br/>卫星波段名称
==================一个又一个又一个又一个又一个又一个又一个又一个又一个又一个、B2、B3、B4、B4、B5、B7、B7、B6、B7、B6、像素_ B1、B2、B3、B4、B4、B4、B4、B5、B7、B7、B6、B7、B6、B7、B6、B7、B6、B7、B7、B7、B7、B7、B7、B7、B7、B7、B7、B7、B7、B6、B7、B7、B6、B7、B6、B7、B6、B6、B7、B6、B6、Bb3,b4,b5,b6,b7,b10,像素qa
=========一、二、二、三、三、三、三、三、三、三、三、四、四、四、四、四、四、四、四、四、四、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六、六ty描述
==========================================
系统:索引-陆地卫星场景id
-必须采用地球引擎格式(例如,lc08_044033_20170716)
系统:自1970年以来以毫秒为单位的开始图像日期时间
卫星ITE-用于确定哪种陆地卫星类型(用于波段重命名)必须是:陆地卫星4号、陆地卫星5号、陆地卫星7号或陆地卫星8号将陆地卫星集合1的类实例化为图像,使用图像。from_sentinel2_toa()方法。


================= =============================================
Sentinel-2A B2, B3, B4, B8, B11, B12, QA60
Sentinel-2B B2, B3, B4, B8, B11, B12, QA60
================= =============================================

================= ====================================
属性描述
====================================
系统:索引-哨兵2场景ID
-必须采用地球引擎格式(例如20180716T183929 T185042_T10SGJ)
系统:时间启动图像DA1970年以来以毫秒为单位的时间
航天器名称-用于确定哨兵2类型
-必须是:哨兵-2a或哨兵-2b
-当前未使用或未选中
=========================================>

=========>示例
----

…代码块::python

=======


ndvi-et模型可以计算以下变量:

归一化差异植被指数[unitless]
etf
参考et[unitless]
etr
参考et(苜蓿)[mm]
et
实际et[mm]

al"calculate"方法,可用于返回多个变量的多波段图像(参见示例…


reference et
==



reference et数据源使用"etr戋source"和"etr戋band"参数进行控制。

ET(ETR),如果使用草地参考ET(ETO),则不会返回有效结果。

/openet/cimis/daily
https://cimis.water.ca.gov/spatialdata.aspx


example notebooks
==



>调用openet ndvi et模型的各种方法的详细jupyter笔记本在"examples"文件夹中提供。

*`computing daily et for a single landsat image<;examples/single_image.ipynb>;`
*`计算"天桥"et图像集合<;examples/collection_overpass.ipynb>;`
*`计算插值每月et图像集合<;examples/collection_interpolate.ipynb>;`

安装
====


基于python openet ndvi的et模块可以通过pip安装:

…代码块::console



dependencies
==

*`earthingine api<;https://github.com/google/earthingine api>;`\br/>*`openet core<;https://github.com/open et/openet core beta>;`\br/>
openet namespace包
===



每个openet模型都存储在"openet"文件夹(命名空间)中。然后可以将模型作为主openet模块的"点"子模块导入。

…代码块::console


|构建图像::https://travis-ci.org/open-et/openet-ndvi-beta.svg?branch=master
:alt:build status
:target:https://travis ci.org/open et/openet ndvi beta
。|版本图像::https://badge.fury.io/py/openet ndvi.svg
:alt:pypi上的最新版本
:目标:https://badge.fury.io/py/openet ndvi
。| codecov图像::https://codecov.io/gh/open et/openet ndvi beta/branch/master/graphs/badge.svg
:alt:coverage status
:目标:https://codecov.io/gh/open-et/openet-ndvi-beta

欢迎加入QQ群-->: 979659372 Python中文网_新手群

推荐PyPI第三方库


热门话题
java本机方法的源代码可用吗?   java如何使父方法抛出异常?   java Android以编程方式设置不同屏幕大小/密度的布局   java如何使用一个变量来管理所有客户端请求   java输入一个txt文件,每行有一组数字   json java从jsonobject获取jsonarray错误   java将一个(WAV)写入一个文件只会说一个单词(最后一个单词)   java Telnet忽略原始字节   proguard java。运行桌面应用程序时出现lang.VerifyError   java用左键移动JLabel?   java如何在jText区域验证选项卡?   文件服务器客户端Javasocket编程中的字符串搜索   java省略了JSTL中的最后一个逗号<c:out>   java如何找到if或else代码已执行的次数?   java JavaScript WebSocket send()方法未执行   浮点数声明上的java标识符预期错误   java这是指二进制搜索算法吗?   编译mod at:reobfJar java时的minecraft问题。util。拉链ZipException:重复条目   java检测特定的震动运动(如图所示:D)