0、下定决心

当当入手数字图像处理一本，从此开此正式跨入数字图像处理大门。以前虽然多多少少接触过这些东西，也做过一些相关的事情，但感觉都不够系统，也不够专业，从今天开始，一步一步地学习下去，相信会有成果的。

1、本章知识点结构

下面的思维导图主要用xmind制作。知识点是对第1章的一些个人认为比较重要的知识概括了一下，然后加了自己的一点点理解。

2、书中例子实践

第一个例子，试一试手。

 function f = twodsin(A, u0, v0, M, N)

 f = zeros(M, N);

 for c = :N

     v0y = v0 * (c - );

     for r = :M

         u0x = u0 * (r - );

         f(r, c) = A * sin(u0x + v0y);

     end

 end

运行速度，只要0.0154秒，但是书上的用了0.0471秒，看来我的电脑比作者的电脑好了不少哈，或者是高版本的MATLAB(2014a)，优化做得更好。

>> timeit(@() twodsin(, /(*pi), /(*pi), , ))

ans =

    0.0154

将第一个例子的第2行注释掉，运行如下，果然速度还是下降了不少，所以还是预分配内存比较好。

>> timeit(@() twodsin(, /(*pi), /(*pi), , ))

ans =

    0.0263

就是一个二维的sin(x,y)函数，它的幅度就表示亮度，如下图所示，黑白正弦起伏，感觉很不错的，MATLAB果然很强大。

>> f=twodsin(, /(*pi), /(*pi), , );

>> imshow(f,[]);

当然此例子的关键不在画出图像，而是要通过向量化的方式将上述函数中的两个for循环去掉，从而提高算法的速度。

 function f = twodsin2(A, u0, v0, M, N)

 r = :M-;

 c = :N-;

 [C, R] = meshgrid(c, r);

 f = A * sin(u0 * R + v0 * C);

同样的参数，只用了0.0072，比上个0.0154，快了50%多。

>> timeit(@() twodsin2(, /(*pi), /(*pi), , ))

ans =

    0.0072

对meshgrid(M,N)函数，我是直接把它理解为生个一个网格，这个网格的就是MxN的面积大小，而向量M和N中的元素，就是这个网格上的点，可以用如下图来表示，更加直观一点，如理解有误，欢迎指正。

上述网格表示的是如下meshgrid，有了图，应该非常容易理解。

>> [X,Y]=meshgrid(:,:)

X =

Y =

总结一下

通过对本章的学习，对图像的理解，更加深入，可以从数学的角度来看一幅图像，同时，也更加清楚以明白图像的信息是如何存储的。这一章中，也讲子许多MATLAB的技巧，很实用，也深深体会到了MATLAB处理矩阵的强大功能，这一优势确实是C/C++/C#等语言无法比拟的。路漫漫，吾必将上下求索。Keep on！！

本博文纯原创，未经博主同意（我的邮箱：endless@139.com），请勿转载。有想交流学习的小伙伴，可以加我的QQ：1449494575，欢迎各位。