C / C++ / MFC

<br />
#include <br />
#include <br />
#include <br />
using namespace std;<br />
int main (void) {<br />
    FILE * in_file;FILE * out_file;<br />
	<br />
	int image[190][245];<br />
<br />
	<br />
    int read_size, write_size;<br />
    int image_size=46175;<br />
  <br />
    in_file = fopen("leaf.raw", "rb");<br />
    if(in_file == NULL) {<br />
               fprintf(stderr,"Could not open file \n");<br />
               exit(1);<br />
    }<br />
    <br />
    read_size = fread(image, 1, image_size, in_file);<br />
    if (read_size != image_size) {<br />
       fprintf(stderr, "Reading error");<br />
       exit(2);<br />
    }<br />
    <br />
    out_file = fopen("newleaf.raw", "wb");<br />
    if(out_file == NULL) {<br />
               fprintf(stderr,"Could not open file \n");<br />
               exit(3);<br />
    }<br />
    <br />
    write_size = fwrite(image, 1, image_size, out_file);<br />
    if (write_size != image_size) {<br />
       fprintf(stderr, "Reading error");<br />
       exit(4);<br />
    }<br />
<br />
	<br />
<br />
 //************************** convolution with Gaussian**********************<br />
	<br />
	int mask [3][3];<br />
	int IconvG [190][245];<br />
	<br />
	mask[0][0]=1;mask[0][1]=2;mask[0][2]=1;<br />
	mask[1][0]=2;mask[1][1]=3;mask[1][2]=2;<br />
	mask[2][0]=1;mask[2][1]=2;mask[2][2]=1;<br />
<br />
	int m,n,u,v,sum,i,j;<br />
	int kCols =3, kRows =3, rows = 190, cols = 245 ,mm,nn,ii,jj;<br />
int kCenterX = kCols / 2;<br />
int kCenterY = kRows / 2;<br />
<br />
for(i=0; i < rows; ++i)              // rows<br />
{<br />
    for(j=0; j < cols; ++j)          // columns<br />
    {<br />
        sum = 0;                     // init to 0 before sum<br />
<br />
        for(m=0; m < kRows; ++m)     // kernel rows<br />
        {<br />
            mm = kRows - 1 - m;      // row index of flipped kernel<br />
<br />
            for(n=0; n < kCols; ++n) // kernel columns<br />
            {<br />
                nn = kCols - 1 - n;  // column index of flipped kernel<br />
<br />
                // index of input signal, used for checking boundary<br />
                ii = i + (m - kCenterY);<br />
                jj = j + (n - kCenterX);<br />
<br />
                // ignore input samples which are out of bound<br />
                if( ii >= 0 && ii < rows && jj >= 0 && jj < cols )<br />
                IconvG[i][j] += image[ii][jj] * mask[mm][nn];<br />
            }<br />
        }<br />
    }<br />
}<br />
<br />
 out_file = fopen("newleaf.raw", "wb");<br />
    if(out_file == NULL) {<br />
               fprintf(stderr,"Could not open file \n");<br />
               exit(3);<br />
    }<br />
    <br />
    write_size = fwrite(IconvG, 1, image_size, out_file);<br />
    if (write_size != image_size) {<br />
       fprintf(stderr, "Reading error");<br />
       exit(4);<br />
    }<br />
	<br />
	<br />
	cout<<"done\n";<br />
    fclose(in_file);<br />
    fclose(out_file);<br />
    return 0;<br />
}<br />