- 在线时间
- 7 小时
- 最后登录
- 2015-2-10
- 注册时间
- 2012-3-22
- 听众数
- 3
- 收听数
- 0
- 能力
- 0 分
- 体力
- 882 点
- 威望
- 0 点
- 阅读权限
- 30
- 积分
- 270
- 相册
- 0
- 日志
- 0
- 记录
- 0
- 帖子
- 9
- 主题
- 0
- 精华
- 0
- 分享
- 0
- 好友
- 1
升级   85% TA的每日心情 | 无聊 2015-2-6 09:42 |
|---|
签到天数: 9 天 [LV.3]偶尔看看II
|
B题
+ B; r, r6 ~& W9 U' l9 v( t8 t主要是应用matlab所提供的灰度共生矩阵
, W" e' M4 W7 ^4 x5 p graycomatrix()
& [9 k( u* R6 d9 i/ A) `5 d9 [ graycoprops();
" M9 H% ~7 H. [ 如果人工干预(例如找到最左边)对于附件1是百分百匹配
0 l0 a7 S e$ C1 A `% ~7 N 附件2是匹配出两个大块1 @% {9 w& D8 l8 d0 w4 v3 m* b
附件3,4,5方法类似* ~! w: L N* P# j3 H5 Y" e9 C
function [d] = readpic3()! p. W# l. g* V+ `; \
%人工看出 008.bmp 第9张图片 为最左边 将008.改为000; t* c' M( n: Z$ m# v7 y/ _
%filename = - h" a; \* I) O" N/ J- N4 ?; s
dirname='E:\Program Files\MATLAB\workfile';
& S" Q+ t" g% afiles=dir([dirname,'\*.bmp']);
7 G5 Z* Q- d6 Z JNum_file=numel(files);
3 e% r# V! U! U( v! wdata=cell(1,Num_file);1 X0 w: n6 z6 O% t5 s& R! X6 b
[M,N]=size(data{1});% M = 1980 N = 728 o0 J* c+ V# L2 U' v
a = zeros(M,N);* H. k/ j& ^& \5 ~3 T# h9 z
b = zeros(M,N);
4 o/ [% }- D e6 zc = zeros(Num_file);
1 [' I, X& H: T8 h( N6 p4 md = zeros(1,Num_file);%存储排序序列6 x! R ?: F+ W
for k = 1:Num_file
/ d* T* G, |: ~% h. J data{k} = imread([dirname '\' files(k).name]);1 K; c! f5 w7 U) S# x6 e) F
end 1 A; L& ~' l7 k, T q
1 ?! @* v* e+ D- O4 n: Dfid = fopen('数据.txt', 'wt');$ I( U! d% L; R/ e
for i = 1:Num_file! V' \6 S8 x2 f" i: D8 Q0 h; O
a = data{i};
1 W+ i; C6 L3 [$ L) p$ r y" d e = ones(1,Num_file);
! y4 ?4 u$ l; y, [9 l for j = 1:Num_file % 要左右匹配两次
# h1 j9 g; y' o, I b = data{j};" X1 | V& k6 }8 ?9 e
img = [uint8(a),uint8(b)];
3 s7 E8 `$ f4 t- k img1 = img(1:1980,70:74); % *****很关键*****
" y5 N+ E! c$ |1 U glcm = graycomatrix(img1,'NumLevels',4,'offset',[0,2]);) e( I3 n* P6 _8 e4 r8 j4 M
% stats = graycoprops(glcm,{'Contrast','Homogeneity','Energy'});%这里仅考虑三个特征值: ^8 `, H* w+ A
% e(j) = abs(stats.Correlation-1)+abs(stats.Homogeneity-1)+abs(stats.Energy-1);%每个特征值减1取绝对值6 w6 F. q; x$ x S' ?: E5 K& h% \
stats = graycoprops(glcm,{'Contrast','Homogeneity','Energy','Correlation'});%这里仅考虑四个特征值
% U+ D7 V& T `$ X4 o2 ]; ^ e(j) = abs(stats.Correlation-1)+abs(stats.Homogeneity-1)+abs(stats.Energy-1)+ stats.Contrast;%每个特征值减1取绝对值$ x4 A6 M- h- q: W! @2 S+ N
fprintf(fid, '%f ',e(j));
" y" U% S# d! ? end( a9 y! Z6 x, k2 D, P! m* Z8 K6 l8 C
fprintf(fid, '\n');" N. @8 [+ G4 I% z
[A,B] = min(e);% A 最小值 B索引 B即为) I$ r, M4 ?9 ~3 f# }! z
d(i)= B-1;
5 W. J6 ?8 `3 q( n4 f( U0 Oend |
-
总评分: 体力 + 25
查看全部评分
|