radon

Преобразование Радона

Свернуть все на странице

Синтаксис

R = radon(I)
R = radon(I,theta)
[R,xp] = radon(___)

Описание

R = radon(I) возвращает преобразование Радона R 2-D полутонового изображения I для углов в области значений [ 0 , 179 ] степеней. Преобразование Радона является проекцией интенсивности изображения вдоль радиальной линии, ориентированной под определенным углом.

пример

R = radon(I,theta) возвращает преобразование Радона для углов, заданных theta.

[R,xp] = radon(___) возвращает вектор xp содержит радиальные координаты, соответствующие каждой строке изображения.

Примеры

свернуть все

Открыть Live Script

Сделайте масштаб осей видимым для этого изображения.

iptsetpref('ImshowAxesVisible','on')

Создайте образец изображения. 

I = zeros(100,100);
I(25:75, 25:75) = 1;

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

theta = 0:180;
[R,xp] = radon(I,theta);

Отобразите преобразование.

imshow(R,[],'Xdata',theta,'Ydata',xp,'InitialMagnification','fit')
xlabel('\theta (degrees)')
ylabel('x''')
colormap(gca,hot), colorbar

Figure contains an axes. The axes contains an object of type image.

Сделать масштаб осей невидимым.

iptsetpref('ImshowAxesVisible','off')

Входные параметры

свернуть все

Полутоновое изображение, заданное как 2-D числовая матрица.

Типы данных: single | double | int8 | int16 | int32 | int64 | uint8 | uint16 | uint32 | uint64 | logical

Углы проекции в степенях, заданные как числовой скаляр или числовой вектор.

Типы данных: double

Выходные аргументы

свернуть все

Радоновое преобразование I изображения, вернулся как одно из следующего.

  • Если theta является скаляром, тогда R - числовой вектор-столбец, содержащее преобразование Радона для theta степени.

  • Если theta является вектором, тогда R - матрица, в которой каждый столбец является преобразованием Радона для одного из углов в theta.

Радиальные координаты, соответствующие каждой строке R, возвращенный как числовой вектор. Радиальные координаты являются значениями вдоль оси x '-, которая ориентирована на theta степени против часовой стрелки от оси X. Источником обеих осей является центральный пиксель изображения, который определяется как

floor((size(I)+1)/2)
Для примера в изображении 20 на 30 центральный пиксель равен (10,15).

Алгоритмы

Преобразование Радона изображения является суммой преобразований Радона каждого отдельного пикселя.

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

Вклад каждого субпикселя пропорционально разделяется на два ближайших интервала в соответствии с расстоянием между проецируемым местоположением и центрами интервалов. Если проекция субпикселя достигает центральной точки интервала, то интервал на осях получает полное значение субпикселя или четверть значение пикселя. Если проекция субпикселя достигает границы между двумя интервалами, значение субпикселя разделяется равномерно между интервалами.

Ссылки

[1] Bracewell, Ronald N., Двумерная визуализация, Englewood Cliffs, NJ, Prentice Hall, 1995, pp. 505-537.

[2] Lim, Jae S., Двумерная обработка сигналов и изображений, Englewood Cliffs, NJ, Prentice Hall, 1990, pp. 42-45.

Расширенные возможности

См. также

| | | | |

Представлено до R2006a

Документация по Image Processing Toolbox

Поддержка

Для просмотра документации необходимо авторизоваться на сайте