图像识别1-Radon变换学习-MATLAB得到热力图
1. fitsread
- 感觉是联网/自带的图片，我甚至没有保存original image程序也可以自动读取/处理图片
I = fitsread('solarspectra.fts');
I = rescale(I);figure
imshow(I)
title('Original Image')
BW = edge(I);
figure
imshow(BW)
title('Edges of Original Image')
theta = 0:179;
[R,xp] = radon(BW,theta);figure
imagesc(theta, xp, R); colormap(hot);
xlabel('\theta (degrees)');
ylabel('x^{\prime} (pixels from center)');
title('R_{\theta} (x^{\prime})');
colorbar
 |  |
|---|
2. imread
%% Radon Transform
function RadonTest()srcImage=imread('radon2.bmp');#原代码：[fileName=input('Input image file name:');srcImage=imread(fileName);]
grayImage=rgb2gray(srcImage);
cannyImage=edge(grayImage,'canny');theta=0:180;
[R,x]=radon(cannyImage,theta);
figure,imagesc(theta,x,R);
title('R_theta X');
xlabel('theta(degree)');
ylabel('X\prime');
colormap(hot);
colorbar;end
对图中直线段进行radon变换：
- 开始报错“MATLAB处理图像时出错：错误使用 rgb2gray＞parse_inputs (line 80)MAP 必须为 m x 3 的数组。”
- 错误原因：读入图像不是GRB三通道的
- 解决办法：修改为24位位图便可
- 右击打开方式，利用画图软件打开，另存为24位图
- radon变换是一种原始灰度图像到（p,o）二维矩阵的映射，映射关系是灰度图像的像素值对参数（p,o）的直线的线积分，或者叫投影，可理解为垂直于这个直线方向的像素值累计求和。
- 这样的话，对原始灰度图像进行二值化处理以后，再进行拉冬变换，得到结果矩阵，找到极值点，就可以确定直线的位置。因为垂直于直线的方向上，线积分会呈现为极值点。个人觉得这里的线积分是不区分方向的，也就是不区分正向负向的概念。
B站关于Radon变换的可视化解释＆CT原理
https://www.bilibili.com/video/BV12J411n7Wb/?spm_id_from=333.1007.top_right_bar_window_history.content.click&vd_source=83f28538cba7cf90472dcefa298c445e
3. 论文内容复现
%% Recurrence code
function Recurrence() srcImage=imread('girl.png');
figure
imshow(srcImage)
title('Original Image')grayImage=rgb2gray(srcImage);
figure
imshow(grayImage)
title('Gray Image')cannyImage=edge(grayImage,'canny');
figure
imshow(cannyImage)
title('Canny Image')theta=0:180;
[R,x]=radon(cannyImage,theta);
figure,imagesc(theta,x,R);
title('R_theta X');
xlabel('theta(degree)');
ylabel('X\prime');
colormap(hot);
colorbar; end
本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！