hsv 源码

颜色直方图的matlab源码，备注非常详细

%%================================

clear
clc
close all
Image = imread（‘image068.jpg‘）;
[MNO] = size（Image）;
[hsv] = 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（ij）<=15||h（ij）>345
H（ij） = 0;
end
if h（ij）<=25&&h（ij）>15
H（ij） = 1;
end
if h（ij）<=45&&h（ij）>25
H（ij） = 2;
end
if h（ij）<=55&&h（ij）>45
H（ij） = 3;
end
if h（ij）<=80&&h（ij）>55
H（ij） = 4;
end
if h（ij）<=108&&h（ij）>80
H（ij） = 5;
end
if h（ij）<=140&&h（ij）>108
H（ij） = 6;
end
if h（ij）<=165&&h（ij）>140
H（ij） = 7;
end
if h（ij）<=190&&h（ij）>165
H（ij） = 8;
end
if h（ij）<=220&&h（ij）>190
H（ij） = 9;
end
if h（ij）<=255&&h（ij）>220
H（ij） = 10;
end
if h（ij）<=275&&h（ij）>255
H（ij） = 11;
end
if h（ij）<=290&&h（ij）>275
H（ij） = 12;
end
if h（ij）<=316&&h（ij）>290
H（ij） = 13;
end
if h（ij）<=330&&h（ij）>316
H（ij） = 14;
end
if h（ij）<=345&&h（ij）>330
H（ij） = 15;
end
end
end
for i = 1:M
for j = 1:N
if s（ij）<=0.15&&s（ij）>0
S（ij） = 1;
end
if s（ij）<=0.4&&s（ij）>0.15
S（ij） = 2;
end
if s（ij）<=0.75&&s（ij）>0.4
S（ij） = 3;
end
if s（ij）<=1&&s（ij）>0.75
S（ij） = 4;
end
end
end
for i = 1:M
for j = 1:N
if v（ij）<=0.15&&v（ij）>0
V（ij） = 1;
end
if v（ij）<=0.4&&v（ij）>0.15
V（ij） = 2;
end
if v（ij）<=0.75&&v（ij）>0.4
V（ij） = 3;
end
if v（ij）<=1&&v（ij）>0.75
V（ij） = 4;
end
end
end
L=16*H+4*S+V;


% 构建4*16二维数组存放H-S数据
Hist = zeros（164）;
for i = 1:M
for j = 1:N
for k = 1:16
for l = 1:4
if l==S（ij）&& k==H（ij）+1
Hist（kl） = Hist（kl）+1;
end
end
end
end
end
for k = 1:16
for l =1:4
His（（k-1）*4+l） = Hist（kl）;%转化为一维数组
end
end
His = His/sum（His）*1000;
size（His）
% 手工绘制彩色图像直方图
% hist_h
m=0;
for j = 1:300
if rem（j16）==1 && m<16
for k = 0:15
for i = 1:200
hist_h（ij+k） = m;
end 
end
m = m+1;
end
end
% hist_s
m=0;
for j = 1:300
if rem（j4） == 1 && m<64
n = rem（m4）;
for k = 0:3 
for i =1:200 
hist_s（ij+k） = n+1; 
end 
end
m = m+1; 
end 
end
% hist_v
for j = 1:256
for i = 1:200
hist_v（ij） = 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（ij） = 0;
end
end
end
end
end

%将h、s、v分量图合并转化为RGB模式

I_H = hsv2rgb（hist_h/16hist_s/4hist_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怎么给出？
%%=======================================

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        2813  2012-11-20 09:41  hsv1.m