三维空间离散点的平面方程拟合

clear clc w=200 e = imread( strcat(int2str(10),'.bmp' )); %分析不同图片将修改此处即可 for i = 0:w-1 for j=0:w-1 z1(i+1,j+1)=e( (i*w+j)+1 +w*w*2);% 1 R, 2 G, 3 B end end I1=im2double(z1'); figure(1); mesh(I1); axis([0 255 0 255 0 1]) %mesh(I(:,:,2)); for i = 0:w-1 for j=0:w-1 x(i*w+j+1)=double(j); y(i*w+j+1)=double(i); z(i*w+j+1)= double( z1(i+1,j+1) );%绿色像素 end end [ a b c d] = zuiXiaoErCheng(x,y,z,w*w); [a0 a1 a2] = sanYuanFangCheng(a,b,c,d) % z = a0*x+a1*y+a2 for i = 0:w-1 for j=0:w-1 z2(i+1,j+1)=int16(a0*j + a1*i+a2); end end figure(2); I2=im2double(z2); mesh(double(double(z2)./255)); axis([0 255 0 255 0 1]) err1 = std2(z1)