matlab 数学建模 图像处理

I = imread('ellipse.BMP');
Ig = rgb2gray(I);
imshow(Ig);
hold on
conicP = ginput(15);
plot(conicP(:,1)', conicP(:,2)', 'xr');     %drawing sample points
a0 = [1 1 1 1 1 1];
f = @(a,x)a(1)*x(:,1).^2+a(2)*x(:,2).^2+a(3)*x(:,1).*x(:,2)+a(4)*x(:,1)+a(5)*x(:,2)+a(6);
p = nlinfit(conicP, zeros(size(conicP, 1), 1), f, a0);
syms x y
conic = p(1)*x^2+p(2)*y^2+p(3)*x*y+p(4)*x+p(5)*y+p(6);
c = ezplot(conic,[0,1152],[0,648]);
set(c, 'Color', 'Green');

代码转载

c(2,1)=[];一般左推，不足为0；
imshow(b)就是直接按b的灰度级画出图
imshow（b，[]),如果[]就是这样，而不是[low,high]那么灰度级最小的将被画成黑色，而最大的灰度级将被画成白色，其他的不变。
如果是[low,high],那么灰度级小于low的将被画成黑色，大于high的将被画成白色，其他的不变。

>> imagesc(pc)
>> colorbar
>> imshow(pc)
>> colorbar
>> imshow(pc,[])
>> figure 
>> imshow(pc)