计算图像的氡变换

将图像读入工作区。将其转换为灰度图像。

我= fitsread（'solarspectra.fts'）;I =重新调节(我);

显示原始图像。

数字imshow（i）标题（'原始图像​​'）

使用该计算二进制边缘图像边缘功能。显示返回的二进制图像边缘功能。

bw =边缘（i）;图imshow（bw）标题（'原始图像​​的边缘'）

计算图像的氡变换，使用氡功能，并显示变换。变换中峰的位置对应于原始图像中的直线的位置。

Theta = 0：179;[r，xp] =氡（bw，theta）;

显示Radon变换的结果。

图ImageC（THETA，XP，R）;Colormap（热）;Xlabel（'\ theta（deacte）'）;ylabel（'x^{\prime}(从中心的像素)'）;标题('r _ {\ theta}（x ^ {\ prime}）'）;彩色栏

最强大的峰值R.对应于 $\theta =1$学位和X'= -80像素来自中心。

解释Radon变换的峰值

要在原始图中可视化此峰值，请找到图像的中心，由覆盖在下面的图像上的蓝色十字架。红色虚线是以一定角度穿过中心的径向线 $\theta =1$程度。如果从中心（向左）沿着该线路-80像素行进，径向线垂直与固体红线相交。这种固体红线是具有氡变换中最强信号的直线。

要进一步解释Radon变换，请检查接下来的四个最强的峰值R.。

两个强大的山峰R.被发现了 $\theta =1$学位，在-84和-87像素的偏移中，来自中心的偏移量。这些峰值对应于最强线左侧的两个红线，覆盖在下面的图像上。

在中心附近发现了另外两个强峰R.。这些峰位位于 $\theta =91.$学位，来自中心的-8和-44像素的偏移。下面的图像中的绿色虚线是以91度的角度通过中心的径向线。如果沿着径向线沿着中心距离-8和-44像素的距离，则径向线垂直与固体绿线相交。这些实体绿线对应于强大的峰值R.。

图像中的微弱线与较弱的峰值相关R.。