12,444,517 members (52,706 online)
alternative version

48K views
23 bookmarked
Posted

# DCT Implementation in C# (DCT of Image)

, 16 Nov 2009 CPOL
 Rate this:
Implementing DCT on Image (2D DCT) in C# and its Inverse

## Introduction

A discrete cosine transform (DCT) expresses a sequence of finitely many data points in terms of a sum of cosine functions oscillating at different frequencies. DCTs are important to numerous applications in science and engineering, from lossy compression of audio and images (where small high-frequency components can be discarded), to spectral methods for the numerical solution of partial differential equations. The use of cosine rather than sine functions is critical in these applications: for compression, it turns out that cosine functions are much more efficient (as explained below, fewer are needed to approximate a typical signal), whereas for differential equations the cosines express a particular choice of boundary conditions.

In particular, a DCT is a Fourier-related transform similar to the discrete Fourier transform (DFT), but using only real numbers. DCTs are equivalent to DFTs of roughly twice the length, operating on real data with even symmetry (since the Fourier transform of a real and even function is real and even), where in some variants the input and/or output data are shifted by half a sample. There are eight standard DCT variants, of which four are common.

## Background

Here I am showing an implementation of 2D DCT. We generate the 2D DCT kernel first and then using matrix operations we perform DCT. I am very lazy to write the theory. Please refer to the proper books for the same.

## Using the Code

Here we perform the following activities:

1. Select the image, scale the image so that it fits the picture box, this scaling is shown in Scaling percentage box.
2. Select a specific part of the image, we select image part of dimensions multiple of 2 (128*128, 256*256).
3. Read the image in 2D Array, we are using pointer arithmetic for this, we have to allow pointer operations in project properties.
{
int i, j;
GreyImage = new int[Width, Height];  //[Row,Column]
Input = new double [Width, Height];  //[Row,Column]
Bitmap image = Obj;
BitmapData bitmapData1 = image.LockBits(new Rectangle
(0, 0, image.Width, image.Height),
PixelFormat.Format32bppArgb);
unsafe
{
byte* imagePointer1 = (byte*)bitmapData1.Scan0;

for (i = 0; i < bitmapData1.Height; i++)
{
for (j = 0; j < bitmapData1.Width; j++)
{
GreyImage[j, i] = (int)((imagePointer1[0] +
imagePointer1[1] + imagePointer1[2]) / 3.0);
Input [j,i]=(double)GreyImage[j,i];
//4 bytes per pixel
imagePointer1 += 4;
}//end for j
//4 bytes per pixel
imagePointer1 += bitmapData1.Stride - (bitmapData1.Width * 4);
}//end for i
}//end unsafe
image.UnlockBits(bitmapData1);
return;
}
4. Create Object of FastDCT2D class, which implements the DCT.
5. To find DCT, first we have to generate 2D DCT kernel and then implement multiplication:
public double[,] GenerateDCTmatrix(int order)
{
int i, j;
int N;
N = order;
double alpha;
double denominator;
double[,] DCTCoeff = new double[N, N];
for (j = 0; j <= N - 1; j++)
{
DCTCoeff[0, j] = Math.Sqrt(1 / (double)N);
}
alpha = Math.Sqrt(2 / (double)N);
denominator = (double)2 * N;
for (j = 0; j <= N - 1; j++)
for (i = 1; i <= N - 1; i++)
{
DCTCoeff[i, j] = alpha * Math.Cos(((2 * j + 1) *
i * 3.14159) / denominator);
}

return (DCTCoeff);
}
public void FastDCT()
{
double[,] temp = new double[Width, Height];
DCTCoefficients = new double[Width, Height];
DCTkernel = new double[Width, Height];
DCTkernel = GenerateDCTmatrix(Order);
temp = multiply(DCTkernel, Input);
DCTCoefficients = multiply(temp, Transpose(DCTkernel));
DCTPlotGenerate();
return;
}

## Points of Interest

DCT plot is also generated using Dynamic Range Compression.

## History

This is my first attempt. I have also added a slow version of DCT in the code, which is given by another author. Any changes or suggestions are welcome.

## You may also be interested in...

 Pro Pro

 First Prev Next
 How get hash of image? Member 112666334-Nov-15 0:45 Member 11266633 4-Nov-15 0:45
 Obviously incomplete John Ktejik5-Feb-14 7:07 John Ktejik 5-Feb-14 7:07
 Error in Displaymap j1s1e123-May-12 23:06 j1s1e1 23-May-12 23:06
 My vote of 2 Member 826082213-Dec-11 10:38 Member 8260822 13-Dec-11 10:38
 My vote of 4 piaoxuele21-Jul-10 16:36 piaoxuele 21-Jul-10 16:36
 Last Visit: 31-Dec-99 18:00     Last Update: 23-Aug-16 15:51 Refresh 1