与Radon radon函数对应的格式为[R，XP]=radon(I，theta )；
其中，r为线积分值，xp为径向坐标，theta为投影角度；
图像投影是指在某个方向上对图像进行线性积分(累积相加)。 如果将图像考虑为二维函数f(x，y )，则该投影是特定方向上的线性积分，例如f ) x，y在垂直方向上的线性积分是在x轴上的投影。 f(x，y )的水平方向的线积分是向y轴的投影。 通过这些投影，可以取得图像向指定方向的突出特性，这有可能用于图像图案识别等处理。
radon(yxdxte )算法通过定向和投影叠加，可以找到最大投影值时的角度，并确定图像的倾斜角度。
Radon变换(yxdxte变换)是指将数字图像矩阵投影变换到特定角度的光线方向上。 这意味着可以沿任意角度的theta进行Radon变换。 radon变换是图像在各个方向上的投影。 下图中的f(x，y )表示图像，r ) x '是图像右下角的投影。 数学上对每个投影方向进行线积分，在图像区域中对每个投影方向加上像素即可。
Radon变换的本质是将原始XY平面内的点映射到AB平面上，原始记录XY平面上的一条直线的所有的点在AB平面上都位于同一点AB平面上的点的积累厚度，可以看出XY平面上线的存在性。 这就是大家公认的Radon变换的本质。
在数学理论中，以图像的中心为原点，用p (直线到原点的距离)和某个特定的方向)代替a、b，也就是说，理解为图像在空间上的投影，用参数表示上述直线，则如下所示。
假设存在表示图像像素的函数f(x，y )，如上图所示。 该函数通过直线l区域的积分如下。
这里，ds是这条直线的微分。 上述积分可以写如下
因此，如果给出一组p和，则得到沿l(p，)的积分值，即指定角度下的像素累积相加。 因此，如图所示，Radon变换是函数f(x，y )的线积分。 假设有许多平行于l的线，他们具有相同的theta，但径向坐标p不同，则matlab所具有的radon变换命令的各theta角度的radon变换结果为两个输出项r (一个) theta角度下的radon变换也就是说，对一张图像进行某个特定角度的Radon变换会生成n个线积分值，每个线积分值与径向坐标xp相对应，将各个角度的Radon变换值汇总构成一张Radon变化图。
使用Radon变换检测直线，检测步骤如下。
(1)使用边缘检测函数edge函数计算二值图像；
)2)计算二值图像的Radon变换
)3)寻找Radon变换的局部极大值，这些极大值的位置是原始图像中直线的位置。
bw=edge(I1 )； imshow(bw ) theta=0:179； [R，XP]=radon(bw，theta )； figure，imagesc(theta，xp，r ) Radon变换用于直线检测，优于Hough变换的是Radon变换可以检测非二值图像，而Radon变换可以检测直线。 灰度值高的线段在p空间中形成亮点，低灰度值的直线在p空间中形成暗点，而Hough在p空间中形成暗点
from:https://blog.csdn.net/Yu 132563/article/details/99228303
3359 blog.csdn.net/kong XP _1/article/details/81448210
Hough转换器直线的表示方式
相对于平面上的一直线，笛卡尔坐标系中常见的有点倾斜式、两点式两种表现方法。 但是，在hough变换中，可以考虑另一种表达方法。 使用(r，theta )表示直线。 在此，r是从该直线到原点的距离，theta是该直线的垂线与x轴所成的角度。 如下图所示。 即，表示lhzdwt变换中的一直线的参数为(r，theta )。 变换同一直线上的n个点，原图像空间的n个点在参数空间中对应得到n条正弦曲线，这些曲线在一点交叉。
2.如何判断多个点是否在同一直线上
当我们的对象变成点时，我们知道一个点可以发射出无数条直线，根据lhzdwt变换的直线表达形式，假设这个点为 i, 则通过这个点的直线我们用（ri,thetai）表示。再假设一个点为 j，则通过点 j 的一系列直线我们用（rj,thetaj）表示。我们知道两点决定一条直线，所以这两个点的直线必定有r_i=r_j,theta_i= theta_j的时候。那如果是三个点呢，假设第三个点是k，则通过k点的一系列直线为（r_k,theta_k）,如果三点在一条直线上，那必定有某个r_i=r_j=r_k = r,theta_i = theta_j= theta_k = theta。
在lhzdwt变换检测直线时我们需要找到这样一样直线.
3.如何检测出直线
假设有N个点，我们要检测其中的直线，也就是我们要找到具体的r和theata。对于上面所说的每个点可以通过无数条直线，这里我们设为n条（通常 n = 180），则我们一起可以找到N*N个（r, theata），对这N*N个（r,theata）,我们可以利用统计学，统计到在theta=某个值theta_i时，多个点的r近似相等于r_i。也就是说这多个点都在直线（r_i,theta_i）上。
4.举例说明
如果空间中有3个点，如何判断这三个点在不在一个直线上，如果在，这条直线的位置为？这个例子中，对于每个点均求过该点的6条直线的（r,theta）坐标，共求了3*6个（r,theta）坐标。可以发现在theta=60时，三个点的r都近似为80.7，由此可判定这三个点都在直线（80.7，60）上。通过 r-o-theta 坐标系可更直观表示这种关系，如下图：图中三个点的（r,theta）曲线汇集在一起，该交点就是同时经过这三个点的直线。
在实际的直线检测情况中，如果超过一定数目的点拥有相同的（r,theta）坐标，那么就可以判定此处有一条直线。在r-O-theta 坐标系图中，明显的交汇点就标示一条检测出的直线。如下图，可以判定出平面上的点共构成了两条直线，即检测出两条直线。
亮点对应笛卡尔坐标系下（图像）的一条直线，越亮，代表在这条直线上的非零像素数越多。（原图像是二值图像，亮度代表像素个数）
lhzdwt变换直线检测的步骤：
如果一幅图像中的像素构成一条直线，那么这些像素坐标值（x, y）在lhzdwt空间对应的曲线一定相交于一个点，所以我们只需要将图像中的所有像素点（坐标值）变换成lhzdwt空间的曲线，并在lhzdwt空间检测曲线交点就可以确定直线了。
1、对边缘二值化图像进行lhzdwt空间变换；
2、在4邻域内找到lhzdwt空间变换的极大值；
3、对这些极大值安装由大到小顺序进行排序，极大值越大，越有可能是直线；
4、输出直线。
MATLAB实现：
close all;clear all;I = imread('scratch.tif');figure;subplot(1,3,1);imshow(I);BW = edge(I,'canny');%Canny提取图像边界[H,T,R] = hough(BW);%计算二值图像的标准lhzdwt变换，H为lhzdwt变换矩阵，I,R为计算lhzdwt变换的角度和半径值subplot(1,3,2);imshow(imadjust(mat2gray(H)),[],'XData',T,'YData',R,'InitialMagnification','fit');%hough变换的图像xlabel('theta'), ylabel('rho');axis on,axis square,hold on;P = houghpeaks(H,5,'threshold',ceil(0.3*max(H(:)))); % 从lhzdwt变换矩阵H中提取5个极值点x = T(P(:,2)); y = R(P(:,1));plot(x,y,'s','color','white');%标出极值点lines=houghlines(BW,T,R,P);%提取线段subplot(1,3,3);imshow(I), hold on;for k = 1:length(lines)xy = [lines(k).point1; lines(k).point2]; plot(xy(:,1),xy(:,2),'LineWidth',2,'Color','green');%画出线段endfrom：https://blog.csdn.net/u010839382/article/details/30810325
下面详细介绍关于hough变换的函数用法：
[H,T,R] = hough(BW,'RhoResolution',0.5,'Theta',-90:0.5:89.5);
计算二值图像BW的标准hough变换。（我在MATLAB2008中，BW不是二值图像，也能处理）H为hough转换矩阵。theta和rho是一维数组，theta记录半径，rho记录角度。 'RhoResolution' 指定一个位于0到 norm(size(BW))间的实散值去作为hough转换的径（rho轴）的步长，默认值为1。 'Theta' 时，用ParameterValue指定 一个矢量去作为hough转换的角度范围，可接受的范围为-90<=theta<90。默认值为-90:89。 0.5表示沿ρ轴的lhzdwt变换的间隔。
左图对应mat矩阵（x-列数，y-行数，角度-（90到0到-90，例如当直线垂直x轴时，角度为0）），右图为生成的H矩阵。
请注意：MATLAB2008和2014的版本中hough的参数 有些许不同
houghpeaks：在lhzdwt变换矩阵里找极值点。线检测和连接用的lhzdwt变换的第一步是用高的计数寻找累加单元(工具箱文本把高的计数单元作为峰值)。因为存在lhzdwt变换参数空间中的量化和典型图像的边缘并不是很完美的直线这样的事实，lhzdwt变换的峰值倾向于相比lhzdwt变换单元更多。函数houghpeaks用任意默认语法来寻找指定的峰值数：
peaks = houghpeaks(H, numpeaks)%提取极值点的个数
peaks = houghpeaks(..., param1, val1,param2, val2)
P = houghpeaks(H,5,'threshold',ceil(0.3*max(H(:))));%'threshold'表示：低于0.3*max(H(:))的像素不会被作为极值点。
houghlines：从lhzdwt变换矩阵中提取线段。一旦一组候选的峰值在lhzdwt变换中被识别出来，如果存在与这些峰值相关的有意义的线段，剩下的就是决定线的起始点和终点。函数houghlines用默认的语法执行这个任务
lines = houghlines(BW, theta, rho,peaks) %peaks为 houghpeaks函数返回的结果
lines = houghlines(..., param1, val1,param2, val2)
lines = houghlines(BW,theta,rho,P,'FillGap',5,'MinLength',7);
FillGap是正的标量，指定了与相同的lhzdwt变换相关的两条线段的距离（距离是怎么算的？若是两条线段在同一直线上，这里的距离是否指的是：两线段的端点之间的距离呢？）。当两条线段之间的距离小于指定的值时，函数houghlines把线段合并为一条线段(默认的距离是20个像素)。
MinLength是正的标量，指定合并的线是保留还是丢弃。如果合并的线比val2指定的值短，就丢弃(默认值是40)
关于hough线检测的演示：https://blog.csdn.net/xsjwangyb/article/details/10917945
from：https://blog.csdn.net/zhanshen112/article/details/79949257?utm_medium=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-2.channel_param&depth_1-utm_source=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-2.channel_param
imadjust(mat2gray(H))：
mat2gray(H功能是实现图像矩阵的归一化操作，使矩阵的每个元素的值都在0和1之间。
imadjust：将灰度图像 I 中的亮度值映射到 J 中的新值并使 1% 的数据是在低高强度和饱和；
from：https://blog.csdn.net/zaishuiyifangxym/article/details/78244308?utm_medium=distribute.pc_aggpage_search_result.none-task-blog-2~all~first_rank_v2~rank_v25-3-78244308.nonecase&utm_term=imadjust%E5%87%BD%E6%95%B0%E4%BD%9C%E7%94%A8
hough线检测OpenCV实现参见：https://blog.csdn.net/qq_30815237/article/details/86750576?ops_request_misc=%257B%2522request%255Fid%2522%253A%2522159849654319724846455228%2522%252C%2522scm%2522%253A%252220140713.130102334.pc%255Fblog.%2522%257D&request_id=159849654319724846455228&biz_id=0&utm_medium=distribute.pc_search_result.none-task-blog-2~blog~first_rank_v2~rank_blog_default-2-86750576.pc_v2_rank_blog_default&utm_term=hough&spm=1018.2118.3001.4187