Imdilate de matlab en C

[LordofDestruction]

Buenas tardes amig@s, tengo unas dudas puntuales, resulta que tengo que realizar un programa que lea un arreglo bidimensional (matriz) con sólo elementos 1 y 0, tal como el tratamiento de una imagen binaria. Y a esta imagen, se le trate como matlab dilatandola con elemento estructurante. Ahí va mi duda.
Tengo lo siguiente:

Código:

#include <stdio.h>
int k,l;
short int matriz[18][42] =
{0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,0,0,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,0,0,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,0,0,0,0,0,1,0,0,0,0,0,0,0,0,0,1,1,1,0,1,0,1,1,1,0,0,0,0,0,0,0,0,
0,0,0,1,1,1,1,1,1,0,0,0,0,0,1,1,1,0,0,0,0,0,0,0,0,1,1,1,0,1,0,1,1,1,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,0,0,0,0,0,0,1,0,1,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,0,0,0,0,0,1,1,1,0,0,0,0,0,0,0,1,1,1,0,0,0,0,0,1,1,1,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,0,0,0,1,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,
0,0,0,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,1,0,1,0,1,0,1,0,1,0,1,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,1,1,1,1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,
0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0};


int a,b;        // Variables que almacenan el la dimensión de la matriz por teclado
int x,y;
int i,j;        // Variables utilizadas para diversas estructuras
int d[18][42];  // Matriz auxiliar usada en Imdilate
int reiniciar;  // Da la opción de reiniciar el programa para usarlo de nuevo.
int aux[20][50];
/*void copiar_x_en_d(){

for(i=0; i<18; i++){
for(j=0; j<42; j++){
d[i][j]=matriz[i][j];
}}}*/
void main(){

   for(i=0;i<18;i++)
        for (j=0;j<42;j++)
           for(k=0;k<3;k++)
                for(l=0;l<3;j++){
                    if(matriz[i][j]==1){
                        d[i][j+1]=1;
                        }
                    }
//copiar_x_en_d();
/*
for(k=1;k<=20;k++)
    for(l=1;l<=20;l++)
        for(i=k-1;i<=k+1;i++)
            for(j=l-1;j<=l+1;j++)
                if(matriz[i][j]==1)
                    aux[k-1][l-1]=1;
for(i=0;i<20;i++)
    for(j=0;j<20;j++)
        d[i][j]=aux[i][j];
*/

printf("\n\n");
system("cls");
printf("La Matriz (Imagen) Dilatada de %d por %d es: \n",18,42);
printf("\n\n");

for (i=0; i<18; i++){
for (j=0; j<42; j++)

printf("%d",d[i][j]);
printf("\n");
}
printf("\n");
printf("\n\n\n");

}

y nada :(
Quien tiene alguna sugerencia?