在matlab中绘制重叠

Question

我有一个奇怪的问题，绘图在图表中重叠，但不在同一轴上的图像中。

我确定我没有在某个地方留下保留，否则它也会在图像本身中重叠。

编辑：我想去掉蓝色的重叠线，我只想在那张照片中出现一条蓝色的线。

这是一个示例：

（注意：黑色图像是 RGB 图像，但我没有绘制那个 atm，所以它意味着在图表上从黑色到白色的过渡。）

alt text http://img541.imageshack.us/img541/3212/parabolaaaaa.png

部分代码：

   for K=1:23

    hold on
    I = fig.img.(['p' num2str(K)]);
    bw=(I);
    imshow(bw)

    ss = bwlabel(bw);
    s = regionprops(ss,'centroid');


    centroids{K} = cat(1,s.Centroid);
    hold(imgca,'on')
    plot(imgca,centroids{K}(:,1), centroids{K}(:,2), 'r*'); hold on;
    x=centroids{K}(:,1);
    y=centroids{K}(:,2);
    points=plot(x,y,'go',x,y,'rx');

    hold on
    axis on
    axis fill
    ccentroids = cat(1,centroids{:});
    C1=ccentroids(:,1);
    C2=ccentroids(:,2);
    set(points,'XData',C1,'YData',C2);

    .
    .  
    .

    p= polyfit(x2,y2,2)
    parabola_x = linspace(-250,640,500);
    parabola_polyval = polyval(p,parabola_x);
    plot(parabola_x,parabola_polyval,'b-');
    .
    .  
    .
    end

有什么想法吗？

Answer 1

你有多条蓝线的原因是因为你用这条线为每次通过循环绘制一条：

plot(parabola_x,parabola_polyval,'b-');

事实上，您在循环中一遍又一遍地绘制所有内容（图像、点和线），而没有清除旧的。

您应该在 for 循环的外部初始化绘图对象，并使用SET 命令在循环内更新它们，而不是仅仅重新绘制它们。我在this answer 中给出了一个例子，你之前问的一个问题是我在哪里讨论使用句柄来绘制对象来修改它们。对于您在此处提供的示例代码，您可以执行以下操作：

hImage = imshow(bw(fig.img.p1));  %# Initialize the image
hold on;                          %# Add to the existing plot
hStar = plot(nan,nan,'r*');       %# Initialize the red star
hPoints = plot(nan,nan,'go',...   %# Initialize the other points
               nan,nan,'rx');
hLine = plot(nan,nan,'b-');       %# Initialize the blue line

for K = 1:23

  I = fig.img.(['p' num2str(K)]);
  bw = (I);
  set(hImage,'CData',bw);  %# Update the image

  ss = bwlabel(bw);
  s = regionprops(ss,'centroid');
  centroids{K} = cat(1,s.Centroid);
  set(hStar,'XData',centroids{K}(:,1),...  %# Update the red star
            'YData',centroids{K}(:,2));
  ccentroids = cat(1,centroids{:});
  C1 = ccentroids(:,1);
  C2 = ccentroids(:,2);
  set(hPoints,'XData',C1,'YData',C2);  %# Update the other points

  ...

  p = polyfit(x2,y2,2);
  parabola_x = linspace(-250,640,500);
  parabola_polyval = polyval(p,parabola_x);
  set(hLine,'XData',parabola_x,...      %# Update the blue line
            'YData',parabola_polyval);

  ...

end