CN104166695A - 一种面向遥感数据内容的剖分、查询及动态显示方法 - Google Patents
一种面向遥感数据内容的剖分、查询及动态显示方法 Download PDFInfo
- Publication number
- CN104166695A CN104166695A CN201410373099.XA CN201410373099A CN104166695A CN 104166695 A CN104166695 A CN 104166695A CN 201410373099 A CN201410373099 A CN 201410373099A CN 104166695 A CN104166695 A CN 104166695A
- Authority
- CN
- China
- Prior art keywords
- subdivision
- dough sheet
- geosot
- remotely
- image data
- 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.)
- Pending
Links
Classifications
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Theoretical Computer Science (AREA)
- Remote Sensing (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明公开了面向遥感数据内容的剖分、查询及动态显示方法,具体为:针对遥感数据,采用设定的剖分层级剖分并建立gsot文件;剖分获得的每一剖分面片,将其GeoSOT编码作为索引号、该剖分面片内的影像数据在遥感数据中的位置信息以及该影像数据的属性内容对应存入gsot文件;获取查询区域以及待查询的属性值,剖分获得查询区域对应的查询面片及查询编码;当一个索引号与一个查询编码由左至右匹配一致且长度大于该查询编码,同时该索引号对应属性值等于待查询的属性值时,提取该索引号对应影像数据。同时本发明还使用上述查询方法进行了基于属性信息的遥感数据的动态显示。
Description
技术领域
本发明属于地球空间信息、遥感数据组织管理技术领域,具体涉及基于地球剖分的遥感数据内容查询方法及动态显示方法。
背景技术
遥感技术为人类研究地球提供了大量的珍贵图像资料,为充分、有效地利用这些资料,实现海量遥感数据的高效共享和发布,需要建立合理的数据查询机制,为遥感数据用户能更便捷地使用遥感数据提供保障。
目前,浏览查询方式是遥感数据应用服务***中最重要的查询方式,其主要功能是向各类用户提供查询检索已编目存档的遥感数据和各级产品数据的有效手段。它根据图像和产品的编目信息提供多种查询内容,如基于时间的查询、基于空间位置的查询和基于地图的查询等,用户可以通过查询界面设置查询条件查询到自己所需要的数据信息。然而,由于这些查询内容都是针对数据文件本身的查询,当用户查询某热点区域内的属性信息时,首先需要查询遥感数据的编目信息,自动、半自动甚至人工地判断该热点是否在此遥感数据文件中,若存在则打开原始数据,分析热点区域对应的属性信息;否则,则放弃。
上述传统的基于编目信息的查询方式,虽然可以很好地适应大部分用户对遥感数据查询的需要,但是在许多特殊的应用领域,需要查询遥感数据的内部某些属性信息时,传统方法较为耗时。下面举一个路径规划的例子
在该应用中,用户不仅仅只是对矢量数据的道路信息进行规划,还需要利用遥感数据,生成沿途的影像条带景观图。以北京到深圳为例,在这过程中,用户需要向遥感数据的服务部门查询覆盖沿途区域的所有影像,以便将这些影像拼成两千多公里长的路径影像条带。覆盖沿途区域的影像中不乏一些局部被云层遮蔽或质量不好的影像,这些影像对于路径规划是没有价值的,在条带拼接前需要将它们剔除,因此沿途区域未被云层覆盖是另一个查询条件。据相关遥感数据部门的统计,在卫星获取的各类遥感数据中,被云层覆盖的遥感影像平均可以达到40%左右,甚至很多区域能达到70%以上。而现有的数据服务方式虽然对遥感数据内云的信息进行了描述,但也仅仅是记录云的含量及百分比。这样的方式,无法对云信息进行定位,即使云量很少,也可能覆盖了感兴趣区域;反之,云量很大,感兴趣区域也未必有云覆盖。
基于现有方式查询时,首先根据影像的编目信息查询覆盖感兴趣区域的影像,然后打开该影像,查看感兴趣区域是否被云层覆盖,若未被覆盖则将整幅影像传输给用户,否则继续判断下一幅影像。如此复杂繁琐的查询过程,对于前面提到的自北京至深圳的影像请求问题,至少需要耗时几天才可完成,且传输给用户的影像往往具有很大的冗余,其覆盖范围大大超出用户关心区域,对数据传输出造成了极大的压力,严重影响查询和利用效率,使得目前的遥感数据查询方法在面对类似问题的情况下,无法有效地给出方案。
因此,采用传统的查询方式,由于在遥感数据编目中缺少内部区域属性信息的描述,并且传统的数据标识方式(例如:坐标表示,元数据表示等)也很难对遥感数据内容进行标识和索引,这使得在遥感数据中查询内容信息的效率极低,从而造成大量遥感数据闲置,给资源应用、发布带来困难。
发明内容
有鉴于此,本发明提供了一种面向遥感数据内容的剖分预处理方法,该方法通过网格统一记录、管理遥感数据的局部属性信息,并通过对查询区域进行相应的剖分,以剖分编码的匹配来实现查询,提高了数据查询的效率,丰富了数据查询的结果。
为达到上述目的,本方法包括如下步骤:
步骤1、针对遥感数据,采用GeoSOT剖分和编码方案中设定的剖分层级剖分该遥感数据并建立剖分索引gsot文件;对于剖分获得的每一剖分面片,将该剖分面片的GeoSOT编码作为索引号、该剖分面片内的影像数据在遥感数据中的位置信息以及该影像数据的属性对应存入剖分索引gsot文件。
步骤2、获取查询区域以及待查询的属性值,使用GeoSOT剖分和编码方案获得查询区域对应的剖分面片和GeoSOT编码,即查询面片及查询编码。
步骤3、当一个索引号与一个查询编码由左至右匹配一致且长度大于该查询编码,同时该索引号对应属性值等于待查询的属性值时,利用剖分索引gsot文件中的位置信息提取该索引号对应影像数据。
进一步地,遥感数据为遥感影像,则步骤1包括以下步骤:
S101、获取遥感影像的经纬度范围;
S102、使用GeoSOT剖分方案,选择设定的剖分层级对遥感影像进行剖分,获得剖分面片的集合;
S103、使用GeoSOT编码方案对S102获得的每个剖分面片进行编码,记为索引号;
S104、确定S102获得的每个剖分面片内的影像数据的四个角点在所述遥感数据中的行列号,依据所述行列号以数据统计或者目视判读的方法确定该行列号范围内的各所需属性类型的属性值;
S105、将属于同一剖分面片的索引号、行列号以及属性值对应写入剖分索引gsot文件中。
进一步地,S102中所述设定的剖分层级的选择方法为:用户指定S102获得的每个剖分面片内的影像数据的大小,使用遥感影像分辨率以及上述影像数据的大小计算理想剖分面片的尺寸,在GeoSOT剖分方案中找到最接近理想剖分面片尺寸的剖分面片,其所属层级即为设定的剖分层级。
进一步地,步骤2中的查询区域由n个经纬度点{(xi,yi)}组成,i=1,2,…n,其中xi和yi分别为第i个经纬度点的经度和纬度,则使用GeoSOT剖分和编码方案的设定剖分层级剖分所述查询区域,获得查询面片的具体步骤为:
S201、计算所述查询区域的最小外包矩形R;
S202、计算外包矩形的长和宽,并取长和宽中的较小值M;
S203、在GeoSOT剖分方案中选取剖分层级j对应的网格尺度θj,使得θj≤M<θj+1;
S204、在GeoSOT剖分方案的第j层级中,找到覆盖所述查询区域的所有面片,从而获得查询面片的集合。
在上述剖分预处理方法的基础上,本发明同时提供了一种动态显示方法,该方法能够通过查询获得属性变化区域,从而使得遥感数据的上球显示更为及时有效,并节省了实时显示占用的资源。
为达到上述目的,本方法包括如下步骤:
步骤一、针对待显示地理区域,其时刻t的遥感影像为Pt,建立Pt的剖分索引gsot文件Wt,并将Pt上图显示。
剖分索引gsot文件Wt是采用GeoSOT剖分和编码方案中的剖分层级A对Pt进行剖分获得的,其中Wt中存有每一剖分面片的GeoSOT编码即索引号、该剖分面片中内的影像数据在Pt中的位置信息以及该影像数据的属性。
步骤二、待显示地理区域在t+i时刻的遥感影像为Pt+i,采用步骤一相同的方式建立Pt+i的剖分索引gsot文件Wt+i。
步骤三、在Wt和Wt+i中,同一索引号对应的属性值若存在不同,则利用Wt+i记录的位置信息,在Pt+i中提取该索引号对应影像数据。
步骤四、在t+i时刻,使用步骤三获得的影像数据对已上图显示的Pt进行更新。
进一步地,步骤一中建立Pt的剖分索引gsot文件Wt的具体过程为:
S1、使用GeoSOT剖分方案中的剖分层级A对Pt进行剖分,获得剖分面片的集合;
S2、使用GeoSOT编码方案对S1获得的每个剖分面片进行编码,记为索引号;
S3、确定S1获得的每个剖分面片内的影像数据的四个角点在Pt中的行列号,依据所述行列号以数据统计或者目视判读的方法确定该行列号范围内的各所需属性类型的各属性值;
S4、将属于同一剖分面片的索引号、行列号以及属性值对应写入Wt中。
进一步地,S1中所述剖分层级A的选择方法为:用户指定S1获得的每个剖分面片内的影像数据的大小,使用遥感影像分辨率以及上述影像数据的大小计算理想剖分面片的尺寸,在GeoSOT剖分方案中找到最接近理想剖分面片尺寸的剖分面片,其所属层级即为剖分层级A。
本发明提供了一种面向遥感数据内容的剖分预处理方法,该方法通过网格统一记录、管理遥感数据的局部属性信息。
为达到上述目的,本方法包括如下步骤:
S101、获取遥感影像的经纬度范围。
S102、使用GeoSOT剖分方案,选择设定的剖分层级对遥感影像进行剖分,获得剖分面片的集合。
S103、使用GeoSOT编码方案对S102获得的每个剖分面片进行编码,记为索引号。
S104、确定S102获得的每个剖分面片内的影像数据的四个角点在所述遥感数据中的行列号,依据所述行列号以数据统计或者目视判读的方法确定该行列号范围内的各所需属性类型的属性值。
S105、将属于同一剖分面片的索引号、行列号以及属性值对应写入剖分索引gsot文件中。
有益效果:
1、本发明所提供的一种面向遥感数据内容的剖分预处理方法,首先利用网格逻辑索引统一记录、管理遥感数据的内容信息,提前将数据内容信息记录下来,从而实现了对遥感数据的同一记录和管理,为快速查询相应属性的遥感数据以开展后续上图显示或者其他工作提供了便利。
2、本发明所提供的基于剖分网格的遥感数据的内容信息查询方法,首先利用网格逻辑索引统一记录、管理遥感数据的内容信息,提前将数据内容信息记录下来,以便于快速获取查询结果,而无须打开原始影像进行判读,特别是针对热点区域的数据内容信息查询效果显著,从而提高遥感数据的查询浏览效率,为后续提高遥感影像数据的内部信息查询效率提供有力的支持。
3、本发明所提供的基于剖分网格的遥感数据的内容信息查询方法,将遥感数据和查询区域均采用具体方法进行了合理化剖分,从而能够确保查询到的用户感兴趣区域数据进行小块传输,而不是包含感兴趣区域的整幅影像传输,为查询后的数据传输减轻了压力,提高了发布共享的效率。
4、本发明通过对于上述属性信息剖分预处理以及查询方法所得效果的应用,提出了一种动态显示方法:针对经过上述预处理后的遥感数据,利用不同时期遥感数据各部分属性信息的变化,查询获得属性变化区域,从而使得遥感数据的上球显示更为及时有效,并节省了实时显示占用的资源。
附图说明
图1是本发明的遥感数据查询方法流程图;
图2是本发明的待处理的遥感影像图;
图3是本发明的气象云图覆盖面片的属性标识图;
图4是本发明的气象云图gsot索引文件结构;
图5是本发明的计算查询区域覆盖面片集的原理图;
图6为本发明的动态显示方法流程图;
图7:(a)为实施例2中的某地2000年土地利用分类图;(b)为实施例2中的某地2010年土地利用分类图;
图8:(a)为实施例2中的某地2000年土地利用分类图覆盖面片的属性值图;(b)为实施例2中的某地2010年土地利用分类图覆盖面片的属性值图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明所使用的GeoSOT编码是基于GeoSOT剖分和编码方案,该方案参见北京大学提出的专利申请:“一种统一现有经纬度剖分网格的方法”(公开号为CN102609525,公开日为2012年7月25日),该专利申请公开了一种GeoSOT地理网格设计方案,用于解决全球地理空间剖分和标识问题。该方案采用全四叉树递归剖分,将地球表面空间从全球至厘米级共进行了32级剖分,每个GeoSOT剖分层级均有其对应大小的GeoSOT网格,GeoSOT网格上下层级之间的面积之比是1/4。该方案对GeoSOT网格进行编码所产生的GeoSOT编码有四进制、二进制一维、二进制二维和十进制四种形式,其中这四种不同进制之间的GeoSOT编码可以两两之间进行转换。
为了能够对带有内容信息的遥感数据进行统一记录和管理,可以采用上述GeoSOT剖分和编码方案,将遥感数据进行剖分并建立剖分索引gsot文件,剖分索引gsot文件应当记录剖分获得的剖分面片的GeoSOT编码以及剖分面片中对应的影像数据和影像数据中的属性值。其中属性值即为遥感数据内容信息。这种gsot文件为以后的查询与其他功能应用提供了极大的方便。例如以下实施例1和实施例2即为该种gsot文件的具体应用。
实施例1、
本实施例是一种基于地球剖分的遥感数据内容查询方法,如图1所示,本方法包括以下步骤:
步骤1、针对遥感数据,采用GeoSOT剖分和编码方案中设定的剖分层级剖分该遥感数据并建立剖分索引gsot文件;对于剖分获得的每一剖分面片,将该剖分面片的GeoSOT编码作为索引号、该剖分面片内的影像数据在遥感数据中的位置信息以及该影像数据的属性对应存入剖分索引gsot文件。
如图2所示为本实施例处理的一幅遥感影像,本步骤所建立的gsot文件可以记录该幅遥感影像的各种属性类型,例如,本实施例中一种属性列dcloud记录了网格中的云覆盖信息,以数字0、1分别记录某网格内无云、有云的信息,即dcloud的属性值为1时代表该网格范围有云,dcloud的属性值为0时代表该网格范围无云。
该步骤中建立gsot文件时采用上述实施例1所提供的步骤,其中具体为:
S101、获取遥感影像的经纬度范围;
S102、使用GeoSOT剖分方案,选择设定的剖分层级对遥感影像进行剖分,获得剖分面片的集合。
设定的剖分层级的选择方法为:用户指定S102获得的每个剖分面片内的影像数据的大小,使用遥感影像分辨率以及上述影像数据的大小计算理想剖分面片的尺寸,在GeoSOT剖分方案中找到最接近理想剖分面片尺寸的剖分面片,其所属层级即为设定的剖分层级。
S103、使用GeoSOT编码方案对S102获得的每个剖分面片进行编码,记为索引号;
S104、确定S102获得的每个剖分面片内的影像数据的四个角点在所述遥感数据中的行列号,依据所述行列号以数据统计或者目视判读的方法确定该行列号范围内的各所需属性类型的属性值。如图2所示为某地区的气象云图,其中只有一种属性类型,即云dcloud,以数字0、1分别记录某区域内无云、有云的信息,即dcloud的属性值为1时代表该行列号范围有云,dcloud的属性值为0时代表该行列号范围无云。
S105、将属于同一剖分面片的索引号、行列号以及属性值对应写入剖分索引gsot文件中;在本实施例中,将索引号、行列号以及dcloud和其他扩展字段组成如图4所示的结构写入剖分索引gsot文件中。
步骤2、获取查询区域以及待查询的属性值,使用GeoSOT剖分和编码方案获得查询区域对应的剖分面片和GeoSOT编码,即查询面片及查询编码。
本实施例中,本步骤中的查询区域由n个经纬度点{(xi,yi)}组成,i=1,2,…n,其中xi和yi分别为第i个经纬度点的经度和纬度,则使用GeoSOT剖分和编码方案的设定剖分层级剖分所述查询区域,则计算查询面片时可以参照图5,分为如下具体步骤:
S201、计算查询区域的最小外包矩形R,R中的经纬度点(X、Y)满足min(xi)≤X≤max(xi),min(yi)≤Y≤max(yi)。
S202、计算外包矩形的长和宽,并取长和宽中的较小值M;M=min{(max(xi)-min(xi)),(max(yi)-min(yi))}。本步骤主要考虑到当所处理的查询区域为条带状区域时,以条带的宽选取剖分层级。
S203、在GeoSOT剖分方案中选取剖分层级j对应的网格尺度θj,使得θj≤M<θj+1;
S204、在GeoSOT剖分方案的第j层级中,找到覆盖上述查询区域的所有面片,从而获得查询面片的集合。这样对于一个条带状区域,当以条带的宽来选取剖分层级时,在该步骤中的查询面片集合组成的区域为覆盖上述查询区域且最小。
步骤3、当一个索引号与一个查询编码由左至右匹配一致且长度大于该查询编码,同时该索引号对应属性值等于待查询的属性值时,利用剖分索引gsot文件中的位置信息提取该索引号对应影像数据。
在本步骤中,每个查询编码均与所述剖分索引gsot文件中的各索引号由左至右进行匹配:
若查询编码与其中一个索引号由左至右完全匹配一致,且查询编码的长度小于该索引号,则判断该索引号对应的属性值是否等于待查询的属性值,若是则提取该索引号对应的行列号,否则不进行提取,所获得的行列号进行汇总,并依据该行列号在遥感影像中获取影像数据。
例如本次查询所计算的查询编码为G0013100121,待查询的属性值dcloud=1,即在图2所示的气象云图中查询编码为G0013100121的查询面片区域内找到有云的区域,则对该气象云图对应的gsot索引文件逐行匹配。G0013100121所包含的区域具有以下性质:前11(10+1)位与G0013100121完全匹配,长度大于或等于11。经查询,得到16个索引号及其对应的dCloud值,如表1所示。
表1
Code | dcloud |
G001310012100 | 1 |
G001310012101 | 1 |
G001310012102 | 1 |
G001310012103 | 1 |
G001310012110 | 1 |
G001310012111 | 1 |
G001310012112 | 1 |
G001310012113 | 1 |
G001310012120 | 1 |
G001310012121 | 1 |
G001310012122 | 1 |
G001310012123 | 1 |
G001310012130 | 1 |
G001310012131 | 1 |
G001310012132 | 1 |
G001310012133 | 1 |
如表1所示,查询区域全部被云层覆盖。
可见,本发明通过基于地球剖分的遥感数据内容查询方法将遥感数据基于其地理范围生成索引gsot文件,并能够利用网格统一记录、管理遥感数据的局部属性信息,在统一遥感数据组织的基础上,能够快速获取查询结果,特别是针对热点区域的内容信息查询效果显著,从而提高遥感数据的查询浏览效率。
实施例3、
针对上述实施例1提供的基于地球剖分的遥感数据内容查询方法,该方法经改进可以适用于对遥感数据进行动态显示。因此本实施例提供了一种间接使用上述查询方法的遥感数据的动态显示方法,该方法流程如图6所示,包括如下步骤:
步骤一、针对待显示地理区域,其时刻t的遥感影像为Pt,建立Pt的剖分索引gsot文件Wt,并将Pt上图显示;
剖分索引gsot文件Wt是采用GeoSOT剖分和编码方案中的剖分层级A对Pt进行剖分获得的,其中Wt中存有每一剖分面片的GeoSOT编码即索引号、该剖分面片中内的影像数据在Pt中的位置信息以及该影像数据的属性;具体的Wt的建立方式可以采用实施例1所提供的方法。
例如,图7(a)和(b)分别为某地2000年、2010年的土地利用分类图,本实施例要对一段时间内土地的植被覆盖变化进行显示,则使用实施例1中步骤1的方法建立第一剖分索引gsot文件时,其属性类型为植被覆盖情况,属性值为L(代表裸地)、G(代表植被)和B(代表建筑)。
步骤二、待显示地理区域在t+i时刻的遥感影像为Pt+i,采用步骤一相同的方式建立Pt+i的剖分索引gsot文件Wt+i。
在本实施例中,仅有2000年以及2010年的影像数据,因此该处起始时刻t设为2000年,i为10年。则对于两幅图像进行剖分并对比如图8(a)和(b)所示。
步骤三、在Wt和Wt+i中,同一索引号对应的属性值若存在不同,则利用Wt+i记录的位置信息,在Pt+i中提取该索引号对应影像数据;。
即对于任一索引号,将Wt+i与Wt中对应该索引号的属性值进行比较,若存在不同,在Wt+i中提取该索引号对应的行列号,否则不进行提取。本实施例中,在图8中(a)和(b)两幅图中不一致的面片会被提取出来。所获得的行列号进行汇总,并依据该行列号在遥感影像中获取影像数据。
步骤四、在t+i时刻,使用步骤三获得的影像数据对已上图显示的Pt进行更新。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种面向遥感数据内容的查询方法,其特征在于,包括如下步骤:
步骤1、针对遥感数据,采用GeoSOT剖分和编码方案中设定的剖分层级剖分该遥感数据并建立剖分索引gsot文件;对于剖分获得的每一剖分面片,将该剖分面片的GeoSOT编码作为索引号、该剖分面片内的影像数据在遥感数据中的位置信息以及该影像数据的属性对应存入剖分索引gsot文件;
步骤2、获取查询区域以及待查询的属性值,使用GeoSOT剖分和编码方案获得查询区域对应的剖分面片和GeoSOT编码,即查询面片及查询编码;
步骤3、当一个索引号与一个查询编码由左至右匹配一致且长度大于该查询编码,同时该索引号对应属性值等于待查询的属性值时,利用剖分索引gsot文件中的位置信息提取该索引号对应影像数据。
2.如权利要求1所述的一种面向遥感数据内容的查询方法,其特征在于,所述遥感数据为遥感影像,则步骤1包括以下步骤:
S101、获取遥感影像的经纬度范围;
S102、使用GeoSOT剖分方案,选择设定的剖分层级对遥感影像进行剖分,获得剖分面片的集合;
S103、使用GeoSOT编码方案对S102获得的每个剖分面片进行编码,记为索引号;
S104、确定S102获得的每个剖分面片内的影像数据的四个角点在所述遥感数据中的行列号,依据所述行列号以数据统计或者目视判读的方法确定该行列号范围内的各所需属性类型的属性值;
S105、将属于同一剖分面片的索引号、行列号以及属性值对应写入剖分索引gsot文件中。
3.如权利要求2所述的一种面向遥感数据内容的查询方法,其特征在于,S102中所述设定的剖分层级的选择方法为:用户指定S102获得的每个剖分面片内的影像数据的大小,使用遥感影像分辨率以及上述影像数据的大小计算理想剖分面片的尺寸,在GeoSOT剖分方案中找到最接近理想剖分面片尺寸的剖分面片,其所属层级即为设定的剖分层级。
4.如权利要求1或2所述的一种面向遥感数据内容的查询方法,其特征在于,所述步骤2中的查询区域由n个经纬度点{(xi,yi)}组成,i=1,2,…n,其中xi和yi分别为第i个经纬度点的经度和纬度,则使用GeoSOT剖分和编码方案的设定剖分层级剖分所述查询区域,获得查询面片的具体步骤为:
S201、计算所述查询区域的最小外包矩形R;
S202、计算外包矩形的长和宽,并取长和宽中的较小值M;
S203、在GeoSOT剖分方案中选取剖分层级j对应的网格尺度θj,使得θj≤M<θj+1;
S204、在GeoSOT剖分方案的第j层级中,找到覆盖所述查询区域的所有面片,从而获得查询面片的集合。
5.一种面向遥感数据内容的动态显示方法,其特征在于,该方法包括如下步骤:
步骤一、针对待显示地理区域,其时刻t的遥感影像为Pt,建立Pt的剖分索引gsot文件Wt,并将Pt上图显示;
所述剖分索引gsot文件Wt是采用GeoSOT剖分和编码方案中的剖分层级A对Pt进行剖分获得的,其中Wt中存有每一剖分面片的GeoSOT编码即索引号、该剖分面片中内的影像数据在Pt中的位置信息以及该影像数据的属性;
步骤二、待显示地理区域在t+i时刻的遥感影像为Pt+i,采用步骤一相同的方式建立Pt+i的剖分索引gsot文件Wt+i;
步骤三、在Wt和Wt+i中,同一索引号对应的属性值若存在不同,则利用Wt+i记录的位置信息,在Pt+i中提取该索引号对应影像数据;
步骤四、在t+i时刻,使用步骤三获得的影像数据对已上图显示的Pt进行更新。
6.如权利要求5所述的一种面向遥感数据内容的动态显示方法,其特征在于,所述步骤一中建立Pt的剖分索引gsot文件Wt的具体过程为:
S1、使用GeoSOT剖分方案中的剖分层级A对Pt进行剖分,获得剖分面片的集合;
S2、使用GeoSOT编码方案对S1获得的每个剖分面片进行编码,记为索引号;
S3、确定S1获得的每个剖分面片内的影像数据的四个角点在Pt中的行列号,依据所述行列号以数据统计或者目视判读的方法确定该行列号范围内的各所需属性类型的各属性值;
S4、将属于同一剖分面片的索引号、行列号以及属性值对应写入Wt中。
7.如权利要求6所述的一种面向遥感数据内容的动态显示方法,其特征在于,S1中所述剖分层级A的选择方法为:用户指定S1获得的每个剖分面片内的影像数据的大小,使用遥感影像分辨率以及上述影像数据的大小计算理想剖分面片的尺寸,在GeoSOT剖分方案中找到最接近理想剖分面片尺寸的剖分面片,其所属层级即为剖分层级A。
8.一种面向遥感数据内容的剖分预处理方法,其特征在于,该方法包括如下步骤:
S101、获取遥感影像的经纬度范围;
S102、使用GeoSOT剖分方案,选择设定的剖分层级对遥感影像进行剖分,获得剖分面片的集合;
S103、使用GeoSOT编码方案对S102获得的每个剖分面片进行编码,记为索引号;
S104、确定S102获得的每个剖分面片内的影像数据的四个角点在所述遥感数据中的行列号,依据所述行列号以数据统计或者目视判读的方法确定该行列号范围内的各所需属性类型的属性值;
S105、将属于同一剖分面片的索引号、行列号以及属性值对应写入剖分索引gsot文件中。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410373099.XA CN104166695A (zh) | 2014-07-31 | 2014-07-31 | 一种面向遥感数据内容的剖分、查询及动态显示方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410373099.XA CN104166695A (zh) | 2014-07-31 | 2014-07-31 | 一种面向遥感数据内容的剖分、查询及动态显示方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104166695A true CN104166695A (zh) | 2014-11-26 |
Family
ID=51910508
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410373099.XA Pending CN104166695A (zh) | 2014-07-31 | 2014-07-31 | 一种面向遥感数据内容的剖分、查询及动态显示方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104166695A (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108984584A (zh) * | 2018-05-23 | 2018-12-11 | 北京都在哪网讯科技有限公司 | 地理信息处理方法及装置 |
CN109241217A (zh) * | 2018-08-14 | 2019-01-18 | 民政部国家减灾中心 | 一种基于地理网格的灾情数据空间统计、查询及展示方法 |
CN109634507A (zh) * | 2018-12-03 | 2019-04-16 | 广东国图勘测地理信息有限公司 | 一种可触控的电子地图控制方法及装置 |
CN110110016A (zh) * | 2018-01-23 | 2019-08-09 | 北京云游九州空间科技有限公司 | 一种遥感影像分析方法及装置 |
CN110245248A (zh) * | 2019-05-31 | 2019-09-17 | 河南大学 | 一种遥感影像关键字查询方法 |
CN110909093A (zh) * | 2019-10-11 | 2020-03-24 | 中国人民解放军战略支援部队信息工程大学 | 智慧地标控制网的构建方法及装置 |
CN111125392A (zh) * | 2019-12-25 | 2020-05-08 | 华中科技大学 | 一种基于矩阵对象存储机制的遥感影像存储、查询方法 |
CN112364188A (zh) * | 2020-11-12 | 2021-02-12 | 中国地质大学(武汉) | 一种遥感影像的索引建立方法、遥感影像检索方法及装置 |
CN113656398A (zh) * | 2021-07-06 | 2021-11-16 | 中铁建港航局集团勘察设计院有限公司 | 面向航道信息数据内容的网格动态显示方法、***及介质 |
CN117033526A (zh) * | 2023-10-09 | 2023-11-10 | 中国地质大学(武汉) | 数据存储方法、数据查询方法、装置、设备及存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103106284A (zh) * | 2013-03-01 | 2013-05-15 | 北京大学 | 用于三维数据球信息关联的剖分中间件和信息关联方法 |
CN103136371A (zh) * | 2013-03-21 | 2013-06-05 | 北京大学 | 多源空间数据的剖分标识生成方法以及数据检索方法 |
CN103150338A (zh) * | 2013-02-04 | 2013-06-12 | 北京大学 | 一种用于影像数据提取的剖分预处理方法及数据提取方法 |
CN103488736A (zh) * | 2013-09-18 | 2014-01-01 | 中国科学技术大学 | 一种多源地理空间信息关联模型的构建方法及*** |
-
2014
- 2014-07-31 CN CN201410373099.XA patent/CN104166695A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103150338A (zh) * | 2013-02-04 | 2013-06-12 | 北京大学 | 一种用于影像数据提取的剖分预处理方法及数据提取方法 |
CN103106284A (zh) * | 2013-03-01 | 2013-05-15 | 北京大学 | 用于三维数据球信息关联的剖分中间件和信息关联方法 |
CN103136371A (zh) * | 2013-03-21 | 2013-06-05 | 北京大学 | 多源空间数据的剖分标识生成方法以及数据检索方法 |
CN103488736A (zh) * | 2013-09-18 | 2014-01-01 | 中国科学技术大学 | 一种多源地理空间信息关联模型的构建方法及*** |
Non-Patent Citations (2)
Title |
---|
吕雪峰,等: "基于GeoSOT区位标识的多源遥感数据组织研究", 《北京大学(自然科学版)》 * |
吕雪峰,等: "基于GeoSOT区位标识的多源遥感数据组织研究", 《北京大学(自然科学版)》, vol. 50, no. 2, 31 March 2014 (2014-03-31) * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110110016B (zh) * | 2018-01-23 | 2021-11-09 | 北京国星宇航科技有限公司 | 一种遥感影像分析方法及装置 |
CN110110016A (zh) * | 2018-01-23 | 2019-08-09 | 北京云游九州空间科技有限公司 | 一种遥感影像分析方法及装置 |
CN108984584A (zh) * | 2018-05-23 | 2018-12-11 | 北京都在哪网讯科技有限公司 | 地理信息处理方法及装置 |
CN109241217A (zh) * | 2018-08-14 | 2019-01-18 | 民政部国家减灾中心 | 一种基于地理网格的灾情数据空间统计、查询及展示方法 |
CN109241217B (zh) * | 2018-08-14 | 2021-08-20 | 民政部国家减灾中心 | 一种基于地理网格的灾情数据空间统计、查询及展示方法 |
CN109634507A (zh) * | 2018-12-03 | 2019-04-16 | 广东国图勘测地理信息有限公司 | 一种可触控的电子地图控制方法及装置 |
CN110245248A (zh) * | 2019-05-31 | 2019-09-17 | 河南大学 | 一种遥感影像关键字查询方法 |
CN110245248B (zh) * | 2019-05-31 | 2022-10-21 | 河南大学 | 一种遥感影像关键字查询方法 |
CN110909093A (zh) * | 2019-10-11 | 2020-03-24 | 中国人民解放军战略支援部队信息工程大学 | 智慧地标控制网的构建方法及装置 |
CN111125392A (zh) * | 2019-12-25 | 2020-05-08 | 华中科技大学 | 一种基于矩阵对象存储机制的遥感影像存储、查询方法 |
CN112364188A (zh) * | 2020-11-12 | 2021-02-12 | 中国地质大学(武汉) | 一种遥感影像的索引建立方法、遥感影像检索方法及装置 |
CN113656398A (zh) * | 2021-07-06 | 2021-11-16 | 中铁建港航局集团勘察设计院有限公司 | 面向航道信息数据内容的网格动态显示方法、***及介质 |
CN117033526A (zh) * | 2023-10-09 | 2023-11-10 | 中国地质大学(武汉) | 数据存储方法、数据查询方法、装置、设备及存储介质 |
CN117033526B (zh) * | 2023-10-09 | 2023-12-29 | 中国地质大学(武汉) | 数据存储方法、数据查询方法、装置、设备及存储介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104166695A (zh) | 一种面向遥感数据内容的剖分、查询及动态显示方法 | |
CN112115198B (zh) | 一种城市遥感智能服务平台 | |
Salvati et al. | Exploring long-term land cover changes in an urban region of southern Europe | |
CN101388043B (zh) | 一种基于小块图片的ogc高性能遥感图像地图服务方法 | |
CN107085600B (zh) | Poi推荐方法、装置、设备及计算机可读存储介质 | |
US20160283096A1 (en) | Method of generating a link by utilizing a picture and system thereof | |
CN107977407A (zh) | 瓦片地图存储方法、装置、存储介质和计算机设备 | |
CN101339570A (zh) | 一种高效的海量遥感数据分布式组织管理方法 | |
CN106484693A (zh) | 地图显示方法、图片块存储方法以及图片显示终端 | |
CN104252489A (zh) | 一种根据经纬度数据快速获得位置文字描述信息的方法 | |
US20190221023A1 (en) | Determining and Presenting Solar Flux Information | |
CN106933833A (zh) | 一种基于空间索引技术的位置信息快速查询方法 | |
Goerlich Gisbert et al. | Clustering cities through urban metrics analysis | |
CN103150338A (zh) | 一种用于影像数据提取的剖分预处理方法及数据提取方法 | |
CN104182475B (zh) | 一种基于掩码技术和剖分编码的位置信息快速检索方法 | |
CN104182472A (zh) | 一种基于内存数据库Redis的土地利用矢量数据存储方法 | |
CN103425770A (zh) | 事件多维度信息显示装置和方法 | |
CN104217414A (zh) | 用于影像拼接的拼接线提取方法和装置 | |
CN102831605B (zh) | 用于大型影像的多尺度图像分割结果快速可视化方法 | |
CN105023089B (zh) | 基于gis平台的城市管理数据监测***及其方法 | |
CN112199528B (zh) | 一种大规模遥感数据在线获取方法 | |
CN109241217A (zh) | 一种基于地理网格的灾情数据空间统计、查询及展示方法 | |
Salvati | “A Chronicle of a Death Foretold”: Urban Expansion and Land Consumption in Rome, Italy | |
CN111797083A (zh) | 一种规范化的基础地理数据制作方法及*** | |
CN105654137A (zh) | 一种海量遥感数据分类组织方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20141126 |