Фильтры шестеренчатой системы
psi = filterbank(sls)sls как 3-D массив с реальным значением. Размер массива M -by- N -by- K, где M и N являются значениями двухэлементного вектора-строки ImageSize sls. K - количество срезов, включая lowpass, K = numshears(sls) + 1.
Загрузка изображения. Создайте систему shearlet, которая может использоваться с изображением.
load clown sls = shearletSystem('ImageSize',size(X))
sls =
shearletSystem with properties:
ImageSize: [200 320]
NumScales: 4
PreserveEnergy: 0
TransformType: 'real'
FilterBoundary: 'periodic'
Precision: 'double'
Получите фильтры резьбы, заданные системой резьбы.
psi = filterbank(sls);
Получите преобразование shearlet изображения.
cfs = sheart2(sls,X);
Подтвердите, что размер третьей размерности psi равен размеру третьей размерности cfs.
[size(psi,3) size(cfs,3)]
ans = 1×2
41 41
Создайте комплексную систему шестеренок для 256 на 256 изображений с усечёнными контурами и четырьмя шкалами.
sls = shearletSystem('TransformType','complex',... 'ImageSize',[256 256],... 'FilterBoundary','truncated',... 'NumScales',4);
Получите фильтры резьбы и их геометрические интерпретации.
[psi,scale,shear,cone] = filterbank(sls);
Различные шкалы могут иметь разное количество сдвигов. Шкалы в системе shearlet варьируются от 0 до 3 включительно. Найдите область значений сдвигов на шкалу.
for k=0:3 ind = find(scale==k); fprintf('scale: %d shear range: [%d:%d]\n',... k,min(shear(ind)),max(shear(ind))) end
scale: 0 shear range: [-1:1] scale: 1 shear range: [-2:2] scale: 2 shear range: [-2:2] scale: 3 shear range: [-3:3]
Отображение уникальных меток конуса частоты.
unique(cone)
ans = 5x1 cell
{'B'}
{'L'}
{'R'}
{'T'}
{'X'}
Найдите резьбовой фильтр в шкале 2 с параметром сдвига, равным 1 и поддержка которого находится в 'R' частотный конус.
vScale = 2;
vShear = 1;
vCone = 'R';
ind = (scale'==vScale)&(shear'==vShear)&([cone{:}]==vCone);
shFilter = psi(:,:,ind);Постройте график резьбового фильтра.
omegax = -1/2:1/256:1/2-1/256; omegay = omegax; surf(omegax,flip(omegay),shFilter,'EdgeColor','none') view(0,90) xlabel('\omega_x') ylabel('\omega_y') str = sprintf('Shearlet Filter: Scale %d / Shear %d / Cone %c',... vScale,vShear,vCone); title(str) axis equal axis tight
Постройте графики поддержек всех резьбовых фильтров шкалы 2 с параметрами сдвига, равными 2. Области, где фильтр поддерживает перекрытие, имеют максимальную яркость.
vScale = 2; vShear = 2; ind = find((scale==vScale).*(abs(shear)==vShear)); tmp = zeros(size(psi,1),size(psi,2)); for k=1:length(ind) tmp = tmp+(psi(:,:,ind(k))~=0); end surf(omegax,flip(omegay),double(tmp),'EdgeColor','none') view(0,90) xlabel('\omega_x') ylabel('\omega_y') str = sprintf('Filter Supports: Scale %d / Abs(Shear) %d',... vScale,vShear); title(str) axis equal axis tight colormap gray colorbar
