Matlab實(shí)現(xiàn)四種HSV色輪圖繪制的示例代碼

更新時(shí)間：2022年07月06日 09:33:22 作者：slandarer

色輪圖就是色彩相位圖，它完整表現(xiàn)了色相環(huán)360度的全部顏色。本文將利用Matlab語(yǔ)言繪制四種不同的HSV色輪圖，感興趣的可以動(dòng)手嘗試一下

前言

有粉絲問(wèn)我圖(d)上的色盤圖咋畫：

明度(V)漸變版

注意，此處是為了還原論文中圖片所以X加了。

%?生成網(wǎng)格
tList=linspace(0,2.*pi,300);
rList=linspace(0,1,100);
[theta,R]=meshgrid(tList,rList);

%?角度及半徑轉(zhuǎn)換為坐標(biāo)
X=cos(theta+pi).*R;
Y=sin(theta).*R;
Z=zeros(size(X));

%?構(gòu)造hsv網(wǎng)格并轉(zhuǎn)換為rgb網(wǎng)格
hsvMesh=cat(3,theta./2./pi,ones(size(R)),R);
rgbMesh=hsv2rgb(hsvMesh);

%?surf繪圖
surf(X,Y,Z,'EdgeColor','none','CData',rgbMesh)
axis?equal;axis([-1,1,-1,1]);
view(0,90);

飽和度(S)漸變版

%?生成網(wǎng)格
tList=linspace(0,2.*pi,300);
rList=linspace(0,1,100);
[theta,R]=meshgrid(tList,rList);

%?角度及半徑轉(zhuǎn)換為坐標(biāo)
X=cos(theta).*R;
Y=sin(theta).*R;
Z=zeros(size(X));

%?構(gòu)造hsv網(wǎng)格并轉(zhuǎn)換為rgb網(wǎng)格
hsvMesh=cat(3,theta./2./pi,R,ones(size(R)));
rgbMesh=hsv2rgb(hsvMesh);

%?surf繪圖
surf(X,Y,Z,'EdgeColor','none','CData',rgbMesh)
axis?equal;axis([-1,1,-1,1]);
view(0,90);

明度飽和度不變花瓣版

%?構(gòu)造花瓣形狀
t=(0:.1:(2.01/3*pi))+pi;
CX=cos(t)+1;
CY=sin(t);
[NCX,NCY]=rotateData(CX,CY,pi/6);
NCX=[NCX,NCX(end:-1:1)];
NCY=[NCY,-NCY(end:-1:1)];

%?循環(huán)繪圖
hold?on;
tt=linspace(0,2*pi,33);tt(end)=[];
for?i=tt
????[tNCX,tNCY]=rotateData(NCX,NCY,i);
????fill(tNCX,tNCY,hsv2rgb([i/2/pi,1,1]),'FaceAlpha',0.2,...
????????'LineWidth',1.5,'EdgeColor',[1,1,1],'EdgeAlpha',0.5)
end
axis?equal;axis([-1.8,1.8,-1.8,1.8]);

%?數(shù)據(jù)旋轉(zhuǎn)角度
function?[X,Y]=rotateData(X,Y,theta)
????rotateMat=[cos(theta),-sin(theta);sin(theta),cos(theta)];
????XY=rotateMat*[X;Y];
????X=XY(1,:);Y=XY(2,:);
end

HSV全空間

%?線性劃分
h=linspace(0,2*pi,31);h=h(1:26);
s=linspace(0,1,9);
v=linspace(0,1,11);

%?繪制上面
hold?on;
[H,S]=meshgrid(h,s);
surf(S.*cos(H),S.*sin(H),H-H+1,hsv2rgb(cat(3,H/2/pi,S,H-H+1)));
%?繪制側(cè)面
[H,V]=meshgrid(h,v);
surf(V.*cos(H),V.*sin(H),V,hsv2rgb(cat(3,H/2/pi,H-H+1,V)));
%?繪制兩個(gè)截面
[S,V]=meshgrid(s,v);
surf(S.*V.*cos(h(1)),S.*V.*sin(h(1)),V,hsv2rgb(cat(3,S-S+h(1)/2/pi,S,V)));
surf(S.*V.*cos(h(end)),S.*V.*sin(h(end)),V,hsv2rgb(cat(3,S-S+h(end)/2/pi,S,V)));

shading?flat;view(60,45);axis?off;
ax=gca;ax.Position=[-1/6,-2/6,1+2/6,1+2/6];