OpenCV实现Matlab的circshift、fftshift、ifftshift函数

最新推荐文章于 2024-07-03 10:52:17 发布

TerryBryant

最新推荐文章于 2024-07-03 10:52:17 发布

阅读量6.1k

点赞数 4

分类专栏：图像处理文章标签： opencv matlab

本文链接：https://blog.csdn.net/dageda1991/article/details/78129205

版权

本文介绍了如何使用OpenCV实现Matlab中的circshift、fftshift和ifftshift函数。通过分析Matlab的源码，作者发现这三个函数都可以归结为circshift的特例。文章详细讲解了circshift的实现，并提供了fftshift和ifftshift的代码示例。最后，作者指出，虽然实现与Matlab的输入格式有所不同，但经过验证，结果与Matlab一致，且适用于一维Mat。欢迎大家在GitHub上交流改进意见。

摘要由CSDN通过智能技术生成

今天挺高兴的，折腾了一天终于搞定了这三个函数，而且在我搜索相关资料的时候，发现不管是百度还是谷歌，都没有人将这三个函数用OpenCV同时实现。感觉挺惊讶的，因为在Matlab代码转OpenCV的时候，特别当用到傅里叶变换时，这类函数应该会很频繁的用到。先不管这些了，这里斗胆来一个全网首发，哈哈。

要实现Matlab的这三个函数，首先可以参考matlab是怎么做的，当在Matlab的命令行输入type circshift.m，可以看到该函数的实现被隐藏了，再看看fftshift和ifftshift，可以发现两者都是调用的circshift。其实理解这三个函数的原理后，就知道fftshift和ifftshift都是circshift的两种特殊情况，所以我们首先实现circshift。

一、circshift的实现

circshift表示的是循环移动，详细说明可以自行百度。它在matlab中的用法是circshift(A,[1 2])，其中A表示待移动矩阵，[1 2]表示矩阵循环向下移动一行，再循环向右移动两行，如下图
这里写图片描述

看起来很简单，但是写起来并不方便，要是用指针操作这些数据，还是很容易写错的。考虑到既然是用OpenCV，那就用Mat自带的函数来实现。这里用到的是copyTo函数，具体思想看代码就能明白：

void circshift(Mat &out, const Point &delta)
{
    Size sz = out.size();

    // 错误检查
    assert(sz.height > 0 && sz.width > 0);
    // 此种情况不需要移动
    if ((sz.height == 1