clear;clc;
%?水平（左右）拼接
%?dataA=imread（‘D:\YlinkSample\IDT\testpicture\test2\1003.JPG‘）;
%?dataB=imread（‘D:\YlinkSample\IDT\testpicture\test2\1002.JPG‘）;
%?rdata1?=?-26;???%重合區域寬度

dataA=imread（‘D:\YlinkSample\IDT\testpicture\1.JPG‘）;
dataB=imread（‘D:\YlinkSample\IDT\testpicture\2.JPG‘）;
rdata1?=?-76;???%重合區域寬度


A=rgb2gray（dataA）;
B=rgb2gray（dataB）;


%《圖像拼接的改進算法》最佳縫合線算法?圖像融合
%先根據之前得到的H矩陣計算重疊區域Rect
[HW~]=size（A）;
L?=?W+1+rdata1;
R?=?W;
n?=?-rdata1;

%最終融合圖的大小
Y=2*W+rdata1+1;?%融合后圖像寬度
D=zeros（HY-13）;???%融合后的圖像
path=zeros（Hn）;%放路徑的矩陣，存放列數序號（1-n）
%放強度值?每條路徑的強度值strength=color^2+geometry
color=zeros（1n）;
geometry=zeros（1n）;

strength1=zeros（1n）;
strength2=zeros（1n）;

Axdao?=?0;
Aydao?=?0;
Bxdao?=?0;
Bydao?=?0;

%計算第一行即初始化的強度值
i?=?1;
for?j=L:R
????y=j-L+1;?%對應B圖的列數
????color（y）=A（ij）-B（iy）;?%灰度值之差
????if（j==L）
????????Bxdao=B（iy+1）+2*B（i+1y+1）;
????????Bydao=B（i+1y）+2*B（i+1y+1）;
????????Axdao=A（ij+1）+2*A（i+1j+1）-A（ij-1）-2*A（i+1j-1）;
????????Aydao=2*A（i+1j-1）+A（i+1j）+2*A（i+1j+1）;
????elseif（j==R）
????????Axdao=?-A（ij-1）-2*A（i+1j-1）;?
????????Aydao=2*A（i+1j-1）+A（i+1j）;
????????Bxdao=B（iy+1）+2*B（i+1y+1）-B（iy-1）-2*B（i+1y-1）;
????????Bydao=2*B（i+1y-1）+B（i+1y）+2*B（i+1y+1）;
????else
????????Axdao=A（ij+1）+2*A（i+1j+1）-A（ij-1）-2*A（i+1j-1）;
????????Bxdao=B（iy+1）+2*B（i+1y+1）-B（iy-1）-2*B（i+1y-1）;
????????Aydao=2*A（i+1j-1）+A（i+1j）+2*A（i+1j+1）;
????????Bydao=2*B（i+1y-1）+B（i+1y）+2*B（i+1y+1）;
????end
????geometry（y）=（Axdao-Bxdao）*（Aydao-Bydao）;
????strength1（y）=color（y）^2+geometry（y）;
????path（1y）=y;
end
color?=?color?.*?0;
geometry?=?geometry?.*?0;

%開始擴展?向下一行?從第二行到倒數第二行?最后一行單獨拿出來?像第一行一樣?因為它的結構差值geometry不好算
for?i=2:H-1
????%先把下一行的強度值全部計算出來?到時候需要比較哪三個就拿出哪三個
????for?j=L:R
????????x?=?i;
????????y?=?j-L+1;
????????color（y）?=?A（ij）?-?B（xy）;
????????if（j==L）
????????????Axdao=2*A（i-1j+1）+A（ij+1）+2*A（i+1j+1）-2*A（i-1j-1）-A（ij-1）-2*A（i+1j-1）;
????????????Bxdao=2*B（x-1y+1）+B（xy+1）+2*B（x+1y+1）;
????????????Aydao=-2*A（i-1j-1）-A（i-1j）-2*A（i-1j+1）+2*A（i+1j-1）+A（i+1j）+2*A（i+1j+1）;
????????????Bydao=-B（x-1y）-2*B（x-1y+1）+B（x+1y）+2*B（x+1y+1）;
????????elseif（j==R）
????????????Axdao=-2*A（i-1j-1）-A（ij-1）-2*A（i+1j-1）;
????????????Bxdao=2*B（x-1y+1）+B（xy+1）+2*B（x+1y+1）-2*B（x-1y-1）-B（xy-1）-2*B（x+1y-1）;
????????????Aydao=-2*A（i-1j-1）-A（i-1j）+2*A（i+1j-1）+A（i+1j）;
????????????Bydao=-2*B（x-1y-1）-B（x-1y）-2*B（x-1y+1）+2*B（x+1y-1）+B（x+1y）+2*B（x+1y+1）;
????????else
????????????Axdao=2*A（i-1j+1）+A（ij+1）+2*A（i+1j+1）-2*A（i-1j-1）-A（ij-1）-2*A（i+1j-1）;
????????????Bxdao=2*B（x-1y+1）+B（xy+1）+2*B（x+1y+1）-2*B（x-1y-1）-B（xy-1）-2*B（x+1y-1）;
????????????Aydao=-2*A（i-1j-1）-A（i-1j）-2*A（i-1j+1）+2*A（i+1j-1）+A（i+1j）+2*A（i+1j+1）;
????????????Bydao=-2*B（x-1y-1）-B（x-1y）-2*B（x-1y+1）+2*B（x+1y-1）+B（x+1y）+2*B（x+1y+1）;
????????end
????????geometry（y）=（Axdao-Bxdao）*（Aydao-Bydao）;
????????strength2（y）=color（y）^2+geometry（y）;
????e

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

?????文件??????18509??2018-04-08?14:29??圖像最佳縫合線拼接算法\1.jpg

?????文件???????5045??2018-04-08?14:29??圖像最佳縫合線拼接算法\2.jpg

?????文件???????8723??2018-04-13?13:48??圖像最佳縫合線拼接算法\bestlinefusionhorizontal.m

?????文件?????????59??2018-04-13?13:50??圖像最佳縫合線拼接算法\參考博客.txt

?????文件?????189438??2018-03-28?15:59??圖像最佳縫合線拼接算法\圖像拼接的改進算法_方賢勇.caj

?????目錄??????????0??2018-04-13?13:50??圖像最佳縫合線拼接算法

-----------?---------??----------?-----??----

???????????????221774????????????????????6