Максимальная частота, присутствующая на изображении в MATLAB

Question

Я предполагаю, что в stackoverflow есть различные формы этого вопроса peresnt. Но я не мог понять, как я могу решить свою проблему.

У меня есть изображение, и я хочу найти частоту содержимого изображения.

img = imread('test.tif'); 
img = rgb2gray(img); 
[N M] = size(img); 
%% Compute power spectrum 
imgf = fftshift(fft2(img)); 
imgfp = (abs(imgf)/(N*M)).^2; 

Я знаю, что для этой цели я должен использовать fft. Но мне было интересно, могу ли я найти максимальную частоту в изображениях с точки зрения определенного значения, например «x cycle/mm» или «x cycles/inch».

Что было бы лучшим способом сделать это?

спасибо.

Answer 1

FFT возвращает данные в массиве, где каждый элемент массива несколько связан с циклами на общую ширину данных (или высоту и т. Д.). Таким образом, вы можете разделить каждый номер бункера FFT на размер изображения в некоторой размерной единице (например, «дюймы»), чтобы получить циклы на единицу измерения (например, на дюйм).

Заметим, что за исключением некоторых очень узко определенных типов изображений (скажем, постоянная амплитуда точно апертурных периодических синусоидальных градиентов), любое содержимое изображения будет забрызгано во всем частотном спектре и диапазоне результатов БПФ. Таким образом, вам, вероятно, придется установить некоторый ненулевой порог для частотного содержимого, прежде чем вы сможете ограничить свою «максимальную» частотную характеристику.