/* ==================================== */
int32_t lreconsdilat(
struct xvimage *g,
struct xvimage *f,
int32_t connex)
/* reconstruction de g sous f */
/* g : image marqueur */
/* f : image masque */
/* resultat dans g */
/* ==================================== */
{
return lgeodilat(g, f, connex, -1);
}
Image<T> Geodilation(Image<T> &G, Image<T> &R, int connex, int niter)
{
Image<T> Geodilat(G.Dimx(), G.Dimy(), G.Dimz());
// Pink Images
struct xvimage* imageG;
struct xvimage* imageR;
struct xvimage* temp;
int32_t typepixel;
if (sizeof(T)==1)
typepixel = VFF_TYP_1_BYTE;
else if (sizeof(T)==2)
typepixel = VFF_TYP_2_BYTE;
else if (sizeof(T)==4)
typepixel = VFF_TYP_4_BYTE;
else
std::cerr<<"Error in Geodilation : ImageType not known"<<std::endl;
imageG=allocheader(NULL,G.Dimx(),G.Dimy(),G.Dimz(),typepixel);
imageG->image_data= G.GetPointer();
imageR=allocheader(NULL,G.Dimx(),G.Dimy(),G.Dimz(),typepixel);
imageR->image_data= R.GetPointer();
temp=copyimage(imageG);
lgeodilat(temp,imageR,connex,niter);
for (int z = 0; z<G.Dimz() ; ++z){
for (int y = 0; y<G.Dimy() ; ++y){
for (int x = 0; x<G.Dimx(); ++x){
Geodilat(x, y, z) = ((T *)(temp->image_data))[x + y * G.Dimx() + z * G.Dimx() * G.Dimy()];
}
}
}
free(imageR);
free(imageG);
free(temp);
return Geodilat;
}
