I'm using Fiji (ImageJ variant) to see how this script aligns the images. Stackreg isn't good because it moves the image up and down, I only want left and right alignment.
How to align two images together automatically: Round 2
3 次查看(过去 30 天)
显示 更早的评论
Hello!
I have my images closely aligned, but not close enough. Is there anyway I can improve object selection in an image or getting the centroid more accurately/ consistently? I have been using this script to locate the centroid of both images and Sean De's imtranslate to modify the image.
Thanks
-Frank
Image

Code
%First Image
i = imread('343TESTA.tif');
j = size(i);
imtool(i);
for k = 1:j(1)
for l = 1:j(2)
if i(k,l) > 160 & i(k,l) > 160
d(k,l) = 1;
else
d(k,l) = 0;
end
end
end
imtool(d);
labela = bwlabel(d);
imagesc(labela);
[labela,num] = bwlabel(d,4);
stats = regionprops(labela,'basic');
stats(1).Area
stats(1).Centroid
max_area = max([stats.Area]);
biggrain = find([stats.Area]== max_area);
p = stats(biggrain).Centroid
%Second Image
a = imread('343TESTB.tif');
b = size(a);
imtool(a);
for c = 1:b(1)
for l = 1:b(2)
if i(c,l) > 160 & a(c,l) > 160
e(c,l) = 1;
else
e(c,l) = 0;
end
end
end
imtool(e);
labelb = bwlabel(e);
imagesc(labelb);
[labelb,num] = bwlabel(e,4);
stats = regionprops(labelb,'basic');
stats(1).Area
stats(1).Centroid
max_area = max([stats.Area]);
biggrain = find([stats.Area]== max_area);
o = stats(biggrain).Centroid
%Calculates the differences in the centroid
displacement = p -o
%IM Translate
v = imread('343TESTA.tif');
u = imread('343TESTB.tif');
u2 = imtranslate(u, [0 -1.0376]); %Note: [y x] for imtranslate
subplot(1,2,1), imshow(v)
subplot(1,2,2), imshow(u2)
imwrite(u2, '343TESTC','tif');
回答(1 个)
Sean de Wolski
2011-5-18
You could try this, though I don't know how well it will work on binary images...
Also, how were these images acquired? Is there a chance that there is a change in thickness/density etc.?
2 个评论
Sean de Wolski
2011-5-18
Maybe. Perhaps try the subpixel phase correlation on the original images, not the images converted to binary.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Image Processing Toolbox 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!