使用活动轮廓时出错

Question

使用活动轮廓时出现错误：

？？？输入的未定义函数或方法“activecontour” 'single' 类型的参数。

==> Segmentasi 中的错误>threshold_Callback 在 114 final = activecontour(image2,mask, 100);

image1 = handles.citra1;
    level=0.008;
    bw = edge(image1,'Canny');
    axes(handles.axes2);
    imshow(bw,[]);

    %active contour
    image2 = bwdist(~bw);
    mask = zeros(size(image2));
    mask(25:end-25,25:end-25) = 1;
    final = activecontour(image2,mask , 100);
    axes(handles.axes5)
    imshow(final,[]);
    handles.data3 = final;
    guidata(hObject,handles);

我正在使用 dicom 文件中的图像（矢状图像）

Answer 1

问题是bwdist函数的输出是一个矩阵，activecontour函数的输入是一个灰度图像。所以在使用之前需要将矩阵转换为灰度图像。这是使用名为mat2gray 的函数完成的。为此，在使用bwdist 后应用mat2gray 如下：

image2 = mat2gray(bwdist(~bw));

然后剩下的代码就可以工作了。检查我的简单示例：

bw = zeros(200,200);
bw(50,50) = 1;
bw(50,150) = 1;
bw(150,50) = 1;
D1 = bwdist(bw);

D2 = mat2gray(D1);
mask = zeros(size(D2));
mask(25:end-25,25:end-25) = 1;
final = activecontour(D2,mask,100);

subplot(1,3,1);imshow(bw)
subplot(1,3,2);imshow(mat2gray(D1))
subplot(1,3,3);imshow(final,[]);