如何提出圆环内边界？

fsnow

我有一组数据需要进行环拟合，现在已完成环的外圆拟合，内圆边界点提取一直不理想（如图），麻烦大家指点一下。
图中空白部分数据设定为NaN。

我的程序：
[filename, pathname] = uigetfile('*.txt','Open Data');
filepath = strcat(pathname,filename);
if length(filepath) == 0
    return;
end
Data = textread(filepath,'','delimiter',';','emptyvalue',NaN);
pcolor(Data), axis equal,axis tight, shading interp;
colorbar; hold on;
% 找出圆环的外边界，并进行最小二乘圆拟合
BoundIdx = []; BoundIdy = [];
for k1=1:size(Data,1)
    NaNIdx=find(isnan(Data(k1,:)));
    JumpSize=diff(NaNIdx);
    BI=find(JumpSize>=2);
    BINum = length(BI);
    if ~(isempty(BI))
         BoundIdx = [BoundIdx NaNIdx(BI(1)) NaNIdx(BI(BINum)+1)];  
         BoundIdy = [BoundIdy k1 k1];
    end
end
[xc,yc,R] = circfit(BoundIdx,BoundIdy);
t = linspace(0,pi*2,100);
z = xc+yc*i+R*exp(i*t);
%plot(BoundIdx,BoundIdy,'kx');
plot(z); hold on;
R = round(R); xc = round(xc); yc = round(yc);
% 找出圆环的内边界，并进行最小二乘圆拟合
BoundIdx = []; BoundIdy = [];
for k1=yc-R:yc+R
    if ~(isnan(Data(k1,xc))&isnan(Data(2*yc-k1,xc)))   
        continue;
    end
    ValueIdx=find(~isnan(Data(k1,xc:xc+R)));
    if ~(isempty(ValueIdx))
         BoundIdx = [BoundIdx ValueIdx(1)-1+xc];  
         BoundIdy = [BoundIdy k1];
    end
    ValueIdx=find(~isnan(Data(k1,xc:-1:xc-R)));
    if ~(isempty(ValueIdx))
         BoundIdx = [BoundIdx -ValueIdx(1)+1+xc];  
         BoundIdy = [BoundIdy k1];
    end
end
[xc2,yc2,R2] = circfit(BoundIdx,BoundIdy);
t = linspace(0,pi*2,100);
z = xc2+yc2*i+R2*exp(i*t);
plot(BoundIdx,BoundIdy,'bx');
plot(z); hold off;

花如月

可以考虑使用最临近插值，然后使用形态学或者边缘检测的方法求出内边界。

另外论坛有一个圆的拟合程序，数据点数足够了，你可以找找看

fsnow

圆拟合程序已经有了。现在想尽量找出环的内圆所有边界点
使用最临近插值不会破坏原有的有效数据吗？

花如月

我没表述准确，这个问题不适合用拟合或者插值
用邻域元素的一些统计特征来填补缺失的数据似乎更科学

[ 本帖最后由 eight 于 2007-11-26 22:49 编辑 ]

心灯

把图形从中间沿最大最小值连线分成两组数据，然后用最优化理论中的Pareto解判断的思想分类就可以把内边界给找出来了。具体内容自己找书看吧。

另外，你外边界都找出来了，内边界不也一样么？只是相对位置关系变了。

[ 本帖最后由 心灯 于 2007-11-26 23:28 编辑 ]