这不是一个真正的答案。我要看着它,但是我发现它是容易重写它为了算出它在做什么。这是它在的地方:
%自相关
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - %
B = abs (fftshift (ifft2 (fft2 (A)。*连词(fft2 (A)))))。/元素个数(一个);
imshow (imcomplement (B))
% %
%阈值矩阵,找到重心
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - %
checkstart = max (B (:));
checkstart =圆(checkstart, 2);
为检查= checkstart: -0.01:0
C = B >检查;
CC = bwconncomp (C);%这应该会更快
如果CC.NumObjects > = 4 & & CC.NumObjects < = 10
打破
结束
结束
s = regionprops (CC、“重心”);%这样做外
imshow (imcomplement (C))
% %
%从中间找到距离和角形心
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - %
%找到实际geometrically-central blob
imcenter = fliplr(大小(C, 1:2)) / 2;
distfromcenter =√总和(vertcat (s.Centroid) -imcenter)。^ 2, 2));
[~,refidx] = min (distfromcenter);%的索引参考峰
%找到所有其他斑点
nonrefidx = 1:元素个数(年代);
nonrefidx (refidx) = [];
%得到距离和角度关于参考
δ= s (refidx)。重心——vertcat (s (nonrefidx) .Centroid);
distfromref =√和三角洲。^ 2,2))
anglefromref = atand(δ(:,2)。/δ(:1));
%减少距离,得到角有关
[mindistance, idx] = min (distfromref);
minangle = anglefromref (idx);
% 45包角(-45)
%这一样东西dif和+ -45测试
anglechange = 45 - mod (45 - minangle, 90);
%旋转图像
D = imrotate (anglechange);
clc
mindistance
anglechange
imshow (imcomplement (D))
% %
% detile图像
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - %
%得到瓷砖
深圳=大小(D, 1:2);
mindistance =圆(mindistance);
瓷砖地板= (sz. / mindistance);% (x, y)
%使用mat2cell ()
M = mindistance *的(瓷砖(1),1);%瓷砖尺寸向量
N = mindistance *的(瓷砖(2),1);
Dt = D (1: mindistance *瓷砖(1),1:mindistance *花砖(2));%的作物
瓷砖= mat2cell (Dt, M, N);% detile
%显示瓷砖
蒙太奇(瓷砖,“边界”10“写成backgroundcolor”,“米”)
提供的图片都是相似的,他们没有一个完整的重复模式的至少在一个方向。虽然似乎有足够的信息来辨别正确的间距,给定的方法为notile找到mindistance给了错误的结果。垫和一个。垫(这些副本),它似乎sqrt(2)倍时失败。
好吧。垫文件似乎detile好,尽管这些模式显然不精确的复制。
我必须想想更多的notile.mat看能做些什么。
编辑:
我取代了第一部分:
%自相关
% - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - %
B = normxcorr2(一个);%这似乎是更可靠
B = imtophat (B, (21));%顶环过滤器
B = mat2gray (B);%正常化
…这似乎为notile工作。垫了。detiling常规仍然不做任何有用的在这种情况下,由于图像中只有一个完整的瓷砖,但起码它可以正确的距离。
