91av视频/亚洲h视频/操亚洲美女/外国一级黄色毛片 - 国产三级三级三级三级

資源簡介

首先進行Harris角點特征提取,而后利用NCC算法進行粗匹配,剔除誤匹配和不匹配向量,基于灰度相關系數,計算配準誤差, 從而得到配準后疊加圖像。適用于輸入圖像存在平移變換和旋轉變換;能夠實現可見光區圖像配準,可實現時間間隔不是很久的多時相遙感影像配準

資源截圖

代碼片段和文件信息

function?[v1_gcpu1_gcpv2_gcpu2_gcp]=GCP_Select2(v1u1v2u2I1I2)
[height1width1]=size(I1);
[height2width2]=size(I2);
%采用NCC(歸一化互相關)算法對兩幅圖像的角點進行粗匹配%
k=15;??????????????????????????????????????????????????????????????????????%鄰域窗口大小(2k+1)*(2k+1))
p=1;
R=zeros(1length(u2))?*?(1e-10);
[v1_new?u1_newv2_newu2_new]?=?deal(zeros(1length(v1)));

for?i=1:length(v1)
????if(((u1(i)+k)0)&&((v1(i)-k)>0))
????????G1_window?=?I1(?v1(i)-k:v1(i)+k?u1(i)-k:u1(i)+k?);
????????G1_avg?=?mean(G1_window(:));
????????for?j=1:length(v2)
????????????G_sum=0;G1_sub=0;G2_sub=0;
????????????if(((u2(j)+k)0)&&((v2(j)-k)>0))&&(abs(v1(i)-v2(j))<80)&&(abs(u1(i)-u2(j))<120)
????????????????G2_window?=?I2(?v2(j)-k:v2(j)+k?u2(j)-k:u2(j)+k?);?????
????????????????G2_avg?=?mean(G2_window(:));
????????????????
????????????????G_sum?=?sum((G1_window(:)?-?G1_avg).*(G2_window(:)?-?G2_avg));
????????????????G1_sub?=?sum((G1_window(:)?-?G1_avg).^2);
????????????????G2_sub?=?sum((G2_window(:)?-?G2_avg).^2);
????????????end
????????????R(j)?=?G_sum?/?(sqrt(G1_sub*G2_sub)?+?1e-10);
????????end
????????R_max=max(R);
????????if?(R_max>0.6)???????????????????????????????????????????????????%閾值0.5,可調,視需要粗配準點數而定
????????????v1_new(p)=v1(i);
????????????u1_new(p)=u1(i);
????????????[~y]=find(R==R_max);????????????????????????????????????????%構建新的角點位置向量,該向量中I1I2角點一一對應
????????????v2_new(p)=v2(y);
????????????u2_new(p)=u2(y);
????????????p=p+1;
????????end
????end
end

p?=?p-1;
[v1_new?u1_newv2_newu2_new]?=?deal(v1_new(1:p)?u1_new(1:p)v2_new(1:p)u2_new(1:p));

%通過簡單方法對匹配點對進行篩選,剔除錯誤匹配點對(剔除距離出現次數最多之外的點),統計兩圖像坐標差值最多的坐標為新的匹配點集%
x_sub=zeros(1length(u1_new));
y_sub=zeros(1length(u1_new));
for?i=1:length(u1_new)
????x_sub(i)=u1_new(i)-u2_new(i);
????y_sub(i)=v1_new(i)-v2_new(i);
end
x=unique(x_sub);
[mn]=hist(x_subx);
x_sub_most=n(find(m==max(m)));
[mn]=hist(y_suby);
y_sub_most=n(find(m==max(m)));
p=1;
[v1_new2?u1_new2v2_new2u2_new2]?=?deal(zeros(1length(u1_new)));%剔除錯誤匹配點對之后的新點對起始坐標存儲位置
for?i=1:length(u1_new)
????if(?(abs(x_sub(i)-x_sub_most(1))<1)&&(abs(y_sub(i)-y_sub_most(1))<1))
????????v1_new2(p)=v1_new(i);u1_new2(p)=u1_new(i);
????????v2_new2(p)=v2_new(i);u2_new2(p)=u2_new(i);
????????p=p+1;
????end
end
p?=?p-1;
[v1_new2?u1_new2v2_new2u2_new2]?=?deal(v1_new2(1:p)?u1_new2(1:p)v2_new2(1:p)u2_new2(1:p));
[v1_gcp?u1_gcp?v2_gcp?u2_gcp]?=?deal(v1_new2?u1_new2v2_new2u2_new2);

