| 2 | 1/1 | 返回列表 |
| 查看: 433 | 回復: 1 | |||
shawmarkk銀蟲 (初入文壇)
|
[求助]
求助用于電子穩(wěn)像的灰度投影算法源碼 已有1人參與
|
| 求助用于電子穩(wěn)像的灰度投影算法源碼!急急急! |
|
%投影法,投影修正,余弦濾波 %水平方向最大偏移正負20個象素 %垂直方向最大偏移正負30個象素 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單位的相關運算 %垂直方向最大偏移正負30個象素 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('列投影') %水平方向最大偏移正負20個象素 %280單位的相關運算 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相對于referenceframe向上運動了|dy| a=zeros(abs(dy),320); guoduframe=[a;currentframe(1:240-abs(dy), ]; elseif dy>0 %dy>0,圖currentframe相對于referenceframe向下運動了|dy|, a=zeros(dy,320); guoduframe=[currentframe(1+abs(dy):240, ;a]; else guoduframe=currentframe; end if dx<0 %dx<0,圖currentframe相對于referenceframe向左運動了|dx| c=zeros(240,abs(dx)); buchangframe=[c,guoduframe(:,1:320-abs(dx))]; elseif dx>0 %dx>0,圖currentframe相對于referenceframe向右運動了|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ā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿北京理工大學本科211材料工程294求調劑 +4 | mikasa的圍巾 2026-03-28 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 0856材料化工調劑 總分330 +11 | zhubinhao 2026-03-27 | 11/550 |
|
|
[考研] 295材料工程專碩求調劑 +7 | 1428151015 2026-03-27 | 7/350 |
|
|
[考研] 321求調劑 +6 | 材料cailiao 2026-03-21 | 6/300 |
|
|
[考研] 求調劑 +6 | 蘆lty 2026-03-25 | 7/350 |
|
|
[考研] 一志愿哈爾濱工業(yè)大學材料與化工方向336分 +8 | 辰沐5211314 2026-03-26 | 8/400 |
|
|
[考研] 277跪求調劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
[考研]
|
18419759900 2026-03-25 | 8/400 |
|
|
[考研] 287求調劑 +10 | land xuxu 2026-03-26 | 10/500 |
|
|
[考研] 085600,材料與化工321分,求調劑 +9 | 大饞小子 2026-03-27 | 9/450 |
|
|
[考研] 一志愿211,335分,0856,求調劑院校和導師 +4 | 傾____蕭 2026-03-27 | 5/250 |
|
|
[考研] 286求調劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
[考研] 一志愿陜師大生物學071000,298分,求調劑 +5 | SYA! 2026-03-23 | 5/250 |
|
|
[考研] 284求調劑 +11 | junqihahaha 2026-03-26 | 12/600 |
|
|
[考研] 294分080500材料科學與工程求調劑 +4 | 柳溪邊 2026-03-26 | 4/200 |
|
|
[考研] 085601求調劑總分293英一數二 +4 | 鋼鐵大炮 2026-03-24 | 4/200 |
|
|
[考研] 07化學303求調劑 +5 | 睿08 2026-03-25 | 5/250 |
|
|
[考研] 26考研-291分-廈門大學(085601)-柔性電子學院材料工程專業(yè)求調劑 +3 | min3 2026-03-24 | 4/200 |
|
|
[考研] 上海電力大學材料防護與新材料重點實驗室招收調劑研究生(材料、化學、電化學,環(huán)境) +4 | 我愛學電池 2026-03-23 | 4/200 |
|
|
[考研] 292求調劑 +4 | 鵝鵝鵝額額額額?/a> 2026-03-24 | 4/200 |
|