Этот пример показывает, как взять преобразование shearlet изображения и восстановить изображение используя только коэффициенты, соответствующие нулевому сдвигу.

Загрузка и отображение изображения круга.

load circleGS
imagesc(circleGS)
colormap gray
axis equal
axis tight

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

[numRows,numCols] = size(circleGS);
sls = shearletSystem('ImageSize',[numRows numCols],'FilterBoundary','truncated');
[psi,scale,shear,cone] = filterbank(sls);

Получите преобразование shearlet изображения.

cfs = sheart2(sls,circleGS);

Найдите индексы резьбовых фильтров, которые соответствуют нулевому сдвигу. Следует иметь в виду, что фильтр lowpass также соответствует нулевому сдвигу.

ind = find((shear==0).*(scale~=-1))'

ind = 1×10

     3     6    10    15    20    25    31    38    46    55

Постройте график одной из шестеренок в частотной плоскости. Поскольку шестерня соответствует нулевому сдвигу, подтвердите, что частотная характеристика сконцентрирована вдоль горизонтальной или вертикальной оси.

sh = 31;
omegax = -1/2:1/numCols:1/2-1/numCols;
omegay = omegax;
figure
surf(omegax,flip(omegay),psi(:,:,sh),'EdgeColor','none')
view(0,90)
xlabel('\omega_x')
ylabel('\omega_y')
axis equal
axis tight
title(['Zero Shear Shearlet: Scale = ',num2str(scale(sh)),', Cone - ',cone{sh}])

Создайте массив, который содержит только коэффициенты shearlet, которые соответствуют нулям сдвига.

cfsx = zeros(size(cfs));
for k=1:length(ind)
    cfsx(:,:,ind(k)) = cfs(:,:,ind(k));
end

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

rec = isheart2(sls,cfsx);
imagesc(rec)
axis equal
axis tight
colormap gray
title('Reconstruction')