一种误差较小的轮廓面积计算算法

1.背景

        基于微分思想的轮廓面积计算方法之一是将多边形轮廓边与X轴会Y轴进行围合,形成一个个梯形,每个梯形的面积有符号,累计求和即得到多边形轮廓的面积。详见博主之前的文章,

记录导致计算轮廓面积出错的一个坑点-CSDN博客文章浏览阅读377次,点赞4次,收藏9次。计算轮廓面积是常见的几何算法话题,获取轮廓面积、计算轮廓法线等场景会涉及到。计算轮廓面积的方法有很多,一种常用的是微积分思路的分段求和办法,即组成轮廓的每条线段与X轴或Y轴进行有向投影,轮廓边线与X轴或Y轴的投影之和即为轮廓的有向面积。其中第0、1、6的投影为正,第2、3、4、5的投影为负,累积即为轮廓面积(绿色填充部分)。如上所示,如果关心面积的正负时,一般将轮廓移到第一或第四象限再求面积。满足上述条件时,逆时针轮廓的面积为正,顺时针轮廓的面积为负。https://blog.csdn.net/baidu_38621657/article/details/140853728

        如果出现轮廓跨Y轴的情况,则需要将轮廓平移到Y轴右侧,此时面积符号为逆正顺负(当然平移到Y轴左侧也可以,如果那样所计算面积的符号为“顺正逆负”)。

        但是这种方法可能会出现精度误差较大,从而导致面积值和符号计算错误的情况,究其原因是因为对轮廓做了偏置处理,而这种处理会带来误差,偏置的距离越大误差可能越大,详细原因见博主此前的文章:《关于float浮点值二进制存储和运算精度损失的话题》。

关于float浮点值二进制存储和运算精度损失的话题_float的二进制存储和输出-CSDN博客文章浏览阅读1k次,点赞8次,收藏18次。浮点值的存储、运算都可能会带来精度损失,了解精度损失背后的机制原因方便我们更好的了解什么情况下会发生精度损失、什么情况下精度损失较大,以及思考怎么避免或减少精度损失。_float的二进制存储和输出https://blog.csdn.net/baidu_38621657/article/details/141027014        需要尽量减少和简化计算过程,保持算法的简便整洁,注意这里的简化不是接受更大的误差的意思,指的是简化逻辑,使用精简准确的公式或算法过程。

对于此种情况(轮廓跨越Y轴)该怎么处理,避免不必要的操作过程,从而减少精度损失呢?

2. 公式推导

博主比较懒, 推导过程就不详细讲述了,有兴趣可以和博主联系或关注博主公众号。

3. 代码

/// <summary>
/// 求有向面积,顺时针为负,逆时针为正
/// </summary>
/// <param name="polygon"></param>
/// <returns></returns>
double GetArea2(const list<Line>& polygon, const Transform& trsW2L)
{//  目前只支持边为线段的轮廓double dArea = 0;for (auto curve : polygon){Vector3f pt0, pt1;Transform::MultPoint(trsW2L, curve.pt0, pt0);Transform::MultPoint(trsW2L, curve.pt1, pt1);double areaThis = 0.0;if (pt0.X >= 0.f){if (pt1.X >= 0.f)areaThis = (pt0.X + pt1.X) * (pt1.Y - pt0.Y) * 0.5;else{double s1 = (pt1.Y - pt0.Y) * pt0.X * pt0.X * 0.5 / (pt0.X - pt1.X);double s2 = -(pt1.Y - pt0.Y) * pt1.X * pt1.X * 0.5 / (pt0.X - pt1.X);areaThis = s1 + s2;}}else{if (pt1.X >= 0.f){double s1 = -(pt1.Y - pt0.Y) * pt0.X * pt0.X * 0.5 / (pt1.X - pt0.X);double s2 = (pt1.Y - pt0.Y) * pt1.X * pt1.X * 0.5 / (pt1.X - pt0.X);areaThis = s1 + s2;}else{areaThis = (pt0.X + pt1.X) * (pt1.Y - pt0.Y) * 0.5;}}dArea += areaThis;}return dArea;
}

各种象限和跨象限情况的单元测试均通过。

欢迎交流,相互学习,公众号:geometrylib

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/bicheng/53197.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

C++ Qt进程间通信机制之QRO、QRemoteObjectHost

