%mean?shift?tracking
function?rgb_meanshift（）
clear;
rgb=imread（‘0329.jpg‘）;%讀取第一幀圖片
figure（1）imshow（rgb）;
[temprect]=imcrop（rgb）;%鼠標截取需要跟蹤的目標，rect=[xywh]分別是左上角坐標和圖片長寬
[abc]=size（temp）;
%目標框中心點坐標
tic_x=rect（1）+rect（3）/2;
tic_y=rect（2）+rect（4）/2;
m_wei=zeros（ab）;%權值矩陣初始化
%定義h為搜索窗口的大小，也可稱為核函數帶寬，包括長度和寬度，分別為目標長度和寬度的一半，h（xy）。
%這里的h是一個固定值，即跟蹤過程中搜索窗口的大小是不變的，這導致?Mean-Shift跟蹤算法在跟蹤過程中不能適應目標尺寸的變化，對大小有改變的目標產生定位不準確。
y（1）=a/2;
y（2）=b/2;
h=y（1）^2+y（2）^2;%帶寬模的平方
%計算權值矩陣
for?i=1:a
????for?j=1:b
????????dist=（i-y（1））^2+（j-y（2））^2;%樣本點距離中心點模的平方
????????m_wei（ij）=1-dist/h;%Epannechnikov核函數
????end
end
C=1/sum（sum（m_wei））;%歸一化系數

%計算目標權值直方圖
hist1=zeros（14096）;
for?i=1:a
????for?j=1:b
????????%rgb顏色空間量化為16*16*16?bins，降維來簡化彩色顏色直方圖
????????%把R，G，B三個通道0-255區間內的值都分為16個區間，即將0-255轉換為0-15
????????q_r=fix（double（temp（ij1））/16）;%fix為趨近0取整函數
????????q_g=fix（double（temp（ij2））/16）;
????????q_b=fix（double（temp（ij3））/16）;
????????q_temp=q_r*256+q_g*16+q_b;%設置每個像素點紅色、綠色、藍色分量所占比重，將三維向量映射為一維的
????????hist1（q_temp+1）=hist1（q_temp+1）+m_wei（ij）;%計算直方圖統計中每個像素點占的權重
????end
end
hist1=hist1*C;
rect（3）=ceil（rect（3））;
rect（4）=ceil（rect（4））;
jishu=1;
%讀取序列圖像
for?i=0329:0380
???????Im=imread（[‘0‘int2str（i）‘.jpg‘]）;
???????jishu=jishu+1;
???????%Im=medfilt2（Im）;
???????num=0;
???????Y=[22];
???????tic
????%mean?shift迭代
????while（（Y（1）^2+Y（2）^2>0.5）&num<20）%迭代條件
???????num=num+1;
???????temp1=imcrop（Imrect）;?
???????%計算侯選區域直方圖
???????%hist2=C*wei_hist（temp1m_weiab）;%target?candidates?pu???????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????????
???????hist2=zeros（14096）;
???????for?i=1:a
?????????for?j=1:b
????????q_r=fix（double（temp1（ij1））/16）;
????????q_g=fix（double（temp1（ij2））/16）;
????????q_b=fix（double（temp1（ij3））/16）;
????????q_temp1（ij）=q_r*256+q_g*16+q_b;
????????hist2（q_temp1（ij）+1）=?hist2（q_temp1（ij）+1）+m_wei（ij）;
?????????end
???????end
???????hist2=hist2*C;
???????%顯示直方圖
???????figure（2）;
???????subplot（121）;
???????plot（hist2）;
???????hold?on;
???????%bdist1?=?bhattacharrya（hist1?hist2）
???????w=weights（hist1hist2）;
???????%w=sqrt（hist1./hist2）;
???????
???????%變量初始化
???????sum_w=0;
???????xw=[00];
???????for?i=1:a
???????????for?j=1:b
??????????????%bitwei（ij）=w（uint32（q_temp1（ij））+1）;
???????????????sum_w=sum_w+w（uint32（q_temp1（ij））+1）;
???????????????xw=xw+w（uint32（q_temp1（ij））+1）*[i-y（1）-0.5j-y（2）-0.5];
???????????end
???????end
???????Y=xw/sum_w;
???????%中心點位置更新
???????rect（1）=rect（1）+Y（2）;
???????rect（2）=rect（2）+Y（1）;
????end
????%跟蹤軌跡矩陣%
????tic_x=[tic_x;rect（1）+rect（3）/2];
????tic_y=[tic_y;rect（2）+rect（4）/2];
???????v

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2018-09-26?08:45??MeanShiftTracking\
?????文件???????51271??2008-11-01?15:15??MeanShiftTracking\0329.jpg
?????文件???????51303??2008-11-01?15:15??MeanShiftTracking\0330.jpg
?????文件???????51393??2008-11-01?15:15??MeanShiftTracking\0331.jpg
?????文件???????51406??2008-11-01?15:15??MeanShiftTracking\0332.jpg
?????文件???????51413??2008-11-01?15:16??MeanShiftTracking\0333.jpg
?????文件???????51313??2008-11-01?15:16??MeanShiftTracking\0334.jpg
?????文件???????51406??2008-11-01?15:16??MeanShiftTracking\0335.jpg
?????文件???????51294??2008-11-01?15:16??MeanShiftTracking\0336.jpg
?????文件???????51331??2008-11-01?15:16??MeanShiftTracking\0337.jpg
?????文件???????51379??2008-11-01?15:16??MeanShiftTracking\0338.jpg
?????文件???????51210??2008-11-01?15:16??MeanShiftTracking\0339.jpg
?????文件???????51337??2008-11-01?15:16??MeanShiftTracking\0340.jpg
?????文件???????51169??2008-11-01?15:16??MeanShiftTracking\0341.jpg
?????文件???????51281??2008-11-01?15:16??MeanShiftTracking\0342.jpg
?????文件???????51367??2008-11-01?15:16??MeanShiftTracking\0343.jpg
?????文件???????51297??2008-11-01?15:16??MeanShiftTracking\0344.jpg
?????文件???????51305??2008-11-01?15:16??MeanShiftTracking\0345.jpg
?????文件???????51435??2008-11-01?15:16??MeanShiftTracking\0346.jpg
?????文件???????51305??2008-11-01?15:16??MeanShiftTracking\0347.jpg
?????文件???????51378??2008-11-01?15:16??MeanShiftTracking\0348.jpg
?????文件???????51562??2008-11-01?15:16??MeanShiftTracking\0349.jpg
?????文件???????51565??2008-11-01?15:16??MeanShiftTracking\0350.jpg
?????文件???????51663??2008-11-01?15:16??MeanShiftTracking\0351.jpg
?????文件???????51619??2008-11-01?15:16??MeanShiftTracking\0352.jpg
?????文件???????51456??2008-11-01?15:16??MeanShiftTracking\0353.jpg
?????文件???????51647??2008-11-01?15:16??MeanShiftTracking\0354.jpg
?????文件???????51472??2008-11-01?15:16??MeanShiftTracking\0355.jpg
?????文件???????51584??2008-11-01?15:16??MeanShiftTracking\0356.jpg
?????文件???????51871??2008-11-01?15:16??MeanShiftTracking\0357.jpg
?????文件???????51759??2008-11-01?15:16??MeanShiftTracking\0358.jpg
............此處省略95個文件信息