CN113701711B - 基于北斗定位与气压计的高精度定位方法及系统 - Google Patents
基于北斗定位与气压计的高精度定位方法及系统 Download PDFInfo
- Publication number
- CN113701711B CN113701711B CN202111027695.9A CN202111027695A CN113701711B CN 113701711 B CN113701711 B CN 113701711B CN 202111027695 A CN202111027695 A CN 202111027695A CN 113701711 B CN113701711 B CN 113701711B
- Authority
- CN
- China
- Prior art keywords
- time
- information sequence
- weather information
- air pressure
- barometer
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000004140 cleaning Methods 0.000 claims description 11
- 230000005855 radiation Effects 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 10
- 239000004973 liquid crystal related substance Substances 0.000 description 4
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- DMBHHRLKUKUOEG-UHFFFAOYSA-N diphenylamine Chemical compound C=1C=CC=CC=1NC1=CC=CC=C1 DMBHHRLKUKUOEG-UHFFFAOYSA-N 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C5/00—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels
- G01C5/06—Measuring height; Measuring distances transverse to line of sight; Levelling between separated points; Surveyors' levels by using barometric means
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/38—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system
- G01S19/39—Determining a navigation solution using signals transmitted by a satellite radio beacon positioning system the satellite radio beacon positioning system transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/42—Determining position
- G01S19/45—Determining position by combining measurements of signals from the satellite radio beacon positioning system with a supplementary measurement
-
- 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)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及坐标定位技术领域,具体地说,涉及一种基于北斗定位与气压计的高精度定位方法及系统。该方法包括如下步骤:步骤S1、基于北斗定位模块获取当前坐标点的平面坐标(x,y);步骤S2,基于气压计获取当前坐标点的高度坐标z;步骤S3,以(x,y,z)作为当前坐标点的坐标并输出。该系统用于实现上述方法。本发明能够较佳地提升气压计的测量精度。
Description
技术领域
本发明涉及坐标定位技术领域,具体地说,涉及一种基于北斗定位与气压计的高精度定位方法及系统。
背景技术
北斗定位的平面方向上坐标精度很高,但高度方向上坐标精度很差,气压计虽然可以测量高度,但是容易受当地气压日较差的影响,导致一天24h内高度存在很大的漂移,导致高度测量不准,所有导致两者结合起来的实际使用效果不好。
发明内容
本发明提供了一种基于北斗定位与气压计的高精度定位方法,其能够克服现有技术的某种或某些缺陷。
根据本发明的基于北斗定位与气压计的高精度定位方法,其包括如下步骤:
步骤S1、基于北斗定位模块获取当前坐标点的平面坐标(x,y);
步骤S2,基于气压计获取当前坐标点的高度坐标z;
步骤S3,以(x,y,z)作为当前坐标点的坐标并输出。
通过本发明中的方法,能够较佳地以北斗定位模块提供的平面坐标及气压计获取的高度坐标作为当前坐标点并输出,故而能够具备较佳的测量精度。
作为优选,步骤S2具体包括如下步骤,
步骤S21,建立气候模型,气候模型用于表征历史时刻的时刻气象信息序列与海平面气压值的对应关系;
步骤S22,基于北斗定位模块获取当前时刻的时刻气象信息序列,并基于相似性判定自气候模型中获取时刻气象信息序列与当前时刻的最相似的历史时刻下的历史海平面气压值;
步骤S23,将气压计中的当前基准气压与经步骤S22中获取的历史海平面气压值进行比较,若比较结果在设定误差范围内则以当前基准气压值作为气压计的基准气压,若比较结果超出设定误差范围则以历史海平面气压值作为气压计的基准气压;
步骤S24,气压计获取当前高度的气压值并结合基准气压获取当前坐标点的海拔高度并作为高度坐标z。
通过建立气候模型,使得能够将气压计中的当前基准气压与气象条件相同或近似的历史时刻的海平面气压值进行比较,并能够设置误差范围,若比较结果在误差范围内,则判定当前基准气压有效进而对当前坐标点的海拔高度进行计算;若比较结果超出误差范围,则能够以历史海平面气压值替换当前基准气压,并对当前坐标点的海拔高度进行计算;故而能够较佳地对气压计的测量精度进行提升。
作为优选,步骤S21具体包括如下步骤,
步骤S211,采集施工现场所处区域的前U年的气象信息集合W,W={Wα|α=1,2,3,…,U},Wα为第α年的年气象信息序列;Wα={Gαβ|α=1,2,3,…,U;β=1,2,3,…,365},Gαβ为第α年第β日的日气象信息序列;Gαβ={Lαβγ|α=1,2,3,…,U;β=1,2,3,…,365;γ=1,2,3,…,q},Lαβγ为第α年第β日第γ时刻的时刻气象信息序列,q为时间序列的长度;
步骤S212,对气象信息集合W中的数据进行清洗;
步骤S213,建立气候模型Cl,气候模型Cl用于建立第α年第β日第γ时刻的气象信息集合Lαβγ与第α年第β日第γ时刻的海平面气压值ATαβγ间的对应关系,其中,Cl={(Lαβγ,ATαβγ)|α=1,2,3,…,U;β=1,2,3,…,365;γ=1,2,3,…,q}。
故而能够较佳地实现气候模型Cl的建立。
作为优选,步骤S212具体包括如下步骤,
步骤S212a,基于傅里叶级数拟合对气象信息集合W中的数据进行粗清洗,进而剔除每年年气象信息序列Wα中判定为噪音的日气象信息序列Gαβ;
步骤S212b,基于置信水平对气象信息集合W中的数据进行细清洗,进而剔除每日日气象信息序列Gαβ中判定为噪音的时刻气象信息序列Lαβγ。
通过上述,能够较佳地剔除噪音,故而能够较佳地提升测量精度。
作为优选,步骤S212a具体包括如下步骤,
步骤S212a1,对日气象信息序列Gαβ中的每个气象指标的日序列进行傅里叶级数拟合,拟合公式为,
其中,Qαβ(β)表示日气象信息序列Gαβ中的一个具体指标的拟合函数;Al和Bl为傅里叶系数,通过拟合获取;l表示为傅里叶阶数,p为傅里叶阶数的取值;ωl为预设参数,取值为4的倍数;
步骤S212a2,对于任意具体气象指标,其拟合函数的Qαβ(β)中的Al和Bl处于设定阈值内,则将对应日的日气象信息序列Gαβ判定为噪音并剔除。
通过上述,能够较佳地剔除对年参数影响不大的日参数,故而能够较佳地剔除无效的数据,以降低数据的大小。
作为优选,步骤S212b具体包括如下步骤,
步骤S212b1,基于公式对每个日气象信息序列Gαβ的时刻气象信息序列Lαβγ中的具体气象指标的置信水平Ul进行计算;/>为对应具体气象指标在日气象信息序列Gαβ中的均值,δ为对应具体气象指标在日气象信息序列Gαβ中的标准差。
步骤S212b2,对于任一时刻气象信息序列Lαβγ,若其中存在置信水平Ul低于0.95的具体气象指标数值,则将对应的时刻气象信息序列Lαβγ判定为噪音并剔除。
通过上述,能够较佳地对每日的日气象信息序列Gαβ进行处理,进而剔除无效的时刻气象信息序列Lαβγ,故而能够较佳地实现数据清洗。
作为优选,时刻气象信息序列为多个气象指标的数值序列,日气象信息序列为当日所有时刻的时刻气象信息序列的序列,年气象信息序列为当年所有日的日气象信息序列的序列。故而能够构建多指标的时刻气象信息序列,从而能够较佳地提升数据处理的精度。
作为优选,所述多个气象指标包括温度、湿度、风向、风速和日照辐射中的一个或多个。故而能够较佳地将对气压产生影响的多种因素均进行考虑。
作为优选,步骤S22中,基于欧式距离对时刻气象信息序列进行进行相似性判定。故而能够较佳地实现相似性判定。
此外,本发明还提供了一种基于北斗定位与气压计的高精度定位系统,其包括北斗定位模块、本地处理模块及气压计,本地处理模块处设置气候模型,气候模型用于对气压计的基准气压进行校正。故而能够较佳地实现任一上述的方法。
附图说明
图1为实施例1中的高精度定位方法及系统的示意图。
具体实施方式
为进一步了解本发明的内容,结合附图和实施例对本发明作详细描述。应当理解的是,实施例仅仅是对本发明进行解释而并非限定。
实施例1
结合图1所示,本实施例提供了一种基于北斗定位与气压计的高精度定位方法,其具有如下步骤:
步骤S1、基于北斗定位模块获取当前坐标点的平面坐标(x,y);
步骤S2,基于气压计获取当前坐标点的高度坐标z;
步骤S3,以(x,y,z)作为当前坐标点的坐标并输出。
通过本实施例中的方法,能够较佳地以北斗定位模块提供的平面坐标及气压计获取的高度坐标作为当前坐标点并输出,故而能够具备较佳的测量精度。
本实施例的步骤S2具体包括如下步骤:
步骤S21,建立气候模型,气候模型用于表征历史时刻的时刻气象信息序列与海平面气压值的对应关系;
步骤S22,基于北斗定位模块获取当前时刻的时刻气象信息序列,并基于相似性判定自气候模型中获取时刻气象信息序列与当前时刻的最相似的历史时刻下的历史海平面气压值;
步骤S23,将气压计中的当前基准气压与经步骤S22中获取的历史海平面气压值进行比较,若比较结果在设定误差范围内则以当前基准气压值作为气压计的基准气压,若比较结果超出设定误差范围则以历史海平面气压值作为气压计的基准气压;
步骤S24,气压计获取当前高度的气压值并结合基准气压获取当前坐标点的海拔高度并作为高度坐标z。
通过建立气候模型,使得能够将气压计中的当前基准气压与气象条件相同或近似的历史时刻的海平面气压值进行比较,并能够设置误差范围,若比较结果在误差范围内,则判定当前基准气压有效进而对当前坐标点的海拔高度进行计算;若比较结果超出误差范围,则能够以历史海平面气压值替换当前基准气压,并对当前坐标点的海拔高度进行计算;故而能够较佳地对气压计的测量精度进行提升。
可以理解的是,气压计作为一个现有器件,其在测量高度时,通过将所测量的气压值与基准气压值进行比较,即可基于气压-高度公式,获取当前测量位置的高度。其测量误差的最主要来源,在于基准气压值的变化。所谓基准气压值即海平面气压值,其不仅受到当前时刻的气象条件的影响,而且还会受到时刻的影响即日较差的影响。而通过本实施例中的方法,能够较佳地以同一地区的历史气象条件下的海平面气压值作为参考,对气压计的基准气压进行校正,故而能够较佳地提升气压计的输出精度。
本实施例中,能够通过北斗定位模块获取当前时刻的气象信息,并能够基于此构建当前时刻的气象信息序列。而后当前时刻的时刻气象信息序列能够输入至气候模型中,进而匹配与当前时刻的气象信息序列最相似的历史时刻的时刻气象信息序列,进而即可输出对应的历史时刻下的历史气压值,进而即可较佳地实现历史海平面气压值与当前基准气压的比较。
本实施例的步骤S22中,是以当前时刻的时刻气象信息序列与历史时刻的时刻气象信息序列进行匹配,在匹配到最相似的历史时刻气象信息序列后,气候模型能够输出该历史时刻气象信息序列所对应的历史时刻及历史海平面气压值。故本实施例中的误差范围的判断规则能够基于如下:
1、判断历史时刻与当前时刻的时间差异是否达到设定时间阈值,如30min;
2、判断历史海平面气压值与当前基准气压的气压差异是否达到设定气压阈值,如0.05hPa。
若上述规则1未满足(即差异值未达到设定阈值),则说明在历史数据中所匹配的与当前时刻气象信息序列最相似的历史时刻的时刻气象信息序列所属的历史时刻,与当前时刻在允许误差内,也即所匹配的历史时刻的时刻气象信息序列有效,该历史时刻的时刻气象信息序列所对应的历史海平面气压值能够作为参考与当前基准气压进行比较。
在上述规则1未满足时,若规则2也未满足(即差异值未达到设定阈值),则说明当前基准气压是有效的,能够作为气压计的基准气压。
在上述规则1未满足时,若规则2满足(即差异值达到设定阈值),则说明当前基准气压是无效的,则以所对应的历史海平面气压值作为气压计的基准气压。
若上述规则1满足(即差异值达到设定阈值),则说明所匹配的历史时刻的时刻气象信息序列所属的历史时刻,与当前时刻超出允许误差范围;也即未匹配到与当前时刻的气象条件相似的历史时刻。此时根据如下公式计算获取气压计的基准气压Pref:
上式中,U表示历史数据所包含的年度数量,Lτ表示与当前时刻相同的历史时刻下的历史海平面气压值,表示Lτ的权重,ετ表示当前时刻的时刻气象信息序列与相同时刻下的历史时刻的时刻气象信息序列的欧氏距离。
基于上述,能够较佳地通过对所有相同时刻的历史时刻下的历史海平面气压值进行加权计算,并将其结果作为气压计的基准气压Pref。
通过上述,能够充分考虑到气象条件及日较差对基准气压的影响,故而能够较佳地实现对气压计的基准气压的校准,故而使得测量精度能够得到较佳的提升。
本实施例的步骤S21具体包括如下步骤,
步骤S211,采集施工现场所处区域的前U年的气象信息集合W,W={Wα|α=1,2,3,…,U},Wα为第α年的年气象信息序列;Wα={Gαβ|α=1,2,3,…,U;β=1,2,3,…,365},Gαβ为第α年第β日的日气象信息序列;Gαβ={Lαβγ|α=1,2,3,…,U;β=1,2,3,…,365;γ=1,2,3,…,q},Lαβγ为第α年第β日第γ时刻的时刻气象信息序列,q为时间序列的长度;
步骤S212,对气象信息集合W中的数据进行清洗;
步骤S213,建立气候模型Cl,气候模型Cl用于建立第α年第β日第γ时刻的气象信息集合Lαβγ与第α年第β日第γ时刻的海平面气压值ATαβγ间的对应关系,其中,Cl={(Lαβγ,ATαβγ)|α=1,2,3,…,U;β=1,2,3,…,365;γ=1,2,3,…,q}。
故而能够较佳地实现气候模型Cl的建立。
本实施例中,U的取值能够为20年。
本实施例的步骤S212具体包括如下步骤,
步骤S212a,基于傅里叶级数拟合对气象信息集合W中的数据进行粗清洗,进而剔除每年年气象信息序列Wα中判定为噪音的日气象信息序列Gαβ;
步骤S212b,基于置信水平对气象信息集合W中的数据进行细清洗,进而剔除每日日气象信息序列Gαβ中判定为噪音的时刻气象信息序列Lαβγ。
通过上述,能够较佳地剔除噪音,故而能够较佳地提升测量精度。
本实施例的步骤S212a具体包括如下步骤,
步骤S212a1,对日气象信息序列Gαβ中的每个气象指标的日序列进行傅里叶级数拟合,拟合公式为,
其中,Qαβ(β)表示日气象信息序列Gαβ中的一个具体指标的拟合函数;Al和Bl为傅里叶系数,通过拟合获取;l表示为傅里叶阶数,p为傅里叶阶数的取值;ωl为预设参数,取值为4的倍数;
步骤S212a2,对于任意具体气象指标,其拟合函数的Qαβ(β)中的Al和Bl处于设定阈值内,则将对应日的日气象信息序列Gαβ判定为噪音并剔除。
通过上述,能够较佳地剔除对年参数影响不大的日参数,故而能够较佳地剔除无效的数据,以降低数据的大小。
其中,l的取值能够为4或8。
其中,S212a2中所设置的阈值能够设定为如0.05等常数。
通过步骤S212a2,能够剔除Al和Bl趋近于零的日气象信息序列Gαβ,故而能够较佳地剔除对于年参数影响不大的日参数。
本实施例的步骤S212b具体包括如下步骤,
步骤S212b1,基于公式对每个日气象信息序列Gαβ的时刻气象信息序列Lαβγ中的具体气象指标的置信水平Ul进行计算;/>为对应具体气象指标在日气象信息序列Gαβ中的均值,δ为对应具体气象指标在日气象信息序列Gαβ中的标准差。
步骤S212b2,对于任一时刻气象信息序列Lαβγ,若其中存在置信水平Ul低于0.95的具体气象指标数值,则将对应的时刻气象信息序列Lαβγ判定为噪音并剔除。
通过上述,能够较佳地对每日的日气象信息序列Gαβ进行处理,进而剔除无效的时刻气象信息序列Lαβγ,故而能够较佳地实现数据清洗。
本实施例中,时刻气象信息序列为多个气象指标的数值序列,日气象信息序列为当日所有时刻的时刻气象信息序列的序列,年气象信息序列为当年所有日的日气象信息序列的序列。故而能够构建多指标的时刻气象信息序列,从而能够较佳地提升数据处理的精度。
本实施例中,所述多个气象指标包括温度、湿度、风向、风速和日照辐射中的一个或多个。故而能够较佳地将对气压产生影响的多种因素均进行考虑。
本实施例中,步骤S22中,基于欧式距离对时刻气象信息序列进行进行相似性判定。故而能够较佳地实现相似性判定。
为了进一步地对本实施例中的方法进行说明,下述采用一个具体的实施例进行说明。
在该具体的实施例中,选取温度(Pm1)、湿度(Pm2)、风向(Pm3)、风速(Pm4)和日照辐射(Pm5)作为气象指标构建时刻气象信息序列。
对于一个历史时刻的时刻气象信息序列Lαβγ,其即为:
其中,和/>分别表示第α年第β日第γ时刻下的温度(Pm1)、湿度(Pm2)、风向(Pm3)、风速(Pm4)和日照辐射(Pm5)的数值。
当前时刻t的时刻气象信息序列即可表示为:
其中,和/>分别表示当前时刻t下的温度(Pm1)、湿度(Pm2)、风向(Pm3)、风速(Pm4)和日照辐射(Pm5)的数值。
故在进行相似性匹配时,相似度(欧式距离)的计算公式即为:
其中,通过逐个计算当前时刻的时刻气象信息序列与历史时刻的时刻气象信息序列的相似度,并取相似度最小的即为所匹配的历史时刻的时刻气象信息序列。
其中,对于一个日气象信息序列Gαβ,其即为:
故在步骤S212a1中,所述的每个气象指标的日序列,即为温度(Pm1)、湿度(Pm2)、风向(Pm3)、风速(Pm4)和日照辐射(Pm5)的当日数值在时刻上的数列。即:
温度(Pm1)气象指标的日序列为:
湿度(Pm2)气象指标的日序列为:
风向(Pm3)气象指标的日序列为:
风速(Pm4)气象指标的日序列为:
日照辐射(Pm5)气象指标的日序列为:
通过步骤S212a1即可较佳地实现对每个气象指标的日序列的拟合,并能在任一气象指标的日序列的傅里叶系数Al和Bl处于设定阈值内时,均将该日的日气象信息序列Gαβ剔除。
在步骤步骤S212b中,即对每个气象指标的日序列的数值进行置信水平的计算,并在任一气象指标的置信水平达不到0.95时,均将该气象指标所在的整个历史时刻的时刻气象信息序列Lαβγ予以剔除。
通过本实施例中的方法,能够较佳地对历史气象数据进行处理,并建立气候模型,通过考虑当前时刻的时刻值及气象条件与历史时刻的时刻值及气象条件的差异,并对气压计的基准气压进行校正,能够较佳地充分考虑到气象条件及日较差对基准气压的影响,故而能够较佳地提升气压计的定位精度。
基于本实施例中的方法,本实施例还提供了一种基于北斗定位与气压计的高精度定位系统,其具有北斗定位模块、本地处理模块及气压计,气候模型设于本地处理模块处,本地处理模块及气压计均设于终端设备处。
以上示意性的对本发明及其实施方式进行了描述,该描述没有限制性,附图中所示的也只是本发明的实施方式之一,实际的结构并不局限于此。所以,如果本领域的普通技术人员受其启示,在不脱离本发明创造宗旨的情况下,不经创造性的设计出与该技术方案相似的结构方式及实施例,均应属于本发明的保护范围。
Claims (8)
1.基于北斗定位与气压计的高精度定位方法,其具有如下步骤:
步骤S1、基于北斗定位模块获取当前坐标点的平面坐标(x,y);
步骤S2,基于气压计获取当前坐标点的高度坐标z;
步骤S3,以(x,y,z)作为当前坐标点的坐标并输出;
步骤S2具体包括如下步骤,
步骤S21,建立气候模型,气候模型用于表征历史时刻的时刻气象信息序列与海平面气压值的对应关系;
步骤S22,基于北斗定位模块获取当前时刻的时刻气象信息序列,并基于相似性判定自气候模型中获取时刻气象信息序列与当前时刻的最相似的历史时刻下的历史海平面气压值,并基于欧式距离对时刻气象信息序列进行进行相似性判定;
其中,逐个计算当前时刻的时刻气象信息序列与历史时刻的时刻气象信息序列的相似度,以欧氏距离最小的为所匹配的历史时刻的时刻气象信息序列,进而获取历史海平面气压值;
步骤S23,将气压计中的当前基准气压与经步骤S22中获取的历史海平面气压值进行比较,并建立误差判断规则;
其中,误差范围的判断规则包括,
规则1、判断历史时刻与当前时刻的时间差异是否达到设定时间阈值,
规则2、判断历史海平面气压值与当前基准气压的气压差异是否达到设定气压阈值;
在进行比较时,
若规则1满足则以如下公式计算获取气压计的基准气压Pref,
上式中,U表示历史数据所包含的年度数量,Lτ表示与当前时刻相同的历史时刻下的历史海平面气压值,表示Lτ的权重,ετ表示当前时刻的时刻气象信息序列与相同时刻下的历史时刻的时刻气象信息序列的欧氏距离;
若规则1未满足且规则2满足,则以经步骤S22中获取的历史海平面气压值作为气压计的基准气压;
若规则1未满足且规则2未满足,则以当前基准气压值作为气压计的基准气压;
步骤S24,气压计获取当前高度的气压值并结合基准气压获取当前坐标点的海拔高度并作为高度坐标z。
2.根据权利要求1所述的基于北斗定位与气压计的高精度定位方法,其特征在于:步骤S21具体包括如下步骤,
步骤S211,采集施工现场所处区域的前U年的气象信息集合W,W={Wα|α=1,2,3,…,U},Wα为第α年的年气象信息序列;Wα={Gαβ|α=1,2,3,…,U;β=1,2,3,…,365},Gαβ为第α年第β日的日气象信息序列;Gαβ={Lαβγ|α=1,2,3,…,U;β=1,2,3,…,365;γ=1,2,3,…,q},Lαβγ为第α年第β日第γ时刻的时刻气象信息序列,q为时间序列的长度;
步骤S212,对气象信息集合W中的数据进行清洗;
步骤S213,建立气候模型Cl,气候模型Cl用于建立第α年第β日第γ时刻的气象信息集合Lαβγ与第α年第β日第γ时刻的海平面气压值ATαβγ间的对应关系,其中,Cl={(Lαβγ,ATαβγ)|α=1,2,3,…,U;β=1,2,3,…,365;γ=1,2,3,…,q}。
3.根据权利要求2所述的基于北斗定位与气压计的高精度定位方法,其特征在于:步骤S212具体包括如下步骤,
步骤S212a,基于傅里叶级数拟合对气象信息集合W中的数据进行粗清洗,进而剔除每年年气象信息序列Wα中判定为噪音的日气象信息序列Gαβ;
步骤S212b,基于置信水平对气象信息集合W中的数据进行细清洗,进而剔除每日日气象信息序列Gαβ中判定为噪音的时刻气象信息序列Lαβγ。
4.根据权利要求3所述的基于北斗定位与气压计的高精度定位方法,其特征在于:步骤S212a具体包括如下步骤,
步骤S212a1,对日气象信息序列Gαβ中的每个气象指标的日序列进行傅里叶级数拟合,拟合公式为,
其中,Qαβ(β)表示日气象信息序列Gαβ中的一个具体指标的拟合函数;Al和Bl为傅里叶系数,通过拟合获取;l表示为傅里叶阶数,p为傅里叶阶数的取值;ωl为预设参数,取值为4的倍数;
步骤S212a2,对于任意具体气象指标,其拟合函数的Qαβ(β)中的Al和Bl处于设定阈值内,则将对应日的日气象信息序列Gαβ判定为噪音并剔除。
5.根据权利要求3所述的基于北斗定位与气压计的高精度定位方法,其特征在于:步骤S212b具体包括如下步骤,
步骤S212b1,基于公式对每个日气象信息序列Gαβ的时刻气象信息序列Lαβγ中的具体气象指标的置信水平Ul进行计算;/>为对应具体气象指标在日气象信息序列Gαβ中的均值,δ为对应具体气象指标在日气象信息序列Gαβ中的标准差;
步骤S212b2,对于任一时刻气象信息序列Lαβγ,若其中存在置信水平Ul低于0.95的具体气象指标数值,则将对应的时刻气象信息序列Lαβγ判定为噪音并剔除。
6.根据权利要求1-5任一所述的基于北斗定位与气压计的高精度定位方法,其特征在于:时刻气象信息序列为多个气象指标的数值序列,日气象信息序列为当日所有时刻的时刻气象信息序列的序列,年气象信息序列为当年所有日的日气象信息序列的序列。
7.根据权利要求6所述的基于北斗定位与气压计的高精度定位方法,其特征在于:所述多个气象指标包括温度、湿度、风向、风速和日照辐射中的一个或多个。
8.基于北斗定位与气压计的高精度定位系统,用于实现权利要求1-7中任一所述的基于北斗定位与气压计的高精度定位方法,其特征在于:包括北斗定位模块、本地处理模块及气压计,本地处理模块处设置气候模型,气候模型用于对气压计的基准气压进行校正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111027695.9A CN113701711B (zh) | 2021-09-02 | 2021-09-02 | 基于北斗定位与气压计的高精度定位方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111027695.9A CN113701711B (zh) | 2021-09-02 | 2021-09-02 | 基于北斗定位与气压计的高精度定位方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113701711A CN113701711A (zh) | 2021-11-26 |
CN113701711B true CN113701711B (zh) | 2023-11-03 |
Family
ID=78658933
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111027695.9A Active CN113701711B (zh) | 2021-09-02 | 2021-09-02 | 基于北斗定位与气压计的高精度定位方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113701711B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115127514B (zh) * | 2022-03-31 | 2024-03-19 | 广东小天才科技有限公司 | 高度测量方法、装置、电子设备及存储介质 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102735213A (zh) * | 2011-04-02 | 2012-10-17 | 中国科学院国家天文台 | 基于基准站的气压相对测量求得精确高程的方法 |
CN104616075A (zh) * | 2015-01-30 | 2015-05-13 | 广西大学 | 一种适合于台风天气下的短期负荷预测方法 |
JP2016024041A (ja) * | 2014-07-18 | 2016-02-08 | Kddi株式会社 | 自ら位置する地上高を推定する移動端末、プログラム及び方法 |
CN105509770A (zh) * | 2016-01-07 | 2016-04-20 | 中国科学院嘉兴微电子与系统工程中心 | 一种gnss和mems组合导航系统中的气压计在线校正方法 |
JP2016109637A (ja) * | 2014-12-10 | 2016-06-20 | ソニー株式会社 | 装置及び方法 |
CN106851585A (zh) * | 2017-01-12 | 2017-06-13 | 杭州电子科技大学 | 一种基于气压计和WiFi的混合楼层定位方法 |
CN108549961A (zh) * | 2018-05-02 | 2018-09-18 | 河海大学 | 一种基于cmip5预估海浪有效波高的方法 |
CN110108254A (zh) * | 2019-04-30 | 2019-08-09 | 大器物联科技(广州)有限公司 | 一种基于气压精确测量高度、垂直运动速度和垂直运动距离的方法 |
CN110866630A (zh) * | 2019-09-27 | 2020-03-06 | 兰州大方电子有限责任公司 | 一种历史相似性天气分析方法 |
CN111639802A (zh) * | 2020-05-28 | 2020-09-08 | 中电投珠海横琴热电有限公司 | 一种燃机机组运行优化指导方法 |
-
2021
- 2021-09-02 CN CN202111027695.9A patent/CN113701711B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102735213A (zh) * | 2011-04-02 | 2012-10-17 | 中国科学院国家天文台 | 基于基准站的气压相对测量求得精确高程的方法 |
JP2016024041A (ja) * | 2014-07-18 | 2016-02-08 | Kddi株式会社 | 自ら位置する地上高を推定する移動端末、プログラム及び方法 |
JP2016109637A (ja) * | 2014-12-10 | 2016-06-20 | ソニー株式会社 | 装置及び方法 |
CN104616075A (zh) * | 2015-01-30 | 2015-05-13 | 广西大学 | 一种适合于台风天气下的短期负荷预测方法 |
CN105509770A (zh) * | 2016-01-07 | 2016-04-20 | 中国科学院嘉兴微电子与系统工程中心 | 一种gnss和mems组合导航系统中的气压计在线校正方法 |
CN106851585A (zh) * | 2017-01-12 | 2017-06-13 | 杭州电子科技大学 | 一种基于气压计和WiFi的混合楼层定位方法 |
CN108549961A (zh) * | 2018-05-02 | 2018-09-18 | 河海大学 | 一种基于cmip5预估海浪有效波高的方法 |
CN110108254A (zh) * | 2019-04-30 | 2019-08-09 | 大器物联科技(广州)有限公司 | 一种基于气压精确测量高度、垂直运动速度和垂直运动距离的方法 |
CN110866630A (zh) * | 2019-09-27 | 2020-03-06 | 兰州大方电子有限责任公司 | 一种历史相似性天气分析方法 |
CN111639802A (zh) * | 2020-05-28 | 2020-09-08 | 中电投珠海横琴热电有限公司 | 一种燃机机组运行优化指导方法 |
Non-Patent Citations (1)
Title |
---|
基于地面移动通信基站的差分气压测高方法;杜晓辉 等;《北京航空航天大学学报》;第第39卷卷(第第1期期);第83-88页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113701711A (zh) | 2021-11-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106053724A (zh) | 基于云计算的气体传感器精度补偿方法及装置 | |
CN113701711B (zh) | 基于北斗定位与气压计的高精度定位方法及系统 | |
CN108692711B (zh) | 基于低空探空火箭的海洋数据处理的实现方法 | |
CN107271081A (zh) | 基于两阶段最小二乘拟合的硅压阻式压力变送器温度补偿方法及装置 | |
US10989837B2 (en) | Method in connection with a radiosonde and system | |
CN113865750A (zh) | 非接触式设备的测温校准方法及非接触式设备 | |
CN113393525B (zh) | 基于精准称重的重量校准方法 | |
CN117490819A (zh) | 一种用于电子秤的误差补偿系统及方法 | |
CN109758703B (zh) | 一种用于消防火场气压高度传感器的误差修正系统及方法 | |
CN115859211B (zh) | 一种基于三温不确定度估算模型的地表温度产品融合方法 | |
CN113655501A (zh) | 基于北斗定位的现场带电体安全距离区域划分系统及方法 | |
Permana et al. | DS18B20 sensor calibration compared with fluke hart scientific standard sensor | |
Lake | Calibration procedures and instrumental accuracy estimates of ATLAS air temperature and relative humidity measurements | |
CN113639880B (zh) | 一种真空红外定标辐射响应一致性方法 | |
CN113654519B (zh) | 基于北斗定位的终端楼层判定装置及方法 | |
CN113932765A (zh) | 基于温度补偿的高精度静力水准仪及其温度补偿方法 | |
CN113721264B (zh) | 基于北斗定位的终端设备的中继与管理系统及方法 | |
CN111044489B (zh) | 一种基于多波长测量获得大气折射率高度分布剖面的方法 | |
CN113654527B (zh) | 基于北斗定位的施工现场全景管控显示方法及系统 | |
CN116380339B (zh) | 一种薄膜规真空计校准方法及相关设备 | |
CN111123406A (zh) | 手持气象仪温度数据拟合方法 | |
CN111024297A (zh) | 一种基于多压力传感器的内腔压力检测校准方法 | |
Wang et al. | Data processing and experiment of barometric altimeter based on recursive least squares filtering | |
CN115523896A (zh) | 海拔获取方法、系统、码表及计算机可读存储介质 | |
SU1530896A1 (ru) | Способ барометрического нивелировани |
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 |