用C编写的程序 用四阶龙格-库塔方法求解初值问题dy/dx=xcos(x+y)

1 原理  2 检测步骤 将参数空间(ρ,θ) 量化成m*n(m为ρ的等份数，n为θ的等份数)个单元，并设置累加器矩阵，初始值为0； 对图像边界上的每一个点（x,y）带入ρ=xcosθ+ysinθ，求得每个θ对应的ρ值，并在ρ和θ所对应的单元，将累加器加1，即：Q(i,j)=Q(i,j)+1； 检验参数空间中每个累加器的值，累加器最大的单元所对应的ρ和θ即为直角坐标系中直线方程的参数。  3 接口 image：二值图像，canny边缘检测输出。这里是result。 rho: 以像素为单位的距

在图片处理中，霍夫变换主要是用来检测图片中的几何形状，包括直线、圆、椭圆等。 在skimage中，霍夫变换是放在tranform模块内，本篇主要讲解霍夫线变换。 对于平面中的一条直线，在笛卡尔坐标系中，可用y=mx+b来表示，其中m为斜率，b为截距。但是如果直线是一条垂直线，则m为无穷大，所有通常我们在另一坐标系中表示直线，即极坐标系下的r=xcos(theta)+ysin(theta)。即可用（r,theta）来表示一条直线。其中r为该直线到原点的距离，theta为该直线的垂线与x轴的夹角。如