MATLAB提取颜色直方图特征

MATLAB提取颜色直方图特征

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

clear
clc
close all
Image = imread（‘1.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怎么给出？
%%=======================================

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     目录           0  2019-04-24 20:10  颜色直方图\
     文件        2806  2018-06-26 09:56  颜色直方图\hsv1.m
     目录           0  2019-04-24 20:10  颜色直方图\图\
     文件      391535  2017-06-15 17:41  颜色直方图\图\一号公路.PNG
     文件      446142  2017-06-15 17:46  颜色直方图\图\上海外滩.PNG
     文件      330449  2017-06-15 18:08  颜色直方图\图\东京塔.PNG
     文件      539426  2017-06-15 17:38  颜色直方图\图\九曲花街.PNG
     文件      407486  2017-06-16 21:50  颜色直方图\图\北京颐和园.PNG
     文件      407011  2017-06-15 18:00  颜色直方图\图\印度尼西亚海神庙.PNG
     文件      440820  2017-06-15 18:14  颜色直方图\图\地狱谷.PNG
     文件      357477  2017-06-15 17:39  颜色直方图\图\大峡谷国家公园.PNG
     文件      487572  2017-06-15 17:38  颜色直方图\图\大都会艺术博物馆.PNG
     文件      446882  2017-06-15 18:10  颜色直方图\图\大阪城公园.PNG
     文件      224343  2017-06-15 17:45  颜色直方图\图\天安门广场.PNG
     文件      561497  2017-06-15 18:07  颜色直方图\图\奈良公园.PNG
     文件      238700  2017-06-15 18:13  颜色直方图\图\富士山‘.PNG
     文件      402059  2017-06-15 17:49  颜色直方图\图\峨眉金顶.PNG
     文件      433491  2017-06-15 17:40  颜色直方图\图\帝国大厦.PNG
     文件      534902  2017-06-15 17:48  颜色直方图\图\成都锦里.PNG
     文件      388199  2017-06-15 17:43  颜色直方图\图\故宫.PNG
     文件      528828  2017-06-15 17:40  颜色直方图\图\斯坦福大学.PNG
     文件      498368  2017-06-15 17:55  颜色直方图\图\新加坡鱼尾狮公园.PNG
     文件      523584  2017-06-15 17:40  颜色直方图\图\时代广场.PNG
     文件      503423  2017-06-15 17:47  颜色直方图\图\杭州西湖.PNG
     文件      531078  2017-06-15 18:04  颜色直方图\图\泰国湄南河.PNG
     文件      480123  2017-06-15 17:53  颜色直方图\图\泰国白庙.PNG
     文件      409416  2017-06-15 17:51  颜色直方图\图\泰国皮皮岛.PNG
     文件      566589  2017-06-15 18:08  颜色直方图\图\浅草寺.PNG
     文件      460221  2017-06-15 18:09  颜色直方图\图\清水寺.PNG
     文件      386072  2017-06-15 17:51  颜色直方图\图\西藏布达拉宫.PNG
     文件      307145  2017-06-15 17:37  颜色直方图\图\金门大桥.PNG
............此处省略4个文件信息