基于OpenCV的圆周率PI值的估计

好久没有写博客了，今天有空写一下自己做的一个小实验，有兴趣的伙伴一起验证一下。刚刚开始学习计算机视觉，OpenCV没有基础，刚刚接触这门语言，就做几个小实验练习一下。

好了，废话不多说，进入正题。吾理小子估计圆周率的值是通过内接圆完成的，首先简述一下原理。

如上图所示，圆内接于正方形。设圆的半径为R，易得圆的面积和正方形的面积如下：

$S_{a}=Pi*R^{2}$

$S_{s}=4*R^{2}$

由此可得：

$Pi=4*\frac{S^{_{a}}}{S_{s}}$

有上式可知，如果知道圆的面积和正方形的面积，就可以求出圆周率Pi的值。

此时，我们假设有一个未知边长的正方形及其内接圆。我们可以通过向正方形内随机撒入多个噪点，由噪点的个数来近似正方形的面积，用落在圆内的噪点个数来等效圆的面积，由此就可以估算圆周率Pi的值了。

接下来，使用OpenCV来实现上述过程。我将源码贴在文末，供各位小伙伴参考。

下面先给出程序的运行过程，以及估计结果。

上面三幅图分别是估计过程中的三个图，最后一幅图是估计结果3.13908，与真实值3.1415还是比较接近的。可以通过进一步增大正方形的边长来提高估计精度。最后贴上源码。


#include <iostream>
#include <opencv2\opencv.hpp>using namespace std;
using namespace cv;int main()
{const int side = 600;        //定义正方形边长const int npixels = 300000;  //定义随机点个数int i, j;Mat s1 = Mat::zeros(side,side,CV_8UC1);Mat s2 = s1.clone();         //图像s1,s2均为0值，即纯黑circle(s1, Point(side / 2, side / 2), side / 2, 255, -1);//在s1正方形图像中开辟一个内接圆imshow("s1", s1); //显示s1图像for (int k = 0; k < npixels; k++)  //产生随机数，撒在图像s2中{i = rand() % side;j = rand() % side;s2.at<uchar>(i, j) = 255;}Mat r;    //定义输出图像rbitwise_and(s1,s2,r); //图像s1与s2进行与运算imshow("s2", s2);     //显示图像s2imshow("r", r);       //显示图像rint Acircle = countNonZero(r);   //统计图像r中的非零点个数int Square = countNonZero(s2);   //统计图像s2中的非零点个数float pi = 4 * (float)Acircle / Square;  //估计圆周率的值cout << "PI的值：" << pi << endl;  //输出圆周率的值waitKey(0);return 0;}

基于OpenCV的圆周率PI值的估计相关推荐

[Julia语言]使用Chudnovsky 算法快速计算圆周率 Pi (π) 值
测试用的电脑是一台10年老电脑,CPU型号:E3 1230V2,3.3GHZ,4核8线程,8GB内存. 用下面的Julia程序,计算1万位的Pi值,耗时为0.26秒. 作为比较,用Julia实现的另一 ...
[JuliaDelphi]快速计算圆周率 Pi (π) 值
算法参考自: 计算圆周率Pi(π)值,精确到小数点后10000位 - 圆周率10000位 - C++ 爱好者分别以Julia 1.72和Delphi 10.3写的算法如下面所示.在一台10年的老电脑 ...
OpenCV进阶（6）基于OpenCV的深度学习人体姿态估计之多人篇
在之前的文章中,我们使用OpenPose模型对一个人进行人体姿势估计.在这篇文章中,我们将讨论如何进行多人姿势估计. 当一张照片中有多个人时,姿态估计会产生多个独立的关键点.我们需要找出哪组关键点属于 ...
c语言计算pi后1000位,计算圆周率 Pi (π)值, 精确到小数点后 10000 位
大家都知道π=3.1415926--无穷多位, 历史上很多人都在计算这个数, 一直认为是一个非常复杂的问题.现在有了电脑, 这个问题就简单了. 电脑可以利用级数计算出很多高精度的值, 有关级数的问题请 ...
计算圆周率 Pi (π)值
大家都知道π=3.1415926--无穷多位, 历史上很多人都在计算这个数, 一直认为是一个非常复杂的问题.现在有了电脑, 这个问题就简单了. 电脑可以利用级数计算出很多高精度的值, 有关级数的问题请 ...
python 圆周率代码_基于Python计算圆周率pi代码实例
一计算公式: 二实现代码 (1) import math from tqdm import tqdm import time total,s,n,t=0.0,1,1.0,1.0 while(mat ...
查看圆周率 Pi值 Pi的前百位 Pi的前百万位
早上忽然看到PI的热搜,就想着看看PI的前一百位甚至更多位是多少,用spark-sql.presto和doris的pi()函数尝试了下,都只给出十几位,使用hadoop提供的example.jar能算 ...
C语言程序设计入门08——初识格式化输出1：格式化输出PI值
格式化输出函数 printf() 作为 C 程序的标准格式化输出函数,printf()的函数声明包含在头文件 <stdio.h> 中,函数的一般形式为: printf ( <待打印 ...
蒙特卡洛方法的使用:计算pi值与积分计算
参考: https://blog.csdn.net/Cuixinyang19_/article/details/79888981 # -*- coding: utf-8 -*- "" ...

基于OpenCV的圆周率PI值的估计

基于OpenCV的圆周率PI值的估计相关推荐

最新文章

热门文章