//2018.1.19_21:09 by Cooper Liu
//Questions? Contact me: [email protected]
//本程序为C++工程，编译环境需要配置openCV
#include <iostream>
#include <opencv2/core/core.hpp>
#include <opencv2/highgui/highgui.hpp>
//#include<opencv2/opencv.hpp>//如果要包含所有库
using namespace cv;
//将mat类型转换为二维数组的类型
uchar** transform(Mat &img, int row, int col)
{
int i = 0, j = 0;
uchar **ptr = (uchar **)malloc(row*sizeof(uchar *));//二维数组ptr[][]
for (i = 0; i < row; i++)
ptr[i] = (uchar *)malloc(col*sizeof(uchar));
//新图的大小以及指向它的指针
Mat img2 = Mat(row, col, CV_8UC1);//图像img2：row*col大小
uchar *ptmp = NULL;//这是关键的指针！！
for (i = 0; i < row; i++)
{
ptmp = img2.ptr<uchar>(i);//指针指向img2的第i行
for (j = 0; j < col; j++)
{
ptr[i][j] = img.at<uchar>(i, j);//img的矩阵数据传给二维数组ptr[][]
ptmp[j] = ptr[i][j];//二维数组数据传给img2的第i行第j列
}
}
// namedWindow("新图");
// imshow("新图", img2);
// 等待100000 ms后窗口自动关闭
waitKey(100000);
return ptr;
}
uchar** mydct（uchar **ptr2,int M,int N）
{
int p, q, m, n;
double ap, aq;
uchar output[M][N];
for (p = 0; p < M - 1; p++)
{
for (q = 0; q < N - 1; q++)
{
if（p == 0）
{
ap = sqrt(1.0 / M);
}
else
{
ap = sqrt(2.0 / M);
}
if（q == 0）
{
aq = sqrt(1.0 / N);
}
else
{
aq = sqrt(2.0 / N);
}
double tmp = 0.0;
for (m = 0; m < M - 1; m++)
{
for (n = 0; n < N - 1;n++)
{
tmp+= ptr2[m][n] * cos((2 * m + 1)*p*PI / (2.0 * M)) * cos((2 * n + 1)*q*PI / (2.0 * N));
}
}
output[p][q] = ap*aq*temp;
}
}
}
//将二维数组转化成mat类型以图片形式输出
uchar** arraytomat（uchar **output,int M,int N）
{
int i, j;
//定义mat类型存放新图片
Mat img2 = Mat(M, N, CV_8UC1);//图像img2：row*col大小
uchar *ptmp = NULL;//这是关键的指针！！
for (i = 0; i < M; i++)
{
ptmp = img2.ptr<uchar>(i);//指针指向img2的第i行
for (j = 0; j < N; j++)
{
//output[i][j] = img.at<uchar>(i, j);//img的矩阵数据传给二维数组ptr[][] 此步骤已经在transform中存入了二维数组
ptmp[j] = output[i][j];//二维数组数据传给img2的第i行第j列
}
namedWindow("新图");
imshow("新图", img2);
}
}
void main()
{
//这一段程序实现：图像img -> ptr[][] -> 图像img2 的传递过程，即图像->二维数组->图像
Mat img = imread("pic.png", CV_LOAD_IMAGE_GRAYSCALE);//读取图像img。0表示转换为灰度图像读入
namedWindow("原图");
imshow("原图", img);
int row = img.rows;
int col = img.cols;
uchar **ptr2 = transform(img, img.rows, img.cols);
//对转换成的矩阵作dct变换
uchar **result=mydct（ptr2, row, col）;
arraytomat(result, row, col);
}修改版本：
2018/8/2 13:42
//2018.1.19_21:09 by Cooper Liu
//Questions? Contact me: [email protected]
//本程序为C++工程，编译环境需要配置openCV
#include <iostream>
#include <opencv2/core/core.hpp>
#include <opencv2/highgui/highgui.hpp>
//#include<opencv2/opencv.hpp>//如果要包含所有库
using namespace cv;
//将mat类型转换为二维数组的类型
uchar** transform(Mat &img, int row, int col)
{
int i = 0, j = 0;
uchar **ptr = (uchar **)malloc(row*sizeof(uchar *));//二维数组ptr[][]
for (i = 0; i < row; i++)
ptr[i] = (uchar *)malloc(col*sizeof(uchar));
//新图的大小以及指向它的指针
Mat img2 = Mat(row, col, CV_8UC1);//图像img2：row*col大小
uchar *ptmp = NULL;//这是关键的指针！！
for (i = 0; i < row; i++)
{
ptmp = img2.ptr<uchar>(i);//指针指向img2的第i行
for (j = 0; j < col; j++)
{
ptr[i][j] = img.at<uchar>(i, j);//img的矩阵数据传给二维数组ptr[][]
ptmp[j] = ptr[i][j];//二维数组数据传给img2的第i行第j列
}
}
// namedWindow("新图");
// imshow("新图", img2);
// 等待100000 ms后窗口自动关闭
waitKey(100000);
return ptr;
}
uchar** mydct（uchar **ptr2,int M,int N）
{
int p, q, m, n;
double ap, aq;
uchar output[M][N];
for (p = 0; p < M - 1; p++)
{
for (q = 0; q < N - 1; q++)
{
if（p == 0）
{
ap = sqrt(1.0 / M);
}
else
{
ap = sqrt(2.0 / M);
}
if（q == 0）
{
aq = sqrt(1.0 / N);
}
else
{
aq = sqrt(2.0 / N);
}
double tmp = 0.0;
for (m = 0; m < M - 1; m++)
{
for (n = 0; n < N - 1;n++)
{
tmp+= ptr2[m][n] * cos((2 * m + 1)*p*PI / (2.0 * M)) * cos((2 * n + 1)*q*PI / (2.0 * N));
}
}
output[p][q] = ap*aq*temp;
}
}
}
//将二维数组转化成mat类型以图片形式输出
uchar** arraytomat（uchar **output,int M,int N）
{
int i, j;
//定义mat类型存放新图片
Mat img2 = Mat(M, N, CV_8UC1);//图像img2：row*col大小
uchar *ptmp = NULL;//这是关键的指针！！
for (i = 0; i < M; i++)
{
ptmp = img2.ptr<uchar>(i);//指针指向img2的第i行
for (j = 0; j < N; j++)
{
//output[i][j] = img.at<uchar>(i, j);//img的矩阵数据传给二维数组ptr[][] 此步骤已经在transform中存入了二维数组
ptmp[j] = output[i][j];//二维数组数据传给img2的第i行第j列
}
namedWindow("新图");
imshow("新图", img2);
}
}
void main()
{
//这一段程序实现：图像img -> ptr[][] -> 图像img2 的传递过程，即图像->二维数组->图像
Mat img = imread("pic.png", CV_LOAD_IMAGE_GRAYSCALE);//读取图像img。0表示转换为灰度图像读入
namedWindow("原图");
imshow("原图", img);
int row = img.rows;
int col = img.cols;
uchar **ptr2 = transform(img, img.rows, img.cols);
//对转换成的矩阵作dct变换
uchar **result=mydct（ptr2, row, col）;
arraytomat(result, row, col);
}
G
M
T
| Detect languageAfrikaansAlbanianAmharicArabicArmenianAzerbaijaniBasqueBelarusianBengaliBosnianBulgarianCatalanCebuanoChichewaChinese (Simplified)Chinese (Traditional)CorsicanCroatianCzechDanishDutchEnglishEsperantoEstonianFilipinoFinnishFrenchFrisianGalicianGeorgianGermanGreekGujaratiHaitian CreoleHausaHawaiianHebrewHindiHmongHungarianIcelandicIgboIndonesianIrishItalianJapaneseJavaneseKannadaKazakhKhmerKoreanKurdishKyrgyzLaoLatinLatvianLithuanianLuxembourgishMacedonianMalagasyMalayMalayalamMalteseMaoriMarathiMongolianMyanmar (Burmese)NepaliNorwegianPashtoPersianPolishPortuguesePunjabiRomanianRussianSamoanScots GaelicSerbianSesothoShonaSindhiSinhalaSlovakSlovenianSomaliSpanishSundaneseSwahiliSwedishTajikTamilTeluguThaiTurkishUkrainianUrduUzbekVietnameseWelshXhosaYiddishYorubaZulu | AfrikaansAlbanianAmharicArabicArmenianAzerbaijaniBasqueBelarusianBengaliBosnianBulgarianCatalanCebuanoChichewaChinese (Simplified)Chinese (Traditional)CorsicanCroatianCzechDanishDutchEnglishEsperantoEstonianFilipinoFinnishFrenchFrisianGalicianGeorgianGermanGreekGujaratiHaitian CreoleHausaHawaiianHebrewHindiHmongHungarianIcelandicIgboIndonesianIrishItalianJapaneseJavaneseKannadaKazakhKhmerKoreanKurdishKyrgyzLaoLatinLatvianLithuanianLuxembourgishMacedonianMalagasyMalayMalayalamMalteseMaoriMarathiMongolianMyanmar (Burmese)NepaliNorwegianPashtoPersianPolishPortuguesePunjabiRomanianRussianSamoanScots GaelicSerbianSesothoShonaSindhiSinhalaSlovakSlovenianSomaliSpanishSundaneseSwahiliSwedishTajikTamilTeluguThaiTurkishUkrainianUrduUzbekVietnameseWelshXhosaYiddishYorubaZulu |
Text-to-speech function is limited to 200 characters
| Options : History : Feedback : Donate | Close |