QQ登录

只需要一步,快速开始

 注册地址  找回密码
查看: 2707|回复: 1
打印 上一主题 下一主题

[建模教程] 1-2、数字图像处理基础

[复制链接]
字体大小: 正常 放大
杨利霞        

5273

主题

82

听众

17万

积分

  • TA的每日心情
    开心
    2021-8-11 17:59
  • 签到天数: 17 天

    [LV.4]偶尔看看III

    网络挑战赛参赛者

    网络挑战赛参赛者

    自我介绍
    本人女,毕业于内蒙古科技大学,担任文职专业,毕业专业英语。

    群组2018美赛大象算法课程

    群组2018美赛护航培训课程

    群组2019年 数学中国站长建

    群组2019年数据分析师课程

    群组2018年大象老师国赛优

    跳转到指定楼层
    1#
    发表于 2019-3-5 16:20 |只看该作者 |正序浏览
    |招呼Ta 关注Ta
    详细资源请下载附件
    , m4 w8 n; S1 W& Q( E9 |$ `
    / v: }) l: P* m1 u数学建模题目中有时会涉及到与数字图像有关的操作。在这类题目中,往往不会涉及到太多与数字图像处理相关的专业知识,但是要求程序员了解图像存储格式与常用基础操作等。3 }' U' P# [3 x+ e$ {5 n
    7 l$ T- L- e3 e
    一、数字图像常用存储格式。- ~' {8 A3 Y% g1 j2 k

    2 V7 O  ~# l6 B" Y* i- Z" r& g数字图像在计算机中以矩阵形式存储,通过一个或多个数字表示每个点的数值大小来控制每个像素点的颜色变化。常见的存储类型有灰度图像、RGB图像、二值图像、HSI图像等等。在这里我们着重介绍灰度图像与RGB图像。2 {2 e4 y% x! M2 O+ o8 {. ], H
    0 w) \6 U/ X$ c1 s) i& m$ b6 E
    灰度图像,顾名思义,是一种单色图像,只有纯白、纯黑、和254个不同深度大小的介于纯白、纯黑之间的灰色强度。对于灰度图,每个像素的亮度由且仅由一个数值来表示,通常数值范围在0到255之间,0表示黑、255表示白,其它值表示处于黑白之间的灰度,数字越小就越白。由于灰度图的存储方式易于理解,比较简单且涵盖了较多的图像信息,在不需要彩色信息时我们常将图像转化为灰度图像进行进一步操作。   {$ J# n$ f* }' P
    3 [* D+ p( }! D0 D# ^0 f* q5 P
    二值图像,可以作为一种特殊的灰度图像来理解。他只有纯黑和纯白两种颜色,在计算机中每个像素点只有0、1两种取值可能。其中0表示黑色,1表示白色。 3 r3 M, ~+ `# N' d3 E) J
    7 S0 I% l5 g+ ~  x
    在中学学习美术时,我们都知道红、黄、蓝是三原色。通过这三个颜色的颜料画在一起叠加即可组合成任意一种颜色。在计算机里,也类似,但三原色不再是红、黄、蓝,而是红(R)、绿(G)、蓝(B)。在计算机里通过控制这三个颜色的量组合在一起,也可以合成任意一种颜色,于是就有了RGB图。RGB图像里,每个像素点由三个数值控制颜色,分别对应红、绿、蓝的分量大小。范围一般也为0~255之间,0表示这个颜色分量没有,255表示这个颜色分量取到最大值。例如:某RGB图中一个像素点的红、绿、蓝分量均为255,则根据光学叠加的理论可知,该点为纯白色。 8 x! \: P4 H  _8 Z) C0 r0 _& D

    1 M) [9 G5 J3 @2 {& `6 F6 H7 o  oHIS图像每个像素点由色调、饱和度、亮度决定。相比较之下,在数学建模类题目里并不常用。
    $ H4 C2 {8 n' @, f; i8 @1 I. D: y! K( ]$ @% V; k: v
    相关实例与代码:* v) M7 ^3 n3 D% F
    & v/ U: ]4 Z: [8 G8 v: d
    (1)假设matlab中有一张RGB图像存放于变量pic中,将其转化灰度图像。
    8 N3 L& F/ D2 s! n1 @& T5 U9 X' w1 J% r/ f. N, T. {" _5 h% j
    newPic= rgb2gray('pic')" C2 \' A3 F! i4 t) v9 y
    1
    % ~, S& l* C- y8 I: t注:rgb2gray函数中的2,按照英文读作two,与to同音。这个函数读为rgb to gray,同理,其他图像存储类型转化函数也类似此函数。
    ( ?! O8 |) T3 ^' H8 ?, g4 r  x
    6 \/ L) G: Q* S$ U9 [/ A5 d: C4 D(2)生成一幅128*128的RGB图像,该图像左上角为红色、右上角为蓝色、左下角为绿色,右下角为白色。
    & i- M6 ~' l5 i9 @% X( n) U( C: S4 k
    代码:. o' K. I9 {' r* X
    2 @1 W' A! _( `; e* {% q' N( Z
    clear;clc;
      Q2 I6 c. Q1 e/ Q8 ]! W$ \picRed=zeros(128,128);      %初始化RGB红色分量为0矩阵
    - k( R- A* j& P) ipicRed(1:64,1:64)=1;
    7 J& T# ?% Y/ K4 s' v% W+ r4 n# VpicRed(65:128,65:128)=1;    %令左上角、右下角的红色分量值为1
    9 {5 Y* S/ a  I* dpicGreen=zeros(128,128);    %初始化RGB绿色分量为0矩阵
    ( m4 g) G6 C5 V# ~' |3 npicGreen(65:128,1:64)=1;# }+ V/ E. s0 i* M, d0 n5 E
    picGreen(65:128,65:128)=1;  %令左上角、右下角的绿色分量值为1. S0 l" ^, S# l, S& y
    picBlue=zeros(128,128);     %初始化RGB蓝色分量为0矩阵
    7 y; W  K/ _4 t6 S8 o/ n6 h! ypicBlue(1:64,65:128)=1;
    , _3 U; j- o! J5 ?picBlue(65:128,65:128)=1;   %令左上角、右下角的蓝色分量值为1/ W1 M# e* p8 {! b, z) `  J! x! [
    pic=cat(3,picRed,picGreen,picBlue); %合成RGB分量3 }8 s) [6 ~! |" e
    imshow(pic);title('生成图像');      %显示图像# [. |6 P( n2 B) ?
    1
    , g8 g. \  T  Z! I% f& g$ A2
    & r, B$ }% B5 R0 u! D$ V37 b5 I( {7 N, h8 x+ \/ h
    4
    : l0 N0 L( s! L" P* b5/ d) g2 O! r& ]' s0 ?
    6
    & O. H" H3 k3 E/ t6 M9 L. r  j  E7
    4 Z/ t4 Z! G2 i  A' X5 f# a8$ ^. Q0 q0 s, K0 k
    9
      E" v/ w0 d3 F/ J2 \! L( D10. W3 M* N0 e  F% g1 x3 u
    110 y! ~* q& ]% ~2 ?7 X) n( Y
    12
    7 \! z6 G# X% Q5 [: x0 k结果:
    9 q3 ^2 `9 G( F* ~! {8 E6 S, n  \$ m
    : h9 C; Y% V" a, L4 i& _( ^$ r分析:
    + R3 m+ b$ t2 s3 \  j% v! `RGB图像存在红,绿,蓝三个颜色的通道。每个像素点由这三个通道的强度合成决定颜色。根据光学理论,当红色分量为255,其他分量为0时,呈现出红色;当绿色分量为255,其他分量为0时,呈现出绿色;当蓝色分量为255,其他分量为0时,呈现出蓝色。当三通道的颜色分量全为255时,显示出的颜色变为白色。同理,当三通道分量全为0时会显示黑色。" B2 j. f% v" L3 F5 ^5 V8 \

    4 n2 r! ^7 u! X3 o  c" C二、数字图像的读取与显示。- Q% z( A5 Y8 l2 u* _' g" C3 W
    $ O3 A3 Q$ }3 f& n
    (1)图像读取——imread
    : {& Y# P2 L3 ^4 |/ Y5 o* b  K# N& I1 D" N5 B% g
    Matlab自带参考文档里imread的用法介绍:
    % b2 g. E7 [/ w- G
    $ E( d' |' l& @+ R' `, J  U3 ^: S8 Himread - Read image from graphics file
    # ?" h. p7 i2 M6 f
    & a# ^+ c' d$ u+ w    This MATLAB function reads a grayscale or color image from the file specified by
    9 v' ~# n& j: l0 k: Z    the string filename.* O) G$ C6 p; R9 q3 z

    / i3 N3 H) X3 M5 o: o    A = imread(filename, fmt)
    , }" I% l7 E! {7 ]0 M    [X, map] = imread(...)
    ) g7 ]& S/ ^9 v0 H- c    [...] = imread(filename)# Z9 _# ]2 J0 ^+ V5 t2 I+ u
        [...] = imread(URL,...)
    8 ^9 K) k& e: I+ x    [...] = imread(...,Param1,Val1,Param2,Val2...)- u4 l% o( m" ?% c& B' f
    1' F+ ]! K7 r8 x; E+ o, B5 c
    2
    & x$ [! H. J0 h9 d3
    % L: r- G. L# F  F4
    + K% A, ]" z1 ^+ `! W5 A5
    : H$ x8 P$ B3 S! \2 F6
    8 p0 v, N) [; r8 l  X9 e1 a8 u7" {! J* @( I' l. R( x3 Y! c, z% p
    8
    & R8 a8 G% w: r9
    2 l! M! `; }$ p' I2 F10
    7 w& ]' }+ ^9 U, a: [以其中一个常用的用法为例:
    ; Z! |% U, r/ F: `* r, E3 O" }0 w: q% K. a6 k0 z
    [...] = imread(filename)4 A7 B3 ^+ g& k- f, k! Q
    1$ g1 z( O7 U3 s  y" Q# Z
    其中filename为待读取图像文件名。若该图像不在matlab默认路径下并不与当前运行.m文件同目录时,需要手动添加路径名。
    , ?: O' _3 ?3 y# ?# Y! l! z2 }* T6 s9 x4 D2 ~, C
    实例:1 ^5 I7 J% y% Z! e" Q! H# E/ u

    0 f  `, n# h, f3 V. e# j>> pic=imread('E:\1.jpg')
    2 }  d1 ^- \5 E0 p& ^5 ~1( O! v5 t# d/ i# e
    该代码将E盘下的1.jpg文件读入了matlab,用变量pic存储。) E( |$ O8 l* Y' Y; V

    : s! `3 T9 i0 I1 g/ |8 q(2)图像显示——imshow" o! U5 l. l# m0 b- _

    6 \9 Z2 u9 S+ n& N- V) P3 L5 UMatlab自带参考文档里imshow的用法介绍:4 d4 c- Z) j  J  j& W4 k$ a
    ' A. W' d6 @- }4 Z. \2 C! H9 p% P
    imshow - Display image$ Z+ E& O: t* r  C- e% x; K4 o5 K

    % F  P* H  w2 H, P1 J- J2 g: e; Z  T    This MATLAB function displays the image I in a Handle Graphics figure, where I4 D+ T% H7 s: ^' G
        is a grayscale, RGB (truecolor), or binary image.
    . H# t+ f) s9 s. u1 s1 X5 i+ V7 i+ o% T2 O1 k& I6 l- v# s/ r( z
        imshow(I). D0 H' F. l* i. Y' G; ^7 J8 h& r
        imshow(I,RI)
    & @' O" R: h6 E) @8 B    imshow(X,map)
    : a9 f% @8 ?, l# g. q1 v) l    imshow(X,RX,map)+ Z; ?3 E% P) D0 d6 O0 R
        imshow(filename)
    6 [; m* _9 _! T; ]# m    imshow(___,Name,Value...)) B; _' A) d) X& q
        imshow(gpuarrayIM,___)
    3 W% `# Q3 y  C4 w- R    imshow(I,[low high])$ Y) z2 ^9 {1 ^7 W6 L* T3 q: P: }
    himage = imshow(___)
    ' E3 S# X* X8 f- }. `- F1/ e2 w2 V( Z) P& P6 k- L. B3 D
    2+ k* i- O! J4 d1 t/ M5 H& L
    3: M. m' [. \- @" X9 v) ~
    44 @) l" G% F, J) F
    59 H) H! Y& R9 {1 d0 D9 Q
    6
    ) Y* K& n9 J0 V. {8 _* @7+ g  u' ^8 h- x; s* @
    8& }) ?. [4 M' U
    9  I# b6 G; n: k- T
    10
    5 |* m1 C2 p( d) j% Y2 y- P9 v1 A- Q2 d11
    9 G: {9 U2 Y9 @6 s12  E5 K- P* X, R. g) Z9 [8 d3 ^
    13
    2 |# j8 o9 W" t" d6 q: b140 R8 r3 E3 O& N5 I4 G. b' K
    以其中一个常用的用法举例:3 W/ j$ e& n; v+ n) A' c

    * X) ~* |- _7 c  _% fimshow(I)
    , r/ T; A  k4 i& L/ w1% P0 s* }+ d0 `  d7 d) A
    其中,I为一个矩阵。该函数将把该矩阵所代表的图像显示出来。
    3 w- v8 U5 w9 t# u) @% y. n% k2 A1 C& p5 M; p5 `& |. w8 }4 C
    注:并不是只有读取的图片所生成的变量才可以显示出图像,一般而言自己创建的、符合要求的矩阵均可通过该函数显示出其代表的图像。
    9 a, z4 C% r; E  R- w: _1 A6 ]6 A1 u+ H% p+ V  `
    (3)批量获取图像
    3 j( \8 h7 e( o5 Z! n% a7 C- V  K' {; C  f0 p  Z
    在做题的时候,有时会需要读取一个文件夹里的所有图像。输入几十行imread的确是一种方法,不过非常低级。此时应使用批量读取图像的方法。
    : p; K% [7 |" l% I& Z5 ~! b: p9 x3 R4 J8 Z
    首先需学习dir函数,可以通过它获取当前文件夹下的文件名等信息。在此列出其几个常用用法:+ H1 Z8 O1 d$ n0 Y

    # T/ z' b: _8 S9 V    filepath=dir('.')       %列出当前目录下所有子文件夹和文件  F' Q/ I* P/ j5 y( ~2 N4 T
        filepath=dir('E:\matlab')   %列出指定目录下所有子文件夹和文件& _2 G( e! {  m0 e* e6 V
        filepath=dir('*.m') %列出当前目录下符合正则表达式的文件夹和文件2 B" R1 w9 ?. L. D3 R" F
    1. L; B1 D6 Z% i' l$ g
    2
    7 {% {' `/ w! ^  k7 `  t! `. I# n1 G3
    ! S+ V# g8 R- J变量filepath用于存放dir函数返回的文件信息。filepath变量不再是一个单纯的数值矩阵,而是一个面向对象的对象矩阵。每个对象包含name、date、bytes等属性。可通过filepath.name调用文件名。
    0 z& V- I4 L* z& |+ t' `3 r
    / L! [# M2 h: d# ]有关正则表达式的知识在此不详细叙述,上述dir函数第三个用法中的正则表达式*.m表示所有文件名以.m结尾的文件。- n/ p3 @, P" g! x9 F
    9 O+ h# I8 Y4 ^) [
    以本校2017年新锐赛题目为例,将文件夹中20个图像统一格式、分别命名为1.jpeg~20.jpeg后,可用如下代码读取:" K/ |( z# }6 C  `  Q+ ^0 g

    ; {8 g2 m& J; F5 w. ~5 bfor i=1:20! m7 I: B4 o* |% e
            imageName=strcat('附件1:英文纵切\',num2str(i-1),'.jpeg');
    & y2 H, E* v4 U: Q        pic(:,:,i)=rgb2gray(imread(imageName));) `: C$ W' T8 F6 O
    end3 u( s8 A7 U6 W/ D$ ?
    16 ]' ^: f' u: G. l
    2
    - |. K2 O, a% Q5 ?3
    ! v0 }2 @* G8 ?0 T9 I4
    8 r& X0 P/ \; ^3 M, X4 |, q, X' m此代码读取了文件夹中的20个数字图像文件并将其转化成了灰度图像存储于三维变量pic中。最高维数字代表第几张图。7 [; G; N$ q" \' L' G2 y2 I
    / ^9 @5 v9 f1 L
    若不事先手动统一名称,可用如下代码读取:
    . `( q( ^, B0 K+ R, A( E5 H0 m) {, m. I6 Q: Y
    filepath=dir(‘*.jpeg’);
    ' }8 L5 S6 q) Vfor i=1:size(filepath)      / j/ o  Z1 ^" k* T# ]* ^8 i
      filename=strcat('附件1:英文纵切\’,filepath(i).name);     9 ^9 [% T, Y! B' u8 \8 c
        pic(:,:,i)=rgb2gray(imread(filename));% a4 i" I; {( C; x$ S3 f
    end. ]4 Z' I# z8 o: g$ @+ z
    1) T' s4 H* d0 ^- Y, D' ]
    21 A  {- H. Y- `1 f& R) Q
    3
    , k, m3 f! M" c# G41 W3 v% D0 N/ L
    5
    ) d8 P  V3 S. x三、数字图像基础操作。
    8 w: {& x* e/ [1 b, x' i1 \! ~' A/ s
    (1)水平镜像
    $ B) O0 }8 y) T. p9 g4 A4 g& d" z5 e
    & T+ d" S! v& m6 n- ]6 s代码:
    5 O. k8 [, G$ I
    5 R  o2 _& D3 o% t2 [' OI=imread('1.jpg');      % 读取图像1.jpg
    4 S1 C3 z! Q/ Q/ y0 W6 S0 }& qsubplot(1,2,1);         % 画1行2列的子图的第一部分" ?% J% a5 v6 b. H. N9 }
    imshow(I);                  % 显示图像
    * F* C# @. y( @' e" ttitle('原图');            % 添加标题“原图”
    8 c; V/ I. H- D* j# R9 S, s! eI=double(I);                % 将I图像中灰度值转化为double类型
    * c* ^1 M/ h; A8 S2 N  A  cH=size(I);                  % 将I图像大小赋值给H
    9 A# @1 D9 ?$ P* h4 yI2(1:H(1,1),1:H(1,2))=I(1:H(1,1),H(1,2):-1:1);  % 行坐标不变,列坐标反转
    8 Q# G; O# J* Fsubplot(1,2,2);         % 画1行2列的子图的第二部分
    * Z5 F3 F# m; ~imshow(uint8(I2));      % 显示图像
    * n% z" s1 }4 R$ d& Htitle('水平镜像');      % 添加标题“水平镜像”; B+ s2 u0 [$ q4 ?3 h9 L9 o; T+ v
    1+ f  \; q0 M2 W' P3 D# ~) }
    2
    5 B3 P: }* {( N! H34 v+ a0 W" h- Y# Z& u2 n
    4
    . \6 V2 W/ Q0 g2 u, g. `58 n: N4 A! o/ ~/ v) x
    6& L% _& F0 H9 W* X; d! g+ G
    72 I3 K6 b0 l5 a  V+ b% ?; v
    8
    % a- ~8 n" M$ M7 n( [0 b( K* Z9 ?9! C9 O# l+ `- c/ E! T
    10' d& F8 k+ y7 c" y% q
    结果: & I1 ?' Z  |) O; K* J: U. R4 A1 x/ G
    4 g# ]( S5 x' \3 u
    (2)大小变换——imresize
    - \; }. A* A5 R: ^7 u
    ( R8 [0 L/ @7 n) {1 }newPic=imresize(pic,[200,200])
    % V* J% L9 z- `. s9 `4 K1) o# e6 i. Z# h' s
    此代码将存放于matlab中的图像变量转化成了200*200大小像素点的新图像,并存放于变量newPic中。
    $ I( U% y# u8 f4 B* j, l2 U2 e5 r, w: y# Z0 B0 H1 a, M
    (3)旋转——imrotate& `% t) F  _6 M5 V- X, M2 i6 O  F

    ) P% L+ z) u5 y/ q8 ~6 W. x/ g( IB=imrotate(A,angle,method,’loose’);2 Y2 G& ^9 T) c2 X
       B=imrotate(A,angle,method,’crop’);. i% n/ Q& @( U- c4 ]
       解释:angle为旋转的角度(默认逆时针旋转)
    " T9 M9 b1 ~, H( _         method为插值方式1 G1 c) V" W0 l
             loose:当图片旋转一定角度之后,若使得图片超出范围大小,不剪掉超出部分" i9 p4 r# x) A0 M# J
             crop:当图片旋转一定角度之后,若使得图片超出范围大小,剪掉超出部分
    # w5 M/ X! X8 i3 f: y- ?13 d4 ^* v, o. i2 D# j
    2
    + M2 {8 M6 c8 u* V3
    " o. Y3 e' Z8 m: `4/ b$ m, G# E1 H, k
    5
    * Y% F3 [3 m  K, k6 A3 G: {68 r  G3 ~; T6 K7 u1 Z9 `5 T
    举例:
    1 C* ?( J: E$ }: q
    6 P8 a6 l& V* W$ F' sclear;" z7 P( ~0 ?- S& y5 r' L0 q
    pic=imread('1.jpg');
    # v% b8 y( k( tpic_new1=imrotate(pic,45,'nearest','crop');
    & R& \# \& f4 @pic_new2=imrotate(pic,45,'nearest','loose');
    " z/ h5 f( ~" j0 b7 f; y% |. M9 ~subplot(1,3,1);                            ! ]/ U  _, s1 p# ]# \6 W
    imshow(pic);% j+ s0 R- O' v3 r) m
    title('原图');
    3 S! W( J/ _; [! e$ Q' wsubplot(1,3,2);
    : E" c( m6 E3 x2 q1 x- aimshow(pic_new1);   
    ' M5 D8 ~0 Y& Q3 C2 k$ btitle('crop');
    1 \2 h/ f* E, v8 w8 a5 \subplot(1,3,3);
    ' O6 X, m# O4 p6 g5 X* f3 Wimshow(pic_new2);+ ~, v' C" J+ R2 Z) A. Y; T
    title('loose');        
    2 S5 I: ?; n$ M& R1
    6 {  K  e& w" O6 z5 }24 x/ K8 A5 g, X" O; O3 p2 s" U
    3
    ' i7 Q. y  ~: ]0 J$ I  L4) B) K) a3 d8 s+ u+ D: Y; a/ D* s& M
    5
    7 H7 o, B2 z" f/ `; [5 [0 \6
    7 T' T! I7 ^" w7 T7
    9 h( m) \5 G+ T2 s! V% Y. v8) e& l0 I9 S$ B) q
    9
    9 j5 \7 U5 ?7 j10
    : g* m7 E+ C" l9 X* i2 i11
    ; c8 u6 M8 X! E6 ~6 c+ ~% Y# K12% t4 k" }" m  G  f9 d; H0 ^
    13
    ' ?, T; N. b" O9 ]) @  V7 x. ^运行结果: , v' X- ]. d( O* M8 r: H5 \5 q

    # y8 w! v) v( H5 a4 `& \除此之外,数字图像还有多种滤波、频域变换、角点检测等知识,设计数字图像处理专业课知识且在数学建模中用的不多,在此不介绍。; a  v' X0 ]9 U' J, t
    --------------------- . f% Y* [0 o" x9 }6 m- P* z! O
    作者:bit_xzh $ L) [$ z5 ?. u
    来源:CSDN 0 [" o& x) x5 p/ h  }1 D
    原文:https://blog.csdn.net/GregoryHanson/article/details/77802804 % H2 Z' {2 o) T6 n2 i3 o: S

    2 ^6 W0 V# i$ Q4 r$ v
    - R% u0 W2 s8 g! n; {) ]/ Q3 X
    , p/ b5 }2 j8 D

    数字图像处理.docx

    15.82 KB, 下载次数: 3, 下载积分: 体力 -2 点

    售价: 5 点体力  [记录]

    zan
    转播转播0 分享淘帖0 分享分享0 收藏收藏1 支持支持0 反对反对0 微信微信

    47

    主题

    3

    听众

    317

    积分

    升级  5.67%

  • TA的每日心情
    开心
    2020-4-28 18:55
  • 签到天数: 39 天

    [LV.5]常住居民I

    回复

    使用道具 举报

    您需要登录后才可以回帖 登录 | 注册地址

    qq
    收缩
    • 电话咨询

    • 04714969085
    fastpost

    关于我们| 联系我们| 诚征英才| 对外合作| 产品服务| QQ

    手机版|Archiver| |繁體中文 手机客户端  

    蒙公网安备 15010502000194号

    Powered by Discuz! X2.5   © 2001-2013 数学建模网-数学中国 ( 蒙ICP备14002410号-3 蒙BBS备-0002号 )     论坛法律顾问:王兆丰

    GMT+8, 2026-4-15 14:37 , Processed in 1.333803 second(s), 60 queries .

    回顶部