%%================================
clear
clc
close all
Image = imread('1.jpg');
[M,N,O] = size(Image);
[h,s,v] = rgb2hsv(Image);
H = h; S = s; V = v;
h = h*360;
%将hsv空间非等间隔量化：
% h量化成16级；
% s量化成4级；
% v量化成4级；
for i = 1:M
for j = 1:N
if h(i,j)<=15||h(i,j)>345
H(i,j) = 0;
end
if h(i,j)<=25&&h(i,j)>15
H(i,j) = 1;
end
if h(i,j)<=45&&h(i,j)>25
H(i,j) = 2;
end
if h(i,j)<=55&&h(i,j)>45
H(i,j) = 3;
end
if h(i,j)<=80&&h(i,j)>55
H(i,j) = 4;
end
if h(i,j)<=108&&h(i,j)>80
H(i,j) = 5;
end
if h(i,j)<=140&&h(i,j)>108
H(i,j) = 6;
end
if h(i,j)<=165&&h(i,j)>140
H(i,j) = 7;
end
if h(i,j)<=190&&h(i,j)>165
H(i,j) = 8;
end
if h(i,j)<=220&&h(i,j)>190
H(i,j) = 9;
end
if h(i,j)<=255&&h(i,j)>220
H(i,j) = 10;
end
if h(i,j)<=275&&h(i,j)>255
H(i,j) = 11;
end
if h(i,j)<=290&&h(i,j)>275
H(i,j) = 12;
end
if h(i,j)<=316&&h(i,j)>290
H(i,j) = 13;
end
if h(i,j)<=330&&h(i,j)>316
H(i,j) = 14;
end
if h(i,j)<=345&&h(i,j)>330
H(i,j) = 15;
end
end
end
for i = 1:M
for j = 1:N
if s(i,j)<=0.15&&s(i,j)>0
S(i,j) = 1;
end
if s(i,j)<=0.4&&s(i,j)>0.15
S(i,j) = 2;
end
if s(i,j)<=0.75&&s(i,j)>0.4
S(i,j) = 3;
end
if s(i,j)<=1&&s(i,j)>0.75
S(i,j) = 4;
end
end
end
for i = 1:M
for j = 1:N
if v(i,j)<=0.15&&v(i,j)>0
V(i,j) = 1;
end
if v(i,j)<=0.4&&v(i,j)>0.15
V(i,j) = 2;
end
if v(i,j)<=0.75&&v(i,j)>0.4
V(i,j) = 3;
end
if v(i,j)<=1&&v(i,j)>0.75
V(i,j) = 4;
end
end
end
L=16*H+4*S+V;
% 构建4*16二维数组存放H-S数据
Hist = zeros(16,4);
for i = 1:M
for j = 1:N
for k = 1:16
for l = 1:4
if l==S(i,j)&& k==H(i,j)+1
Hist(k,l) = Hist(k,l)+1;
end
end
end
end
end
for k = 1:16
for l =1:4
His((k-1)*4+l) = Hist(k,l);%转化为一维数组
end
end
His = His/sum(His)*1000;
size(His)
% 手工绘制彩色图像直方图
% hist_h
m=0;
for j = 1:300
if rem(j,16)==1 && m<16
for k = 0:15
for i = 1:200
hist_h(i,j+k) = m;
end
end
m = m+1;
end
end
% hist_s
m=0;
for j = 1:300
if rem(j,4) == 1 && m<64
n = rem(m,4);
for k = 0:3
for i =1:200
hist_s(i,j+k) = n+1;
end
end
m = m+1;
end
end
% hist_v
for j = 1:256
for i = 1:200
hist_v(i,j) = 0.98;
end
end
% 把His赋值给hist_v
for k = 1:64
for j = 1:256
if floor((j-1)/4) == k
for i = 1:200
if i<200-His(k+1)%i>His(k+1)%
hist_v(i,j) = 0;
end
end
end
end
end
%将h、s、v分量图合并转化为RGB模式
I_H = hsv2rgb(hist_h/16,hist_s/4,hist_v);
% 画图显示
figure(1);
imshow(Image),title('原图');axis on;
figure(2),imshow(H,[]),title('H分量图');axis on;
figure(3),imshow(S,[]),title('S分量图');axis on;
figure(4),imshow(V,[]),title('V分量图');axis on;
figure(5),imshow(I_H,[]),title('H-S直方图');axis on;
%figure(6),imshow(I_rgb,[]),title('色彩量化后的RGB图像');axis on %I_rgb怎么给出？
%%=======================================
matlab颜色直方图.zip （33个子文件） 
颜色直方图 图 
清水寺.PNG 449KB 杭州西湖.PNG 492KB 故宫.PNG 379KB 天安门广场.PNG 219KB 银座.PNG 449KB 泰国湄南河.PNG 519KB 马来西亚丹绒亚路海滩.PNG 228KB 东京塔.PNG 323KB 金门大桥.PNG 300KB 大峡谷国家公园.PNG 349KB 大都会艺术博物馆.PNG 476KB 新加坡鱼尾狮公园.PNG 487KB 泰国白庙.PNG 469KB 成都锦里.PNG 522KB 西藏布达拉宫.PNG 377KB 地狱谷.PNG 430KB 峨眉金顶.PNG 393KB 浅草寺.PNG 553KB 泰国皮皮岛.PNG 400KB 上海外滩.PNG 436KB 奈良公园.PNG 548KB 富士山‘.PNG 233KB 马来西亚美人鱼岛.PNG 500KB 印度尼西亚海神庙.PNG 397KB 时代广场.PNG 511KB 帝国大厦.PNG 423KB 北京颐和园.PNG 398KB 九曲花街.PNG 527KB 大阪城公园.PNG 436KB 一号公路.PNG 382KB 斯坦福大学.PNG 516KB 马来西亚水上清真寺.PNG 518KB
hsv1.m 3KB
评论0
最新资源