Exponenta Banner
exponenta event banner

littlewoodPaleySum

Сумма Литлвуда-Палей

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

Синтаксис

lpsum = littlewoodPaleySum(sf)
lpsum = littlewoodPaleySum(sf,fb)
[lpsum,f] = littlewoodPaleySum(___)

Описание

lpsum = littlewoodPaleySum(sf) возвращает сумму Литлвуда-Палей для 2D наборов фильтров в 2D вейвлете, рассеивающем разложение sf. lpsum является M-by-N-by-Nfb матрица, где M-by-N является матричным размером заполненных фильтров, и Nfb является количеством наборов фильтров.

Поскольку рассеивающееся преобразование является сжимающимся, суммы Литлвуда-Палей не превышают 1.

lpsum = littlewoodPaleySum(sf,fb) возвращает сумму Литлвуда-Палей для заданных наборов фильтров fb. fb является положительным целым числом или вектором положительных целых чисел между 1 и numfilterbanks(sf) включительно. lpsum является M-by-N-by-L матрица, где L является количеством уникальных элементов в fb.

пример

[lpsum,f] = littlewoodPaleySum(___) возвращает пространственные частоты за сумму Литлвуда-Палей. f является матрицей 2D столбца с первым столбцом, содержащим пространственные частоты в x - направление и второй столбец, содержащий пространственные частоты в y - направление.

Примеры

свернуть все

Открыть скрипт

Этот пример показывает, как получить и отобразить сумму Литлвуда-Палей среды рассеивания изображений.

Создайте рассеивающуюся среду с двумя наборами фильтров и добротностями 2 и 1, соответственно.

sf = waveletScattering2('QualityFactors',[2 1]);

Получите суммы Литлвуда-Палей и пространственные частоты этих двух наборов фильтров. Отобразите максимальное значение сумм. Поскольку рассеивающееся преобразование является сжимающимся, суммы не превышают 1.

[lpsum,f] = littlewoodPaleySum(sf);
max(max(lpsum(:,:,1)))
max(max(lpsum(:,:,2)))

ans =

    1.0000


ans =

    1.0000

Отобразите сумму Литлвуда-Палей второго набора фильтров с нулевой сосредоточенной частотой. Обратите внимание, что 2D набор фильтров Morlet, используемый в рассеивающемся преобразовании, не разработан, чтобы получить самые высокие пространственные частоты совместно в x-и направлениях Y.

f(f>1/2) = f(f>1/2)-1;
surf(fftshift(f(:,1)),fftshift(f(:,2)),fftshift(lpsum(:,:,2)))
shading interp
view(0,90)
xlabel('f_x')
ylabel('f_y')
colorbar
title('Q=1')

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

свернуть все

Рассеивание среды разложения, заданной как объект waveletScattering2.

Индекс набора фильтров в изображении, рассеивающем среду разложения, заданную как положительное целое число или вектор положительных целых чисел между 1 и numfilterbanks(sf) включительно. Количество наборов фильтров в sf равно количеству заданного QualityFactors в sf.

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

свернуть все

Сумма Литлвуда-Палей для наборов фильтров в изображении, рассеивающем среду разложения sf, возвращенный как 3-D матрица с действительным знаком. lpsum является M-by-N-by-L матрица, где M-by-N является матричным размером заполненных фильтров, и L не превышает количество наборов фильтров в sf.

Частоты за сумму Литлвуда-Палей, возвращенную как матрица 2D столбца с действительным знаком. Частоты находятся в циклах на пиксель. Первый столбец f содержит пространственные частоты в x - направление, и второй столбец содержит пространственные частоты в y - направление. В этом соглашении преобразование Фурье является 1-периодическим в обеих переменных Фурье.

Смотрите также

Введенный в R2019a

Документация Wavelet Toolbox
Поддержка