Exponenta 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 - направление.

Примеры

свернуть все

Скрипт Open Live Script

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

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

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

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

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

max(max(lpsum(:,:,2)))
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')

Figure contains an axes object. The axes object with title Q=1 contains an object of type surface.

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

свернуть все

Сеть рассеивания вейвлета изображений в виде 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

Поддержка