CN109712235B - 一种数字化宗地图的构建方法 - Google Patents
一种数字化宗地图的构建方法 Download PDFInfo
- Publication number
- CN109712235B CN109712235B CN201811432643.8A CN201811432643A CN109712235B CN 109712235 B CN109712235 B CN 109712235B CN 201811432643 A CN201811432643 A CN 201811432643A CN 109712235 B CN109712235 B CN 109712235B
- Authority
- CN
- China
- Prior art keywords
- arc
- arc section
- area
- calculated
- section
- 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
- 238000010276 construction Methods 0.000 title claims description 7
- 238000000034 method Methods 0.000 claims abstract description 24
- PCTMTFRHKVHKIS-BMFZQQSSSA-N (1s,3r,4e,6e,8e,10e,12e,14e,16e,18s,19r,20r,21s,25r,27r,30r,31r,33s,35r,37s,38r)-3-[(2r,3s,4s,5s,6r)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-19,25,27,30,31,33,35,37-octahydroxy-18,20,21-trimethyl-23-oxo-22,39-dioxabicyclo[33.3.1]nonatriaconta-4,6,8,10 Chemical compound C1C=C2C[C@@H](OS(O)(=O)=O)CC[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H]([C@H](C)CCCC(C)C)[C@@]1(C)CC2.O[C@H]1[C@@H](N)[C@H](O)[C@@H](C)O[C@H]1O[C@H]1/C=C/C=C/C=C/C=C/C=C/C=C/C=C/[C@H](C)[C@@H](O)[C@@H](C)[C@H](C)OC(=O)C[C@H](O)C[C@H](O)CC[C@@H](O)[C@H](O)C[C@H](O)C[C@](O)(C[C@H](O)[C@H]2C(O)=O)O[C@H]2C1 PCTMTFRHKVHKIS-BMFZQQSSSA-N 0.000 claims description 10
- 230000009466 transformation Effects 0.000 claims description 3
- 230000001131 transforming effect Effects 0.000 claims 1
- 230000007547 defect Effects 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Landscapes
- Image Generation (AREA)
Abstract
本发明的目的在于提供一种包含弧段的宗地图构建方法，包括以下步骤：A、根据界址点构建多边形，获取多边形的面积，求出弧段对应弓形的面积；B、实地测量获得弧段数、弧段的弦长，并根据各个弧段的弦长比作为面积比例系数，获得各个弧段对应弓形的面积；C、设定各个弧段半径和弧段圆心角的取值范围，利用二分法，构建求解弧段半径和弧段圆心角的超越方程模型；D、求解超越方程模型，得到各个弧段圆心角的数值，进而计算得到各个弧段对应的半径；E、根据原始宗地图中的界址围成的多边形尺寸参数、各个待计算弧段对应的圆心角及半径参数，绘制宗地图。该方法克服传统人工模式操作复杂、精度低、效率低的缺陷，具有效率高、成本低的特点。
Description
技术领域
本发明涉及计算机图形学与地理信息科学领域，具体涉及一种数字化宗地图的构建方法。
背景技术
在城镇地籍测量中，宗地图通常包括由界址点构成的多边形部分以及弧段部分，弧段部分采用圆弧，其主要原因在于数学模型简单，所以在道路弯道、十字路口交汇、宗地边界曲线等工程设计中使用非常广泛；特别是为了满足市政设计需求，宗地边界的选取有时必须采用曲线，确保与道路设计红线的一致性。众所周知宗地面积具有严肃的法律效益，一旦确定面积就不允许随意更改；上世纪八九十年代前的宗地图多为纸质版，缺乏数字化总地图信息，再加上部分数字化宗地图损毁，仅有纸质图纸的情况下，如何构建对应曲线保证发证面积与图形面积的一致性就成为困扰地籍制图人员的难题。传统的模式多采用人工估算，重复模拟直到满足精度需求为止，效率低下且具有很大不确定性。
发明内容
本发明的目的在于提供一种数字化宗地图的构建方法，该方法克服传统人工模式操作复杂、精度低、效率低的缺陷，具有效率高、成本低的特点。
本发明的技术方案包括以下步骤：
A、根据界址点构建多边形，获取多边形的面积，将多边形的面积与原始宗地图中的宗地总面积进行对比获取差值，从而求出弧段对应弓形的面积；
B、实地测量获得宗地面积所在位置的弧段数，以及各个弧段的两端点之间的距离,从而获得弧段的弦长，并根据各个弧段的弦长比作为面积比例系数，获得各个弧段对应弓形的面积；
C、设定各个弧段半径和弧段圆心角的取值范围，利用二分法，基于各个弧段的弧段对应弓形的面积、以及各个弧段的弦长，构建求解弧段半径和弧段圆心角的超越方程模型；
D、求解超越方程模型，得到各个弧段圆心角的数值，进而计算得到各个弧段对应的半径；
E、根据原始宗地图中的界址围成的多边形尺寸参数、各个待计算弧段对应的圆心角及半径参数，绘制宗地图。
优选地，所述的步骤B中，各个弧段和弦所围成的面积的获取过程如下：
获取原始宗地图中的宗地总面积以及界址围成的多边形面积计算得到弧段对应的弓形总面积，若只包含一个弧段，则将该弧段对应的弓形总面积作为该待计算弧段的弓形面积；若包含多个待计算弧段，则将弧段弓形总面积根据比例系数分配至各个待计算弧段，作为各个待计算弧段的弓形面积，各个待计算弧段的比例系数与其弦长成正比；
各个待计算弧段的弓形面积计算如下：
si＝st*di/(d1+d2+…+dn)；
其中st为弧段弓形总面积，si为第i个待计算弧段的弓形面积，di为第i个弧段对应的弦长，i＝1，2…，n；n为待计算弧段的数量。
优选地，所述步骤C中超越方程模型的构建过程如下：
如图2，弧段对应扇形面积为：(r2*θ)/2；
辅助三角形Ⅰ为：(r2*sin(θ))/2；
弧段弓形面积为：
s＝r2(θ-sin(θ))/2 (1)；
其中r为弧段对应的半径，θ为弧段对应的圆心角；s为弧段弓形面积；
弧段弦长d为：
d＝2*r*sin(θ/2) (2)；
由式(2)进行三角函数变换得到：
结合式(1)、(3)、(4)得到超越方程：
优选地，所述的步骤C中，设定弧段半径和弧段圆心角的取值范围的方法为：
如图3，实测获得弧段上的中点E，基于弧段的两端点AB和中点E，构建三角形ABE和矩形ABCD分别作为辅助三角形Ⅱ和辅助矩形，基于三角形ABE和矩形ABCD；三角形ABE的面积小于弧段对应弓形的面积，弧段对应弓形的面积小于矩形 ABCD的面积；
三角形ABE的面积为为:(d*r/2-(r2*sinθ)/2)；
矩形ABCD的面积为：2*(d*r/2-(r2*sinθ)/2)；
则：(d*r/2-(r2*sinθ)/2)<s<2*(d*r/2-(r2*sinθ)/2) (6)；
其中：
sinθ＝2*sin(θ/2)*cos(θ/2)；
sin(θ/2)＝d/2r；
将式(6)进行展开，得到弧段对应半径以及圆心角的取值范围：
(d4/16+s2)/(d*s)<r<((d4/4+s2)/(2*d*s) (7)；
2*arcsin(d2*s/d4/8+2*s2)<θ<2*arcsin(d2*s/(d4/4+s2)) (8)。
优选地，超越方程(5)的求解过程如下：
根据式(8)得到θ的最优解所在的范围，在上述范围内采用二分法对超越方程 (5)进行求解，求解次数为：log2((r2-r1)/c)+1次，其中r1为公式(7)中的r的下限值，r2为公式(7)中的r的上限值，c为半径精度。
本发明的数字化宗地图的构建方法优选用于仅有纸质宗地图且弧段参数缺失的的情况下，本发明方案通过采用弧段弦长为标准进行弧段弓形面积的分配，解决各个弧段弓形面积难以合理分配的问题，为后续参数计算打下坚实基础；借助辅助三角形Ⅰ构建关于弧段圆心角的超越方程模型，解决了弧段关键参数计算的工具；借助辅助三角形Ⅱ以及辅助矩形，通过几何面积关系计算弧段对应扇形的面积范围，进而求得弧段对应半径以及圆心角，从而根据该范围结合二分法对弧段圆心角的超越方程模型进行求解，得到关键参数：弧段圆心角，进而实现包含弧段的宗地图的绘制。该发明原理严谨，科学高效，克服了人工方法不确定性和效率低下的缺陷问题，具有很强的社会经济价值。
附图说明
图1为本发明数字化宗地图的构建方法的流程示意图；
图2为本发明数字化宗地图的构建方法的半径、弦长、圆心角示意图；
图3为本发明数字化宗地图的构建方法的辅助三角形Ⅰ、辅助三角形Ⅱ、辅助矩形示意图；
图中各部分名称及序号如下：
1为弧段，2为弧段端点(A,B)，3为弦，4为弧段对应圆心。
具体实施方式
下面结合实施例详细阐述本发明。
实施例1
如图1-3所示，本实施例一种数字化宗地图的构建方法，包括以下步骤：
A、根据界址点构建多边形，获取多边形的面积，将多边形的面积与原始宗地图中的宗地总面积进行对比获取差值，从而求出弧段对应弓形的面积；
B、实地测量获得宗地面积所在位置的弧段数，以及各个弧段的两端点之间的距离,从而获得弧段的弦长，并根据各个弧段的弦长比作为面积比例系数，获得各个弧段对应弓形的面积；本实施例中，弧段弦长d＝20米；该弧段对应的弓形面积s＝36.234 平方米；
所述的步骤B中，各个弧段和弦所围成的面积的获取过程如下：
获取原始宗地图中的宗地总面积以及界址围成的多边形面积计算得到弧段对应的弓形总面积，若只包含一个弧段，则将该弧段对应的弓形总面积作为该待计算弧段的弓形面积；若包含多个待计算弧段，则将弧段弓形总面积根据比例系数分配至各个待计算弧段，作为各个待计算弧段的弓形面积，各个待计算弧段的比例系数与其弦长成正比；
各个待计算弧段的弓形面积计算如下：
si＝st*di/(d1+d2+…+dn)；
其中st为弧段弓形总面积，si为第i个待计算弧段的弓形面积，di为第i个弧段对应的弦长，i＝1，2…，n；n为待计算弧段的数量；
C、设定各个弧段半径和弧段圆心角的取值范围，利用二分法，基于各个弧段的弧段对应弓形的面积、以及各个弧段的弦长，构建求解弧段半径和弧段圆心角的超越方程模型；
所述步骤C中超越方程模型的构建过程如下：
弧段对应扇形面积为：(r2*θ)/2；
辅助三角形Ⅰ为：(r2*sin(θ))/2；
弧段弓形面积为：
s＝r2(θ-sin(θ))/2 (1)；
其中r为弧段对应的半径，θ弧段对应的圆心角；s为弧段弓形面积；
弧段弦长d为：
d＝2*r*sin(θ/2) (2)；
由式(2)进行三角函数变换得到：
结合式(1)、(3)、(4)得到超越方程：
所述的步骤C中，设定弧段半径和弧段圆心角的取值范围的方法为：
实测获得弧段上的中点E，基于弧段的两端点AB和中点E，构建三角形ABE和矩形ABCD分别作为辅助三角形Ⅱ和辅助矩形，基于三角形ABE和矩形ABCD；三角形ABE的面积小于弧段对应弓形的面积，弧段对应弓形的面积小于矩形ABCD的面积；
三角形abe的面积为为:(d*r/2-(r2*sinθ)/2)；
矩形abcd的面积为：2*(d*r/2-(r2*sinθ)/2)；
则：(d*r/2-(r2*sinθ)/2)<s<2*(d*r/2-(r2*sinθ)/2) (6)；
其中：
sinθ＝2*sin(θ/2)*cos(θ/2)；
sin(θ/2)＝d/2r；
将式(6)进行展开，得到弧段对应半径以及圆心角的取值范围：
(d4/16+s2)/(d*s)<r<((d4/4+s2)/(2*d*s) (7)；
本实施例中，求解得到15.610<r<28.505；
2*arcsin(d2*s/d4/8+2*s2)<θ<2*arcsin(d2*s/(d4/4+s2)) (8)；
本实施例中，求解得到θ的取值范围为：
θ>2*arcsin(d2*s/d4/8+2*s2)＝0.796537686；
θ<2*arcsin(d2*s/(d4/4+s2))＝1.545108663；
D、求解超越方程模型，得到各个弧段圆心角的数值，进而计算得到各个弧段对应的半径；
根据式(8)得到θ的最优解所在的范围，在上述范围内采用二分法对超越方程 (5)进行求解，求解次数为：log2((r2-r1)/c)+1次，其中r1为公式(7)中的r的下限值，r2为公式(7)中的r的上限值，c为半径精度，本实施例中，c＝0.001米；
求解得到r＝20米，θ＝1.04719753；
E、根据原始宗地图中的界址围成的多边形尺寸参数、各个待计算弧段对应的圆心角及半径参数，绘制宗地图。
Claims (3)
1.一种包含弧段的宗地图构建方法，其特征在于包括以下步骤：
A、根据界址点构建多边形，获取多边形的面积，将多边形的面积与原始宗地图中的宗地总面积进行对比获取差值，从而求出弧段对应弓形的面积；
B、实地测量获得宗地面积所在位置的弧段数，以及各个弧段的两端点之间的距离,从而获得弧段的弦长，并根据各个弧段的弦长比作为面积比例系数，获得各个弧段对应弓形的面积；
C、设定各个弧段半径和弧段圆心角的取值范围，利用二分法，基于各个弧段的弧段对应弓形的面积、以及各个弧段的弦长，构建求解弧段半径和弧段圆心角的超越方程模型；
所述步骤C中超越方程模型的构建过程如下：
弧段对应扇形面积为：(r2*θ)/2；
辅助三角形Ⅰ为：(r2*sin(θ))/2；
弧段弓形面积为：
s＝r2(θ-sin(θ))/2 (1)；
其中r为弧段对应的半径，θ为弧段对应的圆心角；s为弧段弓形面积；
弧段弦长d为：
d＝2*r*sin(θ/2) (2)；
由式(1)进行变换得到：
由式(2)进行三角函数变换得到：
结合式(1)、(3)、(4)得到超越方程：
所述的步骤C中，设定弧段半径和弧段圆心角的取值范围的方法为：
实测获得弧段上的中点E，基于弧段的两端点AB和中点E，构建三角形ABE和矩形ABCD分别作为辅助三角形Ⅱ和辅助矩形，基于三角形ABE和矩形ABCD；三角形ABE的面积小于弧段对应弓形的面积，弧段对应弓形的面积小于矩形ABCD的面积；
三角形ABE的面积为为:(d*r/2-(r2*sinθ)/2)；
矩形ABCD的面积为：2*(d*r/2-(r2*sinθ)/2)；
则：(d*r/2-(r2*sinθ)/2)<s<2*(d*r/2-(r2*sinθ)/2) (6)；
其中：
sinθ＝2*sin(θ/2)*cos(θ/2)；
sin(θ/2)＝d/2r；
将式(6)进行展开，得到弧段对应半径以及圆心角的取值范围：
(d4/16+s2)/(d*s)<r<((d4/4+s2)/(2*d*s) (7)；
2*arcsin(d2*s/d4/8+2*s2)<θ<2*arcsin(d2*s/(d4/4+s2)) (8)；
D、求解超越方程模型，得到各个弧段圆心角的数值，进而计算得到各个弧段对应的半径；
E、根据原始宗地图中的界址围成的多边形尺寸参数、各个待计算弧段对应的圆心角及半径参数，绘制宗地图。
2.如权利要求1所述的包含弧段的宗地图构建方法，其特征在于：
所述的步骤B中，各个弧段和弦所围成的面积的获取过程如下：
获取原始宗地图中的宗地总面积以及界址围成的多边形面积计算得到弧段对应的弓形总面积，若只包含一个弧段，则将该弧段对应的弓形总面积作为该待计算弧段的弓形面积；若包含多个待计算弧段，则将弧段弓形总面积根据比例系数分配至各个待计算弧段，作为各个待计算弧段的弓形面积，各个待计算弧段的比例系数与其弦长成正比；
各个待计算弧段的弓形面积计算如下：
si＝st*di/(d1+d2+…+dn)；
其中st为弧段弓形总面积，si为第i个待计算弧段的弓形面积，di为第i个待计算弧段对应的弦长，i＝1，2…，n；n为待计算弧段的数量。
3.根据权利要求1所述的包含弧段的宗地图构建方法，其特征在于：
超越方程(5)的求解过程如下：
根据式(8)得到θ的最优解所在的范围，在上述范围内采用二分法对超越方程(5)进行求解，求解次数为：log2((r2-r1)/c)+1次，其中r1为公式(7)中的r的下限值，r2为公式(7)中的r的上限值，c为半径精度。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811432643.8A CN109712235B (zh) | 2018-11-28 | 2018-11-28 | 一种数字化宗地图的构建方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811432643.8A CN109712235B (zh) | 2018-11-28 | 2018-11-28 | 一种数字化宗地图的构建方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109712235A CN109712235A (zh) | 2019-05-03 |
| CN109712235B true CN109712235B (zh) | 2020-02-11 |
Family
ID=66254487
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811432643.8A Active CN109712235B (zh) | 2018-11-28 | 2018-11-28 | 一种数字化宗地图的构建方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109712235B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113781644B (zh) * | 2021-08-04 | 2024-03-12 | 城乡院(广州)有限公司 | 一种宗地草图的自适应调整方法、装置、设备及介质 |
Family Cites Families (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100384992B1 (ko) * | 2001-05-22 | 2003-05-23 | 삼아항업(주) | 자동 지도 제작 방법 |
| CN201045575Y (zh) * | 2007-04-26 | 2008-04-09 | 同济大学 | 土地面积测量的数据处理系统 |
| CN101650890A (zh) * | 2009-09-18 | 2010-02-17 | 西安煤航信息产业有限公司 | 一种道路在地图上的处理表示方法 |
| CN103971586B (zh) * | 2013-01-31 | 2016-08-31 | 北京四维图新科技股份有限公司 | 电子地图生成方法及装置 |
| CN103440680B (zh) * | 2013-08-22 | 2015-12-02 | 浙江大学 | 一种基于一范数优化的Polycube可控生成方法 |
| CN105677843B (zh) * | 2016-01-06 | 2018-08-21 | 中科宇图科技股份有限公司 | 一种自动获取宗地四至属性的方法 |
| CN105893520B (zh) * | 2016-03-30 | 2020-03-10 | 广东中冶地理信息股份有限公司 | 一种城市道路面积估算方法 |
| CN108805885A (zh) * | 2018-06-13 | 2018-11-13 | 广州极飞科技有限公司 | 地块分割方法及终端、航线规划方法及移动装置控制方法 |
-
2018
- 2018-11-28 CN CN201811432643.8A patent/CN109712235B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN109712235A (zh) | 2019-05-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Richardus | Project Surveying: Completely revised 2nd edition-General adjustment and optimization techniques with applications to engineering surveying | |
| KR20050006173A (ko) | 디지털 지형도의 작성방법 및 작성장치 | |
| CN104537606B (zh) | 地理坐标投影变化方法 | |
| CN103389080A (zh) | 一种基于地理信息应用城市独立坐标系参数获取方法 | |
| CN104123695A (zh) | 一种实现坐标转换方法 | |
| CN109712235B (zh) | 一种数字化宗地图的构建方法 | |
| CN108981660A (zh) | 一种三角高程的对边测量方法 | |
| CN101976462A (zh) | 一种三维重建方法 | |
| CN111651539B (zh) | 一种利用近景遥感技术实现平面地图元素快速更新的方法 | |
| Lacquement | Recalculating mound volume at moundville | |
| CN201140629Y (zh) | 地形学现地定点交会尺 | |
| CN101726283A (zh) | 航空摄影测量大比例尺测图标志方法 | |
| CN102121825A (zh) | 一种利用gps测量跨河水准高差的方法 | |
| Baiocchi et al. | Georeferencing the historical maps of Rome between the seventeenth and eighteenth centuries | |
| CN111611412B (zh) | 一种遥感影像地物建档数据的生成方法 | |
| CN102426010A (zh) | 一种变分辨率数字高程模型的生成方法 | |
| CN201503269U (zh) | 航空摄影测量大比例尺测图标志点装置 | |
| CN204854717U (zh) | 水准/InSAR一体化监测装置 | |
| JP7373626B1 (ja) | 変換装置、変換方法及び制御プログラム | |
| JP2004028630A (ja) | 路面形状の計測方法 | |
| de Lange | Geoobjects and Reference Systems | |
| CN209910534U (zh) | 一种土工用等分圆基础标尺 | |
| CN107727074B (zh) | 一种测距相机单片摄影测地块面积的方法 | |
| CN106126682B (zh) | 产业地图区域地斑智能绘制方法及系统 | |
| Jandžiković | Generalization and chartometric analysis of the Sitnica river network system |
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 | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20240821 Address after: 528000, Floor 22-23, Building B, Vanke Financial Center, No. 55 Jihua Fifth Road, Zumiao Street, Chancheng District, Foshan City, Guangdong Province (Residence Declaration) Patentee after: Foshan Surveying and Mapping Geographic Information Research Institute Co.,Ltd. Country or region after: China Address before: No. 1 Huijing 1st Street, Chancheng District, Foshan City, Guangdong Province, 528000 Patentee before: FOSHAN SERVEYING MAPPING AND GEOINFORMATION Research Institute Country or region before: China |
|
| TR01 | Transfer of patent right |