CN106248346A - 一种基于几何拓扑学的高精度遥感卫星对地覆盖分析方法 - Google Patents
一种基于几何拓扑学的高精度遥感卫星对地覆盖分析方法 Download PDFInfo
- Publication number
- CN106248346A CN106248346A CN201510321035.XA CN201510321035A CN106248346A CN 106248346 A CN106248346 A CN 106248346A CN 201510321035 A CN201510321035 A CN 201510321035A CN 106248346 A CN106248346 A CN 106248346A
- Authority
- CN
- China
- Prior art keywords
- time
- satellite
- geometric
- target area
- achieve
- 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
Abstract
一种基于几何拓扑学的高精度遥感卫星对地覆盖分析方法。本发明利用几何拓扑学中的几何拓扑运算方法求得遥感卫星对地覆盖过程中每次过境覆盖区域及覆盖分析完毕后总的过境覆盖区域及百分比,比传统的网格点法复杂度要低且精度要高。本发明通过记录在过境过程中幅宽两点的轨迹数据集,利用轨迹数据集组成的几何图形与目标区域几何图形的交集求得每次过境的覆盖区域,以及历次过境覆盖区域几何图形求并得到总覆盖区域,通过记录的每次过境开始时间、结束时间等求得条带数据量,并根据每次过境覆盖区域与条带区域的百分比与条带数据量的乘积求得有效的数据量。通过本发明提出的新技术方案,势必为以后遥感卫星对地覆盖产生巨大影响。
Description
技术领域
本发明涉及航天遥感卫星对地覆盖分析领域,尤指一种利用几何拓扑关系运算进行遥感卫星对地覆盖效能综合分析的新方法。
背景技术
随着空间技术的发展,遥感卫星的数量越来越多,遥感卫星数据在各领域的应用也越来越广,然而,按照“天地统筹、协同建设”的思想,航天遥感卫星系统的发展势必对地面系统和应用系统的建设提出了更高的要求,建设多大规模与处理能力的地面系统和应用系统成为了亟需论证和解决的问题。航天遥感卫星对地覆盖效能分析正是解决这一问题的有效途径和方法。在卫星发射之前利用计算机仿真的方法对卫星的覆盖能力进行有效的分析和评估可以预测卫星的数据获取质量、规模、品种、规格和时效性,为地面系统和应用系统建设指标提供充分的依据。但是,现有的卫星覆盖能力解析法和数值计算法多存在计算精度不高、计算资源占用较大以及卫星多次过境重叠度无法计算等问题。本发明采用几何拓扑学中的几何拓扑运算方法,有效地解决了上述问题,大大提高了遥感卫星对地覆盖综合分析的效率。
发明内容
针对现有卫星对地覆盖方法精度不高,时间复杂度和空间复杂度较大的问题,本发明提出了一种基于几何拓扑学的高精度遥感卫星对地覆盖效能分析方法。该方法利用JTS(Java Transaction Service)拓扑套件中几何拓扑运算的方法,在全球二维平面直角坐标系中将目标分析区域理解成为一个几何图形A(A中可能存在多个非连通的闭合区域),同时,按照一定的步长,将卫星每一飞行时刻幅宽左右两点组成一条线段B,通过计算线段B到几何图形A的距离distance来判断卫星是否处于过境状态。如下(1)为JTS中建立B到A的距离关系,(2)为求得B到A的距离distance。
DistanceOp distOp=new DistanceOp(A,B); (1)
double distance=distOp.distance(); (2)
如果距离distance大于0,则表示线段B与目标分析区域几何图形A不相交,即处于非过境状态;如果距离distance等于0,则表示线段B与几何图形A相交或在其内部,即处于过境状态,记录此时刻幅宽左右两点的经纬度坐标(lat,lon)。随着卫星的飞行,逐步形成了一个从开始过境目标分析区域A到最后离开目标分析区域A的以幅宽左右两点集合为边界点的封闭区域。在全球二维平面直角坐标系下将封闭区域的经纬度(lat,lon)边界点集合转换为二维直角坐标系下的(x,y)点集合并按照顺时针或者逆时针方向将其进行首尾相接构成了卫星每一次过境轨迹的几何图形,进而通过几何拓扑求交运算得到该轨迹几何图形与目标分析区域几何图形的交集即为本次过境对目标分析区域的有效覆盖范围。
本发明提出了一种遥感卫星对地覆盖分析方法,其特征在于,包括以下步骤:
a.构建目标区域几何图形;
b.卫星过境覆盖分析;
c.覆盖数据处理及统计。
进一步,其特征在于,利用的是几何拓扑关系运算进行遥感卫星对地覆盖效能综合分析。
进一步,其特征在于,步骤a:目标区域可以为规则区域、非规则区域和多个不相连的规则或非规则区域。
进一步,其特征在于,步骤b:采用了外接矩形过境判断和目标分析区域过境判断相结合的方法,目的是为了记录每一次过境的最后结束时间,利用幅宽两点组成线段是否离开目标区域外接矩形来判断每一次的过境结束,利用幅宽两点组成线段进入外接矩形后首次进入目标区域判断为每次过境的开始,当存在过境中断时,利用标识flag记录每次离开目标区域的时间,并将最后一次时间记录为过境结束时间,并从每次过境开始到过境结束记录幅宽两点的轨迹数据集。
进一步,其特征在于,步骤c:利用每次过境的幅宽两点轨迹数据集构成的几何图形与目标区域几何图形求得的交集作为每次过境的覆盖区域,利用每次覆盖区域组成的几何图形求得的并集作为总覆盖 区域。
本发明的基于几何拓扑学的遥感卫星对地覆盖分析方法具有以下优点:
1、利用拓扑几何的线面位置关系进行过境判断和覆盖轨迹数据处理,相比传统的网格覆盖分析方法中通过建立全球平面底图直角坐标系,将幅宽两点投影到坐标系中,建立直线方程,通过判断与直线相交的属于目标行政区域的网格来判断卫星瞬时幅宽所扫过的网格来进行过境判断和过境覆盖数据处理从而简化了算法的处理流程,极大的降低了时间复杂度。
2、利用拓扑几何的求交运算计算过境轨迹几何图形和目标区域几何图形的交集得到每一次对目标分析区域的有效覆盖范围,比传统网格点法(通过计算网格点的个数而导致部分非覆盖区域也被算作覆盖区域)的精度更高。
3、传统网格点法进行覆盖处理运算的过程中,网格覆盖次数、是否覆盖等大量网格数据都将存储在内存中,消耗大量的内存资源,而利用拓扑关系进行过境覆盖处理只需要记录边缘点数据即可,处理的数据也是各种边缘点数据,相对来说处理和存储的数据量大大减少,降低了内存资源的消耗。
4、通过矩形区域过境和目标区域过境相结合的方法,规避了卫星恰巧过境目标分析区域存在凹形部分或不连续部分的特殊情况,这是矩形网格算法中所不能的。
附图说明
图 1幅宽两点计算图
图 2基于几何拓扑的遥感卫星对地覆盖分析方法流程图
图 3目标区域几何图形A构建(以中国行政区域为例)
图 4瞬时过境判断图 (以中国行政区域为例)
图 5卫星对地覆盖图
图 6网格法与几何拓扑对地覆盖法对比图
具体实施方式
本发明的基于几何拓扑学的高精度遥感卫星对地覆盖分析方法包括目标区域几何图形构建、卫星过境覆盖区域计算和覆盖数据处理与分析。
进一步,所述目标区域几何图形构建,包括收集目标区域边界经纬度点(lat,lon)集合数据(按照边界顺时针或逆时针方向收集)和将数据转化到二维地图坐标系下,并根据二维地图坐标系下的点集合数据构建几何图形。
进一步,所述卫星过境覆盖区域计算,包括幅宽两点计算、利用(1)、(2)进行过境状态判断和过境时间段内幅宽两点坐标数据的记录。
其中如上图 a 1=sinicosβ、a2=cosisinβ、星下点经纬度幅宽两点R、L的计算公式如下:
纬度公式:
纬度公式:
经度λL公式:
经度λR公式:
当轨道倾角小于90°,卫星沿Z轴飞行方向大于0时,(5)式取“-”,(6)式取“+”,卫星沿Z轴飞行方向小于0时,(5)式取“+”,(6)式取“-”;当轨道倾角大于90°,卫星沿Z轴飞行方向大于0时,(5)式取“+”,(6)式取“-”,当轨道倾角大于90°,卫星沿Z轴飞 行方向小于0时,(5)式取“-”,(6)式取“+”。
进一步,所述覆盖数据处理及分析,包括利用JTS根据每次过境幅宽两点轨迹数据集合构建的几何图形A与目标区域几何图形B求得交集几何图形C=A.intersects(B)即每次过境的覆盖区域、利用交集几何图形的面积及目标区域几何图形面积求得每次过境覆盖百分比、利用JTS根据历次过境交集C1、C2、C3……几何图形求得并集几何图形D=(C1.union(C2)).union(C3)……即为最终的全部覆盖区域、根据全部覆盖区域D几何图形面积及目标几何图形面积求得总覆盖百分比、根据历次过境的结束时间、开始时间、卫星曝光时间、分辨率、幅宽长度、量化位数和波段数通过公式(7)求得每次过境的条带数据量和根据每次过境的覆盖面积、条带面积、求得的条带数据量通过公式(8)求得每次过境的有效数据量,其中条带数据量和有效数据量的公式如下:
条带数据量:Dband=((Tend-Tstart)/Tsensor)*Lbreadth/resolution*Q*B (7)
Tend过境结束时间、Tstart过境开始时间、Tsensor曝光时间、Lbreadth幅宽长度、Q量化位数、B波段数
有效数据量:Dvalid=(Svaild/Sband)*Dband (8)
Svaild有效过境面积、Sband过境条带面积。
如图 1所示,本发明的方法在实施中需要目标区域几何图形构建、卫星过境覆盖区域计算和覆盖数据处理与分析三个步骤。
如图 2所示,目标区域几何图形构建,以中国行政区域为例,首先分开收集台湾、海南、内陆等互不相连区域的数据,构成List<double[]>的数据集,每一个double数组为一个区域的数据集,这些数据集点顺序为顺时针或逆时针方向,然后将这些数据集利用JTS构建为多个闭合的几何图形Geometry。
如图 3所示,卫星过境判断图,在卫星过境覆盖区域计算的过程中,以中国行政区域为例,在图中点A处是过境中断地方,定义:卫星在绕地球一圈的过程中从幅宽两点组成的线段T首次进入目标区域到最终离开目标区域为一次过境,因此由于目标区域的不规则性,就会出现图中点A处的情况(线段T在一次过境中存在离开了目标区域又再次进入目标区域),点A处的情况如果不加以处理,则经过点 A处的此次过境将会由于算法中设定的从进入目标区域到离开目标区域为一次过境而被默认为两次过境,故本发明中利用目标区域外接矩形过境判断和目标区域过境判断相结合的方法来处理点A处情况,首先由于目标区域外接矩形为规则区域,不会出现类似点A处的情况,故将算法中设定为:如果幅宽两点组成的线段连续经过外接矩形区域且离开了外接矩形区域,则一次过境结束,过境次数会加1,但此时如果记录当前时刻为过境结束时间就是错误的,因此还需要对目标区域在一次过境中何时最后一次出境进行判断,首先在算法中设定标识flag为true,当幅宽两点组成的线段T进入到目标区域时,即变为false,一旦线段T离开目标区域时,判断flag是否为false,如果是则记录此刻的时间Ltime且设置flag为true,通过标识flag的设置,即可有效的判断在一次过境中每次离开目标区域时的时刻,当线段T离开目标区域外接矩形时即可将最后一次离开目标区域的时刻记录为此次过境结束时间,在卫星过境的同时会记录从开始过境时刻到结束过境时刻每经过一个步长卫星幅宽两点的值,形成每次过境的幅宽轨迹数据集,最终形成历次过境的开始时间与结束时间集合和幅宽轨迹数据集。
如图 4所示,卫星对地覆盖图,通过对覆盖数据处理与分析,最后形成了卫星对地覆盖图,首先根据在卫星过境覆盖分析过程中获得的每次过境的幅宽轨迹数据集,将其构建为一个封闭的几何图形,然后与目标区域几何图形求交集,获得的几何图形即为此次过境的覆盖区域,由图 5可以看到几何拓扑法比传统的网格点法求得的覆盖区域更为准确,然后再对每次过境的覆盖区域几何图形求并集所得的几何图形即为总覆盖区域,用总覆盖区域几何图形的面积及目标区域几何图形的面积即为求出总覆盖百分比,然后根据历次过境的结束时间、开始时间、卫星曝光时间、分辨率、幅宽长度、量化位数和波段数通过公式(7)求得每次过境的条带数据量,然后根据每次过境的覆盖面积、条带面积、求得的条带数据量通过公式(8)求得每次过境的有效数据量,从而形成卫星对地的综合覆盖特性。
以上示例,仅为说明本发明的技术特征和可实施例,其目的是在于使该领域的技术人员能够了解本发明的内容并具体实施。由此,凡 根据本发明的构思作出的变换和修饰,均包含在本发明的权利要求范围内,依专利法提出申请。
Claims (5)
1.一种遥感卫星对地覆盖分析方法,其特征在于,包括以下步骤:
a.构建目标区域几何图形;
b.卫星过境覆盖分析;
c.覆盖数据处理及统计。
2.如权利要求1所述的方法,其特征在于利用的是几何拓扑关系进行覆盖分析而非传统的网格。
3.如权利要求1所述的方法,其特征在于,步骤a:目标区域可以为规则区域、非规则区域和多个不相连的规则或非规则区域。
4.如权利要求1所述的方法,其特征在于,步骤b:采用了目标分析区域过境判断和其外接矩形过境判断相结合的方法,目的是为了规避卫星轨迹恰好经过了目标分析区域的凹形区域或不连续区域而判断为卫星本次过境已经离开了目标分析区域的特殊情况,如图3中A处所示,利用卫星每一时刻载荷幅宽两点构成的线段是否离开目标区域外接矩形来判断每一次的过境结束与否,同时,利用该线段进入外接矩形并首次进入目标区域作为该次过境的开始,当存在过境中断时,利用变量LTime记录每次离开目标区域的时间,并将最后一次时间记录为过境结束时间,并从每次过境开始到过境结束记录幅宽两点的轨迹数据集。
5.如权利要求1所述的方法,其特征在于,步骤c:利用步骤b得到的每次过境轨迹数据集构成的几何图形与目标区域几何图形通过几何拓扑运算求交得到本次实际的目标分析区域覆盖范围,进而利用每一次求得的实际目标分析区域覆盖范围,又通过几何拓扑运算求并得到所有分析卫星及载荷在某个时间段内对目标分析区域总的覆盖范围,然后根据历次过境的结束时间、开始时间、卫星曝光时间、分辨率、幅宽长度、量化位数和波段数通过公式(7)求得每次过境的条带数据量和根据每次过境的覆盖面积、条带面积、求得的条带数据量再通过公式(8)求得每次过境的有效数据量 。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510321035.XA CN106248346A (zh) | 2015-06-12 | 2015-06-12 | 一种基于几何拓扑学的高精度遥感卫星对地覆盖分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510321035.XA CN106248346A (zh) | 2015-06-12 | 2015-06-12 | 一种基于几何拓扑学的高精度遥感卫星对地覆盖分析方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106248346A true CN106248346A (zh) | 2016-12-21 |
Family
ID=57626336
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510321035.XA Pending CN106248346A (zh) | 2015-06-12 | 2015-06-12 | 一种基于几何拓扑学的高精度遥感卫星对地覆盖分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106248346A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108254770A (zh) * | 2018-01-24 | 2018-07-06 | 合肥工业大学 | 成像卫星搜索陆地移动目标的方法及装置 |
CN108388687A (zh) * | 2018-01-05 | 2018-08-10 | 合肥工业大学 | 考虑传感器侧摆的多成像卫星区域覆盖优化方法 |
CN111336994A (zh) * | 2020-02-06 | 2020-06-26 | 北京理工大学 | 基于图形学与数值计算相结合的遥感卫星覆盖分析方法 |
CN111639423A (zh) * | 2020-05-21 | 2020-09-08 | 北京天工科仪空间技术有限公司 | 遥感卫星对地覆盖观测的方法及存储介质及设备 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020145606A1 (en) * | 2000-12-27 | 2002-10-10 | Isaac Levanon | Efficient image parcel texture rendering with T-junction crack elimination |
US20030229456A1 (en) * | 2002-03-07 | 2003-12-11 | The Gov. Of The U.S.A. As Represented By The Secretary Of The Dept. Of Health & Human Services | Methods for predicting properties of molecules |
CN103093098A (zh) * | 2013-01-21 | 2013-05-08 | 武汉大学 | 一种卫星光学传感器动态观测能力的定量评价方法 |
CN103870697A (zh) * | 2014-03-20 | 2014-06-18 | 中国资源卫星应用中心 | 基于超越不等式的卫星覆盖预报方法 |
CN106156472A (zh) * | 2015-04-20 | 2016-11-23 | 中国科学院遥感与数字地球研究所 | 一种遥感卫星对地覆盖能力分析方法 |
-
2015
- 2015-06-12 CN CN201510321035.XA patent/CN106248346A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020145606A1 (en) * | 2000-12-27 | 2002-10-10 | Isaac Levanon | Efficient image parcel texture rendering with T-junction crack elimination |
US20030229456A1 (en) * | 2002-03-07 | 2003-12-11 | The Gov. Of The U.S.A. As Represented By The Secretary Of The Dept. Of Health & Human Services | Methods for predicting properties of molecules |
CN103093098A (zh) * | 2013-01-21 | 2013-05-08 | 武汉大学 | 一种卫星光学传感器动态观测能力的定量评价方法 |
CN103870697A (zh) * | 2014-03-20 | 2014-06-18 | 中国资源卫星应用中心 | 基于超越不等式的卫星覆盖预报方法 |
CN106156472A (zh) * | 2015-04-20 | 2016-11-23 | 中国科学院遥感与数字地球研究所 | 一种遥感卫星对地覆盖能力分析方法 |
Non-Patent Citations (2)
Title |
---|
罗伏军 等: "遥感卫星过境仿真系统研究", 《计算机仿真》 * |
翁慧慧: "遥感卫星对地覆盖分析与仿真", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108388687A (zh) * | 2018-01-05 | 2018-08-10 | 合肥工业大学 | 考虑传感器侧摆的多成像卫星区域覆盖优化方法 |
CN108254770A (zh) * | 2018-01-24 | 2018-07-06 | 合肥工业大学 | 成像卫星搜索陆地移动目标的方法及装置 |
CN108254770B (zh) * | 2018-01-24 | 2020-04-21 | 合肥工业大学 | 成像卫星搜索陆地移动目标的方法及装置 |
CN111336994A (zh) * | 2020-02-06 | 2020-06-26 | 北京理工大学 | 基于图形学与数值计算相结合的遥感卫星覆盖分析方法 |
CN111639423A (zh) * | 2020-05-21 | 2020-09-08 | 北京天工科仪空间技术有限公司 | 遥感卫星对地覆盖观测的方法及存储介质及设备 |
CN111639423B (zh) * | 2020-05-21 | 2024-02-09 | 北京天工科仪空间技术有限公司 | 遥感卫星对地覆盖观测的方法及存储介质及设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Di Napoli et al. | Machine learning ensemble modelling as a tool to improve landslide susceptibility mapping reliability | |
Song et al. | Pedoclimatic zone-based three-dimensional soil organic carbon mapping in China | |
Shruthi et al. | Object-based gully feature extraction using high spatial resolution imagery | |
Haapanen et al. | Delineation of forest/nonforest land use classes using nearest neighbor methods | |
Lombardo et al. | Exploiting Maximum Entropy method and ASTER data for assessing debris flow and debris slide susceptibility for the Giampilieri catchment (north‐eastern Sicily, Italy) | |
CN104821013A (zh) | 基于大地坐标系数字高程模型的地表面积提取方法及系统 | |
Lacoste et al. | Evaluating large-extent spatial modeling approaches: A case study for soil depth for France | |
CN106248346A (zh) | 一种基于几何拓扑学的高精度遥感卫星对地覆盖分析方法 | |
Arsanjani et al. | A morphological approach to predicting urban expansion | |
Liu et al. | An improved line-of-sight method for visibility analysis in 3D complex landscapes | |
CN111336994A (zh) | 基于图形学与数值计算相结合的遥感卫星覆盖分析方法 | |
EP3469516B1 (en) | Method and system for improving the resolution of sensor data | |
Nussbaum et al. | Mapping of soil properties at high resolution in Switzerland using boosted geoadditive models | |
Obiahu et al. | Spatiotemporal analysis of land use land cover changes and built-up expansion projection in predominantly dystric nitosol of Ebonyi state, Southeastern, Nigeria | |
Hnatushenko et al. | Satellite monitoring of consequences of illegal extraction of amber in Ukraine | |
Hasan et al. | LULC changes to riverine flooding: A case study on the Jamuna River, Bangladesh using the multilayer perceptron model | |
Wang et al. | Role of reservoir construction in regional land use change in Pengxi River basin upstream of the Three Gorges Reservoir in China | |
Xiang et al. | Mapping potential wetlands by a new framework method using random forest algorithm and big Earth data: A case study in China's Yangtze River Basin | |
Löw et al. | Land cover change detection in the Aralkum with multi-source satellite datasets | |
Raj et al. | Towards evaluating gully erosion volume and erosion rates in the Chambal badlands, Central India | |
Hind et al. | Assessment of land use–land cover changes using GIS, remote sensing, and CA–Markov model: a case study of Algiers, Algeria | |
Cheng et al. | Multi-scale Feature Fusion and Transformer Network for urban green space segmentation from high-resolution remote sensing images | |
Gatto et al. | A simplified semi-quantitative procedure based on the SLIP model for landslide risk assessment: The case study of Gioiosa Marea (Sicily, Italy) | |
Liu et al. | CALC-2020: a new baseline land cover map at 10 m resolution for the circumpolar Arctic | |
Perzl et al. | Geodata requirements for mapping protective functions and effects of forests |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161221 |