CN110411982B - 一种利用地面气象观测资料估算地表反照率的方法 - Google Patents
一种利用地面气象观测资料估算地表反照率的方法 Download PDFInfo
- Publication number
- CN110411982B CN110411982B CN201910599095.6A CN201910599095A CN110411982B CN 110411982 B CN110411982 B CN 110411982B CN 201910599095 A CN201910599095 A CN 201910599095A CN 110411982 B CN110411982 B CN 110411982B
- Authority
- CN
- China
- Prior art keywords
- day
- albedo
- month
- average
- earth surface
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
- G01N21/4738—Diffuse reflection, e.g. also for testing fluids, fibrous materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01W—METEOROLOGY
- G01W1/00—Meteorology
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Atmospheric Sciences (AREA)
- Biodiversity & Conservation Biology (AREA)
- Ecology (AREA)
- Environmental Sciences (AREA)
- Building Environments (AREA)
- Photometry And Measurement Of Optical Pulse Characteristics (AREA)
Abstract
本发明公开了一种利用地面气象观测资料估算地表反照率的方法,首先获取目标站点的气象观测数据,包括该站点所在的地理位置、逐日的辐射作用层的叶面指数LAIk和日照时数nk,根据目标站点的地位位置计算太阳倾角δk;然后根据不同计算需求的时间尺度对数据进行整理,之后根据太阳倾角和地理位置数据计算目标站点的太阳正午时角天顶距,获得的数据计算晴空地表反照率的天顶值;对时角积分得到晴天地表反照率的日平均值;之后采用日照时数考虑云量对反照率的影响,计算地表反照率的日平均值。本发明具有适用范围广、计算精度高、便于应用等优点,是一种计算简便,计算精度高的方法。
Description
技术领域
本发明属于地面观测领域,涉及一种地表反照率计算方法,具体涉及一种利用地面气象观测资料估算地表反照率的方法。
背景技术
地表反照率是指地表反射辐射通量与入射辐射通量的比值,是模拟各种农林气象、生物化学和生物地球化学过程的关键参数。地表反照率主要受辐射作用层特性(例如裸土的颜色、湿度、粗糙度等)和其他外部因素(云量、风向、积雪)的影响。反照率可用两个背靠背水平安装的辐射计直接测量得到。然而,受直接测量有关的经济和技术的限制,仅有极少数的站点观测到地表反照率。例如,在2500多个气象站中,只有18个站点观测逐日的反射辐射,用于估算地表反照率。此外,由于地表特性和各种外部因素随时空变化,从特定站点测量得到的地表反照率不适用于其他未监测的地点。
为此,许多科技工作者致力于地表反照率估算方法的研究,提出了许多地表反照率的计算方法,大致分为基于地面观测的方法和基于遥感的方法两大类。虽然基于地面观测的模型在输入方面相当简单,在应用上也较为方便,但是大多数模型都是针对几个特定站点拟合得到的,能否推广应用于其他地方还有待检验。基于遥感的方法通常用于较大空间尺度的研究中,以往的许多研究都集中在需要遥感测量的地表反照率区域值上。自1965年Conover首次将卫星影像应用于确定云层和地表的反照率以来,科技工作者一直致力于从卫星平台利用中分辨率成像光谱仪(Moderate Resolution ImagingSpectroradiometer,MODIS)等遥感测量资料推算区域表面反照率。尽管MODIS反照率产品因其低成本的可获取性,被广泛用于区域表面反照率的估算。然而,这些产品的最小分辨率为1公里,周期为16天,这限制其在较小空间尺度和假定表面(如参考作物)上的应用。因此需要一种成本低廉,数据准确度高的地表反照率获取方法。
发明内容
本发明的目的是提出一种估算地表反照率的气象学方法,能够简便、快速、精确地利用气象观测资料估算地表反照率。
为了解决上述技术问题,本发明采用技术方案是:
一种利用地面气象观测资料估算地表反照率的方法,其特征在于,包括以下步骤:
步骤1、气象观测数据收集与整理,获取目标站点的气象观测数据,包括该站点所在的地理位置、逐日的辐射作用层的叶面指数LAIk和日照时数nk,其中k为表示时序的下标,根据目标站点的地理位置计算太阳倾角δk;
步骤2、根据太阳倾角和目标站点的地理位置计算目标站点的太阳正午时角天顶距,计算公式如下:
上式中,Zk,表示第k日或第k月平均太阳正午时角天顶距;
步骤3、根据步骤1和步骤2获得的数据计算晴空地表反照率的天顶值αzz,k;
步骤4、计算晴天反照率的日或月平均值,由于太阳角天顶距随太阳东升西落的时角不断变化,第k日或第k月的地表反照率平均值是随时间变化,对时角积分得到晴天地表反照率的日或月平均值;
步骤5、计算地表反照率的日或月平均值,采用日照时数考虑云量对反照率的影响,则地表反照率的日或月平均值为
上式中:C1和C2为常数,Nk为第k日或第k月平均最大可能日照时数。
作为优选,所述步骤1中,根据需要计算表反照率的时间尺度,分为以下四类对输入数据进行整理与分析:
第一类、对于需要计算长系列逐日的日平均地表反照率,采用长系列的逐日日照时数和叶面指数进行计算,对于缺测的日照时数,采用相邻时间的算术平均值代替,此时k为日序数,k=1对应每年的1月1日,k=365或366对应每年的12月31日,则第k日的太阳倾角δk计算公式为:
第二类、对于需要计算多年平均的逐日日平均地表反照率,算术平均得到多年平均的逐日日照时数和叶面指数,此时k为日序数,k=1对应每年的1月1日,k=365或366对应每年的12月31日,第k日的太阳倾角δk采用公式(1)计算得到;
第三类、对于需要计算长系列逐月的月平均地表反照率,将日照时数和叶面指数算术平均得到逐月的日照时数和叶面指数,此时k为月序数,k=1对应1月,k=12对应12月,第k月的太阳倾角的平均值δk计算公式为
δk=23.2sin(29.5k-94) 公式(2)
第四类、对于需要计算多年平均的逐月地表反照率,算术平均得到多年平均的逐月日照时数和叶面指数,此时k为月序数,k=1对应1月,k=12对应12月,第k月的太阳倾角的平均值δk采用公式(2)计算得到。
作为优选,所述步骤3中,所述晴空地表反照率的天顶值计算公式如下:
αzz,k=αzv,k[1-exp(-ζ·LAIk)]+exp(-ζ·LAIk)αzs,k
上式中:ζ为常数;αzs,k,αzv,k分别为太阳正午时裸土和植被反照率的天顶值,与太阳正午时角天顶距、海拔、降水量等有关,其计算公式为
αzs,k=Cs1Pr -0.5αzzsexp(bZk)+Cs2
αzv,k=Cv1Pr -0.5αzzvexp(bZk)+Cv2
上式中:b为与辐射作用层类型有关的常数;αzzs,αzzv,Cs1,Cs2,Cv1和Cv2为常数,采用观测反射辐射的气象站资料进行率定;Pr为该站点大气压强与标准大气压的比值。
作为优选,所述步骤3中,Cv2和Cs2的取值和该站点的多年平均降水量P有关,分为以下三种情况,1)当P≥400mm,即该站点位于湿润半湿润区时,Cs2=0.177,Cv2=0.145;2)当200≤P<400mm,即该站点位于半干旱区时,Cs2=0.216,Cv2=0.155;3)当P<200,即该站点位于干旱区时,Cs2=0.218,Cv2=0.164。
作为优选,所述步骤3中,所述Pr与该站点的海拔有关,其计算公式为
上式中:h为该站点的海拔。
作为优选,所述步骤4中,晴天地表反照率的日或月平均值的计算公式如下:
上式中,α0,k为第k日或第k月晴天地表反照率的平均值。
作为优选,所述步骤5中,最大可能日照时数Nk的计算公式如下:
本发明有益效果是:
本发明与现有技术相比具有下列优点和效果:1)本发明提出方法综合考虑纬度、海拔、日倾角、土壤湿度、植被特性等因素对地表反照率的影响,能够精确地反映地表反照率的时空变异特性;2)本发明提出方法采用我国气象站资料对本发明提出方法中的参数进行率定,能够用于我国不同气候区、干湿区,不同时间尺度下地表反照率的估算,适用范围广;3)本发明提出的方法是解析方法,可用于参考作物等假定表面反照率的计算,应用范围广;4)在辐射平衡分析中,本发明利用气象观测资料估算地表反照率,计算所需的基础数据与辐射平衡分析衔接,便于应用;5)与现有方法相比,本发明具有适用范围广、计算精度高、便于应用等优点,是一种计算简便,计算精度高的方法。
附图说明
图1为本发明实施例利用地面气象观测资料计算地表反照率的流程图。
具体实施方式
本发明提供一种利用地面气象观测资料估算地表反照率的方法,该方法包括以下步骤:
步骤1、目标站点的气象观测数据收集与整理。获取目标站点的气象观测数据,包括该站点所在的地理位置(经度L、纬度和海拔h)、逐日的辐射作用层的叶面指数LAIk和日照时数nk,其中k为表示序号的下标。纬度采用弧度表示,若位于北半球取正值,位于南半球则取负值。根据计算结果的时间尺度,分为以下四类对输入数据进行整理与分析。
1)对于需要计算长系列逐日的日平均地表反照率,采用长系列的逐日日照时数和叶面指数进行计算,对于缺测的日照时数,采用相邻时间的算术平均值代替。此时k为日序数,k=1对应每年的1月1日,k=365(平年)或366(润年)对应每年的12月31日,则第k日的太阳倾角δk计算公式为:
2)对于需要计算多年平均的逐日日平均地表反照率,算术平均得到多年平均的逐日日照时数和叶面指数。此时k为日序数,k=1对应每年的1月1日,k=365(平年)或366(润年)对应每年的12月31日,第k日的太阳倾角δk采用公式(1)计算得到。
3)对于需要计算长系列逐月的月平均地表反照率,将日照时数和叶面指数算术平均得到逐月的日照时数和叶面指数,此时k为月序数,k=1对应1月,k=12对应12月,第k月的太阳倾角的平均值δk计算公式为
δk=23.2sin(29.5k-94) 公式(2)
4)对于需要计算多年平均的逐月地表反照率,算术平均得到多年平均的逐月日照时数和叶面指数,此时k为月序数,k=1对应1月,k=12对应12月,第k月的太阳倾角的平均值δk采用公式(2)计算得到。
步骤2、计算目标站点的太阳正午时角天顶距。采用下公式计算正午太阳天顶角:
公式(3)中:Zk为第k日或第k月平均太阳正午时角天顶距。
步骤3、确定晴空地表反照率的天顶值。采用下公式计算晴空地表反照率的天顶值:
αzz,k=αzv,k[1-exp(-ζ·LAIk)]+exp(-ζ·LAIk)αzs,k 公式(4)
公式(4)中:ζ为常数,对于草地,取ζ=0.50;αzs,k,αzv,k分别为太阳正午时裸土和植被反照率的天顶值,与太阳正午时角天顶距、海拔、降水量等有关,其计算公式为
αzs,k=Cs1Pr -0.5αzzsexp(bZk)+Cs2 公式(5)
αzv,k=Cv1Pr -0.5αzzvexp(bZk)+Cv2 公式(6)
公式(5)和公式(6)中:b为与辐射作用层类型有关的常数,对于草地,取b=0.6090;αzzs,αzzv,Cs1,Cs2,Cv1和Cv2为常数,采用观测反射辐射的气象站资料进行率定。对中国而言,采用额济纳旗、北京、成都、上海、郑州、昆明、武汉、贵阳、三亚、广州、兰州、拉萨、喀什、格尔木等14个观测反射辐射的气象站1993至2015年的逐日观测资料进行率定,确定Cv1=—1.0,Cs1=—1.0,αzzs=0.050,αzzv=0.030,Cv2和Cs2的取值和该站点的多年平均降水量P有关,分为以下三种情况,1)当P≥400mm,即该站点位于湿润半湿润区时,Cs2=0.177,Cv2=0.145;2)当200≤P<400mm,即该站点位于半干旱区时,Cs2=0.216,Cv2=0.155;3)当P<200,即该站点位于干旱区时,Cs2=0.218,Cv2=0.164;Pr为该站点大气压强与标准大气压的比值,与该站点的海拔有关,其计算公式为
公式(7)中,h为该站点的海拔。
步骤4、计算晴天反照率的日或月平均值。由于太阳角天顶距随太阳东升西落的时角不断变化,第k日或第k月的地表反照率平均值是随时间变化的,对时角积分得到晴天地表反照率的日或月平均值为
公式(8)中:α0,k为第k日或第k月晴天地表反照率的平均值。
步骤5、计算地表反照率的日或月平均值。采用日照时数考虑云量对反照率的影响,则地表反照率的日或月平均值为
公式(9)中,C1和C2为常数,C1=0.832,C2=0.032;Nk为第k日或第k月平均最大可能日照时数,其计算公式为
实施例1
某灌区位于热带地区,灌区附近有气象站A(东经109.58°,北纬18.22°,高程419m),观测有1953年1月1日至2015年12月31日的逐日日照时数、全辐射和反射辐射,计算站点A的多年平均的月平均地表反照率。
采用本发明的技术方案推求该站点平均反照率的步骤如下:
按照步骤1对目标站点A的气象观测数据收集与整理。该站点位于北半球,取根据计算结果的时间尺度,将该气象站1993年至2015年观测的日照时数算术平均得到多年平均的月均日照时数,见表1。月均叶面指数采用多年平均值,见表1。同时,k表示月序数,k=1表示1月,k=12表示12月,采用公式(2)计算该站点各月的月均日倾角。
表1是本发明用于算例的计算过程和结果
按照步骤2计算目标站点的太阳正午时角天顶距。采用公式(3)计算该站点各月的太阳正午时角天顶距,见表1。
按照步骤3确定晴空地表反照率的天顶值。该站点高程h=419m,代入公式(6),Pr=0.98;多年平均降水量大于400mm,该站点位于湿润地区,取Cs2=0.177,Cv2=0.145代入公式(5)和公式(6),得到各月的αzs,k和αzv,k,由公式(4)计算得到αzz,k的值,见表1。
按照步骤4计算晴天反照率的日平均值。取b=0.6090,将各月的αzz,k,Zk代入公式(7),得各月的月均地表反照率α0,k,见表1。
按照步骤5计算月均地表反照率。采用公式(10)计算各月的月均最大可能日照时数,代入公式(9)得到各月的月均地表反照率。与该站点的实测值相比,平均相对误差仅为3.3%,表明本发明提出方法的计算精度较高,可用于地表反照率的计算。
上述实施例仅为本发明技术方案说明的算法举例,并不代表对本发明技术方案保护范围的限定,上述为详细展开说明内容均为现有技术中的已有公知常识或者采用已有公知常识即可实现。
Claims (4)
1.一种利用地面气象观测资料估算地表反照率的方法,其特征在于,包括以下步骤:
步骤1、气象观测数据收集与整理,获取目标站点的气象观测数据,包括该站点所在的地理位置、逐日或逐月的辐射作用层的叶面指数LAIk和日照时数nk,其中k为表示时序的下标,根据目标站点的地理位置计算太阳倾角δk;
步骤2、根据太阳倾角和目标站点的地理位置计算目标站点的太阳正午时角天顶距,计算公式如下:
上式中,Zk表示第k日或第k月平均的太阳正午时角天顶距;
步骤3、根据步骤1和步骤2获得的数据计算晴空地表反照率的天顶值αzz,k;
步骤4、计算晴天反照率的日或月平均值,由于太阳角天顶距随太阳东升西落的时角不断变化,第k日或第k月平均的地表反照率是随时间变化,对时角积分得到晴天地表反照率的日或月平均值;
所述步骤4中,晴天地表反照率的日或月平均值的计算公式如下:
上式中,α0,k为第k日或第k月晴天地表反照率的平均值;
步骤5、计算地表反照率的日或月平均值,采用日照时数考虑云量对反照率的影响,则地表反照率的日或月平均值为
上式中:C1和C2为常数,Nk为第k日或第k月平均最大可能日照时数;
所述步骤3中,所述晴空地表反照率的天顶值计算公式如下:
αzz,k=αzv,k[1-exp(-ζ·LAIk)]+exp(-ζ·LAIk)αzs,k
上式中:ζ为常数;αzs,k,αzv,k分别为太阳正午时裸土和植被反照率的天顶值,与太阳正午时角天顶距、海拔、降水量等有关,其计算公式为
αzs,k=Cs1Pr -0.5αzzsexp(bZk)+Cs2
αzv,k=Cv1Pr -0.5αzzvexp(bZk)+Cv2
上式中:b为与辐射作用层类型有关的常数;αzzs,αzzv,Cs1,Cs2,Cv1和Cv2为常数,采用观测反射辐射的气象站资料进行率定;Pr为该站点大气压强与标准大气压的比值。
2.如权利要求1所述的方法,其特征在于:所述步骤3中,Cv2和Cs2的取值和该站点的多年平均降水量P有关,分为以下三种情况,1)当P≥400mm,即该站点位于湿润半湿润区时,Cs2=0.177,Cv2=0.145;2)当200≤P<400mm,即该站点位于半干旱区时,Cs2=0.216,Cv2=0.155;3)当P<200,即该站点位于干旱区时,Cs2=0.218,Cv2=0.164。
4.如权利要求1所述的方法,其特征在于:所述步骤1中,k表示逐月数据;
第三类、对于需要计算长系列逐月的月平均地表反照率,将日照时数和叶面指数算术平均得到逐月的日照时数和叶面指数,此时k为月序数,k=1对应1月,k=12对应12月,第k月的太阳倾角的平均值δk计算公式为
δk=23.2sin(29.5k-94) 公式(2)
第四类、对于需要计算多年平均的逐月地表反照率,算术平均得到多年平均的逐月日照时数和叶面指数,此时k为月序数,k=1对应1月,k=12对应12月,第k月的太阳倾角的平均值δk采用公式(2)计算得到。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910599095.6A CN110411982B (zh) | 2019-07-04 | 2019-07-04 | 一种利用地面气象观测资料估算地表反照率的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910599095.6A CN110411982B (zh) | 2019-07-04 | 2019-07-04 | 一种利用地面气象观测资料估算地表反照率的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110411982A CN110411982A (zh) | 2019-11-05 |
CN110411982B true CN110411982B (zh) | 2021-09-24 |
Family
ID=68360388
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910599095.6A Active CN110411982B (zh) | 2019-07-04 | 2019-07-04 | 一种利用地面气象观测资料估算地表反照率的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110411982B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102103076A (zh) * | 2011-02-01 | 2011-06-22 | 中国科学院遥感应用研究所 | 地表反照率反演方法及*** |
CN102435586A (zh) * | 2011-09-16 | 2012-05-02 | 北京师范大学 | 地表反照率产品的生成方法及*** |
CN104297816A (zh) * | 2014-10-27 | 2015-01-21 | 广东省东莞市气象局 | 紫外辐射表的标定校正方法 |
CN106575450A (zh) * | 2014-05-13 | 2017-04-19 | 河谷控股Ip有限责任公司 | 通过反照率模型、***和方法的增强现实内容渲染 |
KR20180129443A (ko) * | 2017-05-26 | 2018-12-05 | 이화여자대학교 산학협력단 | 생장이 종료된 식생의 구조적 특징을 반영한 눈 알베도 결정 장치 및 그 방법 |
CN109948175A (zh) * | 2019-01-14 | 2019-06-28 | 广州地理研究所 | 基于气象数据的卫星遥感反照率缺失值反演方法 |
-
2019
- 2019-07-04 CN CN201910599095.6A patent/CN110411982B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102103076A (zh) * | 2011-02-01 | 2011-06-22 | 中国科学院遥感应用研究所 | 地表反照率反演方法及*** |
CN102435586A (zh) * | 2011-09-16 | 2012-05-02 | 北京师范大学 | 地表反照率产品的生成方法及*** |
CN106575450A (zh) * | 2014-05-13 | 2017-04-19 | 河谷控股Ip有限责任公司 | 通过反照率模型、***和方法的增强现实内容渲染 |
CN104297816A (zh) * | 2014-10-27 | 2015-01-21 | 广东省东莞市气象局 | 紫外辐射表的标定校正方法 |
KR20180129443A (ko) * | 2017-05-26 | 2018-12-05 | 이화여자대학교 산학협력단 | 생장이 종료된 식생의 구조적 특징을 반영한 눈 알베도 결정 장치 및 그 방법 |
CN109948175A (zh) * | 2019-01-14 | 2019-06-28 | 广州地理研究所 | 基于气象数据的卫星遥感反照率缺失值反演方法 |
Non-Patent Citations (3)
Title |
---|
A numerical analysis of aggregation error in evapotranspiration estimates due to heterogeneity of soil moisture and leaf area index;Qiting Chen et al.;《Agricultural and forest Meteorology》;20190226;第335-350页 * |
Characterization of radiation regimes in nonrandom forest canopies;CHRISTOPHER J.;《Tree Physiology》;19980206;第695--706页 * |
Parameterization of snow-free land surface albedo as a function of vegetation phenology based on MODIS data and applied in climate modelling;Diana Rechid et al.;《Theor Appl Climatol》;20080412;第245-255页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110411982A (zh) | 2019-11-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110716185B (zh) | 一种用于星载微波辐射计的在轨辐射校正方法 | |
CN110567891B (zh) | 一种冬小麦冠层叶绿素估测***及方法 | |
Gao et al. | Multisensor comparisons and validation of MODIS vegetation indices at the semiarid Jornada experimental range | |
Rautiainen et al. | Coupling forest canopy and understory reflectance in the Arctic latitudes of Finland | |
Han et al. | Monitoring model of corn lodging based on Sentinel-1 radar image | |
CN112986158B (zh) | 基于无人机多光谱数据的甜菜氮素营养检测方法及*** | |
CN114202675A (zh) | 一种基于无人机多光谱遥感影像的作物生长参数测定方法及*** | |
Zhang et al. | Development of the direct-estimation albedo algorithm for snow-free Landsat TM albedo retrievals using field flux measurements | |
CN106501816A (zh) | 一种枣树冠层氮含量的卫星遥感监测方法 | |
Colombo et al. | Phenological monitoring of grassland and larch in the Alps from Terra and Aqua MODIS images | |
CN115375036A (zh) | 遥感、光能利用率模型与气象融合的作物成熟期预测方法 | |
CN110210142B (zh) | 一种南方大型灌区水稻实时需水量测算方法 | |
CN102435309A (zh) | 农业田间成像高光谱仪场地反射率定标方法及*** | |
Yamamoto et al. | Uncertainty quantification in land surface temperature retrieved from Himawari-8/AHI data by operational algorithms | |
Shibayama et al. | Continuous monitoring of visible and near-infrared band reflectance from a rice paddy for determining nitrogen uptake using digital cameras | |
CN113063740A (zh) | 一种基于多源遥感数据的小麦冠层氮含量监测方法 | |
CN110411982B (zh) | 一种利用地面气象观测资料估算地表反照率的方法 | |
CN110658325B (zh) | 一种确定sm监测中lst/fvc空间的方法 | |
CN116628405A (zh) | 一种估算草地植被高度的遥感方法 | |
Liu et al. | Estimating wheat coverage using multispectral images collected by unmanned aerial vehicles and a new sensor | |
Liu et al. | Field campaign for quantitative remote sensing in Beijing | |
CN115096854A (zh) | 一种人工靶标支持的塔基高光谱遥感同步大气校正方法 | |
Yadav et al. | Comparison of GNSS and INSAT 3D sounder retrieved precipitable water vapour and validation with GPS sonde data over Indian sub-continent | |
Hasegawa et al. | Seasonal change of bidirectional reflectance distribution function in mature Japanese larch forests and their phenology at the foot of Mt. Yatsugatake, central Japan | |
Zhang et al. | Solar irradiance estimated from FY-2 data at some north China sites |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |