CN104239706B - 一种地面观测气温时空数据集的制备方法 - Google Patents
一种地面观测气温时空数据集的制备方法 Download PDFInfo
- Publication number
- CN104239706B CN104239706B CN201410449297.XA CN201410449297A CN104239706B CN 104239706 B CN104239706 B CN 104239706B CN 201410449297 A CN201410449297 A CN 201410449297A CN 104239706 B CN104239706 B CN 104239706B
- Authority
- CN
- China
- Prior art keywords
- temperature
- data
- model
- time
- grid
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种地面观测气温时空数据集的制备方法,首先确定预制备数据集的时空分辨率与区域范围;对气温观测数据进行预处理,形成规范的驱动数据集;制备研究区域的宏观地理因素与微观地理因素栅格数据集,使其与研究区域的时空分辨率相对应;在时间维度上,通过构建小波神经网络动态组合预测模型,对气温数据的变化特性进行模拟预测,形成气温未来变化的时间序列集;在空间维度上对目标区域逐格点的计算气温数值,生成气温空间变化趋势的空间数据集;基于关键影响因子,对气温时间维和空间维同时进行模拟计算生成气温未来时空变化的情景预测数据集。本发明具有良好的扩展性,可推广到地学台站其他观测要素的时空数据产品制备与生产应用中。
Description
技术领域
本发明涉及地学领域数据挖掘、数据融合与数据制备领域,具体的说,是一种地面观测气温时空数据集的制备方法。
背景技术
气温对人类与自然系统有着重要的影响,其作为反映地球表层系统热量状况的综合环境指数,不仅是参与自然地域系统界限划分的关键指标,也是地表陆面过程模型模拟的重要参数,因此气温被广泛地应用于全球变化研究的各个领域。
根据已知地理空间的特性来探索未知地理空间的特性是许多地理研究的第一步,也是地理学的基本问题。现代社会对环境空间信息的需求日益广泛。高分辨率、栅格化的气温数据不仅是地学模拟的重要驱动数据,也是决策管理的重要信息来源,尤其是气候变化的空间数据信息,其对于区域水文、水资源分析以及区域水资源管理、旱涝灾害管理、生态环境治理都具有重要意义。
气温数据来源于地面观测网络,例如气象站、水文站、生态观测站及一些临时采样点等,尽管这些站点总体上越来越密,然而受观测手段、资源配置、区域特征等方面的限制,我们不可能每个地理格网点(或带)都建立观测站点或者架设观测仪器,特别是偏远地区,地学观测站点分布十分稀疏,如西北五省区国家气象站点不足200个,因此我们一般只能得到离散的、分布不规则的、详尽程度不一的观测点(或采样点)数据。如何利用这些有限的站点去获得时空连续的地学变量并制备数据产品,是地学应用研究首要解决的问题,也是近年来国内外地学研究的重要课题之一。
目前,在地学研究中对气温数据的研究集中于数学统计方法上,这些方法包括在时间序列上的传统模型如回归分析模型,自回归模型、滑动平均模型和一些统计指标,而在空间分布上,一般采用空间插值的方式对其分布特征与规律进行研究,比较常用的方法主要有泰森多边形方法、反距离权重法、径向基函数法、多项式插值法、及各种形式的克里格方法等。这些常规的模型对于非线性气温观测序列存在着局限性。传统的时间序列分析是从纯数学的角度来分析的,但是由于实际应用中时间序列具有不规则、混沌等非线性特征,很难对系统建立理想的模型,预测也显得无能为力,很难对于气温时间序列的变化趋势进行精确的把握,并且在时空分布上忽略了对其有重要影响的环境影响因子,制备的数据集空间变化上严重依赖于观测网格点值。
此外许多研究人员尝试利用GCMs(Global Climate Models)或者RCMs(RegionalClimate Models)模拟结果作为未来气候变化的情景模式,但是这些模式不仅存在着概念上的争论而且计算复杂,目前数据分辨率也很粗糙。
随着社会的发展和科学技术的进步,科研、生产与生活对环境观测数据的精细化的需求越来越迫切,不仅需要高分辨率的栅格数据,而且对地学变量的未来趋势以及过去情景的数据资料也十分关注,然而现阶段的分析方法还难以提供上述所要求的时空尺度分辨率较高的数据产品,因此对地学数据进行时空研究势在必行。
与此同时科研人员也不断提出新理论新方法,以克服气温数据中大量的不确定性和复杂性,其中小波分析与人工神经网络在地学研究中表现最为卓越,人工神经网络由于其具有自组织、自学习以及非线性逼近能力和小波分析的时频域聚焦能力,对时间序列的分析得到了有效的应用。这些方法仅针对环境观测要素的时空特性的某一方面进行了理论验证研究与应用,而未从系统上考虑方法间的有效耦合与集成,构成环境观测要素的时空分析系统,尤其是未来情境预测,至今还没有形成一种有效的地面观测气温时空数据集制备方法。
发明内容
本发明所要解决的技术问题在于弥补现有技术不足,将新理论新方法应用于气温数据的数据挖掘中,其目的在于提供一种地面观测气温时空数据集的制备方法,应用该方法能够真实直观地反映气温的时空变化规律,结构清晰明了、易于实现。
本发明的目的是采用以下技术方案解决上述技术问题。
一种地面观测气温时空数据集的制备方法,其步骤是:
1)数据准备与处理
1.1).首先确定气温数据制备目标区域的边界范围、时间分辨率TEMP tp 与空间分辨率TEMP sp ,利用区域矢量边界数据制作区域的ASCII GRID格式的栅格本底数据,并标定每个网格单元的经度和纬度,即GRIDi,j(λ,)形式,其中i为行号,j为列号,λ为该网格中心点的经度,为该网格中心的纬度;
1.2). 收集整理研究区域内的全部气象观测站点的气温观测记录,并按照站点位置信息以及观测气温数值,制备内容为STATIONid(No,ID,DATE,λ,,z,TEMP)的TXT格式文件,其中:No为观测记录编号,ID为世界气象组织WMO规划的气象站点唯一编码,DATE为观测记录时间,z为站点高程,TEMP为对应时间与地点的气温数值,文件名称为该站点的ID编号;
1.3). 对区域全部站点的观测记录,进行序列长度对比分析,确定序列时间维度的模拟预测的分段划分方案,其时间序列划分方式在模型设置与控制中进行定义;
1.4). 获取目标区域的数字高程模型( Digital Elevation Model,简称DEM)数据,首先利用STATIONid所包含的z值对其进行数据质量检验,对通过验证的数据,利用GIS软件通过拼接、裁剪、投影转换、重采样等预处理,生成与目标区域相对应的空间分辨率为TEMP sp 的ASCII GRID格式的栅格数据,并在此基础上提取目标区域的坡度(Gslope)、坡向(Gaspect)、高程( z )信息;
1.5). 坡度函数计算任一单元和八邻域单元间高程距离的变化的最大变化比率,即坡度值,地面平坦则坡度值较低,地势陡峭,则坡度值较高,其计算公式为:
式中Slope deg 以角度表示的坡度值[0-90°],arctan为反正切函数,表示水平变化率,为垂直变化率,其计算方法为:
式中a~i为网格高程值,x_cell_size和y_cell_siz为网格纵横向大小,在本发明中采用正方形网格,即x_cell_size = y_cell_siz=GRID i,j (TEMP sp ) i = GRID i,j (TEMP sp ) j ;
1.6). 坡向按照顺时针方向从0°(正北方)到360°进行度量,坡向的值代表着斜坡面的方向,平坦的坡面没有方向,赋值为-1 ;
1.7). NDVI(Normalized Difference Vegetation Index,归一化植被指数)数据经拼接、裁剪、坐标变换、重采样与研究区域范围相对应,时间分辨率为TEMP tp 与空间分辨率TEMP sp ;
2) 气温时间序列预测
2.1)首先对气温观测时间序列数据集进行自相关(ACF)和偏相关(PACF)分析,以相关性最大值确定预见期T;
2.2) 对气温时间观测序列进行小波分解,分解为一系列的高低频成分,然后把分解得到的不同频率成分,分解的层次取决于序列的复杂性,考虑到计算效率的问题,一般情况下n取值3层,小波函数采用DB小波系数族;
2.3) 把作为神经网络模型的输入部分,形成小波与神经网络的组合模型,即为
2.4)把时刻通过质量控制的预测结果TEMP t+T 纳入下一次样本集,对[(t+T)+T]时刻的气温值TEMP (t+T)+T做出预测,而样本总数不变,动态循环;
2.5) 模型输出内容为STATIONid_date(No,ID,DATE,λ,,TEMP)的TXT文本文件,并以“站点ID_时间”为文件名称;
3)气温空间化计算模拟
3.1)基于气温影响因子数据集与神经网络模型构建气温的空间化模拟模型,即,式中:f ANN 为气温空间化分布式计算的人工神经网络模型,TEMP i,j 为GRID i,j (λ i,j , i,j )的温度,λ i,j 为网格点GRID i,j 的经度, i,j 为纬度,z i,j 为海拔高度,G i,j 为地形因子(坡度、坡向等),V i,j 为地表覆盖,etc为其他与气温直接相关,并可获取的环境因子;
3.2)模型的输入内容为STATIONid(No, ID, DATE, λ, , z, Gslope, Gaspect,NDVI, etc)的TXT文本文件,按照模型运行设置的数据集划分方案,依据不同的组合方式对气温与环境影响因子之间的关系进行学习模拟,以确定模型最优输入因子组合;
3.3)确定输入因子组合后,并且模拟结果满足预设评价指标时,把区域内全部网格点的环境因子值,按照上述3.2) 确定的模式和格式,逐一输入模型,即得到全区域范围的气温空间化数据,气温空间化计算模拟输出结果为GRID i,j (λ, , DATE, TEMP)的文本文件,可基于GIS软件,对TXT格式数据进行栅格数据转换,并进行多维可视化显示;
4) 气温未来时空变化预测
4.1)气温未来时空变化预测将排除与时间密切相关的环境因子,则须修正气温空间化模型,即气温的未来空间变化预测模型为;
4.2) 模型的输入格式为内容为STATION id (No, ID, DATE, λ, , z, G slope , G aspect ) TXT文本格式,按照不同的组合方式确定模型最优输入,以确定模型最优输入因子组合;
4.3)的预测输出数据集与观测数据整体作为修正后的空间化模型总体样本集,按照最优组合进行空间趋势预测模拟,从而生成气温的未来变化时空数据集;
4.4) 模型输出结果为GRID i,j (No,λ, , DATE, TEMP)格式的文本文件,可基于GIS软件,对这些数据结果进行栅格数据转换,并进行多维可视化显示;
5) 气温时空数据集使用
已完全生成气温的时空数据集,对于任意一网格点的温度数据可通过纬度λ和经度查询,对于站点的查询则可以利用站点ID进行快速查询,按照具体应用情况转换为所需的数据格式文档;
6) 模型设置与控制
6.1) 人工神经网络模型采用三层神经网络模型,各层都有一个连接权重矩阵,其中上标表示层之间的权重矩阵和偏差,通常将输入部分成为输入层,第一层和第二层也称为隐含层,第三层作为网络的输出称为输出,则神经网络的输出矢量为:
其隐含层采用正切双弯曲转移函数(Tan-Sigmoid Transfer Function,TSTF),其定义为
6.2)本发明涉及的人工神经网络模型的数据集按照70%,15%,15%的比例进行划分,其中70%为训练样本集,15%为检验样本集,15%为测试样本集。
对输入数据的进行归一化处理,消除原始数据量纲和取值范围不同的影响,其变换公式为:
式中,S为原始数据集序列;S min ,S max 分别是原始数据集中的最小值和最大值;S’为归一化后的数据序列,S’ min ,S’ max 为归一化数据序列的最大值与最小值,公式中S’ min =﹣1,S’ max = 1。
模型评价指标为用平均绝对误差MAE、平均相对误差MRE与均方根误差RMSE作为检验精度的标准;MRE和MAE是一种比较直观的精度评价方法,它直观地反映出反映样本数据估值的总体误差或精度水平,而RMSE反映出样本数据的估值灵敏度和极值,其分别定义为:
注:O i 为观测值,S i 为模拟值,N为点值个数
当模型通过符合预期模型评价指标时,即可开展该区域气温时空数据产品的自动化制备。
本发明的优点和产生的有效效果:
本发明对气温数据进行了深度信息挖掘,实现了由多点气温观测记录到气温时空数据立方体的转变,由有限个观测点到任一网格点的数据扩展延伸。本发明在气温空间化上充分考虑了影响气温的众多环境因子,空间分辨率与数据制备范围可自定义,相对于常规空间插值模型相比更加合理、精确;在气温时间序列预测上建立的小波神经网络动态组合预测模型可以基于气温数据的变化特征,预测其未来变化趋势;通过两者的耦合,可生成气温未来时空变化的情景数据集。本发明有效地提高了气温时空分析挖掘的能力,为地区热量资源分析、农业生产建设、生态环境治理乃至全球变化深层次的地学研究提供良好的数据产品,此外该方法具有良好的扩展性,可推广到地学台站其他观测要素的时空数据产品制备与生产应用中。
附图说明
图1为本发明一种地面观测气温时空数据集制备方法流程图;
图2为本发明所涉及的数据结构与表达形式;
图3为本发明所采用的三层神经网络模型示意图;
图4为本发明实施例中国西北地区地形图与气象站点的空间分布;
图5为本发明实施例中国西北地区各个年份气象站数目分布;
图6为本发明实施例DEM与站点高程对比图(虚线为趋势线);
图7为本发明实施例制备的中国西北地区气温环境因子数据库;
图 8 本发明实施例中气温与环境影响因子相关性分析;
图 9 本发明实施例生成的2001-2005年中国西北地区年均气温时空数据集,2006年为未来情景数据。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明:
本发明公开了一种地面观测气温时空数据集的制备方法,包括以下部分:首先定义确定研究区域,形成本底栅格数据,对气温数据进行预处理,制备本底栅格数据,制备宏观地理因素(地理位置、地形、自然地理环境)与微观地理因素(坡度、坡向、小地形形态、植被、土壤、遮蔽度)等的驱动数据集,通过比较分析最终确定了经度、纬度、海拔(高程)、坡度、坡向与植被覆盖指数对温度空间分布重要的影响关系;其次构建小波神经网络动态预测模型,对气温在时间纬度上进行变化特性分析与趋势预测,在空间维度上利用多种环境因子驱动的人工神经网络模型对目标区域逐格点的计算气温数值,得出空间变化趋势;在此基础上可开展气温要素的时间数列集与空间数据集制备,并可对其时间维和空间维同时进行预测生成未来变化情景数据集。
本发明在地面观测气温时空数据集制备时,其基本流程如图1所示,具体实施方案由以下步骤来实现:
第一步:确定待制备气温时空数据集的区域范围与时空分辨率。本实施例以中国西北地区为研究区域,东经73°-112°,北纬31°-49°之间 (见图 4 ),气温数据制备时间范围为2001-2005年,并对2006年的年均气温变化情境做出预测,时间分辨率为年,空间分辨率为千米。该区域作为中国气候变化的敏感区和生态脆弱区,与其它地区相比,地貌类型复杂多样,气候变化非常复杂,影响因素众多,选择本例更能验证本发明方法的有效性。
第二步:数据准备与制备
1)中国形成区划图作为边界数据,并制备成ESRI公司的Shapefile格式,利用ArcGIS生成1km×1km的4227893个网格格点的栅格数据(见图7 Grid图层)。
2)气温观测站点数据为中国气象科学数据共享网提供的中国西北地区地面气候资料月多要素数据集,利用自编程序,提取每个站点的年均气温,站点的空间位置信息,以站点ID为文件名生成格式为STATIONid(No,ID,DATE,λ,,z,TEMP)格式的TXT文件,对其整理后共200个站点的数据集,数据集中存在着无效值、缺省值和空值。为了使数据序列具有分析意义,这些年份的数据采用多年气温均值进行代替。
3)对本实施例气象台站的气温记录进行分析,各站点数据集并非等长,其最长记录长度为55年,最短记录长度为23年,数据集规则如下:设序列长度为n,则训练样本集为ROUND (n×0.70),检验样本集为ROUND (n×0.15),测试样本集为(n - ROUND (n×0.70)- ROUND (n×0.15)),ROUND为四舍五入函数。
4)DEM数据采用来源于美国国家海洋和大气局的美国国家地球物理数据中心(NGDC,http://www.ngdc.noaa.gov/mgg/global/global.html) GLOBE项目提供的全球数字高程模型,利用ArcGIS软件从中获取中国西北地区的DEM数据,其空间分辨率为30″(约1km),基于Albers投影坐标系统,把原始数据重采样为1km×1km的格点数据(图7 DEM层)。为了检验DEM数据的精度,基于气象站点坐标提取其对应的DEM高程,并与站点的高程值进行对比,结果见图6,通过比较发现:两者呈较好的相关性,R2=0.99,DEM的高程与气象站点的高程相差50m以内的格点占82%,相差100m的部分站7%,平均差值0.5m,基本可满足使用要求。
5)按照本发明坡度,坡向制备方法,基于获取的DEM数据,利用ArcGIS软件获取该区域范围内的坡度坡向数据(图7 Slope与Aspect图层 )。
6)NDVI采用EOS的极地轨道环境遥感卫星Terra的中分辨率成像光谱仪MODIS(Moderate-Resolution Imaging Spectroradiometer) MOD13A3(Vegetation IndicesMonthly L3 Global 1km)产品,地面分辨率为1km,时间分辨率为月。NDVI数据经ENVI软件拼接、裁剪、坐标变换等与研究区时空尺度对应(图7 NDVI图层 )。
第三步:中国西北地区气温时间序列模拟预测
1)所有站点观测序列存在1年的滞后自相关,因此T = 1。选择2001年到2005年的年均气温,进行数据序列预测验证,并对2006年的西北地区年气温均值进行预测。
2)利用Matlab软件对每个站点观测序列进行3层小波分解,小波函数选择db2,分解为一系列的高低频成分。
3)按照本发明2.3所述,中国西北地区时间序列分析模型即为
4)模型参数设置为最小均方根误差且相关系数大于0.9,最大训练为1000次,对其2001-2006年的中国西北地区200个站点的气温值进行模拟预测,形成中国西北地区2001-2006年年均气温时间序列预测集,进行保存(见表 1)
表 1中国西北地区2001-2006年气温时间序列预测集
5)预报精度见表2,预测结果显示模型对气象站观测序列的整体模拟结果效果较好,绝对误差(MAE)均小于0.5°C,平均相对误差(MRE)在10%左右,均小于15%,均方根误差(RMSE)在0.5°C左右,相关系数R的平方值(R 2 )接近于1,因此该站点预测结果具有较高的可信度。
表2 中国西北地区气象站点2001-2005年均气温时间序列预测结果评价分析
第四步:中国西北地区年均气温空间化模拟
1) 根据环境因子对气象的影响,进行相关分析,如图8所示,在空间位置中,经度主要涉及到离海的远近,而纬度主要考虑太阳辐射的角度对温度分布的影响。根据研究区域气象观测站点的温度分析可以看出,南北 (YZ) 方向存在着温度上升趋势,在东西(XZ)方向存着温度下降-上升的宽“U”型趋势。温度与高程呈显著的负相关(R=0.91),NDVI也具有明显的负相关表现,坡度、坡向与气温的关系较为复杂,难以用线性函数进行描述。按照本发明的与所建立的气温相关的环境影响因子(经度,纬度,高程,坡度,坡向,NDVI),进行不同组合方案的模拟比较,以确定空间化产品制备方案,同时为了验证本发明方法在气温空间化方面的优越性,特与其他经典方法进行了对比。
2) 模型结果详见表3,从表中我们可以看出随着高程对本区域气温空间分布的影响最为显著,其次是NDVI,坡度、坡向对整体的精度表现不明显,通过与常规插值方法的对比,从中可以看出所建立模型的优越性,特别是考虑所有环境因子影响的模型精度最高,本发明构建的空间模型总体上误差较低,按照评价指标,均有明显的优越性,在12种方案中,前四种最优模型均是基于本发明构建的模型。
表3 中国西北地区气象站点2001-2005年均气温不同空间化模型结果评价分析
3)基于上述结果中国西北地区年均气温空间化模型为,可以基于现有数据进行空间数据集制备生产,制备的数据结构如图9,中国西北地区2000-2005年的年均气温空间数据集,并可直接应用科研、生产与生活中。
4) 由于NDVI严重依赖于时间的变化,是与气温类似的时间变化序列,无法直接获取其未来变化数据,并且与从评价指标上相差不大,每个网格节点的空间位置几乎没有任何变化,具备长期预测条件,因此中国西北地区年均气温未来情境模型为,其中按此模式可生成中国西北地区2006年年均气温的时空数据集,如图9 T2006。
Claims (1)
1.一种地面观测气温时空数据集的制备方法,其特征在于包含以下步骤:
数据准备与处理
a.首先确定气温数据制备目标区域的边界范围、时间分辨率TEMPtp与空间分辨率TEMPsp,利用区域矢量边界数据制作区域的ASCII GRID格式的栅格本底数据,并标定每个网格单元的经度和纬度,即GRIDi,j(λ,)形式,其中i为行号,j为列号,λ为该网格中心点的经度,为该网格中心的纬度;
b.收集整理研究区域内的全部气象观测站点的气温观测记录,并按照站点位置信息以及观测气温数值,制备内容为STATIONid(No,ID,DATE,λ,z,TEMP)的TXT格式文件,其中:No为观测记录编号,ID为世界气象组织WMO规划的气象站点唯一编码,DATE为观测记录时间,z为站点高程,TEMP为对应时间与地点的气温数值,文件名称为该站点的ID编号;
c.对区域全部站点的观测记录,进行序列长度对比分析,确定序列时间维度的模拟预测的分段划分方案,其时间序列划分方式在模型设置与控制中进行定义;
获取目标区域的数字高程模型(Digital ElevationModel,简称DEM)数据,首先利用STATIONid所包含的z值对其进行数据质量检验,对通过验证的数据,利用GIS软件通过拼接、裁剪、投影转换、重采样预处理,生成与目标区域相对应的空间分辨率为TEMPsp的ASCIIGRID格式的栅格数据,并在此基础上提取目标区域的坡度、坡向、高程信息;
e.坡度函数计算任一单元和八邻域单元间高程距离的变化的最大变化比率,即坡度值,地面平坦则坡度值较低,地势陡峭,则坡度值较高,其计算公式为:
式中Slopedeg以角度表示的坡度值[0-90°],arctan为反正切函数,表示水平变化率,为垂直变化率,其计算方法为:
式中a~i为网格高程值,x_cell_size和y_cell_siz为网格纵横向大小,所述网格为正方形网格,即
x_cell_size=y_cell_siz=GRIDi,j(TEMPsp)i=GRIDi,j(TEMPsp)j;
f.坡向按照顺时针方向从正北方0°到360°进行度量,坡向的值代表着斜坡面的方向,平坦的坡面没有方向,赋值为-1;
g.NDVI数据经拼接、裁剪、坐标变换、重采样与研究区域范围相对应,时间分辨率为TEMPtp与空间分辨率TEMPsp;
气温时间序列预测
h.首先对气温观测时间序列数据集进行自相关(ACF)和偏相关(PACF)分析,以相关性最大值确定预见期T;
i.对气温时间观测序列进行小波分解,分解为一系列的高低频成分{D1,D2,……,Dn,An},然后把分解得到的不同频率成分,分解的层次取决于序列的复杂性,考虑到计算效率的问题,n取值3层,小波函数采用DB小波系数族;
j.把{D1,D2,……,Dn,An}作为神经网络模型的输入部分,形成小波与神经网络的组合模型,即为
k.把时刻通过质量控制的预测结果TEMPt+T纳入下一次样本集,对[(t+T)+T]时刻的气温值TEMP(t+T)+T做出预测,而样本总数不变,动态循环;
l.模型输出内容为STATIONid_date文本文件,并以“站点ID_时间”为文件名称;
气温空间化计算模拟
m.基于气温影响因子数据集与神经网络模型构建气温的空间化模拟模型,即式中:fANN为气温空间化分布式计算的人工神经网络模型,TEMPi,j为的温度,λi,j为网格点GRIDi,j的经度,为纬度,zi,j为海拔高度,Gi,j为地形因子,包括坡度、坡向,Vi,j为地表覆盖,etc为其他与气温直接相关,并可获取的环境因子;
n.模型的输入内容为的TXT文本文件,按照模型运行设置的数据集划分方案,依据不同的组合方式对气温与环境影响因子之间的关系进行学习模拟,以确定模型最优输入因子组合;
o.确定输入因子组合后,并且模拟结果满足预设评价指标时,把区域内全部网格点的环境因子值,按照上述m.确定的模式,逐一输入模型,即得到全区域范围的气温空间化数据,气温空间化计算模拟输出结果为GRIDi,j(λ,DATE,TEMP)的文本文件,可基于GIS软件,对TXT格式数据进行栅格数据转换,并进行多维可视化显示;
气温未来时空变化预测
p.气温未来时空变化预测将排除与时间密切相关的环境因子,则须修正气温空间化模型,即气温的未来空间变化预测模型为
q.模型的输入内容为STATIONid文本格式,按照不同的组合方式确定模型最优输入,以确定模型最优输入因子组合;
r.的预测输出数据集与观测数据整体作为修正后的空间化模型总体样本集,按照最优组合进行空间趋势预测模拟,从而生成气温的未来变化时空数据集;
q.模型输出结果为GRIDi,j格式的文本文件,可基于GIS软件,对这些数据结果进行栅格数据转换,并进行多维可视化显示;
气温时空数据集使用
t.已完全生成气温的时空数据集,对于任意一网格点的温度数据可通过纬度λ和经度查询,对于站点的查询则可以利用站点ID进行快速查询,按照具体应用情况转换为所需的数据格式文档;
模型设置与控制
w.人工神经网络模型采用三层神经网络模型,各层都有一个连接权重矩阵其中上标表示i1,j1层之间的权重矩阵和偏差bi,通常将输入部分成为输入层,第一层和第二层也称为隐含层,第三层作为网络的输出称为输出,则神经网络的输出矢量为:
a3=f3(LW3,2f2(LW2,1f1(LW1,1+b1)+b2)+b3)
x.所述人工神经网络模型的数据集按照70%,15%,15%的比例进行划分,其中70%为训练样本集,15%为检验样本集,15%为测试样本集;
y.对输入数据的进行归一化处理,消除原始数据量纲和取值范围不同的影响,其变换公式为:
式中,S为原始数据集序列;Smin,Smax分别是原始数据集中的最小值和最大值;S’为归一化后的数据序列,S’min,S’max为归一化数据序列的最大值与最小值,公式中S’min=﹣1,S’max=1;
z.模型评价指标为用平均绝对误差MAE、平均相对误差MRE与均方根误差RMSE为检验精度的标准;MRE和MAE是一种直观的精度评价方法,它直观地反映出样本数据估值的总体误差或精度水平,而RMSE反映出样本数据的估值灵敏度和极值,其分别定义为:
注:Oi为观测值,Si为模拟值,N为点值个数
当模型通过符合预期模型评价指标时,即可开展该区域气温时空数据产品的自动化制备。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410449297.XA CN104239706B (zh) | 2014-09-05 | 2014-09-05 | 一种地面观测气温时空数据集的制备方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410449297.XA CN104239706B (zh) | 2014-09-05 | 2014-09-05 | 一种地面观测气温时空数据集的制备方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104239706A CN104239706A (zh) | 2014-12-24 |
CN104239706B true CN104239706B (zh) | 2017-08-15 |
Family
ID=52227756
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410449297.XA Active CN104239706B (zh) | 2014-09-05 | 2014-09-05 | 一种地面观测气温时空数据集的制备方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104239706B (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104484539B (zh) * | 2015-01-04 | 2017-08-11 | 北京师范大学 | 一种多参数时空气象场模拟方法及装置 |
CN106647355B (zh) * | 2016-11-09 | 2023-11-07 | 中国民用航空飞行学院 | 一种飞行情景环境评价的数据处理方法和系统 |
CN106682415A (zh) * | 2016-12-23 | 2017-05-17 | 中国科学院深圳先进技术研究院 | 一种业务数据显示方法及装置 |
CN107942410B (zh) * | 2017-09-19 | 2020-04-21 | 武汉船用机械有限责任公司 | 一种极地服务温度的预测方法和装置 |
CN108007440A (zh) * | 2017-11-27 | 2018-05-08 | 天津科技大学 | 一种原位测量营养盐垂向通量的方法 |
CN109191408B (zh) * | 2018-04-19 | 2022-03-01 | 中国气象局公共气象服务中心 | 快速循环地面气象融合方法、装置及服务器 |
WO2019224347A1 (en) * | 2018-05-25 | 2019-11-28 | Basf Agro Trademarks Gmbh | Training a machine learning algorithm and predicting a value for a weather data variable, especially at a field or sub-field level |
US10767348B2 (en) * | 2018-07-30 | 2020-09-08 | Deere & Company | Machine stability detection and control |
CN109885808B (zh) * | 2018-12-26 | 2022-09-13 | 西安建筑科技大学 | 一种近地表气象要素计算方法 |
CN109902120B (zh) * | 2019-01-16 | 2023-06-02 | 中山大学 | 一种面向栅格文本降水数据的可视化分析方法 |
CN110503142B (zh) * | 2019-08-14 | 2023-05-09 | 淮阴师范学院 | 基于时空随机效应加法过程模型的多源遥感产品层次贝叶斯融合方法 |
CN110633282A (zh) * | 2019-09-18 | 2019-12-31 | 四川九洲空管科技有限责任公司 | 一种空域资源多级三维网格化方法及工具 |
CN110929808A (zh) * | 2019-12-11 | 2020-03-27 | 国网湖南省电力有限公司 | 舞动温度的多要素智能修正方法及系统 |
CN111460637B (zh) * | 2020-03-20 | 2023-07-21 | 成都市环境保护科学研究院 | 基于数值方法的城市通风潜势定量评估方法 |
CN111750838B (zh) * | 2020-07-07 | 2022-06-28 | 中南民族大学 | 农业用地规划图的生成方法、装置、设备及存储介质 |
CN112115414B (zh) * | 2020-07-16 | 2023-10-24 | 华东师范大学 | 一种广布种分布范围的预测方法 |
CN111912537B (zh) * | 2020-07-31 | 2021-04-06 | 中国科学院地理科学与资源研究所 | 基于格网的高温实时预警发布和显示方法、系统及设备 |
CN112684520B (zh) * | 2020-12-30 | 2022-12-13 | 北京墨迹风云科技股份有限公司 | 一种气象预报的订正方法、装置、计算机设备及存储介质 |
CN113761756B (zh) * | 2021-09-26 | 2022-05-06 | 中国农业科学院农业资源与农业区划研究所 | 一种表面温度高温和低温数据集重构方法 |
CN114997055B (zh) * | 2022-06-06 | 2024-04-05 | 安徽理工大学 | 海平面温度时频域变化特性分析方法 |
CN115690773B (zh) * | 2022-12-26 | 2023-04-07 | 武汉天际航信息科技股份有限公司 | Dem分区及重建方法、计算设备和存储介质 |
CN117313307A (zh) * | 2023-06-26 | 2023-12-29 | 北京航空航天大学 | 一种融合时空环境信息的气候模型模拟温度数据修正方法 |
CN117370714B (zh) * | 2023-12-07 | 2024-03-19 | 南京气象科技创新研究院 | 一种代表站定量确定方法 |
-
2014
- 2014-09-05 CN CN201410449297.XA patent/CN104239706B/zh active Active
Non-Patent Citations (9)
Title |
---|
2001-2010年甘肃省植被覆盖的时空变化;崔丹丹等;《中国沙漠》;20140731;第34卷(第4期);第1161-1166页 * |
DEM结构特征对坡度坡向的影响分析;刘学军等;《地理与地理信息科学》;20041130;第20卷(第6期);第1-5,39页 * |
Spatial Interpolation of meteorology monitoring data for Western China using Back-Propagation Artificial Neural Networks;Yaonan Zhang等;《IEEE》;20101231;第1-5页 * |
基于ArcGIS的坡度分析;杜朝正;《资源开发与市场》;20091231;第25卷(第1期);第17-18页 * |
基于小波过程神经网络的短期风速预测方法研究;姚莹;《中国优秀硕士学位论文全文数据库(电子期刊)》;20140331;第2014年卷(第3期);第C042-365页 * |
多年平均气温空间化BP神经网络模型的模拟分析;张赛;《地球信息科学学报》;20110831;第13卷(第4期);第534-538页 * |
我国地面气温格点化数据集的研制;张强等;《中国科技资源导刊》;20080531;第40卷(第3期);第12-20页 * |
气象数据集说明文档;z23wo;《百度文库》;20140124;第1-4页 * |
气象监测数据的时空特征分析与建模;刘晓晓等;《地理空间信息》;20090831;第7卷(第4期);第104-107页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104239706A (zh) | 2014-12-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104239706B (zh) | 一种地面观测气温时空数据集的制备方法 | |
Lindberg et al. | Urban Multi-scale Environmental Predictor (UMEP): An integrated tool for city-based climate services | |
Agarwal et al. | Delineation of groundwater potential zone: An AHP/ANP approach | |
Camera et al. | Evaluation of interpolation techniques for the creation of gridded daily precipitation (1× 1 km2); Cyprus, 1980–2010 | |
Samanta et al. | Interpolation of climate variables and temperature modeling | |
Moiz et al. | A systematic decision support tool for robust hydropower site selection in poorly gauged basins | |
Long et al. | Reconstruction of historical arable land use patterns using constrained cellular automata: A case study of Jiangsu, China | |
Solaiman et al. | Development of probability based intensity-duration-frequency curves under climate change | |
Mortensen | Wind resource assessment using WAsP software | |
Kale et al. | Land-use and land-cover change in Western Ghats of India | |
CN110427594A (zh) | 适合小型实验室的气象要素数据获取系统 | |
CN101216481A (zh) | 一种反映区域土壤自然侵蚀程度的方法 | |
Chaaban et al. | Accuracy assessment of ESA WorldCover 2020 and ESRI 2020 land cover maps for a Region in Syria | |
Park et al. | Comparison of precipitation distributions in precipitation data sets representing 1km spatial resolution over South Korea produced by PRISM, IDW, and cokriging | |
Zhang et al. | A global dataset of daily near-surface air temperature at 1-km resolution (2003–2020) | |
CN102682335A (zh) | 精确确定区域对流层延迟的神经网络方法 | |
Bell et al. | The MaRIUS‐G2G datasets: Grid‐to‐Grid model estimates of flow and soil moisture for Great Britain using observed and climate model driving data | |
Pangali Sharma et al. | Assimilation of snowmelt runoff model (SRM) using satellite remote sensing data in Budhi Gandaki River Basin, Nepal | |
CN113902580A (zh) | 一种基于随机森林模型的历史耕地分布重建方法 | |
Guo et al. | How the variations of terrain factors affect the optimal interpolation methods for multiple types of climatic elements? | |
Faith et al. | The BioRap biodiversity assessment and planning study for Papua New Guinea | |
Khalil et al. | Time series analysis of land surface temperature and drivers of urban heat island effect based on remotely sensed data to develop a prediction model | |
Liu et al. | Estimation of precipitation induced by tropical cyclones based on machine‐learning‐enhanced analogue identification of numerical prediction | |
Li et al. | Evaluating the influencing factors of urbanization in the Xinjiang Uygur Autonomous Region over the past 27 years based on VIIRS-DNB and DMSP/OLS nightlight imageries | |
Mubiru et al. | Monthly average daily global solar irradiation maps for Uganda: A location in the equatorial region |
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 |
Effective date of registration: 20200817 Address after: 730000 No.318, Donggang West Road, Chengguan District, Lanzhou City, Gansu Province Patentee after: NORTHWEST INSTITUTE OF ECO-ENVIRONMENT AND RESOURCES, CAS Address before: Chengguan District of Gansu city of Lanzhou province Donggang West Road 730000 No. 320 Patentee before: Institute of environment and Engineering in cold and dry areas, Chinese Academy of Sciences |
|
TR01 | Transfer of patent right |