文章目录 QRO_serverQRO_client深入浅出C++ Qt开发技术专栏:https://blog.csdn.net/yao_hou/category_9276099.html Qt高级开发视频教程 https://edu.csdn.net/course/detail/37447?spm=1001.2014.3001.5507 C++零基础入门视频教程 https://edu.csdn.net/course/detail/379…

Rust 学习笔记 3:一般性编程概念

上一篇&#xff1a;Rust 学习笔记 2&#xff1a;猜数字游戏 文章目录 1. 前言2. 背景3. Rust 中的一般性编程概念3.1 变量及其可变性(Mutability)3.1.1 变量定义3.1.2 常量3.1.3 变量隐藏(Shadowing) 3.2 基本类型3.2.1 标量(scalar)类型3.2.1.1 整型(Integer Types)3.2.1.2 浮…

项目拆解:短视频冷门赛道—ai绘画+温馨小屋,引流变现全攻略

在这个快节奏的时代&#xff0c;工作、学习、家庭的重担仿佛三座大山&#xff0c;让人喘不过气&#xff0c;心情时常跌入谷底。就像蜗牛遇到威胁会缩进壳里&#xff0c;我们也会在疲惫和忧虑时&#xff0c;渴望一个属于自己的温暖小窝&#xff0c;来安放疲惫的心灵。而自媒体平…

Flink 1.14.* Flink窗口创建和窗口计算源码

解析Flink如何创建的窗口&#xff0c;和以聚合函数为例&#xff0c;窗口如何计算聚合函数 一、构建不同窗口的build类1、全局窗口2、创建按键分流后的窗口 二、在使用窗口处理数据流时&#xff0c;不同窗口创建的都是窗口算子WindowOperator1、聚合函数实现2、创建全局窗口(入参…

SpringFrameWork学习笔记

本笔记基于【尚硅谷新版SSM框架全套视频教程&#xff0c;Spring6SpringBoot3最新SSM企业级开发】https://www.bilibili.com/video/BV1AP411s7D7?vd_sourcea91dafe0f846ad7bd19625e392cf76d8 总结 资料获取网址&#xff1a;https://www.wolai.com/v5Kuct5ZtPeVBk4NBUGBWF 技术…

Java项目: 基于SpringBoot+mysql房产销售系统 (含源码+数据库+开题报告+答辩PPT+毕业论文)

一、项目简介 本项目是一套基于SpringBootmysql房产销售系统 包含&#xff1a;项目源码、数据库脚本等&#xff0c;该项目附带全部源码可作为毕设使用。 项目都经过严格调试&#xff0c;eclipse或者idea 确保可以运行&#xff01; 该系统功能完善、界面美观、操作简单、功能齐…

数学基础 -- 线性代数之LU分解

LU分解 LU分解&#xff08;LU Decomposition&#xff09;是线性代数中非常重要的一种矩阵分解方法。它将一个方阵分解为一个下三角矩阵&#xff08;L矩阵&#xff09;和一个上三角矩阵&#xff08;U矩阵&#xff09;的乘积。在数值线性代数中&#xff0c;LU分解广泛用于求解线…

Halcon基于灰度值的模板匹配

Halcon基于灰度值的模板匹配 基于灰度值的模板匹配是最经典的模板匹配算法&#xff0c;也是最早提出来的模板匹配算法。这种算法的根本思想是&#xff0c;计算模板图像与检测图像之间的像素灰度差值的绝对值总和&#xff08;SAD方法&#xff09;或者平方差总和&#xff08;SSD…

Spring解析

目录 容器与 bean 1) 容器接口 演示1 - BeanFactory 与 ApplicationContext 的区别 代码参考 收获&#x1f4a1; 演示2 - 国际化 2) 容器实现 演示1 - DefaultListableBeanFactory 代码参考 收获&#x1f4a1; 演示2 - 常见 ApplicationContext 实现 代码参考 收获…

ico格式怎么转换?5个软件让你轻松转换文件格式

ico格式怎么转换&#xff1f;5个软件让你轻松转换文件格式 ICO格式是常用于网站图标和应用程序图标的文件格式&#xff0c;虽然它很常见&#xff0c;但并非所有图像编辑软件都支持直接保存为ICO格式。如果你需要将其他格式的图片&#xff08;如PNG、JPG等&#xff09;转换为IC…

读书学习笔记入门 # Datawhale X 李宏毅苹果书 AI夏令营

