离散正弦变换(Discrete SineTransform, DST):
类似于DCT,离散正弦变换也有八种形式,常用的是DST-Ⅰ,一般我们提到DST而不特殊说明即指DST-Ⅰ。
这里直接给出DST-Ⅰ的定义:
正变换:
逆变换:
因此可得到变换矩阵SN:
变换矩阵SN是一个正交矩阵,而且它还是一个对称矩阵(转置与本身相等),所以正变换与反变换是相同的。
基于OpenCV离散正弦变换:
void MainWindow::on_pushButton_10_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)(sin(pi*i*j/(N+1))*pow(2/(N+1),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();
}
效果图: