CN105184076A - 一种遥感地表温度数据的时空一体化融合方法 - Google Patents
一种遥感地表温度数据的时空一体化融合方法 Download PDFInfo
- Publication number
- CN105184076A CN105184076A CN201510557427.6A CN201510557427A CN105184076A CN 105184076 A CN105184076 A CN 105184076A CN 201510557427 A CN201510557427 A CN 201510557427A CN 105184076 A CN105184076 A CN 105184076A
- Authority
- CN
- China
- Prior art keywords
- pixel
- temperature data
- surface temperature
- time
- window
- 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.)
- Granted
Links
- 238000007500 overflow downdraw method Methods 0.000 title claims abstract description 10
- 230000008859 change Effects 0.000 claims abstract description 29
- 238000000034 method Methods 0.000 claims abstract description 14
- 230000002123 temporal effect Effects 0.000 claims description 19
- 230000004927 fusion Effects 0.000 claims description 14
- 238000012216 screening Methods 0.000 claims description 5
- 238000012952 Resampling Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 4
- 238000005520 cutting process Methods 0.000 claims description 4
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 230000008439 repair process Effects 0.000 claims description 2
- 230000000295 complement effect Effects 0.000 abstract description 3
- 238000012545 processing Methods 0.000 abstract description 2
- 230000009467 reduction Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 241001269238 Data Species 0.000 description 2
- 230000000422 nocturnal effect Effects 0.000 description 2
- VAYOSLLFUXYJDT-RDTXWAMCSA-N Lysergic acid diethylamide Chemical compound C1=CC(C=2[C@H](N(C)C[C@@H](C=2)C(=O)N(CC)CC)C2)=C3C2=CNC3=C1 VAYOSLLFUXYJDT-RDTXWAMCSA-N 0.000 description 1
- 230000002411 adverse Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002689 soil Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Landscapes
- Radiation Pyrometers (AREA)
Abstract
本发明涉及定量遥感数据处理技术领域,公开了一种遥感地表温度数据的时空一体化融合方法,首先对遥感地表温度数据进行预处理,其次,利用多时相温度数据和对应的时间求解地表温度昼夜变化模型参数,实现空间无缝,最后在极轨卫星地表温度影像数据上逐窗口筛选相似像元,结合由多时相静止卫星与上述模型构成的时间差异、尺度差异、相似度差异以及欧氏距离等权重因子,共同构建时空一体化融合模型并求解生成空间无缝的高时空分辨率地表温度数据集。本发明充分利用极轨卫星与多时相静止卫星地表温度数据的时空互补信息和静止卫星地表温度的昼夜变化规律,不需要任何辅助数据,易于实现,实用价值高,可用于时空无缝的遥感地表温度数据业务化运行。
Description
技术领域
本发明涉及定量遥感数据处理技术领域,尤其是涉及一种遥感地表温度数据的时空一体化融合方法。
背景技术
地表温度(LST,LandSurfaceTemperature)是研究地-气间物质与能量交换的重要参数,广泛应用于城市热岛效应、土壤湿度估算、森林火灾监测、气候变化等领域。卫星热红外遥感是获取地表温度的重要手段,然而,遥感观测系统需要在空间分辨率和时间分辨率之间进行权衡与取舍,难以同时获得高时间和高空间分辨率的温度数据,极大限制了其应用的广度与深度。例如:静止卫星地表温度数据的时间分辨率最高达到15分钟,却具有至少3千米的低空间分辨率;而极地卫星提供的地表温度数据虽然具有较高的空间分辨率(约100米-1000米),但其时间分辨率一般都较低。另一方面,现有发布的地表温度数据由于云、气溶胶等不利天气的影响往往出现空值现象,导致空间不连续。因此,如何在不改变现有观测系统的条件下,实现遥感地表温度数据空间无缝的同时提高遥感LST的空间分辨率和时间分辨率,成为热红外遥感应用中的关键问题之一。
为了提高遥感LST数据的时空分辨率,国内外学者已经开展了一些研究。早期的研究主要有遥感LST空间降尺度/时间扩展:在空间降尺度方面,李小文院士提出将传统地学趋势面用于遥感尺度转换研究,其主要原理都是引入更高空间分辨率的地理要素构造趋势面,通过建立趋势面与LST间的关系完成降尺度转换。而对于遥感LST时间尺度扩展,一般是对多时相数据进行时相滤波、谐波分析等方法重建出中间时间的数据。以上针对单一遥感参量来提高时空分辨率的方法,对辅助数据质量或长时序性要求较高,且提高的分辨率程度有限,在某些应用也显得无能为力。因此,需发展一种新的LST时空降尺度方法。
多传感器遥感LST时空融合技术正是为了满足这种需求而发展起来的一种新技术。与上述针对单一遥感参量的空间降尺度、时间扩展方法不同的是,遥感LST时空融合是通过对具有不同空间尺度和时间尺度的LST进行融合处理,生成时间相对密集的高空间分辨率LST数据;其优势在于无需其它辅助数据条件,同时生成的数据具有输入数据中最高空间分辨率和最高时间分辨率。
然而,无论是传统时空降尺度还是时空融合,在算法实施过程中都假设所输入数据空间无缝,这个假设对于地表温度数据(由于受云等影响)经常不能满足;此外,在基于多时相、多尺度温度数据时空融合时,很少顾及参量昼夜变化规律、遥感观测条件等的差别,直接影响了融合影像的求解精度。因此如何利用多时相、多尺度时空互补信息实现遥感地表温度数据空间无缝同时提高其时空分辨率,需要研究新的时空一体化融合方法。
发明内容
本发明所要解决的问题是:针对现有观测系统的技术限制和遥感地表温度数据时空不连续现象,提供一种遥感地表温度数据的时空一体化融合方法。
为解决上述技术问题,本发明提出的技术方案为:一种遥感地表温度数据的时空一体化融合方法,设有一个空间无缝的极轨卫星地表温度数据M和同一地区的多时相静止卫星地表温度数据Gtp(p=0,1,…n且n>10);进一步假设在时间序列t0、t1…tn-1、tn上有一景地表温度数据G与极轨卫星地表温度数据M获取时间上对应,则实现时空一体化融合包括以下步骤:
步骤一、对地表温度数据M和Gtp(p=0,1,…n)进行预处理,包括几何配准、裁剪、温度单位转换、空值标记、重采样;
步骤二、基于昼夜变化模型的空间修复:根据当前时相待修复像元在其它时相对应有效像元值的数量分如下两种情况求解所述昼夜变化模型的参数,并进行空值修复;
情况1)若当前时相的待修复像元在其它时相上有足够非空值像元,直接利用这些有效值与对应时间构成的数据对按需求分白天和晚上进行所述昼夜变化模型参数拟合;
情况2)若待修复像元的其它时相有效像元及其对应时间构成的数据对不足以进行温度昼夜变化模型参数拟合时,则以相邻时相tx对应的非空值像元为中心像元开一个窗口,设窗口大小为wa×wa,wa根据G的特征而定,在窗口内筛选相似像元,筛选条件为:邻域像元与中心像元差的绝对值小于或等于该中心像元相关的函数,用公式表达为:其中G(i,tx)为时间tx时静止卫星地表温度数据G上窗口内邻域像元,G(wa/2,tx)为静止卫星地表温度数据G上窗口内中心像元,da为常量,使得筛选的相似像元满足其它时相上有足够多的有效数据对;利用这些有效数据对进行所述昼夜变化模型参数的拟合数据数据。
步骤三、构建遥感地表温度数据时空一体化融合模型:
在地表温度数据M上逐个以像元为中心像元开一个大小为wb×wb的窗口,窗口大小视M特征而定,按以下公式在窗口内筛选相似像元:
|M(i,t)-M(wb/2,t)|≤db×2M(wb/2,t)
其中M(i,t)为时间t时M上窗口内邻域像元,M(wb/2,t)为M上窗口内中心像元,db为常量,
计算相似像元地表温度数据在当前时间和预测时间间的温度变化值与中心像元所在昼夜变化模型的温度变化值的差异TDi如下,
TDi=|G(i,t)-G(i,tx)-(Temp(wb/2,t)-Temp(wb/2,tx))|
其中,G(i,t)为当前时间t的静止卫星地表温度数据G上窗口内第i个邻域像元,G(i,tx)为时间tx时静止卫星地表温度数据G上窗口内第i个邻域像元,i为正整数,其最大值为窗口大小;
结合TDi及其它参量构建以下相似像元的权重函数表达式:
Wi=(1/(Vi·SDi)/∑1/(Vi·SDi))
其中SDi是窗口内相似像元与中心像元的相似度、Vi是由M上窗口内相似像元与中心像元温度变化差异、欧式距离以及不同分辨率地表温度数据间尺度差异共同构成;假设地表温度数据M和Gtn在时间tn时对应,作为基础时间对,其它时间序列tx∈{t0、t1、…tn-1}作为预测时间,得到时空一体化融合模型:
其中,M(wb/2,tx)为融合的空间无缝且高时空分辨率的地表温度数据,具有静止卫星G的时间分辨率和极轨卫星M的空间分辨率。
进一步的,在所述步骤二中,由于考虑到观测条件、传感器误差、反演误差的原因使得现有地表温度数据存在包含一些噪声、误差较大的离群点,采用M估计理论提高模型参数拟合精度。
所述步骤二中,当待修复像元的其它时相像元都为空值时,选择当前时相邻域有效像元值,包括原有的和建模求解的;利用空间内插技术完成修复。
本发明充分结合多尺度、多时相传感器间的时空互补信息和静止卫星地表温度的昼夜变化规律,利用静止卫星高时间分辨率特性建立地表温度昼夜变化模型,进而实现空值像元的修复,同时与较高分辨率极轨卫星温度数据构建时空融合模型,最终生成空间无缝的高时空分辨率地表温度数据集。而且不需要任何辅助数据,易于实现,实用价值高,可用于业务化运行。本发明与现有技术相比具有空间修复精确、时空权重函数合理、时空一体化融合等优点。主要体现在:首先,根据待修复空值像元在时间维上的有效像元数分情况求解地表温度昼夜变化模型参数,并利用M估计等稳健回归方法使模型参数最优化,从而提高修复精度;其次,通过引入相似像元本身的相似度差异因子、基于温度昼夜变化模型的时间差异因子等构建更加合理的时空权重函数;最后,利用多时相、多尺度地表温度数据进行时空一体化融合可实现空间无缝的同时提高其时空分辨率。总之,新方法可以批量生成空间上连续且时间上相对密集的地表温度数据集,特别适用于具有高时间分辨率的静止卫星温度数据。
附图说明
图1是本发明的实施例流程图。
具体实施方式
以下结合附图和实施例详细说明本发明技术方案。
假设有一个空间无缝的高空间分辨率极轨卫星地表温度数据M和同一地区的多时相静止卫星地表温度数据Gtp(p=0,1,…n,且n>10)。进一步假设在时间序列t0、t1…tn-1、tn上有一景地表温度数据G与极轨卫星地表温度数据M获取时间上对应,通过时空一体化融合就能获取空间无缝且具有极轨卫星地表温度数据空间分辨率的Mtp(p=0,1,…n)。
如图1所示,本发明的实施例流程可以分为三个阶段:
第一阶段:地表温度数据的预处理
对地表温度数据M和Gtp(p=0,1,…n)进行预处理,一般包括几何配准、裁剪、温度单位转换、空值标记、重采样。几何配准是通过地面控制点及相应的经纬度信息完成投影转换与定位,根据所感兴趣的研究区域范围分别对所有温度数据进行裁剪,将所有温度数据转换成相同的单位,对地表温度数据Gtp空值像元进行标记,在完成空间修复后将Gtp(p=0,1,…n)重采样与地表温度数据M相匹配的空间分辨率。
第二阶段、空间无缝处理
步骤2-1、根据预处理中空值像元标记结果,若所有时相Gtp均无空值,直接进入第三阶段;否则,进入步骤2-2;
步骤2-2、确定地表温度昼夜变化模型参数求解的数据对(温度与时间)。根据当前时相待修复像元在其它时相对应有效像元值的数量分如下三种情况确定数据对:
情况一、若当前时相的待修复空值像元在其它时相上有足够数量(大于6)非空值像元,可确定这些有效温度Temp(t)及其观测时间t作为如下地表温度昼夜变化模型参数求解的数据对:
其中,b1和b2表示如下:
b1=T0+Tacos(β(ts-tm))-b2(2)
b2=[-βTasin(β(ts-tm))]/α(3)
式中个参数定义如表1所示
表1:各参数的意义
情况二、当待修复像元的其它时相有效像元数及其观测时间所组成的数据对不足以求解温度昼夜变化模型参数时,则以相邻时相tx对应的非空值像元为中心像元开一个窗口,设窗口大小为wa×wa,wa根据G的特征而定,在窗口内筛选相似像元,筛选条件为:邻域像元与中心像元差的绝对值小于或等于该中心像元相关的函数,用公式表达为:
其中G(i,tx)为时间tx时静止卫星地表温度数据G上窗口内第i个邻域像元,i为正整数,其最大值为窗口大小,G(wa/2,tx)为静止卫星地表温度数据G上窗口内中心像元,da为常量,使得筛选的相似像元满足其它时相上有足够多的有效像元。利用这些有效像元值及其观测时间来求温度昼夜变化模型参数。
情况三、当待修复像元的其它时相像元都为空值时,暂时不予确定数据对。
步骤2-3、温度昼夜变化模型参数求解。依次针对每一个能求其模型参数的空值像元,将静止卫星观测时间与对应的地表温度值带入公式(1)进行参数拟合求解。由于观测条件、传感器误差、反演误差等原因使得现有地表温度数据可能包含一些噪声、误差较大的离群点,为了提高参数求解的精度,必须加以剔除。为此,引入了基于M-估计理论的稳健拟合方法,计算出最优的参数。
令εi,t为(i,t)像素点的地表温度差异,即:
使其满足所有像元的地表温度差异和最小:其中ρ(·)有多种代价函数,本发明中使用经典的Huber代价函数:
其中,k为Huber参数。
为了方便求解,将尺度差异和最小化问题转化为迭代重加权最小二乘问题,即
这里
表示权重,ρ′(εi,t)是Huber函数ρ(εi,t)的导数,通过这种方式依次可以求出所有能建模空值像元的模型参数值,用于后续像元填充。
步骤2-4、空值像元填充。当空值像元所在的模型参数获得后,将上述空值像元观测时间带入模型,即可求得改时间所对应的温度值,完成像元填充。
在情况三下,由于没有建模,部分像元依然会存在空值,此时可选择当前时相邻域有效像元值(包括原有的和建模求解的),利用空间内插技术完成修复。由于静止卫星地表温度数据数量能够保证(目前发布的静止卫星一天最少24景,最多96景),这种特殊情况一般很少见。
第三阶段、时空一体化融合处理
步骤3-1、在M上以逐个以像元为中心像元开一个大小为wb×wb的窗口,wb根据M的特征而定,在窗口内筛选相似像元,筛选条件为:邻域像元与中心像元差的绝对值小于或等于该中心像元相关的函数,用公式表达为:
|M(i,t)-M(wb/2,t)|≤db×2M(wb/2,t)
其中M(i,t)为时间t时地表温度数据M上窗口内邻域像元,M(wb/2,t)为相应窗口内中心像元,db为常量,对不同的传感器会有稍微的变化;窗口内满足筛选条件的邻域像元即为相似像元。
步骤3-2、实施例根据M和Gtp之间的尺度差异、相似像元相似度、不同时相Gtp之间的时间差异以及欧氏距离,构建相似像元的权重函数表达式。对当前窗口内相似像元赋权,首先计算窗口内相似像元到中心像元的距离Di,表达式为:
这样的表达式是为了控制距离因子在1到1+20.5之间。
其次,计算窗口内各相似像元M(i,t)与中心像元M(wb/2,t)的相似度SDi,在计算之前须先求出它们的相似性差异Si如下,
Si=|M(i,t)-M(wb/2,t)|
由于相似性差异越小,表明该相似像元具有与中心像元更高的相似度,计算公式如下,
SDi=exp(-Si)/∑exp(-Si)
再次,计算当前时间不同分辨率地表温度数据M(i,t)与G(i,t)间的尺度差异Ri如下,
Ri=|M(i,t)-G(i,t)|
最后,计算相似像元地表温度数据在当前时间和预测时间间的温度变化值与中心像元所在模型的温度变化值的差异TDi如下,
TDi=|G(i,t)-G(i,tx)-(Temp(wb/2,t)-Temp(wb/2,tx))|
其中,G(i,t)为当前时间t的静止卫星地表温度数据G上窗口内第i个邻域像元;
合并上述尺度因子、时间因子和距离因子,得到Ei如下,
Ei=ln(Ri*100+1)*ln(TDi*100+1)*Di
为了与相似度因子匹配,需要对Ei进行归一化,得到
Vi=Ei/∑Ei(14)
为保证权重总和为1,建立如下归一化的权重函数Wit
Wi=(1/(Vi·SDi)/∑1/(Vi·SDi))
该式表示尺度差异越小、时间差异越小、相似度越高且离中心像元距离越近的邻域像元获得的权重值越大。
步骤3-3、计算融合后的地表温度值。设地表温度数据M和Gtn在时间tn时对应,可作为基础时间对,其它时间序列tx∈{t0、t1、…tn-1}作为预测时间,由地表温度数据M和Gtn和窗口内相似像元及其权重函数加权计算窗口内中心像元的温度值,
其中,M(wb/2,tx)为融合的空间无缝且高时空分辨率的地表温度数据,具有静止卫星G的时间分辨率和极轨卫星M的空间分辨率;
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (3)
1.一种遥感地表温度数据的时空一体化融合方法,其特征在于:设有一个空间无缝的极轨卫星地表温度数据M和同一地区的多时相静止卫星地表温度数据Gtp(p=0,1,…n且n>10);进一步假设在时间序列t0、t1…tn-1、tn上有一景地表温度数据G与极轨卫星地表温度数据M获取时间上对应,则实现时空一体化融合包括以下步骤:
步骤一、对地表温度数据M和Gtp(p=0,1,…n)进行预处理,包括几何配准、裁剪、温度单位转换、空值标记、重采样;
步骤二、基于昼夜变化模型的空间修复:根据当前时相待修复像元在其它时相对应有效像元值的数量分如下两种情况求解所述昼夜变化模型的参数,并进行空值修复;
情况1)若当前时相的待修复像元在其它时相上有足够非空值像元,直接利用这些有效值与对应时间构成的数据对按需求分白天和晚上进行所述昼夜变化模型参数拟合;
情况2)若待修复像元的其它时相有效像元及其对应时间构成的数据对不足以进行温度昼夜变化模型参数拟合时,则以相邻时相tx对应的非空值像元为中心像元开一个窗口,设窗口大小为wa×wa,wa根据G的特征而定,在窗口内筛选相似像元,筛选条件为:邻域像元与中心像元差的绝对值小于或等于该中心像元相关的函数,用公式表达为:其中G(i,tx)为时间tx时静止卫星地表温度数据G上窗口内邻域像元,G(wa/2,tx)为静止卫星地表温度数据G上窗口内中心像元,da为常量,使得筛选的相似像元满足其它时相上有足够多的有效数据对;利用这些有效数据对进行所述昼夜变化模型参数的拟合数据数据。
步骤三、构建遥感地表温度数据时空一体化融合模型:
在地表温度数据M上逐个以像元为中心像元开一个大小为wb×wb的窗口,窗口大小视M特征而定,按以下公式在窗口内筛选相似像元:
|M(i,t)-M(wb/2,t)|≤db×2M(wb/2,t)
其中M(i,t)为时间t时M上窗口内邻域像元,M(wb/2,t)为M上窗口内中心像元,db为常量,
计算相似像元地表温度数据在当前时间和预测时间间的温度变化值与中心像元所在昼夜变化模型的温度变化值的差异TDi如下,
TDi=|G(i,t)-G(i,tx)-(Temp(wb/2,t)-Temp(wb/2,tx))|
其中,G(i,t)为当前时间t的静止卫星地表温度数据G上窗口内第i个邻域像元,G(i,tx)为时间tx时静止卫星地表温度数据G上窗口内第i个邻域像元,i为正整数,其最大值为窗口大小;
结合TDi及其它参量构建以下相似像元的权重函数表达式:
Wi=(1/(Vi·SDi)/∑1/(Vi·SDi))
其中SDi是窗口内相似像元与中心像元的相似度、Vi是由M上窗口内相似像元与中心像元温度变化差异、欧式距离以及不同分辨率地表温度数据间尺度差异共同构成;假设地表温度数据M和Gtn在时间tn时对应,作为基础时间对,其它时间序列tx∈{t0、t1、…tn-1}作为预测时间,得到时空一体化融合模型:
其中,M(wb/2,tx)为融合的空间无缝且高时空分辨率的地表温度数据,具有静止卫星G的时间分辨率和极轨卫星M的空间分辨率。
2.根据权利要求1所述的遥感地表温度数据的时空一体化融合方法,其特征在于:在所述步骤二中,由于考虑到观测条件、传感器误差、反演误差的原因使得现有地表温度数据存在包含一些噪声、误差较大的离群点,采用M估计理论提高模型参数拟合精度。
3.根据权利要求1所述的遥感地表温度数据的时空一体化融合方法,其特征在于:所述步骤二中,当待修复像元的其它时相像元都为空值时,选择当前时相邻域有效像元值,包括原有的和建模求解的;利用空间内插技术完成修复。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510557427.6A CN105184076B (zh) | 2015-09-02 | 2015-09-02 | 一种遥感地表温度数据的时空一体化融合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510557427.6A CN105184076B (zh) | 2015-09-02 | 2015-09-02 | 一种遥感地表温度数据的时空一体化融合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105184076A true CN105184076A (zh) | 2015-12-23 |
CN105184076B CN105184076B (zh) | 2017-11-17 |
Family
ID=54906154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510557427.6A Expired - Fee Related CN105184076B (zh) | 2015-09-02 | 2015-09-02 | 一种遥感地表温度数据的时空一体化融合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105184076B (zh) |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105868529A (zh) * | 2016-03-18 | 2016-08-17 | 北京师范大学 | 一种基于遥感的近地表日均大气温度反演方法 |
CN106169058A (zh) * | 2015-12-25 | 2016-11-30 | 中国石油大学(华东) | 一种基于微波遥感与时空信息的云下像元lst估算方法 |
CN106340018A (zh) * | 2016-08-31 | 2017-01-18 | 中国水利水电科学研究院 | 水文气象要素空间插值最优分辨率的确定方法 |
CN107423537A (zh) * | 2017-01-25 | 2017-12-01 | 河海大学 | 一种基于自适应阈值的地表温度降尺度的方法 |
CN107682074A (zh) * | 2017-11-08 | 2018-02-09 | 南京天际易达通信技术有限公司 | 一种卫星上行信号发射时间补偿方法、装置及通信系统 |
CN107942410A (zh) * | 2017-09-19 | 2018-04-20 | 武汉船用机械有限责任公司 | 一种极地服务温度的预测方法和装置 |
CN109186774A (zh) * | 2018-08-30 | 2019-01-11 | 清华大学 | 地表温度信息获取方法、装置、计算机设备和存储介质 |
CN109903234A (zh) * | 2019-01-18 | 2019-06-18 | 重庆邮电大学 | 一种城市热力景观定量描述及多尺度特征分析方法 |
CN110532340A (zh) * | 2019-09-03 | 2019-12-03 | 华东师范大学 | 空间信息时空元数据构建方法 |
CN111159310A (zh) * | 2019-12-13 | 2020-05-15 | 中国地质大学(武汉) | 可扩展的、带有信息增益策略的图像生成时空融合算法 |
CN111651411A (zh) * | 2020-04-21 | 2020-09-11 | 成都信息工程大学 | 一种复杂地形遥感土壤水分产品降尺度方法 |
CN112464746A (zh) * | 2020-11-10 | 2021-03-09 | 清华苏州环境创新研究院 | 一种卫星影像和机器学习的水质监测方法及系统 |
CN114184280A (zh) * | 2021-12-07 | 2022-03-15 | 自然资源部国土卫星遥感应用中心 | 一种基于热量平衡的地表温度时间归一化方法 |
CN115359369A (zh) * | 2022-10-19 | 2022-11-18 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种基于时相自适应的山地卫星影像融合方法和系统 |
CN115905447A (zh) * | 2022-12-23 | 2023-04-04 | 中国农业科学院农业资源与农业区划研究所 | 基于广义三角帽不确定度模型的地表温度产品融合方法 |
CN116450700A (zh) * | 2023-03-09 | 2023-07-18 | 中国农业科学院农业资源与农业区划研究所 | 极轨卫星地表温度时间归一化方法、装置及电子设备 |
CN116611587A (zh) * | 2023-07-19 | 2023-08-18 | 中国气象局公共气象服务中心(国家预警信息发布中心) | 一种基于极轨-静止卫星融合技术的太阳能资源预测方法 |
CN117216504A (zh) * | 2023-11-08 | 2023-12-12 | 中国农业科学院农业资源与农业区划研究所 | 一种极轨卫星日均地表温度遥感估算方法 |
CN117930298A (zh) * | 2024-03-25 | 2024-04-26 | 中国科学院空天信息创新研究院 | 基于卫星温度和姿态误差建模的静轨卫星定位误差修正方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103034988A (zh) * | 2012-12-18 | 2013-04-10 | 武汉大学 | 一种任意传感器数量的时空定量遥感融合方法 |
CN103399322A (zh) * | 2013-08-02 | 2013-11-20 | 中国科学院寒区旱区环境与工程研究所 | 一种基于modis地表温度的多年冻土图自动更新方法 |
CN103983360A (zh) * | 2014-05-30 | 2014-08-13 | 中国科学院遥感与数字地球研究所 | 一种基于hj-1b irs卫星数据反演lst的方法 |
CN104236716A (zh) * | 2013-06-17 | 2014-12-24 | 中国科学院城市环境研究所 | 基于成对hj-1b影像时空信息的地表温度反演方法 |
CN104360351A (zh) * | 2014-12-04 | 2015-02-18 | 中国科学院南京地理与湖泊研究所 | 一种基于遥感数据的农业区地表温度高精度反演方法 |
CN104637027A (zh) * | 2015-02-26 | 2015-05-20 | 武汉大学 | 顾及非局部特性与时空变化的遥感数据时空定量融合方法 |
-
2015
- 2015-09-02 CN CN201510557427.6A patent/CN105184076B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103034988A (zh) * | 2012-12-18 | 2013-04-10 | 武汉大学 | 一种任意传感器数量的时空定量遥感融合方法 |
CN104236716A (zh) * | 2013-06-17 | 2014-12-24 | 中国科学院城市环境研究所 | 基于成对hj-1b影像时空信息的地表温度反演方法 |
CN103399322A (zh) * | 2013-08-02 | 2013-11-20 | 中国科学院寒区旱区环境与工程研究所 | 一种基于modis地表温度的多年冻土图自动更新方法 |
CN103983360A (zh) * | 2014-05-30 | 2014-08-13 | 中国科学院遥感与数字地球研究所 | 一种基于hj-1b irs卫星数据反演lst的方法 |
CN104360351A (zh) * | 2014-12-04 | 2015-02-18 | 中国科学院南京地理与湖泊研究所 | 一种基于遥感数据的农业区地表温度高精度反演方法 |
CN104637027A (zh) * | 2015-02-26 | 2015-05-20 | 武汉大学 | 顾及非局部特性与时空变化的遥感数据时空定量融合方法 |
Non-Patent Citations (4)
Title |
---|
何连财: "基于遥感对不同下垫面地表温度时空变化分析", 《测绘与空间地理信息》 * |
周义等: "热红外遥感图像中云覆盖像元地表温度估算研究进展", 《光谱学与光谱分析》 * |
杨丽萍等: "西安城区地表温度的遥感反演与时空演变分析", 《兰州大学学报(自然科学版)》 * |
杨贵军等: "黑河流域ASTER 与MODIS 融合生成高分辨率地表温度的验证", 《农业工程学报》 * |
Cited By (34)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106169058A (zh) * | 2015-12-25 | 2016-11-30 | 中国石油大学(华东) | 一种基于微波遥感与时空信息的云下像元lst估算方法 |
CN105868529A (zh) * | 2016-03-18 | 2016-08-17 | 北京师范大学 | 一种基于遥感的近地表日均大气温度反演方法 |
CN106340018B (zh) * | 2016-08-31 | 2019-07-12 | 中国水利水电科学研究院 | 水文气象要素空间插值最优分辨率的确定方法 |
CN106340018A (zh) * | 2016-08-31 | 2017-01-18 | 中国水利水电科学研究院 | 水文气象要素空间插值最优分辨率的确定方法 |
CN107423537A (zh) * | 2017-01-25 | 2017-12-01 | 河海大学 | 一种基于自适应阈值的地表温度降尺度的方法 |
CN107423537B (zh) * | 2017-01-25 | 2020-09-18 | 河海大学 | 一种基于自适应阈值的地表温度降尺度的方法 |
CN107942410A (zh) * | 2017-09-19 | 2018-04-20 | 武汉船用机械有限责任公司 | 一种极地服务温度的预测方法和装置 |
CN107682074A (zh) * | 2017-11-08 | 2018-02-09 | 南京天际易达通信技术有限公司 | 一种卫星上行信号发射时间补偿方法、装置及通信系统 |
CN107682074B (zh) * | 2017-11-08 | 2024-03-29 | 南京天际易达通信技术有限公司 | 一种卫星上行信号发射时间补偿方法、装置及通信系统 |
CN109186774B (zh) * | 2018-08-30 | 2019-11-22 | 清华大学 | 地表温度信息获取方法、装置、计算机设备和存储介质 |
CN109186774A (zh) * | 2018-08-30 | 2019-01-11 | 清华大学 | 地表温度信息获取方法、装置、计算机设备和存储介质 |
CN109903234A (zh) * | 2019-01-18 | 2019-06-18 | 重庆邮电大学 | 一种城市热力景观定量描述及多尺度特征分析方法 |
CN109903234B (zh) * | 2019-01-18 | 2023-04-07 | 重庆邮电大学 | 一种城市热力景观定量描述及多尺度特征分析方法 |
CN110532340A (zh) * | 2019-09-03 | 2019-12-03 | 华东师范大学 | 空间信息时空元数据构建方法 |
CN110532340B (zh) * | 2019-09-03 | 2022-01-28 | 华东师范大学 | 空间信息时空元数据构建方法 |
CN111159310A (zh) * | 2019-12-13 | 2020-05-15 | 中国地质大学(武汉) | 可扩展的、带有信息增益策略的图像生成时空融合算法 |
CN111159310B (zh) * | 2019-12-13 | 2023-09-29 | 中国地质大学(武汉) | 可扩展的、带有信息增益策略的图像生成时空融合方法 |
CN111651411A (zh) * | 2020-04-21 | 2020-09-11 | 成都信息工程大学 | 一种复杂地形遥感土壤水分产品降尺度方法 |
CN111651411B (zh) * | 2020-04-21 | 2023-03-10 | 成都信息工程大学 | 一种复杂地形遥感土壤水分产品降尺度方法 |
CN112464746B (zh) * | 2020-11-10 | 2023-09-12 | 清华苏州环境创新研究院 | 一种卫星影像和机器学习的水质监测方法及系统 |
CN112464746A (zh) * | 2020-11-10 | 2021-03-09 | 清华苏州环境创新研究院 | 一种卫星影像和机器学习的水质监测方法及系统 |
CN114184280A (zh) * | 2021-12-07 | 2022-03-15 | 自然资源部国土卫星遥感应用中心 | 一种基于热量平衡的地表温度时间归一化方法 |
CN114184280B (zh) * | 2021-12-07 | 2024-03-19 | 自然资源部国土卫星遥感应用中心 | 一种基于热量平衡的地表温度时间归一化方法 |
CN115359369A (zh) * | 2022-10-19 | 2022-11-18 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种基于时相自适应的山地卫星影像融合方法和系统 |
CN115905447B (zh) * | 2022-12-23 | 2023-06-30 | 中国农业科学院农业资源与农业区划研究所 | 基于广义三角帽不确定度模型的地表温度产品融合方法 |
CN115905447A (zh) * | 2022-12-23 | 2023-04-04 | 中国农业科学院农业资源与农业区划研究所 | 基于广义三角帽不确定度模型的地表温度产品融合方法 |
CN116450700B (zh) * | 2023-03-09 | 2024-02-13 | 中国农业科学院农业资源与农业区划研究所 | 极轨卫星地表温度时间归一化方法、装置及电子设备 |
CN116450700A (zh) * | 2023-03-09 | 2023-07-18 | 中国农业科学院农业资源与农业区划研究所 | 极轨卫星地表温度时间归一化方法、装置及电子设备 |
CN116611587B (zh) * | 2023-07-19 | 2023-11-07 | 中国气象局公共气象服务中心(国家预警信息发布中心) | 一种基于极轨-静止卫星融合技术的太阳能资源预测方法 |
CN116611587A (zh) * | 2023-07-19 | 2023-08-18 | 中国气象局公共气象服务中心(国家预警信息发布中心) | 一种基于极轨-静止卫星融合技术的太阳能资源预测方法 |
CN117216504B (zh) * | 2023-11-08 | 2024-02-13 | 中国农业科学院农业资源与农业区划研究所 | 一种极轨卫星日均地表温度遥感估算方法 |
CN117216504A (zh) * | 2023-11-08 | 2023-12-12 | 中国农业科学院农业资源与农业区划研究所 | 一种极轨卫星日均地表温度遥感估算方法 |
CN117930298A (zh) * | 2024-03-25 | 2024-04-26 | 中国科学院空天信息创新研究院 | 基于卫星温度和姿态误差建模的静轨卫星定位误差修正方法及装置 |
CN117930298B (zh) * | 2024-03-25 | 2024-05-28 | 中国科学院空天信息创新研究院 | 基于卫星温度和姿态误差建模的静轨卫星定位误差修正方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105184076B (zh) | 2017-11-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105184076A (zh) | 一种遥感地表温度数据的时空一体化融合方法 | |
CN110738252B (zh) | 空间自相关的机器学习卫星降水数据降尺度方法、系统 | |
CN112016052B (zh) | 基于多源数据的近地表日最高气温估算方法、系统及终端 | |
Voyant et al. | Prediction intervals for global solar irradiation forecasting using regression trees methods | |
Zhang et al. | Analysis on spatial structure of landuse change based on remote sensing and geographical information system | |
CN104239706A (zh) | 一种地面观测气温时空数据集的制备方法 | |
Song et al. | An improved surface soil moisture downscaling approach over cloudy areas based on geographically weighted regression | |
CN104156968A (zh) | 大面积地形复杂区域无人机序列影像快速无缝拼接方法 | |
CN107341513B (zh) | 基于稳健的固定阶数滤波模型的多源海洋表面温度遥感产品融合方法 | |
CN103198447A (zh) | 一种基于卫星云图的风矢场实时度量方法 | |
CN103034988B (zh) | 一种任意传感器数量的时空定量遥感融合方法 | |
Harty et al. | Intra-hour cloud index forecasting with data assimilation | |
CN111652404A (zh) | 一种全天气地表温度反演方法和系统 | |
CN115795399A (zh) | 一种多源遥感降水数据自适应融合方法和系统 | |
CN114387531A (zh) | 一种基于改进地理加权回归模型的地表温度降尺度方法 | |
Mandal | Geo-information based spatio-temporal modeling of urban land use and land cover change in Butwal municipality, Nepal | |
CN109840539B (zh) | 一种基于地块图斑的遥感时空数据融合方法 | |
Xiong et al. | Fusing Landsat-7, Landsat-8 and Sentinel-2 surface reflectance to generate dense time series images with 10m spatial resolution | |
CN110766229B (zh) | 一种基于“降尺度-融合”的高分辨率每日降水制图方法 | |
Quan et al. | Generating 60–100 m, hourly, all-weather land surface temperatures based on the Landsat, ECOSTRESS, and reanalysis temperature combination (LERC) | |
CN115795402B (zh) | 一种基于变分法的多源降水数据融合方法和系统 | |
CN111177652A (zh) | 一种遥感降水数据的空间降尺度方法及系统 | |
CN107688712B (zh) | 一种基于dem和ndvi的气温降尺度方法 | |
Li et al. | Prediction of AOD data by geographical and temporal weighted regression with nonlinear principal component analysis | |
CN115563582A (zh) | 油气管道沿线地质风险多尺度评估方法、装置及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20180830 Address after: 530023 5 Jian Zheng Road, Qingxiu District, Nanning, the Guangxi Zhuang Autonomous Region Patentee after: The Guangxi Zhuang Autonomous Region Basic Geographic Information Center Address before: 230601 No. 111 Kowloon Road, Hefei economic and Technological Development Zone, Anhui Patentee before: Anhui University |
|
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171117 Termination date: 20180902 |