CN103714175B - 一种低分辨率卫星数据产品空间查询方法及装置 - Google Patents

一种低分辨率卫星数据产品空间查询方法及装置 Download PDF

Info

Publication number
CN103714175B
CN103714175B CN201410006747.8A CN201410006747A CN103714175B CN 103714175 B CN103714175 B CN 103714175B CN 201410006747 A CN201410006747 A CN 201410006747A CN 103714175 B CN103714175 B CN 103714175B
Authority
CN
China
Prior art keywords
satellite data
lookup table
spatial
target area
space
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.)
Expired - Fee Related
Application number
CN201410006747.8A
Other languages
English (en)
Other versions
CN103714175A (zh
Inventor
邹巨洪
林明森
邹斌
曾韬
崔松雪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NATIONAL SATELLITE OCEAN APPLICATION SERVICE
Original Assignee
NATIONAL SATELLITE OCEAN APPLICATION SERVICE
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by NATIONAL SATELLITE OCEAN APPLICATION SERVICE filed Critical NATIONAL SATELLITE OCEAN APPLICATION SERVICE
Priority to CN201410006747.8A priority Critical patent/CN103714175B/zh
Publication of CN103714175A publication Critical patent/CN103714175A/zh
Application granted granted Critical
Publication of CN103714175B publication Critical patent/CN103714175B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F16/00Information retrieval; Database structures therefor; File system structures therefor
    • G06F16/10File systems; File servers
    • G06F16/14Details of searching files based on file metadata

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Library & Information Science (AREA)
  • Data Mining & Analysis (AREA)
  • Databases & Information Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Information Retrieval, Db Structures And Fs Structures Therefor (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)

Abstract

本发明公开了一种低分辨率卫星数据产品空间查询方法及装置,涉及微波遥感技术领域,能够在不降低查询结果精度的情况下,提高低分辨率卫星数据产品空间查询的效率。本发明的主要方法包括:以天为单位对卫星数据产品制作一个与之对应的二维网格化的空间查找表;并制作一个与所述空间查找表对应、并且大小相等的文件查找表;当查找目标区域内的卫星数据产品时,查询空间查找表的目标区域,确定目标区域内是否存在卫星数据空间覆盖;若确定空间查找表的目标区域内存在卫星数据空间覆盖,则从文件查找表的目标区域中获取对应卫星数据产品的文件名,进而根据该文件名获取卫星数据产品的文件。本发明主要用于低分辨率卫星数据产品空间查询的过程中。

Description

一种低分辨率卫星数据产品空间查询方法及装置
技术领域
本发明涉及卫星通信领域,尤其涉及微波遥感技术领域,特别是一种低分辨率卫星数据产品空间查询方法及装置。
背景技术
地球遥感定义为使用电磁辐射获取海洋、陆地和大气的信息而没有物理接触要观测的物体、表面或者现象。遥感与船基仪器的观测不同,如海温和风速的观测,船基观测时直接在测量点用温度计和风速计测量;遥感则是从目标的反射或发射的辐射特性间接测量获得的。遥感观测的遥感器可从绑在鸽子身上的简单相机到复杂的多光谱卫星扫描仪。
卫星遥感器不与要观测的目标进行物理接触,它通过接收的辐射强度和频率分布中反演得出目标的物理特性。这个分布依赖于目标本身的辐射特性的和大气传输中的变化等。遥感器接收到的辐射主要有三个来源:表面发射的黑体辐射,反射的太阳辐射和雷达发射信号的后向散射。因此,根据遥感器是否主动发射信号,可分为主动遥感器和被动遥感器。不同的目标或则物理现象在不同的电磁频段响应各有差异,在设计遥感器时,需根据研究的现象的不同选择合适的频段。按遥感器的探测频率,目前的海洋卫星遥感器主要采用光学遥感器、红外遥感器和微波遥感器。观测的现象不同,对遥感器空间分辨率的要求也各有差异。例如,对舰船等目标的识别,需要m级的分辨率,而对海面风场的观测,则通常只需要10km量级的空间分辨率。按遥感器的空间分辨率,又可分为高分辨率卫星遥感器和低分辨率卫星遥感器。所谓高分辨率和低分辨率,是一个相对的概念,为便于区分,这里将m级分辨率的遥感器称为高分辨率遥感器,km级分辨率的遥感器称为低分辨率遥感器。例如,海洋二号微波散射计空间分辨率为25km,则为低分辨率遥感器。
每次卫星过境,按照卫星遥感器的幅宽,会形成一个卫星条带观测数据。高分辨率卫星数据量通常非常巨大,为便于管理和后续处理,将整轨按照幅宽的宽度进行切割,形成一个个的“景”,也就是接近等边的平行四边形。比如一景环境卫星是360km幅宽,也就按360km形成一景,即360km×360km。低分辨率卫星数据量相对较小,通常将绕地球飞行一周的数据作为一个文件存储,形成一个个的“轨”,以海洋二号微波散射计为例,其L2B级海面风场数据产品就是按轨存储。
通常情况下,高分辨率卫星数据产品按景存储,通过每景图像四个角点的经纬度坐标可大致确定该景图像所覆盖的地理范围,从而实现比较方便的空间查询。以资源卫星分发系统为例,用户只需要输入经纬度和时间范围,就可方便的查询满足空间查询条件的卫星数据产品。但是按轨存储的低分辨率卫星每轨数据产品在全球范围内均有分布,不能通过角点坐标确定图像所覆盖的地理区域,如果对每个文件进行遍历的方式查询,查询效率显然很低,这给这类数据产品的空间查询带来了一定的困难。
发明内容
本发明提供了一种低分辨率卫星数据产品空间查询方法及装置,能够在不降低查询结果精度的情况下,提高低分辨率卫星数据产品空间查询的效率。
为达到上述目的,本发明的实施例提供以下技术方案,该技术方案包括:
一方面,本发明的实施例提供一种低分辨率卫星数据产品空间查询方法,包括:
以天为单位对卫星数据产品制作一个与之对应的二维网格化的空间查找表,记录当天的卫星数据空间覆盖情况;并制作一个与所述空间查找表对应、并且大小相等的文件查找表,记录有卫星数据空间覆盖的空间查找表结点所对应的卫星数据产品的文件名;
当查找目标区域内的卫星数据产品时,查询所述空间查找表的所述目标区域,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖;
若确定所述空间查找表的目标区域内存在卫星数据空间覆盖,则从所述文件查找表的所述目标区域中获取对应卫星数据产品的文件名,进而根据所述文件名获取所述卫星数据产品的文件。
另一方面,本发明的实施例还提供一种低分辨率卫星数据产品空间查询装置,制作单元,用于以天为单位对卫星数据产品制作一个与之对应的二维网格化的空间查找表,记录当天的卫星数据空间覆盖情况;并制作一个与所述空间查找表对应、并且大小相等的文件查找表,记录有卫星数据空间覆盖的空间查找表结点所对应的卫星数据产品的文件名;
查找单元,用于当查找目标区域内的卫星数据产品时,查询所述空间查找表的所述目标区域,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖;
获取单元,用于在确定所述空间查找表的目标区域内存在卫星数据空间覆盖时,从所述文件查找表的所述目标区域中获取对应卫星数据产品的文件名,进而根据所述文件名获取所述卫星数据产品的文件。
本发明实施例提供的低分辨率卫星数据产品空间查询方法及装置,以天为单位建立二值化的空间查找表,在进行卫星数据产品查询时,只需根据目标区域的坐标信息,查找空间查找表对应的目标区域的进行查找,便可获得对应的卫星数据产品的文件名,进而获得卫星数据产品,不需要对每个文件进行遍历的空间查询,提高了查询效率。
附图说明
图1示出了本发明实施例提供的一种低分辨率卫星数据产品空间查询方法的流程图;
图2示出了本发明实施例提供的一种以天为单位对卫星数据产品制作一个与之对应的二维网格化的空间查找表和文件查找表的方法流程图;
图3示出了本发明实施例提供的一种根据空间分辨率创建两个二维数组,分别用于记录空间查找表和文件查找表的方法流程图;
图4示出了本发明实施例提供的一种将所述卫星数据映射到所述空间查找表中并赋值的方法流程图;
图5示出了本发明实施例提供的一种查询所述空间查找表的所述目标区域是否存在卫星数据空间覆盖的方法流程图;
图6示出了本发明实施例提供的一种确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖的方法流程图;
图7示出了本发明实施例提供的一种低分辨率卫星数据产品空间查询装置的组成框图;
图8示出了本发明实施例提供的一种制作单元的组成框图;
图9示出了本发明实施例提供的一种查找单元的组成框图。
具体实施方式
本发明的实施例提供一种低分辨率卫星数据产品空间查询方法,如图1所示,该方法包括:
101、以天为单位对卫星数据产品制作一个与之对应的二维网格化的空间查找表,记录当天的卫星数据空间覆盖情况;并制作一个与所述空间查找表对应、并且大小相等的文件查找表,记录有卫星数据空间覆盖的空间查找表结点所对应的卫星数据产品的文件名。
需要说明的是,之所以按天为单位制作空间查找表和文件查找表是因为:首先,卫星资料的文件名上可以查到年月日的信息,且大部分文件能提供的时间精度只能到年月日,不能到小时这个层次;其次,对大部分海洋卫星观测资料来说,对单个卫星的单天数据产品进行网格化,不会出现单个网格结点出现重复观测的情况;如果超过一天,对单个网格结点,会出现超过一个观测结果的情况,而本发明实施例的网格只能记录一个观测结果,这样的话会损失相当多的观测数据,影响查询的效果。
102、当查找目标区域内的卫星数据产品时,查询所述空间查找表的所述目标区域,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖。
103、若确定所述空间查找表的目标区域内存在卫星数据空间覆盖,则从所述文件查找表的所述目标区域中获取对应卫星数据产品的文件名,进而根据所述文件名获取所述卫星数据产品的文件。
进一步的,在以天为单位对卫星数据产品制作一个与之对应的二维网格化的空间查找表并制作一个与所述空间查找表对应、并且大小相等的文件查找表时,本发明实施例提供但不局限于以下方法实现,如图2所示,该方法包括:
201、设定包含年月日的时间,以所述时间为索引,根据文件名查询当天的卫星数据产品。
需要说明的是,卫星数据产品在存储的时候,是分处理级存储且有规范的文件命名规则,根据数据文件名可获取数据处理级别、类型(观测要素)、轨道号、数据观测年月日等信息。以HY-2A卫星散射计L2B级海洋风场数据产品为例,其文件命名规则如下为:H2A_SM2Byyyymmdd_NNNNN.h5。其中:H2A代表HY-2A卫星;SM代表载荷为散射计;2B代表数据处理级别为L2B级;yyyy代表观测开始时间的年;mm代表观测开始时间月;dd代表观测开始时间的天;YYYYMMDD代表产品所覆盖的年月日;NNNNN代表从最南端为起始点的轨道号;.h5代表文件存储类型为hdf5格式。因此,要查询卫星数据产品,只要设定待查时间,并且以该时间为索引,根据文件名查询当天的卫星数据产品即可。
其中,海洋号卫星微波散射计(HY2-SCAT)为我国第一个可业务化运行的微波散射计,于2011年8月搭载海洋二号卫星发射升空。HY2-SCAT主要用于全球海面风场观测,测风风速范围为4~24m/s,风速精度为2m/s或10%;风向测量范围为0~360°,风向精度为±20°。HY2-SCAT工作频率为13.256GHz,采用笔形波束圆锥扫描方式,通过笔形波束以固定仰角围绕天底方向旋转,在卫星平台顺轨方向的运动中形成一定的地面覆盖刈幅;散射计系统包括VV和HH两个极化方式,分别以不同入射角进行观测,在平台的运动过程中对同一分辨单元可获取不同极化方式,不同入射角度的多次后向散射系数(σ0)测量结果,以克服海面风场方向反演的多值模糊问题。其中内波束采用HH极化方式,入射角为41°,对应地面刈幅宽度约为1350km。外波束采用VV极化方式,入射角为48°,对应地面刈幅宽度约为1700km。
HY-2A散射计L2B产品由L2B产品处理软件生成,为经过反演得到的海面风矢量产品。L2B产品数据产品以轨道为单位进行组织,即每个轨道的风矢量测量数据构成一个L2B文件。L2B产品中的每个数据元素都可以通过风矢量单元的行、列号进行索引。L2B风矢量单元行的延伸方向与星下线相垂直,列的延伸方向与星下线方向相一致。
202、设定空间分辨率,根据空间分辨率创建两个二维数组,分别用于记录空间查找表和文件查找表,并对所述空间查找表和文件查找表赋初值。
其中,在所述根据空间分辨率创建两个二维数组,分别用于记录空间查找表和文件查找表时,可以采用但不局限于以下方法,如图3所示,该方法包括:
301、根据所述空间分辨率,将全球划分为M°×N°的经纬网格,行方向是纬度方向,列方向是经度方向。
以海洋二号卫星散射计为例,其空间分辨率为25km,则可将全球划分为0.25°×0.25°的经纬网格,行方向是纬度方向,列方向是经度方向,全球的网格包含720行,1440列。其中,每个网格的中心经纬度可由下式求出:
Lat[i]=(180/IGRID)×(i+0.5)-90 (公式1)
Lon[j]=(360/JGRID)×(j+0.5) (公式2)
其中,IGRID、JGRID分别为网格的行数、列数。
其中,在对所述空间查找表和文件查找表赋初值时,一般会对数组的每个元素赋初值为0,这样方便以后填写新的数据值,但是本发明实施例对此不进行限制,也可以赋于其他的初值。
302、创建与所述经纬网格对应的大小为M×N的两个二维数组,分别用于记录空间查找表和文件查找表。
基于步骤301的经纬度网格划分,创建与等经纬网对应的大小为720×1440的两个二维数组,分别用于存储空间查找表和文件查找表。
203、将所述卫星数据产品映射到所述空间查找表中并赋值,并将所述卫星数据产品对应的文件名录入所述文件查找表。
其中,在将所述卫星数据映射到所述空间查找表中并赋值,并将所述卫星数据产品对应的文件名录入所述文件查找表时,可以通过但不局限于以下的方法实现,如图4所示,该方法包括:
401、读入以天为单位的卫星数据产品,并记录输入读入的卫星数据产品的文件名。
以海洋二号卫星散射计为例,读入单轨海洋二号卫星散射计L2级数据产品,并记录输入数据产品的文件名。
需要说明的是,对单天的数据产品,如果分升轨和降轨,除了在极地地区,其它地区单天只有一次观测结果。对海洋遥感卫星来说,除了专门针对极地应用的卫星,大部分的观测区域都到不了两极地区。同时,对除了极地应用之外,其它大部分海洋相关应用都不需要极地地区的观测数据。那么HY-2A散射计L2级数据产品文件是以轨道为单位进行组织,一个轨道文件包括,卫星平台围绕其空间轨道旋转一周所获取的全部数据。当卫星由最南端向最北端运行时,对应的轨道为升轨,由最北端向最南端运行时,对应的轨道为降轨。本发明根据这一约定,将HY-2A散射计L2级数据产品文件的前半部分数据读出,以构成升轨文件,后半部分数据则构成降轨文件。
402、获取读入的卫星数据产品单个数据单元的经纬度信息。
读入L2级产品数据产品中单个风矢量单元的经纬度信息和风速、风向、时间(时分秒)信息,并进行质量控制,若未能通过质量控制,则跳过后面步骤,进行下一个风矢量单元的遍历。
403、根据所述经纬度信息计算所述数据单元在所述空间查找表中对应的行列号。
利用下式,计算风矢量单元在查找表网格空间中对应的行列号;
I_index=ROUND((lat+90.)×IGRID/180-0.5) (公式3)
J_index=ROUND(lon×JGRID/360-0.5) (公式4)
404、根据所述行列号将对应的空间查找表单元赋值为1,并将所述卫星数据产品的文件名计入所述文件查找表对应的单元。
重复上述步骤401-404操作,直到单天所有数据产品均已读入。
重新设定时间,重复上述步骤201-203操作,直到所有的L2B级数据产品文件均已生成与之对应的空间查找表和文件查找表。
进一步的,在查询所述空间查找表的所述目标区域,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖时,可以采用但不局限于以下方法实现,如图5所示,该方法包括:
501、根据所述目标区域设定所述目标区域的角点坐标;将该角点坐标标记为lat_min、lat_max、lon_min和lon_max,分别代表纬度最小值、纬度最大值和经度最小值与经度最大值。
502、根据所述角点坐标获得所述目标区域的行列号范围。
其中,在根据所述角点坐标获得所述目标区域的行列号范围时,可以采用,但不局限于通过以下公式实现,该公式包括:
i_min=ROUND((lat_min+90.)×IGRID/180-0.5) (公式5)
i_max=ROUND((lat_max+90.)×IGRID/180-0.5) (公式6)
j_min=ROUND(lon_min×JGRID/360-0.5) (公式7)
j_max=ROUND(lon_max×JGRID/360-0.5) (公式8)
其中i_min、i_max、j_min和j_max分别代表网格行号最小值、网格行号最大值和网格列号最小值与网格列号最大值。
503、读入单天的空间查找表,根据所述行列号范围,对读入的空间查找表相应行列号范围内的数据进行计算,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖。
其中,在对读入的空间查找表相应行列号范围内的数据进行计算,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖时,可以采用但不局限于以下方法,如图6所示,该方法包括:
601、对读入的空间查找表相应行列号范围内的数据进行求和运算。
602、和运算后所得和是否大于0;若求和后大于0,则执行603;否则执行604。
其中,当有卫星数据产品在相应表格赋值为1时,求和后大于0,则说明在目标区域内当天有观测结果,即有卫星数据产品;否则说明在目标区域内当天没有观测结果,即没有卫星数据产品。
603、确定所述空间查找表的目标区域内存在卫星数据空间覆盖。
604、确定所述空间查找表的目标区域内不存在卫星数据空间覆盖。
进一步的,本发明的实施例还提供一种低分辨率卫星数据产品空间查询装置,如图7所示,该装置包括:
制作单元71,用于以天为单位对卫星数据产品制作一个与之对应的二维网格化的空间查找表,记录当天的卫星数据空间覆盖情况;并制作一个与所述空间查找表对应、并且大小相等的文件查找表,记录有卫星数据空间覆盖的空间查找表结点所对应的卫星数据产品的文件名。
查找单元72,用于当查找目标区域内的卫星数据产品时,查询所述空间查找表的所述目标区域,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖。
获取单元73,用于在确定所述空间查找表的目标区域内存在卫星数据空间覆盖时,从所述文件查找表的所述目标区域中获取对应卫星数据产品的文件名,进而根据所述文件名获取所述卫星数据产品的文件。
进一步的,如图8所示,所述制作单元71包括:
查询模块711,用于设定包含年月日的时间,以所述时间为索引,根据文件名查询当天的数据产品。
创建模块712,用于设定空间分辨率,根据空间分辨率创建两个二维数组,分别用于记录空间查找表和文件查找表,并对所述空间查找表和文件查找表赋初值。其中,所述创建模块712用于根据空间分辨率创建两个二维数组,分别用于记录空间查找表和文件查找表时,可以通过但不局限于以下的方式实现,具体包括:根据所述空间分辨率,将全球划分为M°×N°的经纬网格,行方向是纬度方向,列方向是经度方向;创建与所述经纬网格对应的大小为M×N的两个二维数组,分别用于记录空间查找表和文件查找表。
操作模块713,用于将所述卫星数据产品映射到所述空间查找表中并赋值,并将所述卫星数据产品对应的文件名录入所述文件查找表。其中,所述操作模块713用于将所述卫星数据产品映射到所述空间查找表中并赋值,并将所述卫星数据产品对应的文件名录入所述文件查找表时,可以通过但不局限于以下的方式实现,具体包括:读入以天为单位的卫星数据产品,并记录输入所述读入的卫星数据产品的文件名;获取读入的卫星数据产品单个数据单元的经纬度信息;根据所述经纬度信息计算所述数据单元在所述空间查找表中对应的行列号;根据所述行列号将对应的空间查找表单元赋值为1,并将所述卫星数据产品的文件名计入所述文件查找表对应的单元。
进一步的,如图9所示,所述查找单元72包括:
设置模块721,用于根据所述目标区域设定所述目标区域的角点坐标。
获取模块722,用于根据所述角点坐标获得所述目标区域的行列号范围。
读取模块723,用于读入单天的空间查找表。
查询模块724,用于根据所述行列号范围,对读入的空间查找表相应行列号范围内的数据进行计算,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖。其中,所述查询模块724具体用于,对读入的空间查找表相应行列号范围内的数据进行求和运算;若求和后大于0,则确定所述空间查找表的目标区域内存在卫星数据空间覆盖;否则确定所述空间查找表的目标区域内不存在卫星数据空间覆盖。
需要说明的是,本发明实施例提供的低分辨率卫星数据产品空间查询装置的组成单元及模块的其他描述,请参考方法部分的相关描述,本发明实施例此处将不再赘述。
本发明实施例提供的低分辨率卫星数据产品空间查询方法及装置,与现有技术相比,现有技术对单天数据产品的空间查询,需要对单天的所有产品文件进行遍历(以HY-2A散射计L2B级风场产品空间查询为例,在最坏的情况下,需要对单天所有的HY-2A散射计L2B级数据产品进行遍历。本发明实施例只需要对查找表指定范围内的结点进行一次求和操作,就能够得出查询结果,提高了查询效率。
并且,与现有技术相比,现有技术在计算观测点是否在待查区域内时,需要根据观测点的经纬度以及待查区域的坐标信息进行位置判断,运算比较复杂。本发明提供的实施例采用二值化查找表技术,只需通过待查区域的坐标信息,进行查找表的行列号范围计算一次操作,不需要对每个观测点都进行位置判断操作,提高了查询效率。
采用本发明的技术,可以在不降低空间查询精度的情况下,提高空间查询的效率。由于海洋遥感数据处理量巨大,多为海量数据,本发明的技术可以提高数据查询的效率两个数量级以上。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种低分辨率卫星数据产品空间查询方法,其特征在于,包括:
以天为单位对卫星数据产品制作一个与之对应的二维网格化的空间查找表,记录当天的卫星数据空间覆盖情况;并制作一个与所述空间查找表对应、并且大小相等的文件查找表,记录有卫星数据空间覆盖的空间查找表结点所对应的卫星数据产品的文件名;
当查找目标区域内的卫星数据产品时,查询所述空间查找表的所述目标区域,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖;
若确定所述空间查找表的目标区域内存在卫星数据空间覆盖,则从所述文件查找表的所述目标区域中获取对应卫星数据产品的文件名,进而根据所述文件名获取所述卫星数据产品的文件;
根据空间分辨率创建两个二维数组,分别用于记录所述空间查找表和所述文件查找表,对所述空间查找表和所述文件查找表赋初值;
根据所述空间分辨率,将全球划分为M°×N°的经纬网格,行方向是纬度方向,列方向是经度方向;
创建与所述经纬网格对应的大小为M×N的所述两个二维数组。
2.根据权利要求1所述的低分辨率卫星数据产品空间查询方法,其特征在于,所述以天为单位对卫星数据产品制作一个与之对应的二维网格化的空间查找表,记录当天的卫星数据空间覆盖情况;并制作一个与所述空间查找表对应、并且大小相等的文件查找表,记录有卫星数据空间覆盖的空间查找表结点所对应的卫星数据产品的文件名,包括:
设定包含年月日的时间,以所述时间为索引,根据文件名查询当天的数据产品;
设定空间分辨率,根据空间分辨率创建两个二维数组,分别用于记录空间查找表和文件查找表,并对所述空间查找表和文件查找表赋初值;
将所述卫星数据产品映射到所述空间查找表中并赋值,并将所述卫星数据产品对应的文件名录入所述文件查找表。
3.根据权利要求2所述的低分辨率卫星数据产品空间查询方法,其特征在于,所述根据空间分辨率创建两个二维数组,分别用于记录空间查找表和文件查找表包括:
根据所述空间分辨率,将全球划分为M°×N°的经纬网格,行方向是纬度方向,列方向是经度方向;
创建与所述经纬网格对应的大小为M×N的两个二维数组,分别用于记录空间查找表和文件查找表。
4.根据权利要求3所述的低分辨率卫星数据产品空间查询方法,其特征在于,将所述卫星数据映射到所述空间查找表中并赋值,并将所述卫星数据产品对应的文件名录入所述文件查找表包括:
读入以天为单位的卫星数据产品,并记录输入读入的卫星数据产品的文件名;
获取读入的卫星数据产品单个数据单元的经纬度信息;
根据所述经纬度信息计算所述数据单元在所述空间查找表中对应的行列号;
根据所述行列号将对应的空间查找表单元赋值为1,并将所述卫星数据产品的文件名计入所述文件查找表对应的单元。
5.根据权利要求1-4中任一项所述的低分辨率卫星数据产品空间查询方法,其特征在于,查询所述空间查找表的所述目标区域,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖包括:
根据所述目标区域设定所述目标区域的角点坐标;
根据所述角点坐标获得所述目标区域的行列号范围;
读入单天的空间查找表,根据所述行列号范围,对读入的空间查找表相应行列号范围内的数据进行计算,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖。
6.根据权利要求5所述的低分辨率卫星数据产品空间查询方法,其特征在于,对读入的空间查找表相应行列号范围内的数据进行计算,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖包括:
对读入的空间查找表相应行列号范围内的数据进行求和运算;
若求和后大于0,则确定所述空间查找表的目标区域内存在卫星数据空间覆盖;否则确定所述空间查找表的目标区域内不存在卫星数据空间覆盖。
7.一种低分辨率卫星数据产品空间查询装置,其特征在于,包括:
制作单元,用于以天为单位对卫星数据产品制作一个与之对应的二维网格化的空间查找表,记录当天的卫星数据空间覆盖情况;并制作一个与所述空间查找表对应、并且大小相等的文件查找表,记录有卫星数据空间覆盖的空间查找表结点所对应的卫星数据产品的文件名;
查找单元,用于当查找目标区域内的卫星数据产品时,查询所述空间查找表的所述目标区域,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖;
获取单元,用于在确定所述空间查找表的目标区域内存在卫星数据空间覆盖时,从所述文件查找表的所述目标区域中获取对应卫星数据产品的文件名,进而根据所述文件名获取所述卫星数据产品的文件;
根据空间分辨率创建两个二维数组,分别用于记录所述空间查找表和所述文件查找表,对所述空间查找表和所述文件查找表赋初值;
根据所述空间分辨率,将全球划分为M°×N°的经纬网格,行方向是纬度方向,列方向是经度方向;
创建与所述经纬网格对应的大小为M×N的所述两个二维数组。
8.根据权利要求7所述的低分辨率卫星数据产品空间查询装置,其特征在于,所述制作单元包括:
查询模块,用于设定包含年月日的时间,以所述时间为索引,根据文件名查询当天的数据产品;
创建模块,用于设定空间分辨率,根据空间分辨率创建两个二维数组,分别用于记录空间查找表和文件查找表,并对所述空间查找表和文件查找表赋初值;
操作模块,用于将所述卫星数据产品映射到所述空间查找表中并赋值,并将所述卫星数据产品对应的文件名录入所述文件查找表。
9.根据权利要求7或8所述的低分辨率卫星数据产品空间查询装置,其特征在于,所述查找单元包括:
设置模块,用于根据所述目标区域设定所述目标区域的角点坐标;
获取模块,用于根据所述角点坐标获得所述目标区域的行列号范围;
读取模块,用于读入单天的空间查找表;
查询模块,用于根据所述行列号范围,对读入的空间查找表相应行列号范围内的数据进行计算,确定所述空间查找表的目标区域内是否存在卫星数据空间覆盖。
10.根据权利要求9所述的低分辨率卫星数据产品空间查询装置,其特征在于,所述查询模块具体用于,
对读入的空间查找表相应行列号范围内的数据进行求和运算;
若求和后大于0,则确定所述空间查找表的目标区域内存在卫星数据空间覆盖;否则确定所述空间查找表的目标区域内不存在卫星数据空间覆盖。
CN201410006747.8A 2014-01-07 2014-01-07 一种低分辨率卫星数据产品空间查询方法及装置 Expired - Fee Related CN103714175B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410006747.8A CN103714175B (zh) 2014-01-07 2014-01-07 一种低分辨率卫星数据产品空间查询方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410006747.8A CN103714175B (zh) 2014-01-07 2014-01-07 一种低分辨率卫星数据产品空间查询方法及装置

Publications (2)

Publication Number Publication Date
CN103714175A CN103714175A (zh) 2014-04-09
CN103714175B true CN103714175B (zh) 2018-01-09

Family

ID=50407150

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410006747.8A Expired - Fee Related CN103714175B (zh) 2014-01-07 2014-01-07 一种低分辨率卫星数据产品空间查询方法及装置

Country Status (1)

Country Link
CN (1) CN103714175B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105159895A (zh) 2014-05-28 2015-12-16 国际商业机器公司 用于存储和查询数据的方法和系统
CN105676191B (zh) * 2016-01-08 2016-11-30 国家卫星海洋应用中心 基于海洋微波遥感技术的散射计数据转化方法及装置
CN107656238B (zh) * 2017-09-12 2020-10-30 中国海洋大学 一种基于双查找表方法的高频地波雷达测向新方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6393293B1 (en) * 2000-01-21 2002-05-21 Motorola, Inc. Method for rectangular parameterization of a geographic area using a geolocation algorithm
CN101339570A (zh) * 2008-08-12 2009-01-07 北京航空航天大学 一种高效的海量遥感数据分布式组织管理方法
CN101388043A (zh) * 2008-09-26 2009-03-18 北京航空航天大学 一种基于小块图片的ogc高性能遥感图像地图服务方法
CN101477748A (zh) * 2009-01-19 2009-07-08 深圳市同洲电子股份有限公司 交通状况的管理方法及系统
CN101515411A (zh) * 2009-03-26 2009-08-26 郑经涛 一种地图
CN103390115A (zh) * 2013-08-02 2013-11-13 国家卫星海洋应用中心 一种海洋卫星遥感观测数据匹配方法和系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6393293B1 (en) * 2000-01-21 2002-05-21 Motorola, Inc. Method for rectangular parameterization of a geographic area using a geolocation algorithm
CN101339570A (zh) * 2008-08-12 2009-01-07 北京航空航天大学 一种高效的海量遥感数据分布式组织管理方法
CN101388043A (zh) * 2008-09-26 2009-03-18 北京航空航天大学 一种基于小块图片的ogc高性能遥感图像地图服务方法
CN101477748A (zh) * 2009-01-19 2009-07-08 深圳市同洲电子股份有限公司 交通状况的管理方法及系统
CN101515411A (zh) * 2009-03-26 2009-08-26 郑经涛 一种地图
CN103390115A (zh) * 2013-08-02 2013-11-13 国家卫星海洋应用中心 一种海洋卫星遥感观测数据匹配方法和系统

Also Published As

Publication number Publication date
CN103714175A (zh) 2014-04-09

Similar Documents

Publication Publication Date Title
Unwin et al. Spaceborne GNSS-reflectometry on TechDemoSat-1: Early mission operations and exploitation
Fore et al. Combined active/passive retrievals of ocean vector wind and sea surface salinity with SMAP
Rajmohan et al. Revamping land coverage analysis using aerial satellite image mapping
Jiang et al. The HY-2 satellite and its preliminary assessment
Liu et al. Tropical Rainfall Measuring Mission (TRMM) precipitation data and services for research and applications
Moon et al. RADARSAT-2 and coastal applications: Surface wind, waterline, and intertidal flat roughness
Moreira et al. A drone-borne multiband DInSAR: Results and applications
CN103714175B (zh) 一种低分辨率卫星数据产品空间查询方法及装置
Ermida et al. Inversion of AMSR‐E observations for land surface temperature estimation: 2. Global comparison with infrared satellite temperature
Hong et al. Self‐organizing nonlinear output (SONO): A neural network suitable for cloud patch–based rainfall estimation at small scales
Koma et al. Better together? Assessing different remote sensing products for predicting habitat suitability of wetland birds
Perwitasari et al. Coordinated airglow observations between IMAP/VISI and a ground‐based all‐sky imager on concentric gravity wave in the mesopause
Ramapriyan The role and evolution of NASA's Earth science data systems
Wang et al. All-weather near-surface air temperature estimation based on satellite data over the Tibetan Plateau
Wright et al. Estimating GPS signal loss in a natural deciduous forest using sky photography
Amici et al. Geological classification of Volcano Teide by hyperspectral and multispectral satellite data
Zhao et al. Potential soil moisture product from the Chinese HY-2 scanning microwave radiometer and its initial assessment
Aulicino et al. Remote Sensing Applications in Satellite Oceanography
CN105676191B (zh) 基于海洋微波遥感技术的散射计数据转化方法及装置
Kalluri et al. A hierarchical data archiving and processing system to generate custom tailored products from AVHRR data
Estes et al. Characteristics, sources, and management of remotely-sensed data
Lyu et al. Estimating Geo‐Referenced Cloud‐Base Height With Whole‐Sky Imagers
Duffourg et al. Simulation of satellite infrared radiances for convective‐scale data assimilation over the Mediterranean
Papadopoulos et al. Spatio-temporal evolution of the deformation around the Rio-Patras fault (Greece) observed by synthetic aperture radar interferometry from 1993 to 2017
Li et al. Downscaling GNSS-R based vegetation water content product using random forest model

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20180109

Termination date: 20190107

CF01 Termination of patent right due to non-payment of annual fee