CN116912446B - 基于自闭合正交多边形的遥感建筑物样本快速采集方法 - Google Patents
基于自闭合正交多边形的遥感建筑物样本快速采集方法 Download PDFInfo
- Publication number
- CN116912446B CN116912446B CN202311174785.XA CN202311174785A CN116912446B CN 116912446 B CN116912446 B CN 116912446B CN 202311174785 A CN202311174785 A CN 202311174785A CN 116912446 B CN116912446 B CN 116912446B
- Authority
- CN
- China
- Prior art keywords
- line
- orthogonal
- building
- point
- remote sensing
- 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 70
- 238000002372 labelling Methods 0.000 claims abstract description 33
- 238000012937 correction Methods 0.000 claims abstract description 20
- 230000008569 process Effects 0.000 claims abstract description 9
- 238000004364 calculation method Methods 0.000 claims description 6
- 230000003993 interaction Effects 0.000 claims description 6
- 230000008030 elimination Effects 0.000 claims description 4
- 238000003379 elimination reaction Methods 0.000 claims description 4
- 230000000007 visual effect Effects 0.000 claims description 3
- 230000003213 activating effect Effects 0.000 claims description 2
- 238000012887 quadratic function Methods 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 abstract description 4
- 238000011835 investigation Methods 0.000 abstract description 2
- 238000009877 rendering Methods 0.000 abstract description 2
- 230000000694 effects Effects 0.000 description 7
- 238000010276 construction Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000003321 amplification Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 238000003199 nucleic acid amplification method Methods 0.000 description 3
- 238000009826 distribution Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000007665 sagging Methods 0.000 description 2
- 208000028257 Joubert syndrome with oculorenal defect Diseases 0.000 description 1
- 238000013473 artificial intelligence Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013523 data management Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
- G01C11/04—Interpretation of pictures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/05—Geographic models
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Computer Graphics (AREA)
- Theoretical Computer Science (AREA)
- Processing Or Creating Images (AREA)
Abstract
本发明公开了一种基于自闭合正交多边形的遥感建筑物样本快速采集方法,包括:1、高分辨率遥感影像在线加载;2、激发正交多边形绘制组件;3、实时计算正交辅助线;4、实现建筑物拐点正交修正;5、实现结束时自动闭合;6、坐标串提交并存储。本发明可应用于多种分辨率的建筑物采集,将为构建大规模遥感建筑物标注样本提供快速采集,综合效率提升14%,采集建筑物样本的正交规则化质量大幅度提升,提升遥感建筑物语义识别泛化性和精确度,减少数据标注不足对语义识别的瓶颈问题。本发明可服务遥感建筑物语义识别,可广泛应用于城市扩张AI监测、城市调查AI监测、LOD1级实景三维建筑物面图形采集等。
Description
技术领域
本发明属于地理信息系统、遥感信息应用、人工智能应用方法的交叉领域,尤其涉及一种基于自闭合正交多边形的遥感建筑物样本快速采集方法。
背景技术
高分辨率遥感影像语义识别是当前AI应用的重要方向之一。遥感建筑物语义识别是高分辨率遥感影像语义识别中典型应用场景之一,可广泛应用于城市扩张监测、城市调查监测、LOD1级实景三维生成等。受限于遥感影像的“同物异谱、同谱异物”及“地理空间是巨型复杂系统”,遥感建筑物语义识别的泛化能力、精确度极其依赖于遥感建筑物标注样本的数据规模,因此,构建大规模的遥感建筑物标注样本是提升遥感建筑物语义识别质量的关键。
当前,构建大规模的遥感建筑物标注样本均依赖于交互式人工采集,其面临着采集效率不高、高分辨率遥感影像数据体量大、分发管理协同难度大,极大制约了构建大规模的遥感建筑物标注样本构建的效率,极大增加了构建大规模的遥感建筑物标注样本的生产成本。同时,受限于绝大多数建筑物的形状的正交规则性,传统基于GIS软件绘制多边形形式的标注图形难以满足遥感建筑物标注样本的形状要求,从而极大程度限制了采集的遥感建筑物标注样本规范性。
因此,研制一种基于在线地图服务模式的、可支持自动正交辅助、自动闭合的遥感建筑物标注样本采集方法,亟待突破。
发明内容
发明目的:本发明的目的在于提供一种基于自闭合正交多边形的遥感建筑物样本快速采集方法,支持自动正交辅助、自动闭合的遥感建筑物标注样本采集,大幅度提升人工采集效率,大幅度提升遥感建筑物标注样本的正交规范性。
技术方案:本发明的一种基于自闭合正交多边形的遥感建筑物样本快速采集方法,实现了高分辨率遥感影像在线加载与显示;构建了正交多边形绘制地图组件;支持实时计算正交辅助线生成;实现建筑物拐点正交修正;支持结束时自动正交闭合;支持坐标串提交并存储。具体步骤如下:
步骤1、高分辨率遥感影像在线加载:将影像通过OGC WMTS服务形式,在浏览器前端,通过Leaflet形式加载并显示;
步骤2、激发正交多边形绘制组件:构建正交多边形绘制地图组件,并激活地图组件;
步骤3、实时计算正交辅助线:延房屋屋顶长边,绘制初始建筑物边线,完成后,调取正交辅助线算符,实时生成已有线的正交辅助线;
步骤4、实现建筑物拐点正交修正:绘制建筑物拐点,经正交拐点修正算符,得到正交修正后建筑物拐点,并形成过程建筑物边线,重复所述建筑物拐点正交修正步骤,直至所有需采集的建筑物拐点完成采集;
步骤5、实现结束时自动闭合:在地图组件中调取自动闭合算符,自动闭合形成遥感建筑物正交多边形标注图形,得到建筑物样本标注;
步骤6、坐标串提交并存储:在浏览器端,将遥感建筑物正交多边形标注图形生成坐标串,请求WebService服务,提交并存储至数据库中,即完成建筑物样本的快速采集。
进一步的,步骤1具体包括如下步骤:
步骤1.1、引入Leaflet库:在HTML文件中引入Leaflet库的CSS和JavaScript文件;
步骤1.2、创建Leaflet地图:在JavaScript中,创建一个Leaflet地图并将其显示在指定的<div>元素中;
步骤1.3、添加地图图层:通过Leaflet的L.tileLayer方法添加WMTS图层,L.tileLayer方法的关键参数包括:URL模板及Options;其中,URL模板:必需参数,用于指定瓦片图层的URL,URL模板中使用占位符{z}、{x}和{y}来代表缩放级别、瓦片的列索引和行索引,Leaflet会根据当前地图的缩放级别和视图范围动态加载相应的瓦片;Options为一个包含各种选项的JavaScript对象,可设置瓦片图层的行为和样式。
进一步的,步骤2具体包括如下步骤:
步骤2.1、创建控件类:创建一个控件类,所述控件类继承自L.Control,控件类包含控件的DOM元素和其他逻辑;
步骤2.2、添加控件到地图:在JavaScript中,通过myControl.addTo(map)将正交多边形绘制控件添加到地图上;
步骤2.3、样式和交互:通过在onAdd方法中创建控件的DOM元素,并添加样式和交互,自定义控件类中的onAdd方法在控件被添加到地图上时被调用,而onRemove方法在控件从地图上移除时被调用,在onRemove方法中添加控件从地图上移除时的清理逻辑。
进一步的,步骤3(图2)具体包括如下步骤:
步骤3.1、AB直线方程表示:在平面坐标系内,两点A和B/>之间的直线以由点斜式方程表示,直线的点斜式方程表示为:/>,其中/>是直线斜率,/>;
步骤3.2、计算正交辅助线方程:要找到与直线AB垂直的直线CD,找到直线AB的斜率的相反数,并通过B来确定正交辅助线方程,正交辅助线垂直斜率为,通过点B/>构建与直线AB垂直的直线CD,使用点斜式方程来表示CD:/>;
步骤3.3、计算正交辅助线两端点坐标:已知正交辅助线CD方程,在实际绘制时,需要以线段表示,因此需找到其两个端点坐标;设定量端点距离中心点B为,计算C、D坐标,以D点为例,根据正交辅助线CD方程/>,将点D的坐标代入方程得/>;同时,已知点D与点B之间距离为/>,得距离公式:/> ,得到方程组:
;
通过消元法,变换为含单个未知变量的一元二次函数,以为例,消元后得到:
;
通过牛顿法,求解,得/>两个值,分别代表/>两端点坐标;
步骤3.4、绘制CD辅助线:将线两个端点绘制,然后连线,设置虚线线型和颜色,形成正交辅助线。
进一步的,步骤3.3中,所述牛顿法是一种数值逼近方法,用于求解方程的根,牛顿法是通过选择一个初始近似解,并使用该近似解和函数的导数来构造一个切线,然后通过切线与x轴的交点作为新的近似解,这个过程会不断迭代,直到找到满足精度要求的根(图3);具体包括:
(1)选择一个初始近似解;
(2)计算函数在/>处的值/>;
(3)计算函数在/>处的导数/>;
(4)使用切线方程,将切线与x轴求交点,得到新的近似解;检查新的近似解/>与旧的近似解/>之间的差异,如果满足精度要求,则停止迭代并将/>作为方程的近似解;否则,将/>作为新的近似解,回到牛顿法步骤(2)继续迭代;直到找到满足精度要求的根或达到最大迭代次数。
进一步的,步骤4具体包括如下步骤:
步骤4.1、绘制E’点(图4步骤4.1):样本标注人员采集建筑物拐点,通过鼠标左击,绘制建筑物拐点E’,得坐标E’,此时不满足正交要求,需进行正交拐点修正;
步骤4.2、正交拐点E修正计算(图4步骤4.1):过E’点,构建AB的平行线A’B’,A’B’方程可描述为:,求取A’B’与CD的交点E,实现E’正交改正至E,A’B’与CD两直线方程可表示为:
;
取交点得:;
其中,表示正交辅助线垂直斜率,/>表示直线斜率。
绘制E点(图4步骤4.3),绘制BE直线。
进一步的,步骤5(图5)具体包括如下步骤:
步骤5.1、自动闭合算符首先需明确起始线与终止线的平行和垂直关系:通过已采集线段的奇数偶数关系来确定起始线与终止线的平行和垂直关系,若已采集线段数为奇数,则起始线与终止线为平行关系;若已采集线段数为偶数,则起始线与终止线为垂直关系;
步骤5.2、自动闭合算符中,其次需计算出闭合点,闭合点每次都会存在两种情况,每种情况又包括两种之情况,分别讨论如下:
当起始线与终止线为垂直关系,即自闭合情况a:(1)、子情况a2和a4情况,AB延长线与UV延长线求交,简称“延延求交”,得到候选闭合点M;(2)、子情况a1和a3情况,AB垂直线与UV垂直线求交,简称“垂垂求交”,得到候选闭合点M;
当起始线与终止线为平行关系即,自闭合情况b:(1)、子情况b2和b4情况,AB垂直线与UV延长线求交,简称“垂延求交”,得到候选闭合点M;(2)、子情况b1和b3情况,AB延长线与UV垂直线求交,简称“延垂求交”,得到候选闭合点M;
步骤5.3、计算多边形面积,保留面积最大情况下的最优闭合点;起始线与终止线为垂直关系、起始线与终止线为平行关系都会分别产生两个交点,从两个交点中需要保留一个最优闭合点,保留的原则为:比较两种子情况的多边形面积大小关系,遵循“面积大的情况更符合视觉推理”基本原则;此处给出计算多边形的面积的方法,即简单多边形通过坐标串完成多边形面积计算,假设多边形顶点坐标依次为 ,基于向量叉乘的性质得多边形面积:
;
该公式绝对值旨在确保结果为正数。
进一步的,步骤6中,所述坐标串的格式为GeoJSON格式。
进一步的,所述数据库为POSTGIS数据库。
有益效果:与现有技术相比,本发明具有如下显著优点:
(1)支持在线地图服务模式下的多人协同标注采集,无需进行高分辨遥感影像分发,可有效提升数据管理效率;(2)支持自动正交辅助、自动闭合的遥感建筑物标注样本采集,大幅度提升人工采集效率,大幅度提升遥感建筑物标注样本的正交规范性。(3)本申请可应用于多种分辨率的建筑物采集,将为构建大规模遥感建筑物标注样本提供快速采集,综合效率提升14%,样本正交规则化综合质量大幅度提升,提升遥感建筑物语义识别泛化性和精确度,减少数据标注不足对语义识别的瓶颈问题。
附图说明
图1为本发明流程图;
图2为正交辅助线算符流程图;
图3为牛顿法求解端点伪代码流程图;
图4为正交拐点修正算符流程图;
图5为自动闭合算符图;
图6为研究区与实验数据图;
图7为GeoJSON格式的坐标串;
图8为遥感建筑物正交多边形一级放大标注结果图;
图9为遥感建筑物正交多边形二级放大标注结果图;
图10为遥感建筑物正交多边形三级放大标注结果图;
图11为一级放大标注效果对比图(含影像);
图12为一级放大标注效果对比图(无影像);
图13为二级放大标注效果对比图(无影像);
图14为三级放大标注效果对比图(无影像)。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
如图1所示,本发明的一种基于自闭合正交多边形的遥感建筑物样本快速采集方法,实验步骤如下:
实验前准备:本发明选取南京市栖霞区仙林大学城所在片区作为建筑物采集区域,具体采集范围为118.908°N-118.949°N,32.094°E-32.116°E,建筑物采集面积约为9.5km2,高分辨率遥感影像数据来源于天地图,分辨率为1m,时相为2021年,数据概况如附图6所示。实验系统硬件环境:CPU:Intel(R)Core(TM)i7-7590;内存:16GB;GPU:NVIDIAGeForce GT1080TI;系统软件环境:操作系统:操作系统:Windows10 64 Bit;编程语言Python、GDAL,数据库PostgreSQL、PostGIS插件,GIS服务平台为GeoServer,对比GIS软件选择QGIS桌面软件。
步骤1、高分辨率遥感影像在线加载与显示。将影像以OGC WMTS服务形式,通过Leaflet形式加载,在浏览器前端显示。
步骤1.1、引入Leaflet库。在HTML文件中引入Leaflet库的CSS和JavaScript文件。
步骤1.2、创建Leaflet地图。在JavaScript中,创建一个Leaflet地图并将其显示在指定的<div>元素中。
步骤1.3、添加地图图层。通过Leaflet的L.tileLayer方法添加WMTS图层。L.tileLayer方法的关键参数包括:URL模板及Options(选项参数)。其中,URL模板:必需参数,用于指定瓦片图层的URL,如:http://your_geoserver_url/gwc/service/wmts?REQUEST=GetTile&SERVICE=WMTS&VERSION=1.0.0&LAYER=my_workspace:my_layer&STYLE=&TILEMATRIXSET=EPSG:4326&TILEMATRIX=EPSG:4326:{z}&TILEROW={y}&TILECOL={x}&FORMAT=image/png。Options(选项参数):一个包含各种选项的JavaScript对象,可以设置瓦片图层的行为和样式,如:minZoom: 0, maxZoom: 18, tileSize: 256。
步骤2、激发正交多边形绘制组件,构建正交多边形绘制地图组件,点击激活地图组件。
步骤2.1、创建控件类。首先,创建一个控件类,这个类会继承自L.Control。这个类将包含控件的DOM元素和其他逻辑。
步骤2.2、添加控件到地图。在JavaScript中,将正交多边形绘制控件添加到地图上,可通过myControl.addTo(map)。
步骤2.3、样式和交互。通过在onAdd方法中创建控件的DOM元素,并添加样式和交互。自定义控件类中的onAdd方法在控件被添加到地图上时会被调用,而onRemove方法在控件从地图上移除时会被调用。在onRemove方法中添加控件从地图上移除时的清理逻辑。
步骤3、实时计算正交辅助线,延房屋屋顶长边,左击鼠标,绘制初始建筑物边线,完成后,调取正交辅助线算符,实时生成已有线的正交辅助线(简称:正交辅助线),具体包括4个子步骤,详见附图2,步骤如下:
步骤3.1、AB直线方程表示。在平面坐标系内,两点A和B/>之间的直线可以由点斜式方程表示。直线的点斜式方程可以表示为:/>,其中/>是直线斜率,/>。
步骤3.2、计算正交辅助线方程。要找到与直线AB垂直的直线CD,我们需要找到直线AB的斜率的相反数,并通过B来确定正交辅助线方程,正交辅助线垂直斜率为。通过点B/>构建与直线AB垂直的直线CD,使用点斜式方程来表示CD:/>。
步骤3.3、计算正交辅助线两端点坐标。已知正交辅助线CD方程,在实际绘制时,需要以线段表示,因此需找到其两个端点坐标。设定量端点距离中心点B为,计算C、D坐标,接下来以D点为例。根据正交辅助线CD方程
,将点D的坐标代入方程得。
为简化理论计算过程、提高绘制效率, Leaflet中直接提供了L.GeometryUtil.destination()方法,以快速计算终点坐标D:var endPoint =L.GeometryUtil.destination(startPoint, slope, distance),其中startPoint为起点B坐标,slope为斜率,distance为距离中心点B为/>。同理,可以获得C点坐标。
步骤3.4、绘制CD辅助线。将线两个端点绘制,然后连线,设置虚线线型和颜色,形成正交辅助线。创建BC和BD至地图中:
var latlngs = [startPoint, endPoint];
var polyline = L.polyline(latlngs, { color: 'red' }).addTo(map);
步骤4、实现建筑物拐点正交修正,绘制建筑物拐点,经正交拐点修正算符,如图4所示,得到正交修正后建筑物拐点(简称:正交拐点),并形成过程建筑物边线,重复所述建筑物拐点正交修正步骤,直至所有需采集的建筑物拐点完成采集,具体2个子步骤如下:
步骤4.1、绘制E’点(附图4步骤4.1)。样本标注人员采集建筑物拐点,通过鼠标左击,绘制 建筑物拐点E’,得坐标E’,此时不满足正交要求,需进行正交拐点修正。
步骤4.2、正交拐点E修正计算(附图4步骤4.2)。过E’点,构建AB的平行线A’B’,A’B’方程可描述为:。求取A’B’与CD的交点E,实现E’正交改正至E,A’B’与CD两直线方程可表示为:
取交点得:
绘制E点(附图4步骤4.3),绘制BE直线,具体绘制方式与步骤3.4一致。
重复步骤4.1和4.2,直至待采集建筑物所有拐点完成了采集。
步骤5、实现结束时自动闭合,右击鼠标,调取自动闭合算符,自动闭合形成遥感建筑物正交多边形标注图形,得到建筑物样本标注,具体步骤如下:
步骤5.1、自动闭合算符首先需明确起始线与终止线的平行和垂直关系。通过已采集线段的奇数偶数关系来确定起始线与终止线的平行和垂直关系,若已采集线段数为奇数,则起始线与终止线为平行关系;若已采集线段数为偶数,则起始线与终止线为垂直关系。
步骤5.2、自动闭合算符中,其次需计算出闭合点。闭合点每次都会存在两种情况,每种情况中有两种子情况,如附图5所示,分别讨论如下:
当起始线与终止线为垂直关系,即自闭合情况a:(1)、子情况a2和a4情况,AB延长线与UV延长线求交,简称“延延求交”,得到候选闭合点M;(2)、子情况a1和a3情况,AB垂直线与UV垂直线求交,简称“垂垂求交”,得到候选闭合点M;
当起始线与终止线为平行关系即,自闭合情况b:(1)、子情况b2和b4情况,AB垂直线与UV延长线求交,简称“垂延求交”,得到候选闭合点M;(2)、子情况b1和b3情况,AB延长线与UV垂直线求交,简称“延垂求交”,得到候选闭合点M;
M点求解时,垂线方程构建过程与步骤3.1和3.2一致,交点构建过程与步骤4.1和4.2一致。
步骤5.3、计算多边形面积,保留面积最大情况下的最优闭合点。以上讨论情况中,起始线与终止线为垂直关系、起始线与终止线为平行关系都会分别产生两个交点,从两个交点中需要保留一个最优闭合点,保留的原则为:比较两种子情况的多边形面积大小关系,遵循“面积大的情况更符合视觉推理”基本原则。此处给出计算多边形的面积的方法,即简单多边形(顶点按顺时针或逆时针顺序给出),通过坐标串完成多边形面积计算。假设多边形顶点坐标依次为
基于向量叉乘的性质得多边形面积:
该公式绝对值保证了结果为正数。
步骤6、坐标串提交并存储,在浏览器端,将遥感建筑物正交多边形标注图形生成为GeoJSON格式的坐标串,数据组织格式如图7所示,请求WebService服务,提交并存储至POSTGIS数据库中。
步骤6.1、安装所需的依赖,包括Python 3.8和PostgreSQL 3.2,并安装了以下Python库:Flask、Flask-CORS和GeoAlchemy2。
步骤6.2、初始化Flask应用,创建一个名为app.py的Python文件,并初始化Flask应用。
步骤6.3、创建数据库模型,在app.py文件中创建一个数据库模型,用于定义GeoJSON数据在数据库中的表结构。
步骤6.4、定义路由和处理方法,在app.py文件中定义路由和处理方法,用于接收GeoJSON数据并将其存储到数据库中,最终形成的样本采集数据结果如附图8(一级放大)、附图9(二级放大)、附图10(三级放大)所示。
此外,通过对比实验,分别对同一区域开展采集效率和效果对比(详见详见附图11-附图14),采集建筑物数量为72个(附图11)。方式1、基于本发明的自闭正交多边形方式采集,方式2、基于QGIS软件绘制多边形采集,得到对比效果(附图12)。由附图13及附图14可见,本发明采集样本效果边界规则、正交,有效约束了因人为采集落点不准而导致的多边形变形情况,可有效适用于遥感建筑物样本采集。同时,两区域采集时间分别为12分钟、14分钟,由于本发明的自闭正交多边形方式采集方式可减少绘制结点数量,实现推测闭合,其采集效率也要高于基于QGIS软件绘制多边形采集。综述,本发明方法可有效提升遥感建筑物语义识别泛化性和精确度,减少数据标注不足对语义识别的瓶颈问题。
Claims (6)
1.一种基于自闭合正交多边形的遥感建筑物样本快速采集方法,其特征在于,包括以下步骤:
步骤1、高分辨率遥感影像在线加载:将影像通过OGC WMTS服务形式,在浏览器前端,通过Leaflet形式加载并显示;
步骤2、激发正交多边形绘制组件:构建正交多边形绘制地图组件,并激活地图组件;
步骤3、实时计算正交辅助线:延房屋屋顶长边,绘制初始建筑物边线,完成后,调取正交辅助线算符,实时生成已有线的正交辅助线;
步骤4、实现建筑物拐点正交修正:绘制建筑物拐点,经正交拐点修正算符,得到正交修正后建筑物拐点,并形成过程建筑物边线,重复所述建筑物拐点正交修正步骤,直至所有需采集的建筑物拐点完成采集;
步骤5、实现结束时自动闭合:在地图组件中调取自动闭合算符,自动闭合形成遥感建筑物正交多边形标注图形,得到建筑物样本标注;
步骤6、坐标串提交并存储:在浏览器端,将遥感建筑物正交多边形标注图形生成坐标串,请求WebService服务,提交并存储至数据库中,即完成建筑物样本的快速采集;
步骤3具体包括如下步骤:
步骤3.1、AB直线方程表示:在平面坐标系内,两点A(xA,yA)和B(xB,yB)之间的直线以由点斜式方程表示,直线的点斜式方程表示为:y-yA=m(x-xA),其中m是直线斜率,m=(yB-yA)/(xB-xA);
步骤3.2、计算正交辅助线方程:要找到与直线AB垂直的直线CD,找到直线AB的斜率的相反数,并通过B来确定正交辅助线方程,正交辅助线垂直斜率为通过点B(xB,yB)构建与直线AB垂直的直线CD,使用点斜式方程来表示CD:y-yB=mperpendicular(x-xB);
步骤3.3、计算正交辅助线两端点坐标:已知正交辅助线CD方程,在实际绘制时,需要以线段表示,因此需找到其两个端点坐标;设定量端点距离中心点B为d,计算C、D坐标,以D点为例,根据正交辅助线CD方程y-yB=mperpendicular(x-xB),将点D的坐标代入方程得yD-yB=mperpendicular(xD-xB);同时,已知点D与点B之间距离为d,得距离公式:d2=(xB–xD)2+(yB–yD)2,得到方程组:
通过消元法,变换为含单个未知变量的一元二次函数,以xD为例,消元后得到:
通过牛顿法,求解f(xD)=0,得xD两个值,分别代表xD两端点坐标;
步骤3.4、绘制CD辅助线:将线两个端点绘制,然后连线,设置虚线线型和颜色,形成正交辅助线;
步骤4具体包括如下步骤:
步骤4.1、绘制E′点:样本标注人员采集建筑物拐点,通过鼠标左击,绘制建筑物拐点E′,得坐标E′(xE′,yE′),此时不满足正交要求,需进行正交拐点修正;
步骤4.2、正交拐点E修正计算:过E′点,构建AB的平行线A′B′,A′B′方程可描述为:y-yE′=m(x-xE′),求取A′B′与CD的交点E,实现E′正交改正至E,A′B′与CD两直线方程可表示为:
取交点得:(xE,yE)=((yE′-mxE′-yB+mperpendicularxB)/(mperpendicular-m),mperpendicularxE+yB-mperpendicularxB);
其中,mperpendicular表示正交辅助线垂直斜率,m表示直线斜率;
步骤5具体包括如下步骤:
步骤5.1、自动闭合算符首先需明确起始线与终止线的平行和垂直关系:通过已采集线段的奇数偶数关系来确定起始线与终止线的平行和垂直关系,若已采集线段数为奇数,则起始线与终止线为平行关系;若已采集线段数为偶数,则起始线与终止线为垂直关系;
步骤5.2、自动闭合算符中,其次需计算出闭合点,闭合点每次都会存在两种情况,每种情况又包括两种之情况,分别讨论如下:
当起始线与终止线为垂直关系,即自闭合情况a:(1)、子情况a2和a4情况,AB延长线与UV延长线求交,简称延延求交,得到候选闭合点M;(2)、子情况a1和a3情况,AB垂直线与UV垂直线求交,简称垂垂求交,得到候选闭合点M;
当起始线与终止线为平行关系即,自闭合情况b:(1)、子情况b2和b4情况,AB垂直线与UV延长线求交,垂延求交,得到候选闭合点M;(2)、子情况b1和b3情况,AB延长线与UV垂直线求交,延垂求交,得到候选闭合点M;
步骤5.3、计算多边形面积,保留面积最大情况下的最优闭合点;起始线与终止线为垂直关系、起始线与终止线为平行关系都会分别产生两个交点,从两个交点中需要保留一个最优闭合点,保留的原则为:比较两种子情况的多边形面积大小关系,遵循“面积大的情况更符合视觉推理”基本原则;此处给出计算多边形的面积的方法,即简单多边形通过坐标串完成多边形面积计算,假设多边形顶点坐标依次为(x1,y1),(x2,y2),...,(xn,yn),基于向量叉乘的性质得多边形面积:
Area=0.5*|(x1y2+x2y3+...+xny1)-(y1x2+y2x3+...+ynx1)|
2.根据权利要求1所述的一种基于自闭合正交多边形的遥感建筑物样本快速采集方法,其特征在于,步骤1具体包括如下步骤:
步骤1.1、引入Leaflet库:在HTML文件中引入Leaflet库的CSS和JavaScript文件;
步骤1.2、创建Leaflet地图:在JavaScript中,创建一个Leaflet地图并将其显示在指定的<div>元素中;
步骤1.3、添加地图图层:通过Leaflet的L.tileLayer方法添加WMTS图层,L.tileLayer方法的关键参数包括:URL模板及Options;其中,URL模板:必需参数,用于指定瓦片图层的URL,URL模板中使用占位符{z}、{x}和{y}来代表缩放级别、瓦片的列索引和行索引,Leaflet会根据当前地图的缩放级别和视图范围动态加载相应的瓦片;Options为一个包含各种选项的JavaScript对象,可设置瓦片图层的行为和样式。
3.根据权利要求1所述的一种基于自闭合正交多边形的遥感建筑物样本快速采集方法,其特征在于,步骤2具体包括如下步骤:
步骤2.1、创建控件类:创建一个控件类,所述控件类继承自L.Control,控件类包含控件的DOM元素和其他逻辑;
步骤2.2、添加控件到地图:在JavaScript中,通过myControl.addTo(map)将正交多边形绘制控件添加到地图上;
步骤2.3、样式和交互:通过在onAdd方法中创建控件的DOM元素,并添加样式和交互,自定义控件类中的onAdd方法在控件被添加到地图上时被调用,而onRemove方法在控件从地图上移除时被调用,在onRemove方法中添加控件从地图上移除时的清理逻辑。
4.根据权利要求1所述的一种基于自闭合正交多边形的遥感建筑物样本快速采集方法,其特征在于,步骤3.3中,所述牛顿法是一种数值逼近方法,用于求解方程的根,牛顿法是通过选择一个初始近似解,并使用该近似解和函数的导数来构造一个切线,然后通过切线与x轴的交点作为新的近似解,这个过程会不断迭代,直到找到满足精度要求的根;具体包括:
(1)选择一个初始近似解x0;
(2)计算函数f(xD)在x0处的值f(x0);
(3)计算函数f(xD)在x0处的导数f′(x0);
(4)使用切线方程y=f′(x0)(x-x0)+f(x0),将切线与x轴求交点,得到新的近似解x1;检查新的近似解x1与旧的近似解x0之间的差异,如果满足精度要求,则停止迭代并将x1作为方程的近似解;否则,将x1作为新的近似解,回到牛顿法步骤(2)继续迭代,直到找到满足精度要求的根或达到最大迭代次数。
5.根据权利要求1所述的一种基于自闭合正交多边形的遥感建筑物样本快速采集方法,其特征在于,步骤6中,所述坐标串的格式为GeoJSON格式。
6.根据权利要求1所述的一种基于自闭合正交多边形的遥感建筑物样本快速采集方法,其特征在于,所述数据库为POSTGIS数据库。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311174785.XA CN116912446B (zh) | 2023-09-13 | 2023-09-13 | 基于自闭合正交多边形的遥感建筑物样本快速采集方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311174785.XA CN116912446B (zh) | 2023-09-13 | 2023-09-13 | 基于自闭合正交多边形的遥感建筑物样本快速采集方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116912446A CN116912446A (zh) | 2023-10-20 |
CN116912446B true CN116912446B (zh) | 2023-12-08 |
Family
ID=88360687
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311174785.XA Active CN116912446B (zh) | 2023-09-13 | 2023-09-13 | 基于自闭合正交多边形的遥感建筑物样本快速采集方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116912446B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114898118A (zh) * | 2022-06-02 | 2022-08-12 | 中国能源建设集团江苏省电力设计院有限公司 | 基于多源点云的输电线路房屋拆迁量自动统计方法及系统 |
CN115048014A (zh) * | 2022-05-24 | 2022-09-13 | 长光卫星技术股份有限公司 | 一种基于卫星遥感影像的WebGIS在线目视解译系统、方法、存储介质和电子设备 |
CN116127288A (zh) * | 2023-04-14 | 2023-05-16 | 南京邮电大学 | 基于独立成分分析的纳米孔传感信号噪声去除方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6424933B1 (en) * | 2000-03-17 | 2002-07-23 | Vicinity Corporation | System and method for non-uniform scaled mapping |
US9310674B2 (en) * | 2014-02-20 | 2016-04-12 | International Business Machines Corporation | Mask that provides improved focus control using orthogonal edges |
-
2023
- 2023-09-13 CN CN202311174785.XA patent/CN116912446B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115048014A (zh) * | 2022-05-24 | 2022-09-13 | 长光卫星技术股份有限公司 | 一种基于卫星遥感影像的WebGIS在线目视解译系统、方法、存储介质和电子设备 |
CN114898118A (zh) * | 2022-06-02 | 2022-08-12 | 中国能源建设集团江苏省电力设计院有限公司 | 基于多源点云的输电线路房屋拆迁量自动统计方法及系统 |
CN116127288A (zh) * | 2023-04-14 | 2023-05-16 | 南京邮电大学 | 基于独立成分分析的纳米孔传感信号噪声去除方法及装置 |
Non-Patent Citations (3)
Title |
---|
Linear building pattern recognition in topographical maps combining convex polygon decomposition;Wei, ZW 等;《GEOCARTO INTERNATIONAL》;第37卷(第26期);11365-11389 * |
POLYGON PARTITION AND SHAPE RECTIFICATION FOR AUTOMATIC GENERATION OF 3D BUILDING MODEL;Sugihara, K 等;《IEEE Region 10 Conference (TENCON)》;2311-2314 * |
李亚文 ; .基于规则集的高分辨率遥感影像城市建筑物自动提取.《发展与创新》.2020,(第14期),第263-264页. * |
Also Published As
Publication number | Publication date |
---|---|
CN116912446A (zh) | 2023-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109993827B (zh) | 一种将建筑图纸转换为三维bim模型的立面图识别方法 | |
Croce et al. | From the semantic point cloud to heritage-building information modeling: A semiautomatic approach exploiting machine learning | |
Lopez et al. | A framework for using point cloud data of heritage buildings toward geometry modeling in a BIM context: a case study on Santa Maria La Real De Mave Church | |
Lin et al. | Structure-preserving retargeting of irregular 3d architecture | |
CN113094770B (zh) | 图纸生成方法、装置、计算机设备和存储介质 | |
CN102663800A (zh) | 一种顾及城市意象的城市建筑综合与渲染的方法 | |
Liao et al. | Cluster-based visual abstraction for multivariate scatterplots | |
CN114202223A (zh) | 企业信用风险评分方法、装置、设备及存储介质 | |
CN111968227B (zh) | 三维地质断层网络不确定性分析方法、系统及存储介质 | |
Gao et al. | A practical grid generation procedure for the design of free-form structures | |
CN114239508A (zh) | 表格还原方法和装置、存储介质及电子设备 | |
Li et al. | HiTailor: Interactive transformation and visualization for hierarchical tabular data | |
Demir et al. | Guided proceduralization: Optimizing geometry processing and grammar extraction for architectural models | |
Kao et al. | Value of information in portfolio selection, with a Taiwan stock market application illustration | |
CN116912446B (zh) | 基于自闭合正交多边形的遥感建筑物样本快速采集方法 | |
Wu et al. | Extracting POP: Pairwise orthogonal planes from point cloud using RANSAC | |
CN112132928B (zh) | 矢量图的处理方法、装置、电子设备和计算机可读介质 | |
CN116721218B (zh) | 一种三维不动产模型轻量化方法、系统及设备 | |
Podgorelec | A new constructive approach to constraint-based geometric design | |
Richards-Rissetto et al. | Procedural modeling for ancient Maya cityscapes initial methodological challenges and solutions | |
CN111063003B (zh) | 矿山分布图的制作方法和系统 | |
Bai et al. | VCMM: A visual tool for continuum molecular modeling | |
Eshonkulov | Methodology Of Application Of The Simplex Method In The Optimal Development Of Industrial Enterprises | |
Park et al. | Development of an AI advisor for conceptual land use planning | |
Zhang et al. | Prospect of architectonic grammar reconstruction from dense 3D point clouds: Historical building information modeling (HBIM) of Guangdong cultural heritages |
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 |