forked from USTC-Resource/USTC-Course
-
Notifications
You must be signed in to change notification settings - Fork 1
Commit
This commit does not belong to any branch on this repository, and may belong to a fork outside of the repository.
- Loading branch information
Showing
1 changed file
with
116 additions
and
123 deletions.
There are no files selected for viewing
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| Original file line number | Diff line number | Diff line change |
|---|---|---|
| @@ -1,123 +1,116 @@ | ||
| <!-- style=flat/plastic/social, label force to change the title --> | ||
| # <div align="center">2019 计算机图像学实验 </div> | ||
| [](LICENSE) | ||
| [](.) | ||
|
|
||
| >说明. 最开始我用的 cpp 实现了实验内容(代码在 `cpp`目录下），以及自己实现了 fft, 1d,2d。但是配置 windows上的 opencv 环境失败，代码只经过了静态语法检测，可能还有些地方有 bug。 后来我用的 python 重新实现了除 FFT 的所有算法，并将结果记录如下 | ||
| ## 使用 | ||
| ### 环境 | ||
| - python3.6+ | ||
| - matplotlib | ||
| - numpy | ||
| - cv2 | ||
|
|
||
| ### 运行 | ||
| 当前目录下的所有 python 代码按如下格式执行 | ||
| `python3 lab*.py <IMG_PATH>` | ||
| 如 | ||
| `python3 lab1.py images/lena.bmp` | ||
|
|
||
|
|
||
| # 实验内容 | ||
| ## 图像的点处理 | ||
| ### 灰度的线性变换 | ||
| 输入斜率，截距， 进行一维线性变换 | ||
|
|
||
| ### 灰度拉伸 | ||
| 输入两个转折点(x1,y1),(x2,y2), 进行分段的线性变换 | ||
| ``` | ||
| 当 x<x1, f(x) = y1/x1*x | ||
| 当 x1<=x<=x2, f(x) = (y2-y1)(x-x1)/(x2-x1)+y1 | ||
| 当 x>x2, f(x) = (255-y2)*(x-x2)/(255-x2)+y2 | ||
| ``` | ||
|
|
||
| ### 灰度直方图 | ||
| 输入图像，显示它的灰度直方图, 还可以输入恢复的上限，下限， 显示这个范围内的灰度直方图 | ||
|
|
||
| ### 直方图均衡 | ||
| 扩大灰度范围，减少灰度之间的数量差值 | ||
|
|
||
|
|
||
| 结果如下 | ||
|  | ||
|
|
||
|  | ||
|
|
||
| ## 数字图像的平滑 | ||
| 滤波，去除图像的噪声,均值滤波，中值滤波去除加性噪声；同态滤波去除乘性噪声 | ||
|
|
||
| 先给图像加上 3% 的椒盐噪声，然后分别使用 窗口大小为 3 的均值滤波器和中值滤波器进行滤波 | ||
|
|
||
|
|
||
| 记图像大小 nxm, 窗口大小为 wxw | ||
| 我在实现滤波器时，移动窗口，每次只会更新移进的值，和移出的值。 | ||
| 即 窗口先又移动，每移动一列，就将这列的数据考虑进来，而将移出的那一列剔除。这样在更新窗口的值时只需 O(w) | ||
|
|
||
| 对于均值滤波，时间复杂 O(nmw) | ||
| 对于中值滤波，需要求出中值，则这需要 O(w^2) 的时间才能完成。可以利用 快速选择的算法，在O(lengthOfArray) 时间里找出排任意名次的数，这里找出中值, w\*w/2 | ||
| 总时间复杂度 O(nmw^2) | ||
|
|
||
| 结果如下 | ||
|  | ||
|
|
||
|
|
||
| ## 图像的边缘检测 | ||
| 实验原理： | ||
|
|
||
| 在灰度图像的情况下，所谓的边缘检测可以看成是基于图像像素灰度值在空间的不连续性对图像做出的一种分割。边缘可以用方向和幅度两个特性来描述。一般而言，沿边缘走向方向其幅度值变化较平缓，而沿垂直于边缘走向其幅度值变化较剧烈。 | ||
|
|
||
| 经典的边缘提取方法是考察图像的每个像素在某个邻域内灰度的变化，利用边缘邻近一阶或二阶方向导数变化规律，用简单的方法检测边缘。这种方法称为边缘检测局部算子 | ||
| 法。 | ||
|
|
||
| 边缘检测算子一般有 | ||
| - Roberts 交叉算子 | ||
| - Sobel 模板卷积 | ||
| - Prewitt 同上 | ||
| - Laplace | ||
|
|
||
| 我实现的时 Roberts, Prewitt， 结果如下 | ||
|  | ||
|  | ||
|
|
||
| ## 傅里叶变换 | ||
| 有时间可以专门写一篇介绍 Fourier 变换 | ||
|
|
||
| 实验要求 | ||
| - 对 `images/rect*` 图像作二维Fourier , 显示频谱，然后作幅度变换，将低频移到中心点 | ||
| - Fourier 反变换 幅度，并显示 | ||
| - Fourier 反变换 相位，并显示 | ||
|
|
||
| 结果如下 | ||
|  | ||
|  | ||
|
|
||
|
|
||
| ## 实验总结 | ||
| - 使用 python 的话，img 像素值类型为 `uint8`， 直接相加减可能造成 溢出， 结果在 mod 256 的域中, 所以在可能出现溢出的情况下，我在前面用 `0+` 后面的结果，这样可以将类型提升为 int 的运算而不会出现溢出。 | ||
| - 在实现滤波器的,更新窗口的时候，我只更新变换的列， 只需要 O(w) 的复杂度，如果更新整个窗口，需要 O(w^2) 的复杂度。 | ||
| - 对于 c++ 实现的 快速傅里叶变换，接口定义如下 | ||
| ```c++ | ||
| typedef complex<double> comp ; | ||
|
|
||
| class dft | ||
| { | ||
| public: | ||
| dft(); | ||
| ~dft(); | ||
| bool dft1d(vector<comp>&, vector<comp> const &); | ||
| bool dft2d(vector<comp>&, vector<comp> const &); | ||
| bool idft1d(vector<comp>&, vector<comp> const &); | ||
| bool dft::_dft2d(vector<vector<comp>>& dst, vector<vector<comp>> const &src,bool isInvert=false) | ||
| bool dft::dft2d(vector<vector<comp>>& dst, vector<vector<comp>> const &src) | ||
| bool dft::idft2d(vector<vector<comp>>& dst, vector<vector<comp>> const &src) | ||
| }; | ||
| ``` | ||
| 实现的思路是： | ||
| 首先实现 一维的 dft, idft。 使用 快速傅里叶算法 fft, | ||
| 对每一层， 计算倒序数，进行计算，一个 log(n) 层，每一层计算 n次， 则一维 fft时间复杂度为 `O(nlog(n))` | ||
| 然后利用傅里叶变换的可分离性，计算二维 傅里叶变换2d dft | ||
| 可以先对每行进行 1d dft, 然后对每列进行 1d dft | ||
| 逆变换同理。 | ||
|
|
||
|
|
||
| # <div align="center">2019 计算机图像学实验</div> | ||
| >说明. 最开始我用的 cpp 实现了实验内容(代码在 `cpp`目录下），以及自己实现了 fft, 1d,2d。但是配置 windows上的 opencv 环境失败，代码只经过了静态语法检测，可能还有些地方有 bug。 后来我用的 python 重新实现了除 FFT 的所有算法，并将结果记录如下 | ||
| ## 1.1. 使用 | ||
| ### 1.1.1. 环境 | ||
| - python3.6+ | ||
| - matplotlib | ||
| - numpy | ||
| - cv2 | ||
|
|
||
| ### 1.1.2. 运行 | ||
| 当前目录下的所有 python 代码按如下格式执行 | ||
| `python3 lab*.py <IMG_PATH>` | ||
| 如 | ||
| `python3 lab1.py images/lena.bmp` | ||
|
|
||
|
|
||
| <div align="center"><h1>实验内容</h1> </div> | ||
| ## 2.1. 图像的点处理 | ||
| ### 2.1.1. 灰度的线性变换 | ||
| 输入斜率，截距， 进行一维线性变换 | ||
|
|
||
| ### 2.1.2. 灰度拉伸 | ||
| 输入两个转折点(x1,y1),(x2,y2), 进行分段的线性变换 | ||
| ``` | ||
| 当 x<x1, f(x) = y1/x1*x | ||
| 当 x1<=x<=x2, f(x) = (y2-y1)(x-x1)/(x2-x1)+y1 | ||
| 当 x>x2, f(x) = (255-y2)*(x-x2)/(255-x2)+y2 | ||
| ``` | ||
|
|
||
| ### 2.1.3. 灰度直方图 | ||
| 输入图像，显示它的灰度直方图, 还可以输入恢复的上限，下限， 显示这个范围内的灰度直方图 | ||
|
|
||
| ### 2.1.4. 直方图均衡 | ||
| 扩大灰度范围，减少灰度之间的数量差值 | ||
|
|
||
|
|
||
| 结果如下 | ||
|  | ||
|
|
||
|  | ||
|
|
||
| ## 2.2. 数字图像的平滑 | ||
| 滤波，去除图像的噪声,均值滤波，中值滤波去除加性噪声；同态滤波去除乘性噪声 | ||
|
|
||
| 先给图像加上 3% 的椒盐噪声，然后分别使用 窗口大小为 3 的均值滤波器和中值滤波器进行滤波 | ||
|
|
||
|
|
||
| 记图像大小 nxm, 窗口大小为 wxw | ||
| 我在实现滤波器时，移动窗口，每次只会更新移进的值，和移出的值。 | ||
| 即 窗口先又移动，每移动一列，就将这列的数据考虑进来，而将移出的那一列剔除。**这样在更新窗口的值时只需 O(w)**。 而如果直接更新整个窗口，需要 O(w\*w) | ||
|
|
||
| 对于均值滤波，时间复杂 O(nmw) | ||
| 对于中值滤波，需要求出中值，则这需要 O(w^2) 的时间才能完成。可以利用 快速选择的算法，在O(lengthOfArray) 时间里找出排任意名次的数，这里找出中值, w\*w/2 | ||
| 总时间复杂度 O(nmw^2) | ||
|
|
||
| 结果如下 | ||
|  | ||
|
|
||
|
|
||
| ## 2.3. 图像的边缘检测 | ||
| 实验原理： | ||
|
|
||
| 在灰度图像的情况下，所谓的边缘检测可以看成是基于图像像素灰度值在空间的不连续性对图像做出的一种分割。边缘可以用方向和幅度两个特性来描述。一般而言，沿边缘走向方向其幅度值变化较平缓，而沿垂直于边缘走向其幅度值变化较剧烈。 | ||
|
|
||
| 经典的边缘提取方法是考察图像的每个像素在某个邻域内灰度的变化，利用边缘邻近一阶或二阶方向导数变化规律，用简单的方法检测边缘。这种方法称为边缘检测局部算子 | ||
| 法。 | ||
|
|
||
| 边缘检测算子一般有 | ||
| - Roberts 交叉算子 | ||
| - Sobel 模板卷积 | ||
| - Prewitt 同上 | ||
| - Laplace | ||
|
|
||
| 对于具体实现，需要注意的是 像素值的类型，以及参与的运算。例如 python 中，像素值类型为 `uint8`， 直接相加减可能造成 溢出， 结果在 mod 256 的域中, 所以在可能出现溢出的情况下，我在前面用 `0+` 后面的结果，这样可以将类型提升为 int 的运算而不会出现溢出。 | ||
|
|
||
|
|
||
| 我实现了 Roberts, Prewitt， 结果如下 | ||
|  | ||
|  | ||
|
|
||
| ## 2.4. 傅里叶变换 | ||
|
|
||
| 实验要求 | ||
| - 对 `images/rect*` 图像作二维Fourier , 显示频谱，然后作幅度变换，将低频移到中心点 | ||
| - Fourier 反变换 幅度，并显示 | ||
| - Fourier 反变换 相位，并显示 | ||
|
|
||
| - 对于 c++ 实现的 快速傅里叶变换，接口定义如下 | ||
| ```c++ | ||
| typedef complex<double> comp ; | ||
|
|
||
| class dft | ||
| { | ||
| public: | ||
| dft(); | ||
| ~dft(); | ||
| bool dft1d(vector<comp>&, vector<comp> const &); | ||
| bool dft2d(vector<comp>&, vector<comp> const &); | ||
| bool idft1d(vector<comp>&, vector<comp> const &); | ||
| bool dft::_dft2d(vector<vector<comp>>& dst, vector<vector<comp>> const &src,bool isInvert=false) | ||
| bool dft::dft2d(vector<vector<comp>>& dst, vector<vector<comp>> const &src) | ||
| bool dft::idft2d(vector<vector<comp>>& dst, vector<vector<comp>> const &src) | ||
| }; | ||
| ``` | ||
| 实现的思路是： | ||
| - 首先实现 一维的变换 dft1d, idft1d。 | ||
|
|
||
| - 使用 快速傅里叶算法 fft,对每一层， 计算倒序数，进行计算，一个 log(n) 层，每一层计算 n次， 则一维 fft时间复杂度为 `O(nlog(n))` | ||
|
|
||
| - 然后利用傅里叶变换的可分离性，计算二维 傅里叶变换dft2d, idft2d: 先对每行进行一维变换, 然后对每列进行一维变换。 | ||
|
|
||
| 结果如下 | ||
|  | ||
|  | ||
|
|