基于OpenCV的离散哈特莱变换

最新推荐文章于 2023-12-31 15:10:26 发布

执剑00天涯

最新推荐文章于 2023-12-31 15:10:26 发布

阅读量1.6k

点赞数 1

CC 4.0 BY-SA版权

分类专栏： OpenCV图像处理文章标签： OpenCV 离散哈特莱变换

本文链接：https://blog.youkuaiyun.com/qq_35759050/article/details/72312811

OpenCV图像处理专栏收录该内容

18 篇文章

订阅专栏

本文详细介绍了离散哈特莱变换(DHT)的基本原理及其应用，通过数学公式展示了DHT的正交变换特性，并提供了实现DHT的C++代码示例。此外，还分析了DHT的四种形式及其实现过程。

离散哈特莱变换(Discrete Hartley Transform, DHT)：

DHT公式给成正交变换的形式：

其中：

上面定义式中直接将定标因子写为1/N的开方，这样保证DHT的变换矩阵HN直接是正交矩阵：

观察发现矩阵HN不但是正交矩阵并且还是对称矩阵（转置与本身相等），因此逆矩阵与本身相等（正交矩阵的逆矩阵等于其转置），即DHT的逆变换与正变换形式相同。

DHT也有四种形式，上面给出的是第一种形式，其它参见维基百科【1】。

模仿离散余弦变换，正变换与逆变换相同:

下面是代码：

void MainWindow::on_pushButton_8_clicked()
{
    Mat img = imread("D:/Pic/2.jpg",0);

    Mat src(img.rows,img.cols,CV_64FC1);
    img.convertTo(src,CV_64FC1);

    Mat dst(src.rows,src.cols,CV_64FC1);

    double pi = 3.141592657;
    Mat tmp(200,200,CV_64FC1);
    int N = 200;
    for(int i=0;i<200;i++)
        for(int j=0;j<200;j++)
        {
            tmp.at<double>(j,i) = (double)(cos(2*pi*i*j/N)+sin(2*pi*i*j/N))*pow(1/N,1/2);
        }

    dst = tmp*src*tmp/45000;
    imshow("dst",dst);
    Plot* plot = new Plot(dst);
    //QGridLayout *grid = new QGridLayout(ui->frame);
    //grid->addWidget(plot,0,0);

    //plot->setTitle("3D");
    plot->coordinates()->axes[X1].setLabelString("X(um)");  //只能写在这
    plot->coordinates()->axes[Y1].setLabelString("Y(um)");  //设置坐标轴标签
    plot->coordinates()->axes[Z1].setLabelString("Z(nm)");
    plot->coordinates()->axes[Z2].setLabelString("Z(nm)");
    plot->coordinates()->axes[Z3].setLabelString("Z(nm)");
    plot->coordinates()->axes[Z4].setLabelString("Z(nm)");

    double start,stop;  //设置颜色条
    plot->coordinates()->axes[Z1].limits(start,stop);
    plot->legend()->setLimits(start,stop);
//    plot->legend()->setAutoScale(true);
    plot->resize(700,600);
    plot->show();
    //grid->~QGridLayout();
}

效果图：