I understand you are trying to use Symmetrization mode for 'swt2' and 'iswt2' function. According to MATLAB R2021a documentation, 'swt2' uses periodic extension:
However, you can symmetrically extend the images manually before performing transform using ‘swt2’ function and then remove the redundant coefficients after taking inverse wavelet transform. Assuming the application of performing the transform and reconstruction on a 2-D image, here's an example:
% Load example image from MATLAB documentation
load woman; % This loads the image 'X' into the workspace
portion=0.25; %Portion of image that can be extended symmetrically
% Symmetrically extend the image manually
X_extended = padarray(X, size(X)*portion, 'symmetric');
% Display the original and extended images
figure;
subplot(1, 2, 1);
imshow(X, []);
title('Original Image');
subplot(1, 2, 2);
imshow(X_extended, []);
title('Symmetrically Extended Image');
% Perform stationary wavelet transform on the extended image
wavelet = 'db1';
N = 2; % Level of decomposition
[A, H, V, D] = swt2(X_extended, N, wavelet);
% Perform inverse stationary wavelet transform
X_reconstructed_extended = iswt2(A, H, V, D, wavelet);
% Remove the redundant coefficients to get the original size
X_reconstructed = X_reconstructed_extended(size(X,1)*portion+1:end-size(X,1)*portion,...
size(X,2)*portion+1:end-size(X,2)*portion);
% Display the reconstructed image
figure;
imshow(X_reconstructed, []);
title('Reconstructed Image after Removing Redundant Coefficients', 'FontSize', 8);
% Sanity check if the reconstructed image is close to the original image
tolerance = 1e-12; % Define a small tolerance
mse = mean((X(:) - X_reconstructed(:)).^2);
if mse<tolerance
disp("Input image and Reconstructed image are same.")
else
disp("Reconstructed image differ.")
end
For more information regarding ‘padarray’ function, Please refer to the following MATLAB R2021a documentation
I hope this helps!
