%附錄7：求擬合圖與原切片圖的重合度的matlab程序代碼（baifenbi1.m）
function?baifenbi=baifenbi1（pnjjpn）%輸出擬合圖與原切片圖的重合度
%填充新圖
juzheng=pnjj;?%pnjj為通過dian.m得到的輪廓邊界二值矩陣
%先填充左邊界的右半部分
for?i=1:511
?????for?j=1:511
?????????if（pnjj（ij）==0&pnjj（ij+1）~=0）
????????pnjj（ij+1）=pnjj（ij）;????????
????end
??end
end
you=pnjj;?
%再填充右邊界的左半部分?
for?i=1:512???
?????for?j=512:-1:2??????????
?????????if（juzheng（ij）==0&juzheng（ij-1）~=0）?
?????????????juzheng（ij-1）=juzheng（ij）;??
?????????end
?????end
end
zuo=juzheng;?
shijijuzheng=you|zuo;?%通過矩陣的或運算得到填充后的新圖?
imshow（you|zuo）?
%原圖的黑點的個數?
biaozhunjuzheng=imread（strcat（int2str（pn）‘.bmp‘））;
nbiao=0;?
for?i=1:512?????
????for?j=1:512??????????
????????if（biaozhunjuzheng（ij）==0）????????????
????????????nbiao=nbiao+1;?????????
????????end
????end
end
%新圖與原圖重合部分黑點的個數?
chonghegs=0;
for?i=1:512?????
????for?j=1:512
????????if（biaozhunjuzheng（ij）==0&shijijuzheng（ij）==0）
???????????chonghegs=chonghegs+1;
????????end
????end
end
%求百分比
baifenbi=chonghegs/nbiao;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1139??2015-04-15?09:56??三維重建\baifenbi1.m
?????文件????????1472??2015-04-15?09:49??三維重建\dian.m
?????文件????????1112??2015-04-15?15:50??三維重建\ff.m
?????文件?????????280??2015-04-14?18:32??三維重建\pczx.m
?????文件?????????452??2015-04-15?15:58??三維重建\xueguansanweitu.m
?????文件?????????702??2015-04-15?09:58??三維重建\zhongzhouxiantouying.m
?????文件?????????231??2015-04-15?09:58??三維重建\zhuanhua.m
?????目錄???????????0??2015-04-15?16:03??三維重建\