面积边长固定的四边形对角线长算法及应用.pdf

中国新技术新产品2024 NO.4（下）-102-工 程 技 术许多老旧的宗地图受当时测绘技术条件的限制，只有边长和面积，没有坐标，也没有可参照的地形地物。在新的不动产登记调查作业中，要恢复其形状和界址点坐标的同时，还要保证面积边长不变，手工调整难度很大。任何形状的宗地在大致确定各界址点位置后，均可将该宗地证载面积减去其他不可调整的固定部分面积后的剩余面积，归算到一个四边形上。在保证四边形四个边长一定，面积等于宗地剩余面积的约束条件下，通过迭代计算四边形对角线长，从而将不稳定的四边形变为两个共边的稳定三角形，这样就稳定了该四边形的形状。根据四边形与宗地其他部分的两个衔接点坐标，按照距离交会法可以计算其他两个顶点的坐标，从而恢复老宗地图的界址点。1 四边形对角线长与面积关系四边形是一个不稳定的图形，边长一定的四边形其面积是无法固定的。固定四边形的一条对角线后，四边形就成了两个共边的稳定三角形，四边形也变得稳定且具有唯一的固定的面积。图 1 为典型四边形，根据海伦公式可以求出三角形 123和三角形 134 的面积，这两个三角形的面积之和便是该四边形的面积，于是便得出四边形的面积 S 与对角线长度 d 的函数关系，如公式（1）和公式（2）所示。s ddddddddddddddddd?141412121221343d dddddddd43443?（1）s fddfddfdfddfdddfd?141423232332414d dfdfddfd14114?（2）2 求解角线长度当四边形的面积一定时，公式（1）或公式（2）可求出对角线长度 d 或 f。由于公式（1）和公式（2）不是线性方程，因此无法直接求解，须采用迭代法进行计算。2.1 对角线长度的取值范围及面积极值在图 1 中，设 1、2 顶点是固定的已知点，3、4 顶点为待定点。当压缩或拉伸对角线 d 时，面积S 也将随之变小或变大。2.1.1 对角线长度的最小值在图1的基础上，使顶点3向顶点1靠近，不断缩短d。受4个固定边长的约束，凸四边形将会压缩为凹四边形，面积逐渐变小，较短的边会向里折叠，最终有2条边折叠在一条直线上，此时的图形为一个小三角形加一段回头线，d压缩至最小（图2）。图 2 四边形对角线 d 最小时的示意图d2d2d3d4d13f42d 的最小值 dmin如公式（3）所示。dmin=max（|d1-d2|，|d3-d4|）（3）式中：max（）为取最大值函数。根据海伦公式，可求出 d 最小时的面积，当 dmin=|d1-d2|时，如公式（4）所示。Sdddddddddddddminminminminmin?1434343443 （4）当 dmin=|d3-d4|时，如公式（5）所示。Sdddddddddddddminminminminmin?1412121221 （5）面积边长固定的四边形对角线长算法及应用张锦良1胡隽杉2刘凯程1董建浩2（1.东莞市测绘院，广东 东莞 523000；2.漯河市土地与房屋勘察测绘队，河南 漯河 462000）摘 要：不动产宗地界址的恢复可简化为根据四边形的面积和边长确定其对角线长度的问题。本文通过研究四边形对角线长与面积间的关系，建立了求解四边形对角线长的算法模型。在四边形四个边长和面积固定的约束条件下，经迭代计算求解对角线长度。根据已知界址点的坐标用距离交会法计算四边形的两个待定点坐标，从而恢复宗地界址。关键词：宗地界址恢复；四边形对角线长与面积关系；弦截迭代公式；距离交会中图分类号：P20文献标志码：A注：1、2、3、4为顶点；d1、d2、d3、d4分别为边长；d为1-3对角线长；f为2-4对角线长。图 1 典型的四边形d1d2d3f1432dd4中国新技术新产品2024 NO.4（下）-103-工 程 技 术同理，计算 f 的最小值如公式（6）所示。fmin=max（|d2-d3|，|d4-d1|）（6）当 fmin=|d2-d3|时，计算四边形面积如公式（7）所示。Sfddfddfdfddfdminminminminmin?1441414114 （7）当 fmin=|d4-d1|时，计算四边形面积如公式（8）所示。Sfddfddfdfddfdminminminminmin?1423232332 （8）注意图 2，这个时候的 f 已经跑到了四边形的外面，因此这时的 f 变得没有意义，就不能再用 f 去求解四边形了。对角线 d 或 f 必须在四边形内部，当 d 或 f 跑到四边形外面时，使得公式（1）或公式（2）不成立。在 d 达到图 2 所示的最小值后，若将 f 变小，则会使 d也跑到四边形的外面，而且这时的四边形就会自相交，这在图形的拓扑关系中也是不允许存在的。因此，在后续的研究中，本文仅考虑 d 或 f 在四边形内部时的情形。2.1.2 对角线长度的最大值在图 1 的基础上，将顶点 3 向外拉伸，使其与顶点 1 的距离增加，d 不断变大。受 4 个固定边长的约束，图形将会逐渐撑开变为凸四边形，在该过程中，面积逐渐变大。当变化到 4 个顶点共圆时，四边形的面积达到最大极值。此时继续拉伸 d，四边形将会变得瘦长，最终 2 条较短的边被拉伸为一条直线，图形变为一个大三角形，d 拉伸最大（图 3）。此时得到公式（9）。dmax=min（d1+d2，d3+d4）（9）式中：min（）为取最小值函数。同理，计算 f 的最大值如公式（10）所示。fmax=min（d2+d3，d4+d1）（10）2.1.3 四边形面积的最大值已有研究表明，当四边形的4个顶点共圆时面积最大1，最大面积如公式（11）所示。Sddddddddddddddddmax?142341134212431234?（11）如上所述，四边形面积最大时四个顶点共圆，即为圆内接四边形。根据圆内接四边形对角互补的性质可知，图 1中的 123 与 341 之和为。设 123 为，则 341 为-。设面积最大时 1-3 对角线长度为 dSmax，根据余弦定理并顾及 cos（-）=-cos，得出公式（12）和公式（13）。dS2 max=d12+d22+2d1d2cos （12）dS2 max=d32+d42-2d3d4cos （13）解上述方程组可得公式（14）。dSddddddddddddmax?141232423412221234 （14）同理可得，计算面积最大时 2-4 对角线的长度如公式（15）所示。fSddddddddddddmax?142342124132423441 （15）2.2 构建迭代模型2.2.1 S（d）与 S（f）函数曲线的特性分析由上文可知，当 d 或 f 在四边形内部时，将其作为 2 个三角形的共边，以此计算四边形的面积才是有意义的。在观察对角线 d 变化过程中发现，d 与 f 是联动的关系。在 d 和f 均在四边形内部的区间内，d 变大 f 必然变小。反之，f 变大 d 就变小。在该区间内，d 与 f 是互替的关系，即求解 d与求解 f 可以得到相同的四边形。当 d 拉伸至最大（大三角形）时，f 还有继续变小的可能。当 d 为最大值时（如图 3 所示），顶点 2 继续向顶点 4靠近，图形将变形为凹四边形。f 继续变小，图形最终会稳定在一个带回头线的小三角形上。此时d已经在四边形外且开始变小。在该区间内，d 失去其意义，根据 f 来确定四边形的形状。反之，当 f 拉伸至最大，d 仍继续变小时，f 也失去意义，须根据 d 来确定四边形。这两个区间是 d 与 f 的互补区间。典型的四边形对角线长度与面积关系曲线如图4 所示（为清楚地说明 d 和 f 两个对角线的联动关系，将 f 轴的方向指向了与 d 相反的左边）。结合图 4，总结四边形的 S（d）曲线和 S（f）曲线特性，得出以下结论。1）随着对角线长度的值由小变大，四边形面积会经历一个由小变大、达到极值后再由大到小的过程。在上升区间或下降区间内 S（d）和 S（f）都具有单值属性，因此在各自的上升区间或下降区间内，给定 S 就可以迭代出d 或 f。2）S（d）的下降段（粗虚线部分）与 S（f）对应的上升段（细实线部分）是完全等效的，同一个面积值，二者计算的四边形完全一样。反之，S（f）的下降段（细虚线部分）与 S（d）上升段（粗实线部分）对应部分也是这样。在上述区间内，当面积一定时求解 d 与 f 可以得到相同的四边形，二者可以互替。3）对比图 2 与图 3，Sdmin Sdmax，Sfmin sfmax。即图 4 曲线的上升区间的面积值完全可以涵盖下降区间的面积值。由上述分析可以看出，将 S（d）和 S（f）的上升段（各自的实线部分）合并在一起而忽略其下降段（各自的虚线部分），就能得到一个完整的曲线图。当设计程序时，仅考虑 S（d）和 S（f）的上升段，会使解算模型变得较为简单。另外，由图 4 可以看出，S（d）和 S（f）的下降段比较陡，d 或 f 的误差对 S 的影响也会较大。因此，仅用 S（d）和 S（f）上升段的另一个优势可获得较好的精度。图 3 四边形对角线 d 最大时的示意图d34122d1d2d3d4中国新技术新产品2024 NO.4（下）-104-工 程 技 术若SdminSSmax，则可以用S（d）函数在dminddSmax区间迭代 d 的目标值，以此得到四边形的第一组解。若Sfmin S Smax，则可以用 S（f）函数在 fmin f fSmax区间迭代 f 的目标值，以此得到四边形的第二组解。若 S 在 Sdmin和 Sfmin间，则只能有一组解。2.2.2 迭代公式当四边形边长和面积确定后，要根据公式（1）和公式（2）分别迭代 1-3 对角线长度和 2-4 对角线长度。假定对角线长度为 x，根据公式（1）或公式（2）即可得出四边形面积的计算值 S（x）。设四边形的给定面积为 S，则面积残差为：v（x）=S（x）-S。根据不动产调查技术规程，面积精度为 0.01，因此，将面积残差的容许值 ds（面积容差）设定在比其小 10 倍的 0.001。当残差大于容差时，通过增加或减少 x 值，再次计算残差，直至残差小于或等于容差为止，这样即可迭代出 x 的目标值 d 或 f。实际计算时，采用弦截迭代公式，计算过程如公式（16）所示2。xxxxvvv?110110（16）弦截迭代公式需要两个初值：x0和 x1，这两个初值均不能超过目标值。v0和 v1是根据两个初值计算的面积残差。当迭代 d 时，x0=dmin，当迭代 f 时，x0=fmin。x1=x0+dx，其中 dx为一微小变化量。要保证 dx 引起的面积变化量小于容差。以迭代 d 为例，由图 1 和图 3 可以看出，当 f 最大时，d 的变化对面积的影响最大，可表示为 ds=dxfmax，因此取 dx=ds/fmax。同理，当迭代 f 时，dx=ds/dmax。由于设定面积容差值时已经留有充分冗余，因此求出的初值不会超过目标值。2.2.3 迭代过程以迭代对角线 d 为例。确 定 两 个 初 值：x0=dmin、x1=x0+dx，其 中 dx=ds/fmax，ds=0.001。计算初值的四边形面积残差：将 d=x0、d=x1 代入式（1）,计算 d 取两个初值时的四边形面积 S（x0）、S（x1），然后按照残差公式 v（x）=S（x）-S 计算四边形面积的两个残差：v0=v（x0）、v1=v（x1）。将两个初值 x0、x1和两个残差 v0、v1代入公式（14）计算第一次迭代值 x，并按残差公式计算四边形面积残差 v。若 v=dSmax则结束迭代，取 d=x；否则 x0=x1，v0=v（x0），x1=x，v1=v（x1）重复公式（2）公式（4）。3 计算顶点坐标计算对角线 d 和 f 后，根据顶点 1 和 2 的坐标，用距离交会法可求出顶点 3 和 4 的坐标。距离交会法可以计算两个相交圆的交点坐标。以交会顶点 3（x，y）为例，已知点 1（x1，y1）和点 2（x2，y2）到顶点 3 的距离分别为 r1和 r2，通过求解两个圆的联立方程即可得到交点坐标（x，y），如公式（17）所示。xxxa xxb yyyyya yyb xx?12212112212122/-/?（17）式中：arrr?2212222；brrrrrr?122112222122224；r2=（x2-x1）2+（y2-y1）2。需要说明的是，两个相交圆将会有两个交点，其中一个交点与圆心构成的三角形 1-2-3 为顺时针，另一个则为逆时针。宗地图中的界址点是按顺时针进行编号的，因为当对图1所示的四边形的顶点和边长进行编号时也是采用的顺时针，所以公式（17）所取的是顺时针的那个解（将 b 反符号后是逆时针的那个解）。4 计算机程序实现过程计算机程序实现过程有以下步骤。1）拾取已定的 1、2点坐标，同时反算边长 d1，录入面积和其他 3 个边长。2）根据四边形“任意三边之和大于第四边”，检查录入的边长数据的合理性。若不合理，则终止。3）计算 dmin、Sdmin、fmin、Sfmin、Smax、dSmax、fSmax、dmax、fmax，同时求出面积的最小值：Smin=min（Sdmin，Sfmin）。4）根据面积的上下限，检查录入面积数据的合理性。若不合理，则终止。5）迭代计算求出 d和（或）f。6）根据距离交会法计算待定顶点坐标。7）在CAD 中展绘计算结果。当解算结果有两组时，须人工判别最合适的结果。5 结语宗地界址恢复作业时，首先将能够确定的宗地界址绘制在 CAD 图上，其次将剩余面积归算到一个待定的四边形上，借助计算 边长和宗地剩余面积确定待定的两个顶点，最后恢复所有宗地界址点。参考文献1 王兵权.论边长确定的四边形面积最大值定理 J.中学教研（数学），2010（9）：25-26.2 喻文建.数值分析与算法 M.2版.北京：清华大学出版社，2015.图 4 四边形对角线长度与面积关系曲线SmaxS（d）曲线S（f）曲线SdmaxdsmaxdmaxfmaxfminsfminsfmaxSmaxSdfSfsmaxSdmindmin