本帖最后由 2040849623 于 2014-5-14 10:52 编辑

下面是我的程序，data数据过多，我只取了一部分，运行的结果只有色品图的轮廓线，颜色填充不上去

程序提示错误为

??? Error using ==> TriScatteredInterp

The number of data point locations should equal the number of data point values.

Error in ==> zuhe at 504

FR=TriScatteredInterp(x,y,R,'natural');

论坛内的高手能帮忙看一下吗？不胜感激！我的代码：

data=[360 0.0001299 3.92E-06 0.0006061

361 0.0001458 4.39E-06 0.0006809

362 0.0001638 4.93E-06 0.0007651

363 0.000184 5.53E-06 0.00086

364 0.0002067 6.21E-06 0.0009666

365 0.0002321 6.97E-06 0.001086

366 0.0002607 7.81E-06 0.0012206

367 0.0002931 8.77E-06 0.0013727

368 0.0003294 9.84E-06 0.0015436

369 0.0003699 1.10E-05 0.0017343

370 0.0004149 1.24E-05 0.001946

371 0.0004642 1.39E-05 0.0021778

372 0.000519 1.56E-05 0.0024358

373 0.0005819 1.74E-05 0.002732

374 0.0006552 1.96E-05 0.0030781

375 0.0007416 2.20E-05 0.003486

376 0.000845 2.48E-05 0.0039752

377 0.0009645 2.80E-05 0.0045409

378 0.0010949 3.15E-05 0.0051583

379 0.0012312 3.52E-05 0.0058029

380 0.001368 0.000039 0.00645

381 0.0015021 4.28E-05 0.0070832

382 0.0016423 4.69E-05 0.0077455

383 0.0018024 5.16E-05 0.0085012

384 0.0019958 5.72E-05 0.0094145

385 0.002236 0.000064 0.01055

386 0.0025354 7.23E-05 0.0119658

387 0.0028926 8.22E-05 0.0136559

388 0.0033008 9.35E-05 0.0155881

389 0.0037532 0.0001061 0.0177302

390 0.004243 0.00012 0.02005

391 0.0047624 0.000135 0.0225114

392 0.00533 0.0001515 0.0252029

393 0.0059787 0.0001702 0.0282797

394 0.0067411 0.0001918 0.031897

395 0.00765 0.000217 0.03621】

sumD=data(:,2)+data(:,3)+data(:,4);

x=data(:,2)./sumD;

y=data(:,3)./sumD;

X=data(:,2);%X刺激值

Y=data(:,3);%Y刺激值

Z=data(:,4);%Z刺激值

N=[X,Y,Z];

M=[0.49,0.31,0.2;

0.17697,0.81240,0.01063;

0.00,0.01,0.99];

RGB=0.17697*M^-1*N';

R=RGB(:,1);

G=RGB(:,2);

B=RGB(:,3);

hold on;

grid

plot(x,y);

%获取数组长度

[len,a]=size(x);

%把图形封闭起来

plot([x(1),x(len)],[y(1),y(len)]);

x1=[x;x(1)];

y1=[y;y(1)];

tic

[XX,YY]=meshgrid(0:.01:1,0:.01:1);

in=inpolygon(XX,YY,x1,y1);

index=find(in==1);

FR=TriScatteredInterp(x,y,R,'natural');

FG=TriScatteredInterp(x,y,G,'natural');

FB=TriScatteredInterp(x,y,B,'natural');

for i =1:length(index)

r=FR(XX(index(i)),YY(index(i)));

g=FG(XX(index(i)),YY(index(i)));

b=FB(XX(index(i)),YY(index(i)));

if(r>0.95&&g>0.95&&b>0.95)  %白化处理,防止因插值带来的误差

r=1;g=1;b=1;

end

if(r>1)

r=1;

elseif(r<0)

r=0;

end

if(g>1)

g=1;

elseif(g<0)

g=0;

end

if(b>1)

b=1;

elseif(b<0)

b=0;

end                          %防止rgb的溢出，正确的做法

%M=mapminmax([r g b],0,1);  %归一化处理这样的归一化处理会损失掉rgb颜色的比例

%这就是为什么没有出现白点的原因，还有一个就是插值求得的rgb有的为负

%但是插值求得的值很多接近于0，必须对其进行归一化处理，否则，整体就变成了黑色

%fill(XX(index(i)),YY(index(i)),M); %给一个小的区域填充，会比一个一个地描绘点更快

plot(XX(index(i)),YY(index(i)), 'MarkerSize',10,'Color',[r,g,b]);

end

toc

select=1;   %控制标记

xre=0.333;yre=0.333; %参考白点xy坐标

plot(xre,yre,'ko');

while(select==1)

zoom on

[x0,y0]=ginput(1);

zoom off

if(inpolygon(x0,y0,x1,y1))

if(xre==x0)

x3=xre;

y3=max(y)+0.1;

elseif(yre==y0)

x3=max(x)+0.1;y3=yre;

else

k=(yre-y0)/(xre-x0);%计算参考点和选取点的斜率

end

if(x0>xre)

x3=max(x)+0.1;

y3=yre-k*(xre-x3);

if(x0

x3=-0.1;

y3=yre-k*(xre-x3);

end

end

select=0;

else

errordlg('拾取点在色品图外部，无法计算色度和饱和度.请重新选取','Waring!','on');

pause(1)

end

end

xline=[xre x0 x3];

yline=[yre y0 y3];

[xcross,ycross]=polyxpoly(xline,yline,x1,y1,'unique');

xxline=[xre x0 xcross];

yyline=[yre y0 ycross];

plot(xxline,yyline,'r');

plot(x0,y0,'ko');

plot(xcross,ycross,'ko');

xlabel('x')

ylabel('y')

title('CIE-xyz')

text(xre,yre,' referene point')

text(x0,y0,' input point')

text(xcross,ycross,' cross point')

%色调H，饱和度S的计算

if( (ycross-y(1))*(x(end)-xcross) == (y(end)-ycross)*(xcross-x(1)) )%交点在纯紫线上

H=780;

else

wavelength=data(:,1);

H=griddata(x,y,wavelength,xcross,ycross,'v4');

end

D=sqrt((xre-xcross)^2+(yre-ycross)^2);%参考白点到交点的距离

d=sqrt((xre-x0)^2+(yre-y0)^2);        %参考白点到拾取点的距离

S=d/D;

2014-5-14 10:46 上传

我的结果

2014-5-14 10:49 上传