文章目录 学习目标&#xff1a;学习内容&#xff1a;Task 1 通过案例了解机器学习机器学习&#xff08;Machine Learning&#xff0c;ML&#xff09;和深度学习&#xff08;Deep Learning&#xff0c;DL&#xff09;的基本概念什么是回归&#xff08;regression&#xff09;什么…

深入解析Linux轻量级进程:线程的概念、原理、优缺点及其与进程的关系与区别

&#x1f351;个人主页&#xff1a;Jupiter. &#x1f680; 所属专栏&#xff1a;Linux从入门到进阶 欢迎大家点赞收藏评论&#x1f60a; 目录 &#x1f4da;Linux线程&#x1f4d5;什么是线程*可以使用多进程去并发的执行一个进程的代码&#xff0c;那为什么要由线程呢&#x…

基于CloudflareSpeedTest项目实现git clone加速

1.网络测速 「自选优选 IP」测试 Cloudflare CDN 延迟和速度&#xff0c;获取最快 IP 更多内容参考项目&#xff1a;https://github.com/XIU2/CloudflareSpeedTest 国外很多网站都在使用 Cloudflare CDN&#xff0c;但分配给中国内地访客的 IP 并不友好&#xff08;延迟高、丢…

Pixelmator Pro for Mac 专业图像处理软件【媲美PS的修图软件】

Mac分享吧 文章目录 效果一、下载软件二、开始安装1、双击运行软件&#xff0c;将其从左侧拖入右侧文件夹中&#xff0c;等待安装完毕2、应用程序显示软件图标&#xff0c;表示安装成功 三、运行测试安装完成&#xff01;&#xff01;&#xff01; 效果 一、下载软件 下载软件…

【STM32+HAL库】---- 通用定时器输入捕获PWM信号

硬件开发板&#xff1a;STM32G0B1RET6 软件平台&#xff1a;cubemaxkeilVScode1 新建cubemax工程 1.1 配置系统时钟RCC 1.2 配置定时器 1.2.1 配置输入捕获 选择通用定时器TIM2-Channel 1为输入捕获引脚&#xff0c;对应IO口是PA0,时钟源选择内部时钟源Internal clock,工作模…

Unity实战案例 2D小游戏HappyGlass(模拟水珠)

本案例素材和教程都来自Siki学院&#xff0c;十分感谢教程中的老师 本文仅作学习笔记分享交流&#xff0c;不作任何商业用途 预制体 在这个小案例中&#xff0c;水可以做成圆形但是带碰撞体&#xff0c;碰撞体比图形小一圈&#xff0c;顺便加上Trail renderer组件 材质 将碰撞…

PIL库图像压缩教程

图像压缩教程 本教程展示如何使用Python中的PIL&#xff08;Pillow&#xff09;库进行图像压缩。详细介绍如何加载图像、压缩图像&#xff0c;以及保存压缩后的图像。 步骤包括&#xff1a; 安装所需库加载图像压缩图像保存压缩后的图像 确保环境中安装了Pillow库。如果没有…

Win11 / Win10 系统极化工具,降低游戏延迟效果明显

Win11 / Win10 系统优化工具,降低游戏延迟效果明显 Windows 系统优化就是精简系统一些功能组件、对一些系统功能进行设置等&#xff0c;这样可以减少不必要的硬件资源占用。 全面的系统优化功能外&#xff0c;据不少网友表示通过优化后 CS GO 游戏降低输入延迟效果明显。 免费…

shell脚本示例(用于学习shell语法)

#!/bin/bash # 定义变量 # 设置x264的安装位置 PREFIX"./x264build" # 交叉编译器bin路径,使用aarch64-linux-gnu-gcc -v命令查看 TOOLCHAIN_DIR"/opt/gcc-aarch64-linux-gnu/bin" # 工具链前缀&#xff0c;用于指定交叉编译工具的路径 TOOLCHAIN"…

沃飞长空联合极氪亮相2024世界动力电池大会

9月1日至2日&#xff0c;2024世界动力电池大会在四川宜宾举办&#xff0c;沃飞长空与同属吉利控股集团旗下的新时代豪华科技品牌极氪汽车一同亮相。 现场&#xff0c;双方携手展出了AE200电动垂直起降航空器、极氪009光辉版、极氪001&#xff0c;以及极氪能源、金砖电池、威睿…