Spline smoothing in images

10 vues (au cours des 30 derniers jours)
sudesh le 8 Fév 2013
Commenté : Walter Roberson le 13 Juin 2022
Can anyone guide me how we use splines for interpolation in images? I mean how we relate value of x,y of spline function in an image that has its spatial coordinates and intensity values.
1 commentaireAfficher -1 commentaires plus anciensMasquer -1 commentaires plus anciens
Mohammed Kahya le 31 Août 2016

Connectez-vous pour commenter.

Réponses (3)

Matt J le 8 Fév 2013
2 commentairesAfficher AucuneMasquer Aucune
Matt J le 8 Fév 2013
You would have to change 1 line to convert Example1D.m from an interpolating spline to a smoothing spline
sFit = BasisFine*(BasisFine\sFine);
and similarly for Example2D.m
sudesh le 8 Fév 2013
how can we smooth a local window(like 3X3) of an image using this spline? As I got till now that x and y parameters of spline function are row/column of image/window_of_image.

Connectez-vous pour commenter.

Image Analyst le 8 Fév 2013
A 2D Savitzky Golay filter is similar, and I already have code for that. Basically it replaces each pixel with a polynomial fit over a sliding window. It gives the effect of smoothing. Let me know if you have the Signal Processing Toolbox and would like my demo.
2 commentairesAfficher AucuneMasquer Aucune
sudesh le 8 Fév 2013
I would like to see the demo. I have the signal processing toolbox but can't it be shown just using image processing toolbox. Its like I am purely from computer science background and it takes some time to me to understand the signal part. But I will try to grasp the point.
Image Analyst le 8 Fév 2013
The sgolay function is not contained in the Image Processing Toolbox. What is does is fit a 25 element 1D vector of gray levels to a polynomial (I used 1 but you can use order 3 if you want) and then it replaced the center value with the value from the fitted curve. It does this column by column (sliding the 25 element window down the rows in each column), and then it does it row by row (blurring across columns) to get a 2D blur.
Try this:
% Filter using Savitzky-Golay filtering.
clc; % Clear command window.
clear; % Delete all variables.
close all; % Close all figure windows except those created by imtool.
imtool close all; % Close all figure windows created by imtool.
workspace; % Make sure the workspace panel is showing.
fontSize = 14;
% Check that user has the Image Processing Toolbox installed.
if ~hasSPT
% User does not have the toolbox installed.
message = sprintf('Sorry, but you do not seem to have the Signal Processing Toolbox.\nDo you want to try to continue anyway?');
reply = questdlg(message, 'Toolbox missing', 'Yes', 'No', 'Yes');
% User said No, so exit.
return;
end
end
% Read in standard MATLAB gray scale demo image.
imageArray = double(imageArray);
[rows columns numberOfColorBands] = size(imageArray);
subplot(2, 2, 1);
imshow(imageArray, [0 255]);
title('Original Grayscale Image', 'FontSize', fontSize);
set(gcf, 'Position', get(0,'Screensize')); % Enlarge figure to full screen.
set(gcf,'name','Savitzky-Golay Filter Image Analysis Demo','numbertitle','off')
% Apply the Savitzky-Golay filter.
% First apply it in the vertical (row) direction.
k = 1; % Order of the polynomial
windowSize = 25;
verticallySmoothedImage = sgolayfilt(imageArray, k, windowSize, [], 1);
subplot(2, 2, 2);
imshow(verticallySmoothedImage, [0 255]);
title('Savitzky-Golay filtered in the vertical direction only', 'FontSize', fontSize);
% Apply the Savitzky-Golay filter.
% First apply it in the vertical (row) direction.
k = 1; % Order of the polynomial
windowSize = 25;
horizontallySmoothedImage = sgolayfilt(imageArray, k, windowSize, [], 2);
subplot(2, 2, 3);
imshow(horizontallySmoothedImage, [0 255]);
title('Savitzky-Golay filtered in the horizontal direction only', 'FontSize', fontSize);
doublySmoothedImage = sgolayfilt(verticallySmoothedImage, k, windowSize, [], 2);
subplot(2, 2, 4);
imshow(doublySmoothedImage, [0 255]);
title('Savitzky-Golay filtered in both directions', 'FontSize', fontSize);

Connectez-vous pour commenter.

Mehri Mehrnia le 7 Mai 2022
I have the same question. In all spline questins, we map coordinates of points and then interpolate.
But, for imaging, how I can map the matrix arrays and interpolate????
4 commentairesAfficher 2 commentaires plus anciensMasquer 2 commentaires plus anciens
Mehri Mehrnia le 13 Juin 2022
You mentioned valuable points but I didn't get all of them. what is "Image Analyst posted code for"?
Walter Roberson le 13 Juin 2022

Connectez-vous pour commenter.

Catégories

En savoir plus sur Geometric Transformation and Image Registration dans Help Center et File Exchange

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by