CN115877420B - 一种基于静止卫星的定位方法、系统、电子设备及介质 - Google Patents
一种基于静止卫星的定位方法、系统、电子设备及介质 Download PDFInfo
- Publication number
- CN115877420B CN115877420B CN202211267794.9A CN202211267794A CN115877420B CN 115877420 B CN115877420 B CN 115877420B CN 202211267794 A CN202211267794 A CN 202211267794A CN 115877420 B CN115877420 B CN 115877420B
- Authority
- CN
- China
- Prior art keywords
- longitude
- latitude
- data
- standard
- satellite
- 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
-
- 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
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Image Processing (AREA)
Abstract
本发明提供一种基于静止卫星的定位方法、系统、电子设备及介质,属于卫星定位技术领域,定位方法包括:获取目标卫星的待定位数据;待定位数据中包括对应区域范围内各点的行号和列号;将待定位数据的格式转换为标准格式,得到待定位标准数据;基于预先建立的等间隔经纬度投影查找表,根据待定位标准数据,确定对应区域范围内各点的最优经度及最优纬度;等间隔经纬度投影查找表为预先根据目标卫星的历史数据及标准卫星数据建立的;等间隔经纬度投影查找表中包括对应区域范围内各点的行号、列号、最优经度及最优纬度。能够基于不符合国际静止卫星定位标准的卫星数据进行精准定位。
Description
技术领域
本发明涉及卫星定位领域,特别是涉及一种基于印度静止卫星(INSAT-3DR)的定位方法、系统、电子设备及介质。
背景技术
图像定位是地球静止气象卫星的核心技术之一,也是地球静止气象卫星产品开发的基础和预处理工作。目前对静止卫星图像定位多采用定位网格矩阵的简单方法,该方法对网格内点插值较容易实现,精确定位方法利用星上各种实时参数,误差相对较小。
但是,目前一些卫星的设计并不完全符合国际标准,如印度的INSAT-3DR系列卫星,其定位方式与国际标准并不一致,其圆盘相邻经度的间隔不同,这就使得中国风云四号系列、日本的葵花系列、美国的GOS系列等符合国际静止卫星定位标准的卫星定位方法无法使用,因此,亟需一种能够适用于不符合国际静止卫星定位标准的卫星定位方法。
发明内容
本发明的目的是提供一种基于静止卫星的定位方法、系统、电子设备及介质,可基于不符合国际静止卫星定位标准的卫星数据(如印度静止卫星INSAT-3DR)进行精准定位。
为实现上述目的,本发明提供了如下方案:
一种基于静止卫星的定位方法,包括:
获取目标卫星的待定位数据;所述待定位数据中包括对应区域范围内各点的行号和列号;所述目标卫星为不符合国际静止卫星定位标准的卫星;
将所述待定位数据的格式转换为标准格式,得到待定位标准数据;
基于预先建立的等间隔经纬度投影查找表,根据所述待定位标准数据,确定对应区域范围内各点的最优经度及最优纬度;所述等间隔经纬度投影查找表为预先根据目标卫星的历史数据及标准卫星数据建立的;所述等间隔经纬度投影查找表中包括对应区域范围内各点的行号、列号、最优经度及最优纬度。
可选地,所述目标卫星为印度静止卫星INSAT-3DR。
可选地,所述获取目标卫星的待定位数据,具体包括:
获取目标卫星的原始数据;所述原始数据中包括多个通道数据;
对所述原始数据进行解析,确定各通道数据对应的数组,以得到待定位数据;各通道数据对应的数组中包括对应区域范围内各点的行号和列号。
可选地,所述将所述待定位数据的格式转换为标准格式,得到待定位标准数据,具体包括:
根据标准卫星的数据格式,采用双线性插值法对所述待定位数据进行重采样,得到待定位标准数据。
可选地,所述基于静止卫星的定位方法还包括:
获取目标卫星的历史数据及标准卫星数据;所述历史数据的区域范围与所述标准卫星数据的区域范围相同;所述历史数据中包括对应区域范围内各点的行号和列号;所述标准卫星数据中包括对应区域范围内各点的经度和纬度;
将所述历史数据的格式转换为标准格式,得到历史标准数据;
根据所述历史标准数据及所述目标卫星的轨道信息,生成初始等经纬图像;所述初始等经纬图像中包括对应区域范围内各点的粗略经度及粗略纬度;
基于所述标准卫星数据,根据所述初始等经纬图像,确定等间隔经纬图像;所述等间隔经纬图像中包括对应区域范围内各点的最优经度、最优纬度、行号及列号;
基于所述等间隔经纬图像生成等间隔经纬度投影查找表。
可选地,所述根据所述历史标准数据及所述目标卫星的轨道信息,生成初始等经纬图像,具体包括:
根据所述历史标准数据及所述目标卫星的轨道信息,生成地理网格;所述地理网格中包括多个网格及各网格对应的粗略经度及粗略纬度;各网格代表对应区域范围内的各点;
根据所述地理网格中各网格对应的粗略经度及粗略纬度,采用双线性插值法,对所述地理网格进行插值,以生成初始等经纬图像。
可选地,所述基于所述标准卫星数据,根据所述初始等经纬图像,确定等间隔经纬图像,具体包括:
从所述初始等经纬图像及所述标准卫星数据中选取多个同名点;
确定各同名点在所述初始等经纬图像中的粗略经度及粗略纬度,得到粗略经纬集;
确定各同名点在所述标准卫星数据中的经度及纬度,得到标准经纬集;
对所述粗略经纬集中各同名点的粗略经度、粗略纬度、所述标准经纬集中各同名点的经度和纬度进行线性拟合,确定缩放系数及偏移系数;
根据所述缩放系数及所述偏移系数,计算所述初始等经纬图像中各点的最优经度和最优纬度;
根据所述初始等经纬图像中各点的最优经度和最优纬度,生成等间隔经纬图像。
为实现上述目的,本发明还提供了如下方案:
一种基于静止卫星的定位系统,包括:
数据获取单元,用于获取目标卫星的待定位数据;所述待定位数据中包括对应区域范围内各点的行号和列号;所述目标卫星为不符合国际静止卫星定位标准的卫星;
格式转换单元,与所述数据获取单元连接,用于将所述待定位数据的格式转换为标准格式,得到待定位标准数据;
经纬确定单元,与所述格式转换单元连接,用于基于预先建立的等间隔经纬度投影查找表,根据所述待定位标准数据,确定对应区域范围内各点的最优经度及最优纬度;所述等间隔经纬度投影查找表为预先根据目标卫星的历史数据及标准卫星数据建立的;所述等间隔经纬度投影查找表中包括对应区域范围内各点的行号、列号、最优经度及最优纬度。
为实现上述目的,本发明还提供了如下方案:
一种电子设备,包括存储器及处理器,所述存储器用于存储计算机程序,所述处理器运行所述计算机程序以使所述电子设备执行上述的基于静止卫星的定位方法。
为实现上述目的,本发明还提供了如下方案:
一种计算机可读存储介质,其存储有计算机程序,所述计算机程序被处理器执行时实现上述的基于静止卫星的定位方法。
根据本发明提供的具体实施例,本发明公开了以下技术效果:首先获取目标卫星的待定位数据;其中,待定位数据中包括对应区域范围内各点的行号和列号;然后将待定位数据的格式转换为标准格式,得到待定位标准数据,使其与标准卫星的数据格式一致;再基于预先建立的等间隔经纬度投影查找表,根据待定位标准数据,确定对应区域范围内各点的经纬度,其中,等间隔经纬度投影查找表为预先根据目标卫星的历史数据及标准卫星数据建立的;等间隔经纬度投影查找表中包括对应区域范围内各点的行号、列号、经度及纬度。通过将不符合国际静止卫星定位标准的卫星数据进行格式转换,并查找对应区域内各点的最优经纬度,能够基于不符合国际静止卫星定位标准的卫星数据进行精准定位。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明基于静止卫星的定位方法的流程图;
图2为本发明基于静止卫星的定位系统的模块示意图。
符号说明:
数据获取单元-1,格式转换单元-2,经纬确定单元-3,历史数据获取单元-4,历史格式转换单元-5,初始图像生成单元-6,经纬图像生成单元-7,查找表生成单元-8。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种基于静止卫星的定位方法、系统、电子设备及介质,通过将不符合国际静止卫星定位标准的卫星数据进行格式转换,并查找对应区域内各点的最优经纬度,能够基于不符合国际静止卫星定位标准的卫星数据进行精准定位。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例一
如图1所示,本发明基于静止卫星的定位方法包括:
S1:获取目标卫星的待定位数据。所述待定位数据中包括对应区域范围内各点的行号和列号。在本实施例中,目标卫星为不符合国际静止卫星定位标准的卫星。优选为印度静止卫星INSAT-3DR。
具体地,首先获取目标卫星的原始数据。所述原始数据中包括多个通道数据。再对所述原始数据进行解析,确定各通道数据对应的数组,以得到待定位数据。各通道数据对应的数组中包括对应区域范围内各点的行号和列号。
INSAT-3DR原始数据为整型的DN值,本发明按照无符号16位整数格式读取,其行数、列数分别为文件名中使用下划线分隔的含H后缀项、含W后缀项。其中,可见光通道(VIS)和短波红外通道(SWIR)的行数为11264,列数为19636。中波红外通道(MWIR)、红外分裂窗通道(IR1、IR2)的行数为2816,列数为4909。水汽通道(WV)的行数为1408,列数为4909。
然后将各通道的数据处理为相应行数列数的数组。最终生成每个通道的文件,文件名称中包含通道、时间等信息。所有通道的文件共同构成待定位数据。
S2:将所述待定位数据的格式转换为标准格式,得到待定位标准数据。
具体地,根据标准卫星的数据格式,采用双线性插值法对所述待定位数据进行重采样,得到待定位标准数据。以将待定位数据转换为标准产品相应通道数组的格式。
标准卫星的数据格式中,可见光通道(VIS)和短波红外通道(SWIR)的行数为11264,列数为11220。中波红外通道(MIR)、红外分裂窗通道(TIR1、TIR2)的行数为2816,列数为2805。水汽通道(WV)的行数为1408,列数为1402。
S3:基于预先建立的等间隔经纬度投影查找表,根据所述待定位标准数据,确定对应区域范围内各点的最优经度及最优纬度。所述等间隔经纬度投影查找表为预先根据目标卫星的历史数据及标准卫星数据建立的。所述等间隔经纬度投影查找表中包括对应区域范围内各点的行号、列号、最优经度及最优纬度。
对于不符合国际静止卫星定位标准的卫星,除了标准的线性关系计算之外还要进行重新校正。因此,本发明通过建立等间隔经纬度投影查找表对卫星数据的经纬度进行校正,以提高定位精度。
进一步地,在建立等间隔经纬度投影查找表方面,本发明基于静止卫星的定位方法还包括:
S101:获取目标卫星的历史数据及标准卫星数据。所述历史数据的区域范围与所述标准卫星数据的区域范围相同。所述历史数据中包括对应区域范围内各点的行号和列号。所述标准卫星数据中包括对应区域范围内各点的经度和纬度。
在本实施例中,标准卫星数据的来源为葵花8卫星。
S102:将所述历史数据的格式转换为标准格式,得到历史标准数据。
S103:根据所述历史标准数据及所述目标卫星的轨道信息,生成初始等经纬图像。所述初始等经纬图像中包括对应区域范围内各点的粗略经度及粗略纬度。
具体地,首先根据所述历史标准数据及所述目标卫星的轨道信息,生成地理网格。所述地理网格中包括多个网格及各网格对应的粗略经度及粗略纬度。各网格代表对应区域范围内的各点。
然后根据所述地理网格中各网格对应的粗略经度及粗略纬度,采用双线性插值法,对所述地理网格进行插值,以生成初始等经纬图像。在本实施例中,初始等经纬图像的步长为0.05°。
S104:基于所述标准卫星数据,根据所述初始等经纬图像,确定等间隔经纬图像。所述等间隔经纬图像中包括对应区域范围内各点的最优经度、最优纬度、行号及列号。
具体地,步骤S104包括:
(1)从所述初始等经纬图像及所述标准卫星数据中选取多个同名点。
(2)确定各同名点在所述初始等经纬图像中的粗略经度及粗略纬度,得到粗略经纬集。
具体地,从初始等经纬图像上选取与葵花8卫星同波长通道的同名点,并找对应位置的经纬度,形成INSAT-3DR同名点位置(行列号)——粗略经纬集。
(3)确定各同名点在所述标准卫星数据中的经度及纬度,得到标准经纬集。
(4)对所述粗略经纬集中各同名点的粗略经度、粗略纬度、所述标准经纬集中各同名点的经度和纬度进行线性拟合,确定缩放系数及偏移系数。具体地,将粗略经纬集及标准经纬集中各同名点位置(行列号)与经纬度进行线性拟合,得到缩放系数和偏移系数。
(5)根据所述缩放系数及所述偏移系数,计算所述初始等经纬图像中各点的最优经度和最优纬度。利用缩放系数和偏移系数计算INSAT-3DR各通道数据点对应的最优经纬度。
(6)根据所述初始等经纬图像中各点的最优经度和最优纬度,生成等间隔经纬图像。具体地,利用各通道数据点对应的最优经纬度,生成等间隔经纬图像的标准文件。
S105:基于所述等间隔经纬图像生成等间隔经纬度投影查找表。即通过标准文件中的经纬数据和行列数据,制作等间隔经纬度投影查找表。
此外,由于静止卫星全圆盘成像的特点,边缘数据受载荷扫描角度以及地球曲率影响,越到边缘数据形变越大,数据越失真。所以本发明将初始等经纬图像中边缘的20行/列的数据剔除,使用剩下数据的经纬度范围和分辨率制作等间隔经纬度数组和行列数组。通过标准文件中的经纬数据(最优经度和最优纬度)及其行列数据和等间隔经纬度数组及其行列数组,将相同经纬度的标准文件的行列数据和等经纬的行列数据对应起来,制作所需区域的等间隔经纬度投影查找表。最后使用等间隔经纬度投影查找表将待定位数据转换为等间隔经纬投影数据。
进一步地,等间隔经纬度投影查找表可根据最终成像结果的陆地边界进行调整,不同分辨率通道不同时次卫星数据会有不同程度的偏移,同一通道相近日期相同时次的偏移量一致。
本发明为印度静止卫星(INSAT-3DR)提供了精确定位的方法,可有效填补当前卫星在观测定位技术上的空白。采用本发明的定位方法,对INSAT-3DR卫星的通道数据进行投影定位,将定位结果与海陆边界线进行比对,像元定位平均偏差优于2个像元,证明了本发明提供的定位方法进行印度静止卫星精确定位的可行性。
为了更好的理解本发明的方案,下面结合印度静止卫星(INSAT-3DR)的IR1通道数据进一步进行说明。
(1)数据解析
选择IR1通道文件名为IR1202009170012239408_4909W_2816H_2Depth.raw的原始数据,使用Python环境下的numpy库的fromfile函数按照16位整数进行读取,得到行数为2816,列数为4909的数组。
(2)数据重采样
将上一步中读取到的数组通过OpenCV中的resize函数,使用双线性插值法重采样为行数为2816,列数为2805的标准数据。
(3)基于插值的等经纬数据生成
1)地理网格的计算
取地球赤道半径6378137m计算,赤道地区1°相当于111km,5°就相当于555km,也就是555km/5km=111个像元,即每隔5°有111个像元,所以应取d=5°/111≈0.05°,即像元间的步长为0.05°。最终生成高宽为(25-1)×5°/0.05°+1=2401,步长为0.05°的等经纬地理网格数据,记为g2401×2401×2。具体生成过程如下:
首先生成一个大小为25×25的阵列A,元素为网格的行列号,再将阵列A插值成2401×2401大小的阵列B,则阵列B记录了每个经纬度在标准数据中的行列位置。为提高插值的平滑度,插值过程采用双三次插值(双立方插值或立方卷积插值),由于双三次插值需要插值元素的16邻域,阵列A的边界部分不足16邻域的由边界元素值补齐。
2)等经纬图像像素值的计算
初始等经纬图像即等经纬地理网格数据对应的图像数据,记为图像H,再记原始图像(等经纬地理网格数据)为F,则初始等经纬图像的生成过程如下:
H(i,j)=F[B(i,j,0),B(i,j,1)];
其中,i∈[0,2816],j∈[0,4909]且为整数,F的定义域为具有行列意义的整数,B(i,j,0)为等经纬地理网格数据中i行j列对应的经度位置,B(i,j,1)为等经纬地理网格数据中i行j列对应的纬度位置,二者为float型,要想得到F[B(i,j,0),B(i,j,1)]的值,必须对F进行插值。考虑到卫星分辨率低的特性,像元间的相关性小,因而采用双线性插值法,插值公式为
f(i+u,j+v)=(1-u)(1-v)f(i,j)+(1-u)vf(i,j+1)+u(1-v)f(i+1,j)+uvf(i+1,j+1);u,v∈[0,1];
由此,生成了一幅步长为0.05°的初始等经纬图像。
(4)原始图像定位
首先引入集合的概念,记等间隔经纬图像g为集合M,对初始等经纬图像G′进行稀疏化处理,记为集合N,则M→N,即集合M中的每一个元素总能在集合N中找到唯一一个元素与之对应。记初始等经纬图像G′为集合N′,等间隔经纬图像g稀疏化处理为M′,则N′→M′,即集合N′中的每一个元素总能在集合M′中找到唯一一个元素与之对应。记
B(i,j,0)→B1,B(i,j,1)→B2。
由经纬度与行列的关系可知,对于B1的每列而言,行元素是递增的,对于B2的每行而言,列元素是递增的。对B1的列进行行搜索,对B2的行进行列搜索。简化网格的间距为5°,计算可知,每隔5°约有111个像元,星下点附近的四边形变形小,边缘部分的四边形变形大。边缘部分是插值,而中间部分则是抽取,组合出简化网格内部所有点的行列号以及简化网格内部所有点的地理经纬度。
由于地球曲面的影响,对于非等间距的经纬度数据,找出矢量数据中记录的点、线、面的地理坐标对应点的行列号,再对这些点按矢量数据中的记录次序画线连接起来,形成初始等经纬图像。
投影转换过程中需要缩放和偏移系数,因此从初始等经纬图像图上选点,然后找对应位置的经纬度,进行线性拟合得到缩放和偏移系数,根据缩放系数及偏移系数,计算初始等经纬图像中各点的最优经度和最优纬度。
根据定位计算,原始数据的图像中心列的位置,不是INSAT-3DR轨道所在的东经74度,而是位于东经73.6-73.7度。即本发明确定的INSAT-3DR经纬度优于印度官方发布的定位结果。
(5)制作投影查找表并投影
选择区域为北纬60度到南纬60度,东经10度到东经135度。IR1通道的分辨率为4KM。创建等间隔经纬度数组,读取等间隔经纬图像中4KM的经纬度数组,将相同经纬度的等间隔经纬图像的行列数据和等经纬的行列数据对应起来,制作所需区域的等间隔经纬度投影查找表。然后将标准数据的数组通过等间隔经纬度投影查找表投影为等间隔经纬度数据。
本发明依据静止轨道遥感卫星数据成像的特点,通过行列展宽到再配准的标称投影图完成一个通道数据的定位功能。需根据气象遥感卫星数据判断星下点位置,根据星下点位置,在缺失太阳天顶角和轨道及姿态数据等定位过程依赖的主要数据情况下,生成定位文件,计算卫星天顶角。再根据标称图的纵横坐标通过索引找到展宽图上对应的经纬度,利用定位信息找到对应的亚像元,进行数据插值完成标称投影。
实施例二
为了执行上述实施例一对应的方法,以实现相应的功能和技术效果,下面提供一种基于静止卫星的定位系统。
如图2所示,本实施例提供的基于静止卫星的定位系统包括:数据获取单元1、格式转换单元2及经纬确定单元3。
其中,数据获取单元1用于获取目标卫星的待定位数据所述待定位数据中包括对应区域范围内各点的行号和列号。所述目标卫星为不符合国际静止卫星定位标准的卫星。
格式转换单元2与所述数据获取单元1连接,格式转换单元2用于将所述待定位数据的格式转换为标准格式,得到待定位标准数据。
经纬确定单元3与所述格式转换单元2连接,经纬确定单元3用于基于预先建立的等间隔经纬度投影查找表,根据所述待定位标准数据,确定对应区域范围内各点的最优经度及最优纬度。所述等间隔经纬度投影查找表为预先根据目标卫星的历史数据及标准卫星数据建立的。所述等间隔经纬度投影查找表中包括对应区域范围内各点的行号、列号、最优经度及最优纬度。
进一步地,在建立等间隔经纬度投影查找表方面,本发明基于静止卫星的定位系统还包括:历史数据获取单元4、历史格式转换单元5、初始图像生成单元6、经纬图像生成单元7及查找表生成单元8。
其中,历史数据获取单元4用于获取目标卫星的历史数据及标准卫星数据。所述历史数据的区域范围与所述标准卫星数据的区域范围相同。所述历史数据中包括对应区域范围内各点的行号和列号。所述标准卫星数据中包括对应区域范围内各点的经度和纬度。
历史格式转换单元5与所述历史数据获取单元4连接,历史格式转换单元5用于将所述历史数据的格式转换为标准格式,得到历史标准数据。
初始图像生成单元6与所述历史格式转换单元5连接,初始图像生成单元6用于根据所述历史标准数据及所述目标卫星的轨道信息,生成初始等经纬图像。所述初始等经纬图像中包括对应区域范围内各点的粗略经度及粗略纬度。
经纬图像生成单元7与所述初始图像生成单元6连接,经纬图像生成单元7用于基于所述标准卫星数据,根据所述初始等经纬图像,确定等间隔经纬图像。所述等间隔经纬图像中包括对应区域范围内各点的最优经度、最优纬度、行号及列号。
查找表生成单元8分别与所述经纬图像生成单元7及所述经纬确定单元3连接,查找表生成单元8用于基于所述等间隔经纬图像生成等间隔经纬度投影查找表。
实施例三
本实施例提供一种电子设备,包括存储器及处理器,存储器用于存储计算机程序,处理器运行计算机程序以使电子设备执行实施例一的基于静止卫星的定位方法。
可选地,上述电子设备可以是服务器。
另外,本发明实施例还提供一种计算机可读存储介质,其存储有计算机程序,该计算机程序被处理器执行时实现实施例一的基于静止卫星的定位方法。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种基于静止卫星的定位方法,其特征在于,所述基于静止卫星的定位方法包括:
获取目标卫星的待定位数据;所述待定位数据中包括对应区域范围内各点的行号和列号;所述目标卫星为不符合国际静止卫星定位标准的卫星;
将所述待定位数据的格式转换为标准格式,得到待定位标准数据;
获取目标卫星的历史数据及标准卫星数据;所述历史数据的区域范围与所述标准卫星数据的区域范围相同;所述历史数据中包括对应区域范围内各点的行号和列号;所述标准卫星数据中包括对应区域范围内各点的经度和纬度;
将所述历史数据的格式转换为标准格式,得到历史标准数据;
根据所述历史标准数据及所述目标卫星的轨道信息,生成初始等经纬图像;所述初始等经纬图像中包括对应区域范围内各点的粗略经度及粗略纬度;
基于所述标准卫星数据,根据所述初始等经纬图像,确定等间隔经纬图像,具体包括:从所述初始等经纬图像及所述标准卫星数据中选取多个同名点;确定各同名点在所述初始等经纬图像中的粗略经度及粗略纬度,得到粗略经纬集;确定各同名点在所述标准卫星数据中的经度及纬度,得到标准经纬集;对所述粗略经纬集中各同名点的粗略经度、粗略纬度、所述标准经纬集中各同名点的经度和纬度进行线性拟合,确定缩放系数及偏移系数;根据所述缩放系数及所述偏移系数,计算所述初始等经纬图像中各点的最优经度和最优纬度;根据所述初始等经纬图像中各点的最优经度和最优纬度,生成等间隔经纬图像;所述等间隔经纬图像中包括对应区域范围内各点的最优经度、最优纬度、行号及列号;
基于所述等间隔经纬图像生成等间隔经纬度投影查找表;
基于预先建立的等间隔经纬度投影查找表,根据所述待定位标准数据,确定对应区域范围内各点的最优经度及最优纬度;所述等间隔经纬度投影查找表为预先根据目标卫星的历史数据及标准卫星数据建立的;所述等间隔经纬度投影查找表中包括对应区域范围内各点的行号、列号、最优经度及最优纬度。
2.根据权利要求1所述的基于静止卫星的定位方法,其特征在于,所述目标卫星为印度静止卫星INSAT-3DR。
3.根据权利要求1所述的基于静止卫星的定位方法,其特征在于,所述获取目标卫星的待定位数据,具体包括:
获取目标卫星的原始数据;所述原始数据中包括多个通道数据;
对所述原始数据进行解析,确定各通道数据对应的数组,以得到待定位数据;各通道数据对应的数组中包括对应区域范围内各点的行号和列号。
4.根据权利要求1所述的基于静止卫星的定位方法,其特征在于,所述将所述待定位数据的格式转换为标准格式,得到待定位标准数据,具体包括:
根据标准卫星的数据格式,采用双线性插值法对所述待定位数据进行重采样,得到待定位标准数据。
5.根据权利要求1所述的基于静止卫星的定位方法,其特征在于,所述根据所述历史标准数据及所述目标卫星的轨道信息,生成初始等经纬图像,具体包括:
根据所述历史标准数据及所述目标卫星的轨道信息,生成地理网格;所述地理网格中包括多个网格及各网格对应的粗略经度及粗略纬度;各网格代表对应区域范围内的各点;
根据所述地理网格中各网格对应的粗略经度及粗略纬度,采用双线性插值法,对所述地理网格进行插值,以生成初始等经纬图像。
6.一种基于静止卫星的定位系统,其特征在于,所述基于静止卫星的定位系统包括:
数据获取单元,用于获取目标卫星的待定位数据;所述待定位数据中包括对应区域范围内各点的行号和列号;所述目标卫星为不符合国际静止卫星定位标准的卫星;
格式转换单元,与所述数据获取单元连接,用于将所述待定位数据的格式转换为标准格式,得到待定位标准数据;
经纬确定单元,与所述格式转换单元连接,用于基于预先建立的等间隔经纬度投影查找表,根据所述待定位标准数据,确定对应区域范围内各点的最优经度及最优纬度;所述等间隔经纬度投影查找表为预先根据目标卫星的历史数据及标准卫星数据建立的,具体地,获取目标卫星的历史数据及标准卫星数据;所述历史数据的区域范围与所述标准卫星数据的区域范围相同;所述历史数据中包括对应区域范围内各点的行号和列号;所述标准卫星数据中包括对应区域范围内各点的经度和纬度;将所述历史数据的格式转换为标准格式,得到历史标准数据;根据所述历史标准数据及所述目标卫星的轨道信息,生成初始等经纬图像;所述初始等经纬图像中包括对应区域范围内各点的粗略经度及粗略纬度;基于所述标准卫星数据,根据所述初始等经纬图像,确定等间隔经纬图像,具体包括:从所述初始等经纬图像及所述标准卫星数据中选取多个同名点;确定各同名点在所述初始等经纬图像中的粗略经度及粗略纬度,得到粗略经纬集;确定各同名点在所述标准卫星数据中的经度及纬度,得到标准经纬集;对所述粗略经纬集中各同名点的粗略经度、粗略纬度、所述标准经纬集中各同名点的经度和纬度进行线性拟合,确定缩放系数及偏移系数;根据所述缩放系数及所述偏移系数,计算所述初始等经纬图像中各点的最优经度和最优纬度;根据所述初始等经纬图像中各点的最优经度和最优纬度,生成等间隔经纬图像;所述等间隔经纬图像中包括对应区域范围内各点的最优经度、最优纬度、行号及列号;基于所述等间隔经纬图像生成等间隔经纬度投影查找表;所述等间隔经纬度投影查找表中包括对应区域范围内各点的行号、列号、最优经度及最优纬度。
7.一种电子设备,其特征在于,包括存储器及处理器,所述存储器用于存储计算机程序,所述处理器运行所述计算机程序以使所述电子设备执行权利要求1至5中任一项所述的基于静止卫星的定位方法。
8.一种计算机可读存储介质,其特征在于,其存储有计算机程序,所述计算机程序被处理器执行时实现如权利要求1至5中任一项所述的基于静止卫星的定位方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211267794.9A CN115877420B (zh) | 2022-10-17 | 2022-10-17 | 一种基于静止卫星的定位方法、系统、电子设备及介质 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211267794.9A CN115877420B (zh) | 2022-10-17 | 2022-10-17 | 一种基于静止卫星的定位方法、系统、电子设备及介质 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115877420A CN115877420A (zh) | 2023-03-31 |
CN115877420B true CN115877420B (zh) | 2023-10-10 |
Family
ID=85770460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211267794.9A Active CN115877420B (zh) | 2022-10-17 | 2022-10-17 | 一种基于静止卫星的定位方法、系统、电子设备及介质 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115877420B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101980294A (zh) * | 2010-09-25 | 2011-02-23 | 西北工业大学 | 基于遥感图像的黄河凌汛检测方法 |
CN103413272A (zh) * | 2013-07-22 | 2013-11-27 | 中国科学院遥感与数字地球研究所 | 低空间分辨率多源遥感图像空间一致性校正方法 |
CN106600589A (zh) * | 2016-12-09 | 2017-04-26 | 中国科学院合肥物质科学研究院 | 一种基于海岸线区域遥感图的星载光谱仪指向配准方法 |
CN112419350A (zh) * | 2020-11-20 | 2021-02-26 | 武汉大学 | 基于地物边界信息的遥感影像自动化几何配准方法及系统 |
CN113469896A (zh) * | 2021-05-21 | 2021-10-01 | 贵州师范学院 | 一种提高地球同步轨道卫星对地观测影像几何校正精度的方法 |
CN114972059A (zh) * | 2022-03-18 | 2022-08-30 | 武汉大学 | 基于风云三号卫星成像仪的影像质量提升方法及设备 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI405029B (zh) * | 2009-08-10 | 2013-08-11 | Altek Corp | 影像定位值決定方法、興趣點標記方法與其裝置 |
-
2022
- 2022-10-17 CN CN202211267794.9A patent/CN115877420B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101980294A (zh) * | 2010-09-25 | 2011-02-23 | 西北工业大学 | 基于遥感图像的黄河凌汛检测方法 |
CN103413272A (zh) * | 2013-07-22 | 2013-11-27 | 中国科学院遥感与数字地球研究所 | 低空间分辨率多源遥感图像空间一致性校正方法 |
CN106600589A (zh) * | 2016-12-09 | 2017-04-26 | 中国科学院合肥物质科学研究院 | 一种基于海岸线区域遥感图的星载光谱仪指向配准方法 |
CN112419350A (zh) * | 2020-11-20 | 2021-02-26 | 武汉大学 | 基于地物边界信息的遥感影像自动化几何配准方法及系统 |
CN113469896A (zh) * | 2021-05-21 | 2021-10-01 | 贵州师范学院 | 一种提高地球同步轨道卫星对地观测影像几何校正精度的方法 |
CN114972059A (zh) * | 2022-03-18 | 2022-08-30 | 武汉大学 | 基于风云三号卫星成像仪的影像质量提升方法及设备 |
Non-Patent Citations (4)
Title |
---|
李志乾.基于插值的FY-2静止气象卫星地理定位数据处理方法研究.遥感信息.2012,第120卷(第2期),第1-4节. * |
沈焕锋.ENVI遥感影像处理方法.2009,436-439. * |
王宏博.应用线形插值对风云二数据进行图像定位与投影计算.农业网络信息.2007,全文. * |
钱辉.静止卫星标称投影解析方法及其在FY2-C中的应用.遥感学报.2009,全文. * |
Also Published As
Publication number | Publication date |
---|---|
CN115877420A (zh) | 2023-03-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Toutin | Geometric processing of remote sensing images: models, algorithms and methods | |
Surazakov et al. | Positional accuracy evaluation of declassified Hexagon KH-9 mapping camera imagery | |
CN110111274B (zh) | 一种星载推扫式光学传感器外方位元素定标方法 | |
CN102735216B (zh) | Ccd立体相机三线阵影像数据平差处理方法 | |
CN103413272A (zh) | 低空间分辨率多源遥感图像空间一致性校正方法 | |
CN113900125B (zh) | 星地联合的线阵成像遥感卫星全自主几何定标方法及系统 | |
CN103697864B (zh) | 一种基于大虚拟相机的窄视场双相机影像拼接方法 | |
CN107610054A (zh) | 一种遥感影像数据的预处理方法 | |
CN101900817B (zh) | 一种普适的遥感数据规则格网化方法 | |
CN105571598B (zh) | 一种卫星激光高度计足印相机姿态的测定方法 | |
CN107516291B (zh) | 夜景影像正射纠正处理方法 | |
CN115877420B (zh) | 一种基于静止卫星的定位方法、系统、电子设备及介质 | |
CN111060991B (zh) | 风云静止卫星的晴空辐射产品的生成方法 | |
Dong et al. | Automatic on-orbit geometric calibration framework for geostationary optical satellite imagery using open access data | |
CN114936202B (zh) | 一种极地反照率遥感数据的重构方法、装置及计算机设备 | |
CN111899183A (zh) | 一种基于地理定位数据的遥感图像几何精校正方法 | |
CN113592744B (zh) | 一种适用于高分辨率遥感影像的几何精校正方法 | |
CN109579796A (zh) | 一种投影后影像的区域网平差方法 | |
CN112767454B (zh) | 基于多视角观测sar数据采样分析的叠掩信息补偿方法 | |
Oh et al. | Automated RPCs Bias Compensation for KOMPSAT Imagery Using Orthoimage GCP Chips in Korea | |
CN112213750B (zh) | 高光谱卫片全谱段逐像元成像角度参数处理方法及介质 | |
CN114862723A (zh) | 一种基于测量密集星场的天文望远镜像场畸变标定方法 | |
CN113379648A (zh) | 一种高分七号和资源三号立体影像联合平差方法 | |
CN113029332A (zh) | 卫星云图的预测方法、装置和处理器 | |
CN110307858B (zh) | 一种测绘卫星交汇误差的自适应校正方法 |
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 |