CN107179145B - 海洋表面温度漂流浮标点观测空间代表性评价算法 - Google Patents
海洋表面温度漂流浮标点观测空间代表性评价算法 Download PDFInfo
- Publication number
- CN107179145B CN107179145B CN201710548544.5A CN201710548544A CN107179145B CN 107179145 B CN107179145 B CN 107179145B CN 201710548544 A CN201710548544 A CN 201710548544A CN 107179145 B CN107179145 B CN 107179145B
- Authority
- CN
- China
- Prior art keywords
- error
- observation
- sst
- spatial
- remote sensing
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K15/00—Testing or calibrating of thermometers
- G01K15/007—Testing
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
- G01S7/497—Means for monitoring or calibrating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Radar, Positioning & Navigation (AREA)
- Bioinformatics & Computational Biology (AREA)
- Remote Sensing (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Evolutionary Biology (AREA)
- Algebra (AREA)
- Software Systems (AREA)
- Databases & Information Systems (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
本发明公开了一种海洋表面温度漂流浮标点观测空间代表性评价算法,在进行SST遥感产品验证时,首先需要根据抽样误差借助高分辨率SST遥感产品计算不同尺度的MSD、RMSD和RMD来定量评价SST多点位观测空间代表性误差;然后以漂流浮标点位观测为参考值,计算SST遥感产品的平均偏差、误差标准差及均方根误差并与漂流浮标多点位观测的空间代表性误差进行比较。本发明借助较高空间分辨率的红外SST遥感产品,从细尺度上公里级的抽样误差和粗尺度上栅格级的空间代表性误差两方面定量评价漂流浮标SST多点位观测在不同空间尺度下的空间代表性。
Description
技术领域
本发明涉及海洋遥感监测领域,具体涉及一种海洋表面温度漂流浮标多点观测空间代表性评价算法。
背景技术
海洋表面温度(sea surface temperature,SST)的浮标观测经常用于验证卫星遥感定量反演的SST产品,包括微波和红外遥感产品。由于浮标观测属于点观测,而遥感产品属于面观测,因此就会有这样两个问题:点观测能够在多大程度上很好地代表遥感面观测?在SST遥感产品验证的过程中,在什么样的尺度下应该考虑空间代表性误差?关于空间代表性问题的研究,目前更多地集中在单个观测站点的地表反射率和LAI两个参数的评价上,并且多数针对抽样误差进行评价,而针对SST浮标多点观测的空间代表性评价尚未见到。
发明内容
为解决上述问题,本发明提供了一种海洋表面温度漂流浮标多点位观测空间代表性评价算法,借助较高空间分辨率的红外SST遥感产品,从细尺度上的抽样误差和粗尺度上的空间代表性误差两方面定量评价漂流浮标SST多点位观测在不同空间尺度下的空间代表性。
为实现上述目的,本发明采取的技术方案为:
海洋表面温度漂流浮标多点观测空间代表性评价算法,在进行SST遥感产品进行验证时,首先需要根据抽样误差借助高分辨率SST遥感产品计算不同尺度(细尺度公里级和粗尺度栅格级)MSD、RMSD和RMD来定量评价SST多点位观测的空间代表性误差(其中,粗尺度栅格级的MSD、RMSD、RMD分别是是细尺度公里级高分辨率遥感观测与其多尺度邻域遥感观测的复合);然后以漂流浮标点位观测为参考值,计算SST遥感产品的平均偏差、误差标准差及均方根误差并与漂流浮标点位观测的空间代表性误差进行比较,如果SST遥感产品误差小于漂流浮标点位观测的空间代表性误差,漂流浮标点位观测的空间代表性误差可以忽略,说明在此尺度上漂流浮标的点位观测的空间代表性强;否则,空间代表性误差不可以忽略。
以高分辨率SST遥感产品的空间分辨率为基本单元,选择漂流浮标聚集的单元,计算细尺度公里级的抽样误差(公式1):
其中,SEkilo是细尺度公里级每一个基本单元的抽样误差,Xi为在一个基本单元范围内的n个漂流浮标观测,为基本单元范围内n个漂流浮标观测的平均值。
在抽样误差的基础上通过以下算式计算均方差(MSD)、均方根误差(RMSD)和相对平均偏差(RMD),评价细尺度公里级漂流浮标多点位观测的空间代表性误差:
通过以下方法完成粗尺度栅格级(>20km)空间代表性误差计算:
首先,以高分辨率SST遥感产品的像元为中心,根据公式(1)分别计算其5×5邻域、13×13邻域、25×25邻域的空间抽样误差,记为SEmac;然后,根据公式(5)复合细尺度公里级抽样误差:
SRgrid=SEkilo+SEmac (5)。
其次,根据不同尺度空间抽样误差,以及以下公式定量评价粗尺度栅格级漂流浮标多点位观测的空间代表性误差。
MSDptog=MSDpot4km+MSD4kmtog, (6)
再次,以以漂流浮标点位观测为参考值,计算SST遥感产品的平均偏差、误差标准差及均方根误差并与漂流浮标多点位观测的空间代表性误差进行比较,如果SST遥感产品误差小于漂流浮标点位观测的空间代表性误差,漂流浮标点位观测的空间代表性误差可以忽略,说明在此尺度上SST点位观测具有较强的空间代表性;否则,空间代表性误差不可以忽略。
本发明具有以下有益效果:
本发明借助较高空间分辨率的红外SST遥感产品,定量评价漂流浮标SST多点位观测在不同空间尺度下的空间代表性。本算法在细尺度空间代表性评价中通过抽样误差的计算,并在三维空间内可视化多点漂流浮标空间代表性误差的空间分布特征;在粗尺度上,以漂流浮标观测的点位位置为中心,在其多个邻域尺度上分别复合细尺度上的抽样误差来定量评价粗尺度栅格级漂流浮标多点位观测的空间代表性。
附图说明
图1为本发明实施例中4km尺度下空间代表性误差空间分布图
图2为本发明实施例中4km尺度下空间代表性误差时间特征。
图3为本发明实施例中20km、50km、100km尺度下空间代表性误差空间分布图
图4为本发明实施例中20km、50km、100km尺度下空间代表性误差时间特征。
图5为本发明实施例中20km、50km、100km尺度下MODIS SST空间变异分布(a:20km;b:50km;c:100km)
图6为本发明实施例中20km、50km、100km尺度下MODIS SST空间变异时间特征。
具体实施方式
为了使本发明的目的及优点更加清楚明白,以下结合实施例对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
实施例
此算法已经以亚印太交汇区为研究区,运用于2010年MODIS反演SST产品验证中:数据说明:2010年漂流浮标SST来源于http://www.meds-sdmm.dfo-mpo.gc.ca/isdm-gdsi/drib-bder/svp-vcs/index-eng.asp;MODIS Aqua星Ievel-3map产品,4km空间分辨率,8天合成数据,来源于http://gcmd.nasa.gov。
结果分析
1)4km尺度下空间代表性误差时空分布如图1-图2所示(仅展示部分结果)。
表1 4km尺度下空间代表性误差统计
2)20km,50km,100km尺度下空间代表性时空分布如图3-图4所示。
表2 20km、50km、100km尺度下空间代表性误差统计(第16周)
表3 20km、50km、100km尺度下空间代表性误差统计(第24周)
3)MODIS SST空间变异(尺度:20km,50km,100km)如图5所示。
4)MODIS SST验证
以漂流浮标点位观测为参考值,4km尺度下MODIS SST验证结果:
表4 4km MODIS SST验证结果统计
week | 1 | 2 | 3 | … | 44 | 45 | 46 |
B | -0.1620 | -0.1179 | -0.0854 | … | -0.3154 | -0.2273 | -0.1494 |
RMSE | 0.5530 | 0.5771 | 0.5935 | … | 0.8095 | 0.8209 | 0.5191 |
SD | 0.5290 | 0.5652 | 0.5876 | … | 0.7458 | 0.7891 | 0.4975 |
MaxB | 3.7863 | 3.3995 | 2.5071 | … | 8.9394 | 7.9626 | 4.3301 |
MinB | 6.1000e-07 | 4.9425e-04 | 4.1141e-04 | … | 2.7848e-04 | 5.7405e-05 | 2.2338e-04 |
表5 4km尺度下漂流浮标多点位观测空间代表性误差统计
week | 1 | 2 | 3 | … | 44 | 45 | 46 |
MSD | 0.0361 | 0.0198 | 0.0461 | … | 0.0240 | 0.0139 | 0.0217 |
RMSD | 0.1266 | 0.0978 | 0.1228 | … | 0.0983 | 0.0987 | 0.1074 |
RMD | 0.6482 | 0.4759 | 0.5587 | … | 0.4003 | 0.4112 | 0.4133 |
漂流浮标空间代表性误差小于MODIS 4km SST的验证误差,在4km尺度下直接验证时,漂流浮标点位观测的空间代表性误差可以忽略,说明在4km尺度上漂流浮标点位观测的空间代表性强。
表6-8分别显示了20km、50km、100km尺度下复合空间代表性误差特征
表6 20km尺度下空间代表性误差统计
week | 1 | 2 | 3 | … | 44 | 45 | 46 |
MSD | 0.1238 | 0.0674 | 0.1432 | … | 0.0626 | 0.0635 | 0.0872 |
RMSD | 0.2820 | 0.2282 | 0.2894 | … | 0.2272 | 0.2318 | 0.2578 |
RMD | 0.6000 | 0.4553 | 0.4553 | … | 0.4333 | 0.4344 | 0.4391 |
表7 50km尺度下空间代表性误差统计
week | 1 | 2 | 3 | … | 44 | 45 | 46 |
MSD | 0.1369 | 0.1061 | 0.1854 | … | 0.1559 | 0.0942 | 0.1491 |
RMSD | 0.3257 | 0.3004 | 0.3601 | … | 0.3569 | 0.2867 | 0.3329 |
RMD | 0.6624 | 0.4468 | 0.5722 | … | 0.4132 | 0.4223 | 0.4249 |
表8 100km尺度下空间代表性误差统计
week | 1 | 2 | 3 | … | 44 | 45 | 46 |
M5D | 0.1944 | 0.2506 | 0.2431 | … | 0.2572 | 0.1944 | 0.2144 |
RMSD | 0.3937 | 0.4472 | 0.4331 | … | 0.4631 | 0.4139 | 0.4389 |
RMD | 0.6672 | 0.4998 | 0.5755 | … | 0.4179 | 0.4277 | 0.4296 |
在20km,50km,100km尺度下,漂流浮标点位观测的空间代表性误差明显偏大,因此在卫星产品直接验证时,这三种尺度下漂流浮标的空间代表性误差不能被忽略。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (4)
1.海洋表面温度漂流浮标点观测空间代表性评价算法,其特征在于,在对海洋表面温度SST(sea surface temperature)遥感产品进行验证时,首先需要根据抽样误差借助高分辨率SST遥感产品计算不同尺度的均方差MSD、均方根误差RMSD和相对平均偏差RMD来定量评价SST多点位观测空间代表性误差,其中,粗尺度栅格级的MSD、RMSD、RMD分别是细尺度公里级较高分辨率遥感观测与其多尺度邻域遥感观测的复合;然后以漂流浮标点位观测为参考值,计算SST遥感产品的平均偏差、误差标准差及均方根误差并与漂流浮标多点位观测的空间代表性误差进行比较,如果SST遥感产品误差接近漂流浮标多点位观测的空间代表性误差,漂流浮标多点位观测的空间代表性误差可以忽略,说明在此尺度上SST点位观测具有较强的空间代表性;否则,空间代表性误差不可以忽略。
2.如权利要求1所述的海洋表面温度漂流浮标点观测空间代表性评价算法,其特征在于,通过以下算式计算均方差MSD、均方根误差RMSD和相对平均偏差RMD:
其中,式(1)中n代表:以高空间分辨率SST遥感产品像元为基本单元,在一个基本单元范围内聚集的漂流浮标点位观察的个数。
3.如权利要求2所述的海洋表面温度漂流浮标点观测空间代表性评价算法,其特征在于,通过以下方法完成细尺度公里级<=20km的抽样误差计算:
以高空间分辨率SST遥感产品像元为基本单元,选择漂流浮标聚集的单元,计算细尺度公里级的抽样误差公式(4),并在此基础上通过公式(1)、(2)、(3)计算空间代表性误差,
其中,SEkilo是细尺度公里级每一个基本单元的抽样误差,Xi为在一个基本单元范围内的第i个漂流浮标观测值,为基本单元范围内n个漂流浮标观测的平均值。
4.如权利要求1所述的海洋表面温度漂流浮标点观测空间代表性评价算法,其特征在于,以漂流浮标点位观测为参考值,计算SST遥感产品的平均偏差、误差标准差及均方根误差并与漂流浮标多点位观测的空间代表性误差进行比较,如果SST遥感产品误差接近漂流浮标点位观测的空间代表性误差,漂流浮标多点位观测的空间代表性误差可以忽略,说明在此尺度上SST点位观测具有较强的空间代表性;否则,空间代表性误差不可以忽略,在20km,50km,100km尺度下,漂流浮标点位观测的空间代表性误差明显偏大,因此在卫星产品直接验证时,这三种尺度下漂流浮标的空间代表性误差不能被忽略。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710548544.5A CN107179145B (zh) | 2017-07-01 | 2017-07-01 | 海洋表面温度漂流浮标点观测空间代表性评价算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710548544.5A CN107179145B (zh) | 2017-07-01 | 2017-07-01 | 海洋表面温度漂流浮标点观测空间代表性评价算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107179145A CN107179145A (zh) | 2017-09-19 |
CN107179145B true CN107179145B (zh) | 2019-12-13 |
Family
ID=59845611
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710548544.5A Active CN107179145B (zh) | 2017-07-01 | 2017-07-01 | 海洋表面温度漂流浮标点观测空间代表性评价算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107179145B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103400022A (zh) * | 2013-06-08 | 2013-11-20 | 杭州师范大学 | 一种海表面温度遥感数据集等纬度重构方法 |
CN103886213A (zh) * | 2014-04-03 | 2014-06-25 | 北京大学 | 一种农作物总初级生产力的遥感估算方法及系统 |
CN105300554A (zh) * | 2015-09-14 | 2016-02-03 | 中国人民解放军国防科学技术大学 | 基于分布式光纤传感的多功能海洋环境监测装置和方法 |
CN106404222A (zh) * | 2016-09-26 | 2017-02-15 | 西北工业大学 | 基于组合式高精度测温电缆的海洋温深剖面探测系统 |
CN106598917A (zh) * | 2016-12-07 | 2017-04-26 | 国家海洋局第二海洋研究所 | 一种基于深度信念网络的上层海洋热结构预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4719893B2 (ja) * | 2005-05-02 | 2011-07-06 | 国立大学法人佐賀大学 | 制御装置、制御方法、及びそのプログラム |
-
2017
- 2017-07-01 CN CN201710548544.5A patent/CN107179145B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103400022A (zh) * | 2013-06-08 | 2013-11-20 | 杭州师范大学 | 一种海表面温度遥感数据集等纬度重构方法 |
CN103886213A (zh) * | 2014-04-03 | 2014-06-25 | 北京大学 | 一种农作物总初级生产力的遥感估算方法及系统 |
CN105300554A (zh) * | 2015-09-14 | 2016-02-03 | 中国人民解放军国防科学技术大学 | 基于分布式光纤传感的多功能海洋环境监测装置和方法 |
CN106404222A (zh) * | 2016-09-26 | 2017-02-15 | 西北工业大学 | 基于组合式高精度测温电缆的海洋温深剖面探测系统 |
CN106598917A (zh) * | 2016-12-07 | 2017-04-26 | 国家海洋局第二海洋研究所 | 一种基于深度信念网络的上层海洋热结构预测方法 |
Non-Patent Citations (4)
Title |
---|
Comparison of WindSat and buoy-measured ocean products from 2004 to 2013;张雷 等;《Acta Oceanol.Sin》;20160131;第35卷(第1期);第67-78页 * |
Sea surface effects on the sea surface temperature" estimation by remote sensing;R. YOKOYAMA等;《INT.J.REMOTE SENSING》;19950228;第16卷(第2期);第227-238页 * |
中国近海锋面时空特征研究及现场观测分析;刘泽;《遥感技术与应用》;20130215;第32卷(第6期);全文 * |
台湾海峡海洋表面温度遥感监测模型研究;张春桂 等;《热带气象学报》;20090228(第01期);第291-301页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107179145A (zh) | 2017-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Cordoba et al. | Diagnosing atmospheric motion vector observation errors for an operational high‐resolution data assimilation system | |
Zhang et al. | Sensitivity of vegetation phenology detection to the temporal resolution of satellite data | |
CN104537222B (zh) | 基于遥感的区域植被覆盖对地表气温影响的估算方法 | |
Ariza-Villaverde et al. | Multifractal analysis applied to the study of the accuracy of DEM-based stream derivation | |
CN111191673B (zh) | 一种地表温度降尺度方法及系统 | |
Dong | Automated measurement of sand dune migration using multi-temporal lidar data and GIS | |
CN112819846B (zh) | 面向多云雨地区的基于多载荷遥感图像的水稻估产方法 | |
CN103927743B (zh) | 一种遥感成像中人造目标的检测方法 | |
Debella-Gilo | Bare-earth extraction and DTM generation from photogrammetric point clouds including the use of an existing lower-resolution DTM | |
Lu et al. | Investigating spatial patterns of persistent scatterer interferometry point targets and landslide occurrences in the Arno River Basin | |
CN109191408B (zh) | 快速循环地面气象融合方法、装置及服务器 | |
Montesano et al. | Boreal canopy surfaces from spaceborne stereogrammetry | |
CN115343226A (zh) | 一种基于无人机的多尺度植被覆盖度遥感计算方法 | |
Mukhopadhyay et al. | Spatio-temporal analysis of land use-land cover changes in Delhi using remote sensing and GIS techniques | |
CN103593852A (zh) | 基于同质图斑的高光谱影像异常探测方法 | |
JP2016118502A (ja) | 点群解析処理装置、方法、及びプログラム | |
CN107179145B (zh) | 海洋表面温度漂流浮标点观测空间代表性评价算法 | |
CN109726679B (zh) | 一种遥感分类误差空间分布制图方法 | |
CN112200845A (zh) | 一种图像配准方法和装置 | |
CN115859211A (zh) | 一种基于三温不确定度估算模型的地表温度产品融合方法 | |
CN105389466A (zh) | 一种校正尺度效应的中低分辨率遥感产品真值获取方法 | |
CN113689481A (zh) | 一种基于中分辨率影像的同名云点快速匹配算法 | |
CN112733906A (zh) | 一种全球高分辨率每日地表动量粗糙度长度估算方法 | |
Ke et al. | Reconstructing complete MODIS LST based on temperature gradients in northeastern Qinghai-Tibet Plateau | |
Yoon et al. | Measurement of the water level in reservoirs from TerraSAR-X SAR interferometry and amplitude images |
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 |