CN110502801A - 海洋温度锋自动追踪和特征参数信息提取方法 - Google Patents
海洋温度锋自动追踪和特征参数信息提取方法 Download PDFInfo
- Publication number
- CN110502801A CN110502801A CN201910675669.3A CN201910675669A CN110502801A CN 110502801 A CN110502801 A CN 110502801A CN 201910675669 A CN201910675669 A CN 201910675669A CN 110502801 A CN110502801 A CN 110502801A
- Authority
- CN
- China
- Prior art keywords
- point
- temperature
- gradient
- along
- fronts
- 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
- 238000000034 method Methods 0.000 title claims abstract description 36
- 239000000284 extract Substances 0.000 claims abstract description 11
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000000605 extraction Methods 0.000 claims description 8
- 230000000750 progressive effect Effects 0.000 claims description 5
- 238000002360 preparation method Methods 0.000 claims description 2
- 238000004458 analytical method Methods 0.000 abstract description 8
- 238000012512 characterization method Methods 0.000 abstract description 2
- 238000013480 data collection Methods 0.000 description 5
- 238000001514 detection method Methods 0.000 description 5
- 230000007613 environmental effect Effects 0.000 description 4
- 238000003012 network analysis Methods 0.000 description 4
- 238000012544 monitoring process Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013075 data extraction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000035558 fertility Effects 0.000 description 1
- 230000036039 immunity Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000013535 sea water Substances 0.000 description 1
- 210000003462 vein Anatomy 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C13/00—Surveying specially adapted to open water, e.g. sea, lake, river or canal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
- G01K13/02—Thermometers specially adapted for specific purposes for measuring temperature of moving fluids or granular materials capable of flow
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K13/00—Thermometers specially adapted for specific purposes
- G01K13/02—Thermometers specially adapted for specific purposes for measuring temperature of moving fluids or granular materials capable of flow
- G01K13/026—Thermometers specially adapted for specific purposes for measuring temperature of moving fluids or granular materials capable of flow of moving liquids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Hydrology & Water Resources (AREA)
- Radar, Positioning & Navigation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
- Radiation Pyrometers (AREA)
Abstract
本发明公开了一种海洋温度锋自动追踪和特征参数信息提取方法:(1)温度梯度场的建立;(2)温度锋横断面的确定;(3)温度锋特征参数的确定;(4)锋区中心位置的追踪;(5)遍历海区内所有网格点,得到海区内所有的海洋温度锋,并提取相应的温度锋特征参数信息,用于海洋中尺度现象研究并对未来海洋温度锋的分布做出可能性预测。本发明在温度水平梯度法原理的基础上,以网格遍历和物体质心权重分配的方式自动追踪海洋温度锋面中心位置,同时实现对海洋温度锋的强度、宽度、边界等特征参数信息的协同一致性提取,解决了传统的海洋温度锋提取方法仅能提供锋面的中心位置而不能同时给出锋区的强度、宽度、边界等关键特征参数信息的问题。
Description
技术领域
本发明涉及海洋信息技术,特别涉及一种基于锋面中心线的海洋温度锋自动 追踪和特征参数信息提取方法。
背景技术
海洋锋是海洋中水文要素水平剧烈变化的区域,可以用要素的水平梯度来描 述。锋区是海水辐聚区域,存在次级环流和强烈的混合作用,海洋生产力高,对 于海洋混合与海洋环流、海-气相互作用、渔业、水声传播等诸多领域的研究有 着非常重要的影响。开展海洋锋这一类海洋中尺度特征现象的研究、分析与预报 的首要环节就是对海洋锋区的特征参数,如锋区的中心位置、强度、宽度、边界 等进行信息提取。特别是在利用卫星遥感海表面温度开展海洋温度锋实时监测和 海洋三维温盐流数值分析和预报的过程中,人们需要基于海洋温度场网格分析和 预报数据集,对诸如海洋中尺度涡和温度锋这类所关心的海洋中尺度特征现象进 行快速解释和应用,以便可以高效地为海洋环境安全保障提供海洋特征现象信息 产品服务。
对于海洋锋面的检测方法主要可以归结为两类:一是梯度法(Gradientalgorithm),二是直方图法(Histogram algorithm)。前者主要是根据给定的温度 空间梯度的阈值来进行海洋锋的检测,适用于任意强度的锋面,而后者则可适用 于高背景噪音下的弱海洋锋检测。但无论是梯度法,还是直方图法,通常只是给 出锋面的中心位置,并不能同时给出海洋锋横断面的宽度,即锋面的侧向边界。 而在实际海洋环境安全保障中,往往需要同时给出海洋锋面的中心位置、强度、 横断面的宽度等信息。尤其是在海洋水声的应用中,海洋锋面边缘的确定,对水 声传播起着非常重要的作用。美国海军自2002年起在其业务化运行的涡可分辨 的全球海洋锋和中尺度涡特征现象预报***中,海洋温度锋面边界的预报是其发 布的非常重要的一项海洋环境安全保障产品。因此,在涡可分辨的全球海洋数值 预报***所产生的海洋温度场网格分析和预报数据集的基础上,研究开发一种精度较高的海洋温度锋面自动检测算法,尤其是可同时给出海洋锋面的中心位置、 强度、宽度和边界等特征参数,并且具有广泛适用性,从而达到对海洋温度锋这 类中尺度特征现象进行快速解释和应用的目的,势在必行。
发明内容
本发明针对海洋温度锋特征参数信息提取,主要目的是在利用卫星遥感海表 面温度开展海洋温度锋实时监测和海洋三维温盐流数值分析和预报的过程中,解 决通常所采用的海洋温度锋面的梯度检测方法仅能提供锋面的中心位置而不能 同时给出锋区的强度、宽度、边界等关键特征参数信息这一不足之处,提供一种 基于锋面中心线的海洋温度锋自动追踪方法,可同时提取锋面的中心位置、强度、 宽度、边界等特征参数信息,实现基于海洋温度场网格分析和预报数据集的中尺 度特征现象的快速解释和应用。
本发明所采用的技术方案是:一种海洋温度锋自动追踪和特征参数信息提取 方法,包括以下步骤:
步骤1,温度梯度场的建立:利用给定海域的温度场网格数据,计算温度梯 度场,并对温度梯度场进行内插,获得较温度场网格数据更高分辨率的温度梯度 场网格数据集;
步骤2,温度锋横断面的确定:建立温度锋判别标准,以某一满足温度锋判 别标准的点为温度锋中心位置初猜点,并以温度锋中心位置初猜点为中心,分别 沿温度梯度正负两个方向寻找温度锋临时边界点;
步骤3,温度锋特征参数的确定:采用逐步迭代的方式,确定锋区的中心位 置,提取特征参数信息;
步骤4,锋区中心位置的追踪:根据步骤3确定的锋区中心位置,沿着与锋 区中心位置的温度梯度方向的逆时针垂直方向,追踪下一个温度锋横断面的锋区 中心位置,提取相应特征参数信息;
步骤5,遍历海区内所有网格点,以搜索得到海区内所有的海洋温度锋,并 完成相应的温度锋特征参数信息的提取,用于海洋中尺度现象研究并对未来海洋 温度锋的分布做出可能性预测。
进一步地,步骤1中,所述的温度梯度场的计算方法为:
利用给定的温度场网格数据,分别计算网格中每个网格点X轴方向,Y轴 方向的温度梯度,其中,在计算X轴方向的温度梯度时,纬度不变,经度按照 固定步长变化;在计算Y轴方向的温度梯度时,经度不变,纬度按照固定步长 变化;
任意网格点处的温度梯度GF为:
式中,GX表示所求点沿着X轴方向的温度梯度值,GY表示所求点沿着Y轴方 向的温度梯度值;
其中,所求点沿着X轴方向的温度梯度值GX的计算方法如下:
当所求点为左边界点时,沿X轴方向,取所求点之后的两个点依次命名为 点l2和点l3;设所求点、点l2和点l3沿着X轴方向的温度值分别为Tlspx、Tl2x、 Tl3x,所求点与点l2在X轴方向上的距离为Llsp-l2x,点l2与点l3在X轴方向上 的距离为Ll2-l3x,则:
式中,W1spx为所求点为左边界点时、所求点沿着X轴方向占有的梯度权重、Wl2x为所求点为左边界点时、点l2沿着X轴方向占有的梯度权重、Wl3x为所求点为 左边界点时、点l3沿着X轴方向占有的梯度权重;
当所求点为右边界点时,沿X轴反方向,取所求点之前的两个点依次命名 为点r2和点r1;设点r1、点r2和所求点沿着X轴方向的温度值分别为Tr1x、Tr2x、 Trspx,点r1与点r2在X轴方向上的距离为Lr1-r2x,点r2与所求点在X轴方向上 的距离为Lr2-rspx,则:
式中,Wr1x为所求点为右边界点时、点r1沿着X轴方向占有的梯度权重、Wr2x为所求点为右边界点时、点r2沿着X轴方向占有的梯度权重、Wrspx为所求点为 右边界点时、所求点沿着X轴方向占有的梯度权重;
当所求点既不是左边界点也不是右边界点时,在X轴上,取与所求点左、 右相邻的两个点分别命名为点m1和点m3;设点m1、所求点和点m3沿着X轴 方向的温度值分别为Tm1x、Tmspx、Tm3x,点m1与所求点在X轴方向上的距离为 Lm1-mspx,所求点与点m3在X轴方向上的距离为Lmsp-m3x,则:
式中,Wm1x为所求点即不是左边界点也不是右边界点时、点m1沿着X轴方向占 有的梯度权重、Wmspx为所求点即不是左边界点也不是右边界点时、所求点沿着 X轴方向占有的梯度权重、Wm3x为所求点即不是左边界点也不是右边界点时、 点m3沿着X轴方向占有的梯度权重;
其中,所求点沿着Y轴方向的温度梯度值GY的计算方法如下:
当所求点为下边界点时,沿Y轴方向,取所求点之后的两个点依次命名为 点d2和点d3;设所求点、点d2和点d3沿着Y轴方向的温度值分别为Tdspy、Td2y、 Td3y,所求点与点d2在Y轴方向上的距离为Ldsp-d2y,点d2与点d3在Y轴方向 上的距离为Ld2-d3y,则:
式中,Wdspy为所求点为下边界点时、所求点沿着Y轴方向占有的梯度权重、Wd2y为所求点为下边界点时、点d2沿着Y轴方向占有的梯度权重、Wd3y为所求点为 下边界点时、点d3沿着Y轴方向占有的梯度权重;
当所求点为上边界点时,沿Y轴反方向,取所求点之前的两个点依次命名 为点h2和点h1;设点h1、点h2和所求点沿着Y轴方向的温度值分别为Th1y、Th2y、 Thspy,点h1与点h2在Y轴方向上的距离为Lh1-h2y,点h2与所求点在Y轴方向 上的距离为Lh2-hspy,则:
式中,Wh1y为所求点为上边界点时、点h1沿着Y轴方向占有的梯度权重、Wh2y为所求点为上边界点时、点h2沿着Y轴方向占有的梯度权重、Whspy为所求点为 上边界点时、所求点沿着Y轴方向占有的梯度权重;
当所求点既不是下边界点也不是上边界点时,在Y轴上,取与所求点下、 上相邻的两个点分别命名为点b1和点b3;设点b1、所求点和点b3沿着Y轴方 向的温度值分别为Tb1y、Tbspy、Tb3y,点b1与所求点在Y轴方向上的距离为Lb1-bspy, 所求点与点b3在Y轴方向上的距离为Lbsp-b3y,则:
式中,Wb1y为所求点即不是下边界点也不是上边界点时、点b1沿着Y轴方向 占有的梯度权重、Wbspy为所求点即不是下边界点也不是上边界点时、所求点沿 着Y轴方向占有的梯度权重、Wb3y为所求点即不是下边界点也不是上边界点时、 点b3沿着Y轴方向占有的梯度权重;
所述的较温度场网格数据更高分辨率的温度梯度场网格数据集的获得方法 为:
采取双线性插值的原理将各点梯度值插值到使用的较温度场网格更高分辨 率的网格上:
GF(x,y)=W11×GF(x1,y1)+W21×GF(x2,y1)+W12×GF(x1,y2)+W22×GF(x2,y2) (8)
式中,(x,y)表示较温度场网格数据更高分辨率的温度梯度场网格数据点对应坐标,GF(x,y)表示(x,y)处的温度梯度值;(x1,y1)、(x2,y1)、(x1,y2)、(x2,y2)分 别为(x,y)点所邻近的四个温度场网格数据点对应坐标;GF(x1,y1)、GF(x2,y1)、 GF(x1,y2)、GF(x2,y2)分别表示(x1,y1)、(x2,y1)、(x1,y2)、(x2,y2)处的温度 梯度值;W11、W21、W12、W22分别为(x1,y1)、(x2,y1)、(x1,y2)、(x2,y2)点的 温度梯度值所对应的权重,分别按如下公式计算:
进一步地,步骤2中,所述的建立温度锋判别标准,以某一满足温度锋判别 标准的点为温度锋中心位置初猜点,并以温度锋中心位置初猜点为中心,分别沿 温度梯度正负两个方向寻找温度锋临时边界点,具体为:
基于步骤1建立的较温度场网格数据更高分辨率的温度梯度场网格数据集, 建立温度锋判别标准,温度锋判别标准为:设定温度锋阈值,将梯度值高于设定 温度锋阈值的网格点视为满足温度锋判别标准;搜索满足温度锋判别标准的网格 点,并将满足温度锋判别标准的网格点设为锋区中心位置的初猜点;
以该锋区中心位置的初猜点为中心,分别沿着温度梯度的正负两个方向,以 固定步长搜索并判断这一横断面上点的温度梯度是否满足温度锋判别标准,将在 正负两个方向上搜索到的第一不满足温度锋判别标准的点暂定为该横断面的临 时边界点,将两个临时边界点之间、分别沿着初猜点温度梯度的正负两个方向、 以固定步长搜索到的满足温度锋判别标准的点构成了一条温度锋横断面。
进一步地,步骤3中,所述的采用逐步迭代的方式,确定锋区的中心位置, 提取特征参数信息,包括:
步骤3-1,对步骤2获得的温度锋横断面,按照计算物体质心的方法,以分 配权重的方式计算该温度锋横断面上的锋区中心以及平均强度;
步骤3-2,如果计算得到的锋区中心并非步骤2中的温度锋中心位置初猜点, 则以计算得到的锋区中心为新的初猜点,并以新的初猜点为中心,按照步骤2 所示的方法重新设定温度锋横断面并按照步骤3-1计算获得重新设定的温度锋横 断面上的锋区中心以及平均强度,判断重新设定的温度锋横断面上的锋区中心是 否与新的初猜点重合,直至计算得到的锋区中心与初猜点重合,获得最终的锋区 中心;
步骤3-3,记录最终的锋区中心位置、平均强度、边界点位置;根据两边界 点的位置信息,计算得到锋区宽度:
式中,FWD表示锋区宽度,(xL,yL)和(xR,yR)分别表示锋区左右两个边界点的 经纬度坐标,R表示地球半径。
进一步地,步骤4中,所述的根据步骤3确定的锋区中心位置,沿着与锋区 中心位置的温度梯度方向的逆时针垂直方向,追踪下一个温度锋横断面的锋区中 心位置,提取相应特征参数信息,具体为:
以与锋区中心位置的温度梯度方向的逆时针垂直方向为方向搜索一个设定 步长,根据步骤2获得该条锋上的下一个初猜点,继续进行步骤2,直至触到岸 界或温度梯度值小于判别标准中所设定的温度锋阈值,获得下一个温度锋横断 面;根据步骤3获得下一个温度锋横断面的锋区中心位置,提取该锋区中心位置 的特征参数信息。
本发明的有益效果是:一种基于锋面中心线的海洋温度锋自动追踪方法,在 温度水平梯度法原理的基础上,以网格遍历和物体质心权重分配的方式自动追踪 海洋温度锋面中心位置,同时实现对海洋温度锋的强度、宽度、边界等特征参数 信息的协同一致性提取,从而克服传统的海洋温度锋提取方法仅能提供锋面中心 位置这一缺陷,尤为重要的是可以提供对海洋水声传播有着重要影响的锋面边界 信息,是海洋环境安全保障中不可或缺的一项产品。同时,该方法一脉秉承了传 统的海洋温度锋提取梯度算法的优势,即适用性强,可以对任意海区、任意强度 的海洋温度锋面特征信息提取,并具有较强的抗噪性。此外,该方法既可以用于 卫星遥感海表面温度监测海洋锋,也可以用于海洋三维温盐流数值分析和预报中 海洋中尺度特征现象的快速解释与应用。
附图说明
图1为本发明海洋温度锋自动追踪和特征参数信息提取方法的流程示意图。
图2为基于2012年1月1日卫星遥感海表面温度数据提取的南海北部海域 温度锋分布图。
图3为基于1993年4月30日海洋温度再分析数据提取的日本以南海域水下 200m黑潮温度锋分布图。
具体实施方式
为能进一步了解本发明的发明内容、特点及功效,兹例举以下实施例,并配 合附图详细说明如下:
本发明为一种基于锋面中心线、以网格遍历和物体质心权重分配的方式自动 追踪海洋温度锋面中心位置,并同时提取温度锋区的强度、宽度、边界等特征参 数信息的方法。该方法主要是应用于在开展海洋温度锋这一类海洋中尺度特征现 象的分析和预报时,基于海洋温度场网格分析和预报数据集,提供一种精度较高 且适用性广的海洋温度锋解释分析应用方法,并可同时适用于海洋盐度锋、密度 锋和声速锋等锋的特征参数信息的提取。
如附图1所示,一种海洋温度锋自动追踪和特征参数信息提取方法,包括以 下步骤:
步骤1,温度梯度场的建立
对于拟提取海洋温度锋的海域,利用给定海域的温度场网格数据,建立温度 梯度场,包括:计算温度梯度场,并将其插值到更细分辨率网格上,获得较温度 场网格数据更高分辨率的温度梯度场网格数据集。
第一步:温度梯度场的计算
利用给定的温度场网格数据,分别计算网格中每个网格点X轴方向,Y轴 方向的温度梯度,其中,在计算X轴方向的温度梯度时,纬度不变,经度按照 固定步长变化;在计算Y轴方向的温度梯度时,经度不变,纬度按照固定步长 变化。
任意网格点处的温度梯度GF为:
式中,GX表示所求点沿着X轴方向的温度梯度值,GY表示所求点沿着Y轴方 向的温度梯度值。
其中,所求点沿着X轴方向的温度梯度值GX的计算方法如下:
当所求点为左边界点时,沿X轴方向,取所求点之后的两个点依次命名为 点l2和点l3;设所求点、点l2和点l3沿着X轴方向的温度值分别为Tlspx、Tl2x、 Tl3x,所求点与点l2在X轴方向上的距离为Llsp-l2x,点l2与点l3在X轴方向上 的距离为Ll2-l3x,则:
式中,W1spx为所求点为左边界点时、所求点沿着X轴方向占有的梯度权重、Wl2x为所求点为左边界点时、点l2沿着X轴方向占有的梯度权重、Wl3x为所求点为 左边界点时、点l3沿着X轴方向占有的梯度权重;
当所求点为右边界点时,沿X轴反方向,取所求点之前的两个点依次命名 为点r2和点r1;设点r1、点r2和所求点沿着X轴方向的温度值分别为Tr1x、Tr2x、 Trspx,点r1与点r2在X轴方向上的距离为Lr1-r2x,点r2与所求点在X轴方向上 的距离为Lr2-rspx,则:
式中,Wr1x为所求点为右边界点时、点r1沿着X轴方向占有的梯度权重、Wr2x为所求点为右边界点时、点r2沿着X轴方向占有的梯度权重、Wrspx为所求点为 右边界点时、所求点沿着X轴方向占有的梯度权重;
当所求点既不是左边界点也不是右边界点时,在X轴上,取与所求点左、 右相邻的两个点分别命名为点m1和点m3;设点m1、所求点和点m3沿着X轴 方向的温度值分别为Tm1x、Tmspx、Tm3x,点m1与所求点在X轴方向上的距离为 Lm1-mspx,所求点与点m3在X轴方向上的距离为Lmsp-m3x,则:
式中,Wm1x为所求点即不是左边界点也不是右边界点时、点m1沿着X轴方向占 有的梯度权重、Wmspx为所求点即不是左边界点也不是右边界点时、所求点沿着 X轴方向占有的梯度权重、Wm3x为所求点即不是左边界点也不是右边界点时、 点m3沿着X轴方向占有的梯度权重。
其中,所求点沿着Y轴方向的温度梯度值GY的计算方法如下:
当所求点为下边界点时,沿Y轴方向,取所求点之后的两个点依次命名为 点d2和点d3;设所求点、点d2和点d3沿着Y轴方向的温度值分别为Tdspy、Td2y、 Td3y,所求点与点d2在Y轴方向上的距离为Ldsp-d2y,点d2与点d3在Y轴方向 上的距离为Ld2-d3y,则:
式中,Wdspy为所求点为下边界点时、所求点沿着Y轴方向占有的梯度权重、Wd2y为所求点为下边界点时、点d2沿着Y轴方向占有的梯度权重、Wd3y为所求点为 下边界点时、点d3沿着Y轴方向占有的梯度权重;
当所求点为上边界点时,沿Y轴反方向,取所求点之前的两个点依次命名 为点h2和点h1;设点h1、点h2和所求点沿着Y轴方向的温度值分别为Th1y、Th2y、 Thspy,点h1与点h2在Y轴方向上的距离为Lh1-h2y,点h2与所求点在Y轴方向 上的距离为Lh2-hspy,则:
式中,Wh1y为所求点为上边界点时、点h1沿着Y轴方向占有的梯度权重、Wh2y为所求点为上边界点时、点h2沿着Y轴方向占有的梯度权重、Whspy为所求点为 上边界点时、所求点沿着Y轴方向占有的梯度权重;
当所求点既不是下边界点也不是上边界点时,在Y轴上,取与所求点下、 上相邻的两个点分别命名为点b1和点b3;设点b1、所求点和点b3沿着Y轴方 向的温度值分别为Tb1y、Tbspy、Tb3y,点b1与所求点在Y轴方向上的距离为Lb1-bspy, 所求点与点b3在Y轴方向上的距离为Lbsp-b3y,则:
式中,Wb1y为所求点即不是下边界点也不是上边界点时、点b1沿着Y轴方向 占有的梯度权重、Wbspy为所求点即不是下边界点也不是上边界点时、所求点沿 着Y轴方向占有的梯度权重、Wb3y为所求点即不是下边界点也不是上边界点时、 点b3沿着Y轴方向占有的梯度权重。
第二步:较高分辨率的温度梯度场网格数据集的获得
为尽量减少后续温度锋搜索过程中的噪音,将上述温度梯度场进行内插,以 获得较高分辨率温度梯度场网格数据集。
采取双线性插值的原理将各点梯度值插值到使用的较温度场网格更高分辨 率的网格上:
GF(x,y)=W11×GF(x1,y1)+W21×GF(x2,y1)+W12×GF(x1,y2)+W22×GF(x2,y2) (8)
式中,(x,y)表示较温度场网格数据更高分辨率的温度梯度场网格数据点对应坐标,GF(x,y)表示(x,y)处的温度梯度值;(x1,y1)、(x2,y1)、(x1,y2)、(x2,y2)分 别为(x,y)点所邻近的四个温度场网格数据点对应坐标;GF(x1,y1)、GF(x2,y1)、 GF(x1,y2)、GF(x2,y2)分别表示(x1,y1)、(x2,y1)、(x1,y2)、(x2,y2)处的温度 梯度值;W11、W21、W12、W22分别为(x1,y1)、(x2,y1)、(x1,y2)、(x2,y2)点的 温度梯度值所对应的权重,分别按如下公式计算:
步骤2,温度锋横断面的确定
建立温度锋判别标准,以某一满足温度锋判别标准的点为温度锋中心位置初 猜点,并以温度锋中心位置初猜点为中心,分别沿温度梯度正负两个方向寻找温 度锋临时边界点,具体为:
基于步骤1建立的较温度场网格数据更高分辨率的温度梯度场网格数据集, 建立温度锋判别标准,温度锋判别标准为:设定温度锋阈值,将梯度值高于设定 温度锋阈值的网格点视为满足温度锋判别标准;搜索满足温度锋判别标准的网格 点,并将满足温度锋判别标准的网格点设为锋区中心位置的初猜点;
以该锋区中心位置的初猜点为中心,分别沿着温度梯度的正负两个方向,以 固定步长搜索并判断这一横断面上点的温度梯度是否满足温度锋判别标准,将在 正负两个方向上搜索到的第一不满足温度锋判别标准的点暂定为该横断面的临 时边界点,将两个临时边界点之间、分别沿着初猜点温度梯度的正负两个方向、 以固定步长搜索到的满足温度锋判别标准的点构成了一条温度锋横断面。
步骤3,温度锋特征参数的确定
采用逐步迭代的方式,确定锋区的中心位置,提取温度锋平均强度、宽度和 边界位置等特征参数信息,具体包括:
步骤3-1,对步骤2获得的温度锋横断面,按照计算物体质心的方法,以分 配权重的方式(即梯度绝对值大的点被分配以较大的权重)计算该温度锋横断面 上的锋区中心以及平均强度;
步骤3-2,如果计算得到的锋区中心并非步骤2中的温度锋中心位置初猜点, 则以计算得到的锋区中心为新的初猜点,并以新的初猜点为中心,按照步骤2 所示的方法重新设定温度锋横断面并按照步骤3-1计算获得重新设定的温度锋横 断面上的锋区中心以及平均强度,判断重新设定的温度锋横断面上的锋区中心是 否与新的初猜点重合,直至计算得到的锋区中心与初猜点重合,获得最终的锋区 中心;
步骤3-3,记录最终的锋区中心位置、平均强度、边界点位置;根据两边界 点的位置信息,计算得到锋区宽度:
式中,FWD表示锋区宽度,(xL,yL)和(xR,yR)分别表示锋区左右两个边界点的 经纬度坐标,R表示地球半径。
步骤4,锋区中心位置的追踪
根据步骤3确定的锋区中心位置,以与锋区中心位置的温度梯度方向的逆时 针垂直方向(即将锋区梯度方向逆时针旋转90°,保证温度高值位于右侧)为方 向搜索一个设定步长,根据步骤2获得该条锋上的下一个初猜点,继续进行步骤 2,直至触到岸界或温度梯度值小于判别标准中所设定的温度锋阈值,获得下一 个温度锋横断面;然后,根据步骤3获得下一个温度锋横断面的锋区中心位置, 提取该锋区中心位置的特征参数信息。
步骤5,重复步骤2至步骤4,遍历海区内所有网格点,以搜索得到海区内 所有的海洋温度锋,并完成相应的温度锋特征参数信息的提取,用于海洋中尺度 现象研究(例如中尺度涡),此外也可对未来海洋温度锋的分布做出可能性预测。
图2为基于2012年1月1日卫星遥感海表面温度数据提取的南海北部海域 温度锋分布图,图3为基于1993年4月30日海洋温度再分析数据提取的日本以 南海域水下200m黑潮温度锋分布图。图中黑细实线为温度等值线(单位:℃), 黑粗实线为温度锋的边界,黑粗虚线为温度锋区的中心位置。由图可以看出,采 用上述温度锋提取方法可以较好地将温度锋区的中心位置,以及温度锋的宽度、 边界等特征参数信息提取出来,从而保证了这些信息的协同一致性。同时也表明, 上述温度锋提取方法具有较强的适用性:即对任意海区、任意强度的海洋温度锋 特征参数信息提取,既适用于卫星遥感海表面温度数据,也既适用于海洋三维温 盐流数值分析和预报数据集。
尽管上面结合附图对本发明的优选实施例进行了描述,但是本发明并不局限 于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,并不是限制性的, 本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保 护的范围情况下,还可以做出很多形式,这些均属于本发明的保护范围之内。
Claims (5)
1.一种海洋温度锋自动追踪和特征参数信息提取方法,其特征在于,包括以下步骤:
步骤1,温度梯度场的建立:利用给定海域的温度场网格数据,计算温度梯度场,并对温度梯度场进行内插,获得较温度场网格数据更高分辨率的温度梯度场网格数据集;
步骤2,温度锋横断面的确定:建立温度锋判别标准,以某一满足温度锋判别标准的点为温度锋中心位置初猜点,并以温度锋中心位置初猜点为中心,分别沿温度梯度正负两个方向寻找温度锋临时边界点;
步骤3,温度锋特征参数的确定:采用逐步迭代的方式,确定锋区的中心位置,提取特征参数信息;
步骤4,锋区中心位置的追踪:根据步骤3确定的锋区中心位置,沿着与锋区中心位置的温度梯度方向的逆时针垂直方向,追踪下一个温度锋横断面的锋区中心位置,提取相应特征参数信息;
步骤5,遍历海区内所有网格点,以搜索得到海区内所有的海洋温度锋,并完成相应的温度锋特征参数信息的提取,用于海洋中尺度现象研究并对未来海洋温度锋的分布做出可能性预测。
2.根据权利要求1所述的海洋温度锋自动追踪和特征参数信息提取方法,其特征在于,步骤1中,所述的温度梯度场的计算方法为:
利用给定的温度场网格数据,分别计算网格中每个网格点X轴方向,Y轴方向的温度梯度,其中,在计算X轴方向的温度梯度时,纬度不变,经度按照固定步长变化;在计算Y轴方向的温度梯度时,经度不变,纬度按照固定步长变化;
任意网格点处的温度梯度GF为:
式中,GX表示所求点沿着X轴方向的温度梯度值,GY表示所求点沿着Y轴方向的温度梯度值;
其中,所求点沿着X轴方向的温度梯度值GX的计算方法如下:
当所求点为左边界点时,沿X轴方向,取所求点之后的两个点依次命名为点l2和点l3;设所求点、点l2和点l3沿着X轴方向的温度值分别为Tlspx、Tl2x、Tl3x,所求点与点l2在X轴方向上的距离为Llsp-l2x,点l2与点l3在X轴方向上的距离为Ll2-l3x,则:
式中,W1spx为所求点为左边界点时、所求点沿着X轴方向占有的梯度权重、Wl2x为所求点为左边界点时、点l2沿着X轴方向占有的梯度权重、Wl3x为所求点为左边界点时、点l3沿着X轴方向占有的梯度权重;
当所求点为右边界点时,沿X轴反方向,取所求点之前的两个点依次命名为点r2和点r1;设点r1、点r2和所求点沿着X轴方向的温度值分别为Tr1x、Tr2x、Trspx,点r1与点r2在X轴方向上的距离为Lr1-r2x,点r2与所求点在X轴方向上的距离为Lr2-rspx,则:
式中,Wr1x为所求点为右边界点时、点r1沿着X轴方向占有的梯度权重、Wr2x为所求点为右边界点时、点r2沿着X轴方向占有的梯度权重、Wrspx为所求点为右边界点时、所求点沿着X轴方向占有的梯度权重;
当所求点既不是左边界点也不是右边界点时,在X轴上,取与所求点左、右相邻的两个点分别命名为点m1和点m3;设点m1、所求点和点m3沿着X轴方向的温度值分别为Tm1x、Tmspx、Tm3x,点m1与所求点在X轴方向上的距离为Lm1-mspx,所求点与点m3在X轴方向上的距离为Lmsp-m3x,则:
式中,Wm1x为所求点即不是左边界点也不是右边界点时、点m1沿着X轴方向占有的梯度权重、Wmspx为所求点即不是左边界点也不是右边界点时、所求点沿着X轴方向占有的梯度权重、Wm3x为所求点即不是左边界点也不是右边界点时、点m3沿着X轴方向占有的梯度权重;
其中,所求点沿着Y轴方向的温度梯度值GY的计算方法如下:
当所求点为下边界点时,沿Y轴方向,取所求点之后的两个点依次命名为点d2和点d3;设所求点、点d2和点d3沿着Y轴方向的温度值分别为Tdspy、Td2y、Td3y,所求点与点d2在Y轴方向上的距离为Ldsp-d2y,点d2与点d3在Y轴方向上的距离为Ld2-d3y,则:
式中,Wdspy为所求点为下边界点时、所求点沿着Y轴方向占有的梯度权重、Wd2y为所求点为下边界点时、点d2沿着Y轴方向占有的梯度权重、Wd3y为所求点为下边界点时、点d3沿着Y轴方向占有的梯度权重;
当所求点为上边界点时,沿Y轴反方向,取所求点之前的两个点依次命名为点h2和点h1;设点h1、点h2和所求点沿着Y轴方向的温度值分别为Th1y、Th2y、Thspy,点h1与点h2在Y轴方向上的距离为Lh1-h2y,点h2与所求点在Y轴方向上的距离为Lh2-hspy,则:
式中,Wh1y为所求点为上边界点时、点h1沿着Y轴方向占有的梯度权重、Wh2y为所求点为上边界点时、点h2沿着Y轴方向占有的梯度权重、Whspy为所求点为上边界点时、所求点沿着Y轴方向占有的梯度权重;
当所求点既不是下边界点也不是上边界点时,在Y轴上,取与所求点下、上相邻的两个点分别命名为点b1和点b3;设点b1、所求点和点b3沿着Y轴方向的温度值分别为Tb1y、Tbspy、Tb3y,点b1与所求点在Y轴方向上的距离为Lb1-bspy,所求点与点b3在Y轴方向上的距离为Lbsp-b3y,则:
式中,Wb1y为所求点即不是下边界点也不是上边界点时、点b1沿着Y轴方向占有的梯度权重、Wbspy为所求点即不是下边界点也不是上边界点时、所求点沿着Y轴方向占有的梯度权重、Wb3y为所求点即不是下边界点也不是上边界点时、点b3沿着Y轴方向占有的梯度权重;
所述的较温度场网格数据更高分辨率的温度梯度场网格数据集的获得方法为:
采取双线性插值的原理将各点梯度值插值到使用的较温度场网格更高分辨率的网格上:
GF(x,y)=W11×GF(x1,y1)+W21×GF(x2,y1)+W12×GF(x1,y2)+W22×GF(x2,y2) (8)
式中,(x,y)表示较温度场网格数据更高分辨率的温度梯度场网格数据点对应坐标,GF(x,y)表示(x,y)处的温度梯度值;(x1,y1)、(x2,y1)、(x1,y2)、(x2,y2)分别为(x,y)点所邻近的四个温度场网格数据点对应坐标;GF(x1,y1)、GF(x2,y1)、GF(x1,y2)、GF(x2,y2)分别表示(x1,y1)、(x2,y1)、(x1,y2)、(x2,y2)处的温度梯度值;W11、W21、W12、W22分别为(x1,y1)、(x2,y1)、(x1,y2)、(x2,y2)点的温度梯度值所对应的权重,分别按如下公式计算:
3.根据权利要求1所述的海洋温度锋自动追踪和特征参数信息提取方法,其特征在于,步骤2中,所述的建立温度锋判别标准,以某一满足温度锋判别标准的点为温度锋中心位置初猜点,并以温度锋中心位置初猜点为中心,分别沿温度梯度正负两个方向寻找温度锋临时边界点,具体为:
基于步骤1建立的较温度场网格数据更高分辨率的温度梯度场网格数据集,建立温度锋判别标准,温度锋判别标准为:设定温度锋阈值,将梯度值高于设定温度锋阈值的网格点视为满足温度锋判别标准;搜索满足温度锋判别标准的网格点,并将满足温度锋判别标准的网格点设为锋区中心位置的初猜点;
以该锋区中心位置的初猜点为中心,分别沿着温度梯度的正负两个方向,以固定步长搜索并判断这一横断面上点的温度梯度是否满足温度锋判别标准,将在正负两个方向上搜索到的第一不满足温度锋判别标准的点暂定为该横断面的临时边界点,将两个临时边界点之间、分别沿着初猜点温度梯度的正负两个方向、以固定步长搜索到的满足温度锋判别标准的点构成了一条温度锋横断面。
4.根据权利要求1所述的海洋温度锋自动追踪和特征参数信息提取方法,其特征在于,步骤3中,所述的采用逐步迭代的方式,确定锋区的中心位置,提取特征参数信息,包括:
步骤3-1,对步骤2获得的温度锋横断面,按照计算物体质心的方法,以分配权重的方式计算该温度锋横断面上的锋区中心以及平均强度;
步骤3-2,如果计算得到的锋区中心并非步骤2中的温度锋中心位置初猜点,则以计算得到的锋区中心为新的初猜点,并以新的初猜点为中心,按照步骤2所示的方法重新设定温度锋横断面并按照步骤3-1计算获得重新设定的温度锋横断面上的锋区中心以及平均强度,判断重新设定的温度锋横断面上的锋区中心是否与新的初猜点重合,直至计算得到的锋区中心与初猜点重合,获得最终的锋区中心;
步骤3-3,记录最终的锋区中心位置、平均强度、边界点位置;根据两边界点的位置信息,计算得到锋区宽度:
式中,FWD表示锋区宽度,(xL,yL)和(xR,yR)分别表示锋区左右两个边界点的经纬度坐标,R表示地球半径。
5.根据权利要求1所述的海洋温度锋自动追踪和特征参数信息提取方法,其特征在于,步骤4中,所述的根据步骤3确定的锋区中心位置,沿着与锋区中心位置的温度梯度方向的逆时针垂直方向,追踪下一个温度锋横断面的锋区中心位置,提取相应特征参数信息,具体为:
以与锋区中心位置的温度梯度方向的逆时针垂直方向为方向搜索一个设定步长,根据步骤2获得该条锋上的下一个初猜点,继续进行步骤2,直至触到岸界或温度梯度值小于判别标准中所设定的温度锋阈值,获得下一个温度锋横断面;根据步骤3获得下一个温度锋横断面的锋区中心位置,提取该锋区中心位置的特征参数信息。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910675669.3A CN110502801B (zh) | 2019-07-25 | 2019-07-25 | 海洋温度锋自动追踪和特征参数信息提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910675669.3A CN110502801B (zh) | 2019-07-25 | 2019-07-25 | 海洋温度锋自动追踪和特征参数信息提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110502801A true CN110502801A (zh) | 2019-11-26 |
CN110502801B CN110502801B (zh) | 2023-03-24 |
Family
ID=68586929
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910675669.3A Expired - Fee Related CN110502801B (zh) | 2019-07-25 | 2019-07-25 | 海洋温度锋自动追踪和特征参数信息提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110502801B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110908404A (zh) * | 2019-12-09 | 2020-03-24 | 自然资源部第二海洋研究所 | 一种基于数据驱动的auv智能观测运动方法 |
CN113063906A (zh) * | 2021-02-09 | 2021-07-02 | 国家***南海预报中心(国家***广州海洋预报台) | 一种叶绿素a锋面的检测方法及装置 |
CN114372528A (zh) * | 2022-01-10 | 2022-04-19 | 中国海洋大学 | 一种黑潮断流现象的识别方法及装置 |
CN111860146B (zh) * | 2020-06-11 | 2023-06-09 | 中山大学 | 海洋锋区域的获取方法、装置、计算机设备和存储介质 |
CN116577844A (zh) * | 2023-03-28 | 2023-08-11 | 南京信息工程大学 | 一种东亚冷锋降水自动识别方法、*** |
CN117213448A (zh) * | 2023-11-07 | 2023-12-12 | 中国人民解放军国防科技大学 | 一种海洋次中尺度锋面调查方法 |
CN117851845A (zh) * | 2024-03-07 | 2024-04-09 | 中国海洋大学 | 一种基于聚类算法的海洋亚中尺度锋面提取方法 |
CN117975296A (zh) * | 2024-04-01 | 2024-05-03 | 自然资源部第一海洋研究所 | 一种异常海洋锋面的卫星遥感探测提取识别方法及*** |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110075903A1 (en) * | 2007-10-26 | 2011-03-31 | Consejo Superior De Investigaciones Cientificas | Method and system for the analysis of peculiarities in fingerprints |
CN103020959A (zh) * | 2012-11-24 | 2013-04-03 | 中国科学院地理科学与资源研究所 | 一种基于引力模型的海洋锋信息提取方法 |
CN107274420A (zh) * | 2017-06-15 | 2017-10-20 | 中国水产科学研究院东海水产研究所 | 基于图像分割的海洋锋面提取方法 |
CN109086818A (zh) * | 2018-07-25 | 2018-12-25 | 中国海洋大学 | 海洋锋面识别方法和装置 |
CN109668635A (zh) * | 2019-01-16 | 2019-04-23 | 中国人民解放军61741部队 | 海表温度融合方法与*** |
-
2019
- 2019-07-25 CN CN201910675669.3A patent/CN110502801B/zh not_active Expired - Fee Related
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110075903A1 (en) * | 2007-10-26 | 2011-03-31 | Consejo Superior De Investigaciones Cientificas | Method and system for the analysis of peculiarities in fingerprints |
CN103020959A (zh) * | 2012-11-24 | 2013-04-03 | 中国科学院地理科学与资源研究所 | 一种基于引力模型的海洋锋信息提取方法 |
CN107274420A (zh) * | 2017-06-15 | 2017-10-20 | 中国水产科学研究院东海水产研究所 | 基于图像分割的海洋锋面提取方法 |
CN109086818A (zh) * | 2018-07-25 | 2018-12-25 | 中国海洋大学 | 海洋锋面识别方法和装置 |
CN109668635A (zh) * | 2019-01-16 | 2019-04-23 | 中国人民解放军61741部队 | 海表温度融合方法与*** |
Non-Patent Citations (3)
Title |
---|
NEETU ET AL.: ""Data-Adaptive Prediction of Sea-Surface Temperature in the Arabian Sea"", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
李威: ""台湾以东黑潮锋的数值研究"", 《中国博士学位论文全文数据库 基础科学辑》 * |
王磊 等: ""利用卫星遥感资料对南海北部陆架海洋表层温度锋的分析"", 《中国海洋大学学报》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110908404B (zh) * | 2019-12-09 | 2023-04-18 | 自然资源部第二海洋研究所 | 一种基于数据驱动的auv智能观测运动方法 |
CN110908404A (zh) * | 2019-12-09 | 2020-03-24 | 自然资源部第二海洋研究所 | 一种基于数据驱动的auv智能观测运动方法 |
CN111860146B (zh) * | 2020-06-11 | 2023-06-09 | 中山大学 | 海洋锋区域的获取方法、装置、计算机设备和存储介质 |
CN113063906A (zh) * | 2021-02-09 | 2021-07-02 | 国家***南海预报中心(国家***广州海洋预报台) | 一种叶绿素a锋面的检测方法及装置 |
CN114372528A (zh) * | 2022-01-10 | 2022-04-19 | 中国海洋大学 | 一种黑潮断流现象的识别方法及装置 |
CN116577844B (zh) * | 2023-03-28 | 2024-02-09 | 南京信息工程大学 | 一种东亚冷锋降水自动识别方法、*** |
CN116577844A (zh) * | 2023-03-28 | 2023-08-11 | 南京信息工程大学 | 一种东亚冷锋降水自动识别方法、*** |
CN117213448A (zh) * | 2023-11-07 | 2023-12-12 | 中国人民解放军国防科技大学 | 一种海洋次中尺度锋面调查方法 |
CN117213448B (zh) * | 2023-11-07 | 2024-01-30 | 中国人民解放军国防科技大学 | 一种海洋次中尺度锋面调查方法 |
CN117851845A (zh) * | 2024-03-07 | 2024-04-09 | 中国海洋大学 | 一种基于聚类算法的海洋亚中尺度锋面提取方法 |
CN117851845B (zh) * | 2024-03-07 | 2024-05-17 | 中国海洋大学 | 一种基于聚类算法的海洋亚中尺度锋面提取方法 |
CN117975296A (zh) * | 2024-04-01 | 2024-05-03 | 自然资源部第一海洋研究所 | 一种异常海洋锋面的卫星遥感探测提取识别方法及*** |
CN117975296B (zh) * | 2024-04-01 | 2024-06-11 | 自然资源部第一海洋研究所 | 一种异常海洋锋面的卫星遥感探测提取识别方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN110502801B (zh) | 2023-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110502801A (zh) | 海洋温度锋自动追踪和特征参数信息提取方法 | |
Piola et al. | Circumpolar properties of Antarctic intermediate water and Subantarctic Mode Water | |
Pareta et al. | Quantitative geomorphological analysis of a watershed of Ravi River Basin, HP India | |
Wu et al. | A localized contour tree method for deriving geometric and topological properties of complex surface depressions based on high-resolution topographical data | |
Raj et al. | Quantifying mesoscale eddies in the Lofoten Basin | |
CN106022259B (zh) | 一种基于激光点云三维特征描述模型的山区道路提取方法 | |
Tsuji et al. | Field survey of the East Java earthquake and tsunami of June 3, 1994 | |
Aguiar et al. | A census of Meddies in a long-term high-resolution simulation | |
Barton et al. | Canary and Portugal currents | |
CN109145376A (zh) | 一种基于cad地形图的水库淹没范围确定方法 | |
Wagensommer et al. | New dinosaur tracksites from the Middle Jurassic of Madagascar: ichnotaxonomical, behavioural and palaeoenvironmental implications | |
CN106021290A (zh) | 一种基于多尺度地理信息的社交网络关联挖掘方法 | |
Sharma et al. | Use of geographical information system in hypsometric analysis of Kanhiya Nala watershed | |
CN110516856A (zh) | 基于卷积神经网络的估算海洋次表层温度的方法 | |
CN115358311A (zh) | 地表变形监测多源数据融合处理方法 | |
Rasmussen et al. | The Faroe shelf circulation and its potential impact on the primary production | |
Mustaffa et al. | A study of flood evacuation center using GIS and remote sensing technique | |
CN111523223A (zh) | 超长重现期极值风速的推算方法 | |
Geyh et al. | Groundwater analysis of environmental carbon and other isotopes from the Jakarta basin aquifer, Indonesia | |
Robertson et al. | Calving retreat and proglacial lake growth at H ooker G lacier, S outhern A lps, N ew Z ealand | |
Qi et al. | Topographical mapping of a bare tidal flat outside a mangrove area based on the waterline method and an iterative hydrodynamic model: a case study of Yingluo Bay, South China | |
Romshoo et al. | Glacier thickness and volume estimation in the Upper Indus Basin using modeling and ground penetrating radar measurements | |
Telling et al. | Levee Vertical Land Motion Changes in the Sacramento-San Joaquin Delta | |
Rodgers | Deposition and fracture patterns within the Herkimer quadrangle, NY: implications for the location and character of the Herkimer fault | |
Grant | Workshop on Sea-Level Rise and Subsidence in the Maritimes |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20230324 |
|
CF01 | Termination of patent right due to non-payment of annual fee |