|
|
EDA365欢迎您登录!
您需要 登录 才可以下载或查看,没有帐号?注册
x
5 f0 \5 _5 M+ u* y: C一、源代码
1 Z3 C" O1 u7 S
7 W- n' B& w( N% B) t- function chengxu()
- %第1步
- close all
- I=imread('1.jpg'); %读取图像
- I=rgb2gray(I); %彩色图像转换成灰度图
- % I=im2bw(I); %二值化
- I=edge(double(I)); %检测图像的边缘
- figure
- imshow(I) %显示边缘检测的结果
- %第2步
- [m,n]=size(I); %计算图像的尺寸
- M=3; %定义X方向分割的块数
- N=3; %定义Y方向分割的块数
- mm=floor(m/M); %子块行的长度
- nn=floor(n/N); %子块列的长度
- count=1; %计数器
- figure
- for i=1:M
- for j=1:N
- A=I((i-1)*mm+1:i*mm,(j-1)*nn+1:j*nn); %分割原图像,得到一个子块
- subplot(M,N,count)
- imshow(A) %显示一个子块
- zuoshangjiao=[(i-1)*mm+1 (j-1)*nn+1]; %子块左上角的坐标
- [x,y,k,b]=zikuai(A,zuoshangjiao); %得到子块里白色像素点拟合得到的直线的斜率k和截距b(调用zikuai函数)
- X{count}=x; %保存子块里所有白色像素的x坐标
- Y{count}=y; %保存子块里所有白色像素的y坐标
- K(count)=k; %保存子块里拟合得到的直线的斜率k
- B(count)=b; %保存子块里拟合得到的直线的截距b
- count=count+1; %计数器加1,进行下一个子块的计算
- end
- end' v& _0 s+ N/ |; V8 J- }6 F
6 O% A m/ O7 m- a+ o$ e5 l, U: ?- T, K8 A& v
二、运行结果" L5 W. Q" x: m: W7 o* L
7 l& ?* Y7 k0 [! Z6 i% J5 ]# G& Q+ v6 D2 u8 Y
( _* b6 P! E6 c% u" f- n2 I
. w0 D1 b* X( X W- m0 A. F3 b3 ]& a
, g; H* |( m/ W/ Z% b5 X4 e
_& z. l4 c+ V
3 }! m' A7 N2 ^9 l; j7 U1 F' x* m8 L' ]' z% O! e) v8 b! m
( ?9 c* o* V X; T6 C+ K6 O. {+ V# x- E
: b1 G5 i" Y$ @# ^/ p8 w
|
|