图像处理中常用的正交变换除了傅里叶变换外,还有其他一些有用的正交变换,其中离散余弦就是一种。离散余弦变换表示为DCT(DiscreteCosineTransformation),常用于图像处理和图像识别等。
一维离散余弦变换
正变换
(1)
(2)
式中F(u)是第u个余弦变换系数,u是广义频率变量,u=1,2,3......N-1;f(x)是时域N点序列,x=0,1,2......N-1
反变换
(3)
显然,式(1)式(2)和式(3)构成了一维离散余弦变换对。
二维离散余弦变换
正变换
(4)
式(4)是正变换公式。其中f(x,y)是空间域二维向量之元素,x,y=0,1,2,......N-1;F(u,v)是变换系数阵列之元素。式中表示的阵列为N×N
反变换
(5)
式中的符号意义同正变换式一样。式(4)和式(5)是离散余弦变换的解析式定义。
矩阵表示法
更为简洁的定义方法是采用矩阵式定义。根据以上公式定义可知,离散余弦变换的系数矩阵可以写成如下:
如果令N=4,那么由一维解析式定义可得如下展开式。
写成矩阵式
若定义F(u)为变换矩阵,A为变换系数矩阵,f(x)为时域数据矩阵,则一维离散余弦变换的矩阵定义式可写成如下形式
[F(u)]=[A][f(x)](6)
同理,可得到反变换展开式
写成矩阵式即
[f(x)]=[A]T[F(u)] (7)
二维离散余弦变换也可以写成矩阵式:
[F(u,v)]=[A][f(x,y)][A]T (8)
[f(x,y)]=[A]T[F(u,v)][A]
式中[f(x,y)]是空间数据阵列,A是变换系数阵列,[F(u,v)]是变换矩阵,[A]T是[A]的转置。
对二维图像进行离散余弦变换
由以上对二维离散余弦变换的定义及公式(7)可知,求二维图像的离散余弦变换要进行以下步骤:
1.获得图像的二维数据矩阵f(x,y);
2.求离散余弦变换的系数矩阵[A];
3.求系数矩阵对应的转置矩阵[A]T;
4.根据公式(7)[F(u,v)]=[A][f(x,y)][A]T计算离散余弦变换;
源代码:
package cn.edu.jxau.image;
import java.awt.image.BufferedImage;
/**
* 图像的变换
* @author luoweifu
*
*/
public class Transformation {
/**
* 要进行DCT变换的图片的宽或高
*/
public static final int N = 256;
/**
* 傅里叶变换
* @return
*/
public int[] FFT() {
return null;
}
/**
* 离散余弦变换
* @param pix 原图像的数据矩阵
* @param n 原图像(n*n)的高或宽
* @return 变换后的矩阵数组
*/
public int[] DCT(int[] pix, int n) {
double[][] iMatrix = new double[n][n];
for(int i=0; i<n; i++) {
for(int j=0; j<n; j++) {
iMatrix[i][j] = (double)(pix[i*n + j]);
}
}
double[][] quotient = coefficient(n); //求系数矩阵
double[][] quotientT = transposingMatrix(quotient, n); //转置系数矩阵
double[][] temp = new double[n][n];
temp = matrixMultiply(quotient, iMatrix, n);
iMatrix = matrixMultiply(temp, quotientT, n);
int newpix[] = new int[n*n];
for(int i=0; i<n; i++) {
for(int j=0; j<n; j++) {
newpix[i*n + j] = (int)iMatrix[i][j];
}
}
return newpix;
}
/**
* 矩阵转置
* @param matrix 原矩阵
* @param n 矩阵(n*n)的高或宽
* @return 转置后的矩阵
*/
private double[][] transposingMatrix(double[][] matrix, int n) {
double nMatrix[][] = new double[n][n];
for(int i=0; i<n; i++) {
for(int j=0; j<n; j++) {
nMatrix[i][j] = matrix[j][i];
}
}
return nMatrix;
}
/**
* 求离散余弦变换的系数矩阵
* @param n n*n矩阵的大小
* @return 系数矩阵
*/
private double[][] coefficient(int n) {
double[][] coeff = new double[n][n];
double sqrt = 1.0/Math.sqrt(n);
for(int i=0; i<n; i++) {
coeff[0][i] = sqrt;
}
for(int i=1; i<n; i++) {
for(int j=0; j<n; j++) {
coeff[i][j] = Math.sqrt(2.0/n) * Math.cos(i*Math.PI*(j+0.5)/(double)n);
}
}
return coeff;
}
/**
* 矩阵相乘
* @param A 矩阵A
* @param B 矩阵B
* @param n 矩阵的大小n*n
* @return 结果矩阵
*/
private double[][] matrixMultiply(double[][] A, double[][] B, int n) {
double nMatrix[][] = new double[n][n];
double t = 0.0;
for(int i=0; i<n; i++) {
for(int j=0; j<n; j++) {
t = 0;
for(int k=0; k<n; k++) {
t += A[i][k]*B[k][j];
}
nMatrix[i][j] = t; }
}
return nMatrix;
}
}
分享到:
相关推荐
图像的正交变换源码 傅里叶变换 离散余弦变换 沃尔什变换
不用多说了,经典算法。初学者和需要交作业的学生可以参考。完整的源代码和可执行程序。
visual c++数字图像处理典型算法及实例,离散余弦变换源码,已经通过测试
【图像加密】基于DCT(离散余弦变换)与SVD(奇异值分解)域自适应嵌入水印matlab源码.md
dct傅里叶离散余弦变换图片水印+java+opencv+idea源码包
【信号隐藏】基于小波变换算法DWT和离散余弦变换DCT变换音频数字水印嵌入提取matlab源码.zip
DCT离散余弦变换的图像融合的MATLAB仿真_源码
DCT C# 源代码 对媒体技术 离散余弦变换
【信号隐藏】基于离散余弦变换DCT-SVD变换音频数字水印嵌入提取matlab源码.zip
1、资源内容:基于Matlab使用离散余弦变换进行图片压缩(源码+图片).rar 2、适用人群:计算机,电子信息工程、数学等专业的大学生课程设计、期末大作业或毕业设计,作为“参考资料”使用。 3、解压说明:本资源需要...
多媒体技术原理应用于MATLAB的DCT变化代码示例
拥有多通道图显示功能,可显示傅里叶变换图,直方图,离散余弦变换图。 1、该资源内项目代码经过严格调试,下载即用确保可以运行! 2、该资源适合计算机相关专业(如计科、人工智能、大数据、数学、电子信息等)正在做...
【图像加密】基于DCT(离散余弦变换)与SVD(奇异值分解)域自适应嵌入水印matlab源码.zip
压缩传感重构算法:稀疏基采用离散余弦变换的omp重构算法(Compressed sensing reconstruction algorithm: sparse basis using the discrete cosine transform reconstruction algorithm omp)
1.版本:matlab2019a,不会运行可私信 2.领域:【图像重构】 3.内容:基于3D 离散余弦变换实现图像的快速压缩重构附matlab代码.zip 4.适合人群:本科,硕士等教研学习使用
离散余弦算法的数字水印matlab源码
资源名:PCD 压缩传感_使用并行坐标下降算法_冗余字典由各种小波变换和离散余弦变换被用于图像修复_matlab 资源类型:matlab项目全套源码 源码说明: 全部项目源码都是经过测试校正后百分百成功运行的,如果您下载后...
Matlab领域上传的视频均有对应的完整代码,皆可运行,亲测可用,适合小白; 1、代码压缩包内容 主函数:main.m; 调用函数:其他m文件;无需运行 运行结果效果图; 2、代码运行版本 ...若运行有误,根据提示修改;...
Matlab研究室上传的视频均有对应的完整代码,皆可运行,亲测可用,适合小白; 1、代码压缩包内容 主函数:main.m; 调用函数:其他m文件;无需运行 运行结果效果图; 2、代码运行版本 ...若运行有误,根据提示修改;...
CSDN海神之光上传的代码均可运行,亲测可用,直接替换数据即可,适合小白; 1、代码压缩包内容 ...DSIFT多聚焦图像融合、加权平均法图像融合、泊松彩色图像融合、主成分结合小波离散变换PCA-DWT图像融合