Вычислите преобразование радона изображения

Считайте изображение в рабочую область. Преобразуйте его в полутоновое изображение.

I = fitsread('solarspectra.fts');
I = rescale(I);

Отобразите оригинальное изображение.

figure
imshow(I)
title('Original Image')

Вычислите двоичное изображение краев с помощью edge функция. Отобразите бинарное изображение, возвращенное edge функция.

BW = edge(I);
figure
imshow(BW)
title('Edges of Original Image')

Вычислите преобразование радона изображения, с помощью radon функция и отображение преобразование. Местоположения peaks в преобразовании соответствуют местоположениям прямых линий в оригинальном изображении.

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

Самый сильный пик в R соответствует $\theta =1$ степень и x' = в-80 пикселях от центра.

Интерпретация Peaks радона преобразовывает

Чтобы визуализировать этот пик на исходном рисунке, найдите центр изображения, обозначенного синим крестом наложенным на изображении ниже. Красная пунктирная линия является радиальной линией, которая проходит через центр под углом $\theta =1$ градус. Если вы перемещаетесь вдоль этой линии в-80 пикселях от центра (влево), радиальная линия перпендикулярно пересекает твердую красную линию. Эта твердая красная линия является прямой линией с самым сильным сигналом в Радоне, преобразовывают.

Чтобы интерпретировать Радон преобразовывают далее, исследуют следующие четыре самых сильных peaks в R.

Два сильных peaks в R найдены в $\theta =1$ степень, при смещениях-84 и в-87 пикселях от центра. Этот peaks соответствует этим двум красным линиям слева от самой сильной линии, наложенной на изображении ниже.

Два других сильных peaks найдены около центра R. Этот peaks расположен в $\theta =91$ степени, со смещениями-8 и в-44 пикселях от центра. Зеленая пунктирная линия в изображении ниже является радиальной линией, проходящей через центр под углом 91 градуса. Если вы путешествуете вдоль радиальной линии на расстояние-8 и в-44 пикселях от центра, то радиальная линия перпендикулярно пересекает чисто зеленые линии. Эти чисто зеленые линии соответствуют сильному peaks в R.

Более слабые линии в изображении относятся к более слабому peaks в R.