class point(int x, int y){
	int x;
	int y;
	}

std:vector <point> v;

class FlatSE()
{

FlatSE(int r){new FlatSE; this.makeBall(r);}

void makeBall(int r){
	for (int x=-r;x=<r;x++)
	{	
		for (int y=-r;y=<r;y++)
		{	
			if (x*x+y*y<r*r) v.push_back(new point(x,y);
		}
	}
}

}

/*
 * creation element structurant joe voisinage 4-connexe
 */

joe = new FlatSE(0);
joe.push_back((-1,0));
joe.push_back((1,0));
joe.push_back((0,-1));
joe.push_back((0,1));

CImg<unsigned char> erode(CImg<unsigned char> im, FlatSE se)
{
/*
 * erode doit pour chaque point calculer la valeur min parmi les points
 * correspondant a FlatSE se.
 */
}

CImg<unsigned char> dilate(CImg<unsigned char> im, FlatSE se)