%此程序的功能是提取網格交叉點并在原圖上標記出來

clear
close?all
clc

img=imread（‘qp1.jpg‘）;
imshow（img）title（‘原圖‘）

%水平方向投影
for?i=1:105
a（i）=sum（img（1:100i）1）;
end

%垂直方向投影
for?i=1:105
b（i）=sum（img（i1:100）2）;
end


maxx=max（a[]2）;mdx=maxx/2;
maxy=max（b[]2）;mdy=maxy/2;
c=find（a>mdx）;
d=find（b>mdy）;

e=repmat（c‘110）;
f=reshape（e‘1001）;
g=repmat（d101）;
h=reshape（g‘1001）;

for?i=1:100
k（i1）=f（i1）;
k（i2）=h（i1）;
end

m=k（:1）;
n=k（:2）;
m=m‘;n=n‘;

figureimshow（img）
title（‘red?+?is?the?crossing?point‘）
hold?on
plot（mn‘r.‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????5106??2011-09-09?21:36??網格交叉點的提取\qp1.jpg
?????文件???????58773??2011-10-25?17:15??網格交叉點的提取\qpg52.jpg
?????文件?????????607??2011-10-25?17:31??網格交叉點的提取\wgjcd.m
?????目錄???????????0??2016-01-11?16:47??網格交叉點的提取\