%?%?將I1作為參考圖像,在I2中的對應角點附近尋找I1角點的同名點?%
[ghz]?=?deal(1);??%搜索窗口大小為2*k+1??%灰度窗口大小為2*h+1????%重新篩選之后的地面控制點起始存儲位置
R_G1G2=zeros(33);
for?i=1:length(u1_new2)
????G1_gary?=?I1(v1_new2(i)-h:v1_new2(i)+h?u1_new2(i)-h:u1_new2(i)+h);
????G1_avg?=?mean(G1_gary(:));????????????????????????????????????????????????????????????????%生成I1中灰度窗口中的灰度向量
????for?k=-g:1:g
????????for?l=-g:1:g
????????????G2_gary?=?I2(v2_new2(i)+k-h:v2_new2(i)+k+h?u2_new2(i)+l-h:

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????4035??2019-05-29?20:39??registration\GCP_Select2.m

?????文件???????2107??2018-01-12?11:01??registration\Harris.m

?????文件?????786624??2018-01-12?11:01??registration\I1.tif

?????文件?????786624??2018-01-12?11:01??registration\I2.tif

?????文件?????184572??2018-01-12?11:01??registration\I3.tif

?????文件?????186756??2018-01-12?11:01??registration\I4.tif

?????文件??????89584??2018-01-12?11:01??registration\I5.jpg

?????文件??????86935??2018-01-12?11:01??registration\I6.jpg

?????文件???????2452??2019-05-29?19:59??registration\ImgRegistration.m

?????文件???????3296??2018-01-12?11:01??registration\Registration.m

?????文件???????1758??2019-05-29?15:23??registration\registration_manual.m

?????文件????????673??2018-09-07?14:17??registration\數據\1\1.txt

?????文件?????254517??2018-09-24?21:07??registration\數據\1\Sar(從).png

?????文件?????233601??2018-09-24?19:46??registration\數據\1\ZY(主).png

?????文件????????679??2018-09-07?14:25??registration\數據\2\2.txt

?????文件?????190726??2018-05-18?16:12??registration\數據\2\主.png

?????文件?????258692??2018-09-26?10:13??registration\數據\2\從.png

?????文件???????1281??2018-05-18?20:20??registration\數據\3\lt5tosar.txt

?????文件??????40831??2018-09-25?20:12??registration\數據\3\主.png

?????文件??????60592??2018-05-18?18:14??registration\數據\3\從.png

?????文件??????63481??2018-05-18?20:38??registration\數據\4\dataenvisat.png

?????文件??????47437??2018-09-25?20:09??registration\數據\4\datalt5.png

?????文件???????1298??2018-05-18?21:29??registration\數據\4\lt5tosar.txt

?????文件????????390??2018-09-07?14:29??registration\數據\5\5.txt

?????文件??????36428??2018-09-07?09:07??registration\數據\5\主.png

?????文件??????83199??2018-09-07?09:07??registration\數據\5\從.png

?????目錄??????????0??2019-09-05?08:24??registration\數據\1

?????目錄??????????0??2019-09-05?08:24??registration\數據\2

?????目錄??????????0??2019-09-05?08:24??registration\數據\3

?????目錄??????????0??2019-09-05?08:24??registration\數據\4

............此處省略6個文件信息

評論

共有 條評論