How does imflatfield work?
14 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
What formula does the imflatfield function use when used with a matrix as input?
For example, I have an image of surface temperature (attached) = bad_img. If I apply
corr_img = imflatfield(bad_img,100);
How can I reproduce corr_img without using imflatfield? I know that the gaussian filter can be generated using
iblur=imgaussfilt(bad_img, 100);
But what are the next steps?
1 commentaire
Réponses (1)
DGM
le 2 Oct 2021
Modifié(e) : DGM
le 2 Oct 2021
William's comment is close. Probably close enough to be useful, but not quite exactly the same as what imflatfield() does.
A = imread('printedtext.png');
sigma = 30;
% this is the reference
output1 = imflatfield(A,sigma);
% this is a replication of the behavior
A = im2single(A); % non-double images are converted to single
% this is the default filter size, but it can be set explicitly too
blurred = imgaussfilt(A,sigma,'padding','symmetric','filtersize',2*ceil(2*sigma)+1);
output2 = A*mean(A(:),'omitnan')./blurred; % mean of A, not of blurred image
output2 = im2uint8(output2);
immse(output1,output2) % show that the results match
Bear in mind that if you want to exactly replicate the behavior, there are other things that happen (e.g. dealing with NaNs and Inf) that need to also be considered. You can always just open imflatfield() and read it.
EDIT: oof. This wouldn't be the first time I absentmindedly answered a dead question...
1 commentaire
Arthur Redgate
le 1 Sep 2023
I appreciate your answer as it's helped me 5 years after the asking of the original question!
Voir également
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!