| 2 | 1/1 | 返回列表 |
| 查看: 432 | 回復(fù): 1 | ||
shawmarkk銀蟲(chóng) (初入文壇)
|
[求助]
求助用于電子穩(wěn)像的灰度投影算法源碼 已有1人參與
|
| 求助用于電子穩(wěn)像的灰度投影算法源碼!急急急! |
鐵蟲(chóng) (小有名氣)
|
%投影法,投影修正,余弦濾波 %水平方向最大偏移正負(fù)20個(gè)象素 %垂直方向最大偏移正負(fù)30個(gè)象素 clear all mov=aviread('v2.avi'); g=[]; for t=50:80 mov(t).cdata=rgb2gray(mov(t).cdata); g=cat(3,g,mov(t).cdata); end gr=[]; gr=cat(3,gr,g(:,:,1)); frow=30; fcol=20; referenceframe=g(:,:,1); refprojrow=zeros(240,1); refprojcol=zeros(1,320); %行投影 refprojrowsum=0; for i=1:240 for j=1:320 refprojrow(i)=refprojrow(i)+double(referenceframe(i,j)); end refprojrowsum=refprojrowsum+refprojrow(i); end refprojrowmean=refprojrowsum/240; refprojrow=refprojrow-refprojrowmean; figure;plot(refprojrow); for i=1:240 if (i240-frow) refprojrow(i)=refprojrow(i)*(1+cos(pi*(frow-1-i)/frow))/2; end end hold on plot(refprojrow,'r'); %列投影 refprojcolsum=0; for j=1:320 for i=1:240 refprojcol(j)=refprojcol(j)+double(referenceframe(i,j)); end refprojcolsum=refprojcolsum+refprojcol(j); end refprojcolmean=refprojcolsum/320 refprojcol=refprojcol-refprojcolmean; figure;plot(refprojcol); for j=1:320 if (j240-fcol) refprojcol(j)=refprojcol(j)*(1+cos(pi*(fcol-1-j)/fcol))/2; end end hold on ; plot(refprojcol,'r'); for t=2:20 t currentframe=g(:,:,t); %行投影 curprojrow=zeros(240,1); curprojrowsum=0; for i=1:240 for j=1:320 %refprojrow(i)=refprojrow(i)+double(referenceframe(i,j)); curprojrow(i)=curprojrow(i)+double( currentframe(i,j)); end curprojrowsum=curprojrowsum+curprojrow(i); end curprojrowmean=curprojrowsum/240; curprojrow=curprojrow-curprojrowmean; for i=1:240 if (i240-frow) curprojrow(i)=curprojrow(i)*(1+cos(pi*(frow-1-i)/frow))/2; end end %figure;plot(refprojrow);hold on %plot(curprojrow,'r');xlabel('行投影') %180單位的相關(guān)運(yùn)算 %垂直方向最大偏移正負(fù)30個(gè)象素 cr=zeros(1,61); for j=1:61 for i=1:180 cr(j)=cr(j)+(refprojrow(j+i-1)-curprojrow(30+i))^2; end end %figure;plot(cr,'g'); xlabel('垂直方向') [b,jmin]=min(cr); dy=31-jmin %列投影 curprojcol=zeros(1,320); curprojcolsum=0; for j=1:320 for i=1:240 %refprojcol(j)=refprojcol(j)+double(referenceframe(i,j)); curprojcol(j)=curprojcol(j)+double(currentframe(i,j)); end curprojcolsum=curprojcolsum+curprojcol(j); end curprojcolmean=curprojcolsum/320; curprojcol=curprojcol-curprojcolmean; for j=1:320 if (j240-fcol) curprojcol(j)=curprojcol(j)*(1+cos(pi*(fcol-1-j)/fcol))/2; end end %figure;plot(refprojcol);hold on; %plot(curprojcol,'r'); xlabel('列投影') %水平方向最大偏移正負(fù)20個(gè)象素 %280單位的相關(guān)運(yùn)算 cc=zeros(1,41); for i=1:41 for j=1:280 cc(i)=cc(i)+(refprojcol(j+i-1)-curprojcol(20+j))^2; end end %figure;plot(cc,'g'); xlabel('水平方向') [a,imin]=min(cc); dx=21-imin if dy<0 %dy<0,圖currentframe相對(duì)于referenceframe向上運(yùn)動(dòng)了|dy| a=zeros(abs(dy),320); guoduframe=[a;currentframe(1:240-abs(dy), ]; elseif dy>0 %dy>0,圖currentframe相對(duì)于referenceframe向下運(yùn)動(dòng)了|dy|, a=zeros(dy,320); guoduframe=[currentframe(1+abs(dy):240, ;a]; else guoduframe=currentframe; end if dx<0 %dx<0,圖currentframe相對(duì)于referenceframe向左運(yùn)動(dòng)了|dx| c=zeros(240,abs(dx)); buchangframe=[c,guoduframe(:,1:320-abs(dx))]; elseif dx>0 %dx>0,圖currentframe相對(duì)于referenceframe向右運(yùn)動(dòng)了|dx| c=zeros(240,dx); buchangframe=[guoduframe(:,abs(dx)+1:320),c]; else buchangframe=guoduframe; end %figure;imshow(buchangframe); gr=cat(3,gr,buchangframe); %referenceframe=buchangframe; %refprojrow=curprojrow; %refprojcol=curprojcol; end |
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 320分,材料與化工專(zhuān)業(yè),求調(diào)劑 +9 | 一定上岸aaa 2026-03-27 | 13/650 |
|
|---|---|---|---|---|
|
[考研]
|
y7czhao 2026-03-26 | 10/500 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見(jiàn)α 2026-03-25 | 3/150 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機(jī)化學(xué) +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
[考研] 藥學(xué)105500求調(diào)劑 +3 | Ssun。。 2026-03-28 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué)324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
[考研] 283求調(diào)劑(080500) +4 | A child 2026-03-27 | 4/200 |
|
|
[考研] 085600材料與化工調(diào)劑 +10 | A-哆啦Z夢(mèng) 2026-03-23 | 16/800 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫(kù) 2026-03-26 | 5/250 |
|
|
[考研] 一志愿211,335分,0856,求調(diào)劑院校和導(dǎo)師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[考研] 333求調(diào)劑 +6 | wfh030413@ 2026-03-23 | 6/300 |
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
[考研] 環(huán)境專(zhuān)碩324分求調(diào)劑推薦 +5 | 軒小寧—— 2026-03-26 | 5/250 |
|
|
[考研] 一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研經(jīng)歷 +6 | cccchenso 2026-03-23 | 6/300 |
|
|
[考研] 296求調(diào)劑 +4 | 汪?! 2026-03-25 | 7/350 |
|
|
[考研] 0703化學(xué)調(diào)劑,求導(dǎo)師收 +7 | 天天好運(yùn)來(lái)上岸?/a> 2026-03-24 | 7/350 |
|
|
[考研] 300求調(diào)劑,材料科學(xué)英一數(shù)二 +5 | leaflight 2026-03-24 | 5/250 |
|
|
[考研] 315分,誠(chéng)求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 306求調(diào)劑 +5 | 來(lái)好運(yùn)來(lái)來(lái)來(lái) 2026-03-22 | 5/250 |
|