Image3Dlinear.hxx

Go to the documentation of this file.

/* ImLib3D
 * Copyright (c) 2001, ULP-IPB Strasbourg.
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or (at
 * your option) any later version.
 * 
 * This program is distributed in the hope that it will be useful, but
 * WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
 * General Public License for more details.
 * 
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
 */

#ifndef _Image3Dlinear_hxx
#define _Image3Dlinear_hxx

#include<ImLib3D/Interpolator3D.hpp>

// *******************************************************
// ********************** Image3Dlinear ********************
// *******************************************************
template<class Im3DValue>
Image3Dlinear<Im3DValue>::~Image3Dlinear()
{
    if(destroyInterpolator && interpolator){delete interpolator;}
    this->DeAllocate(); 
}

// *************************** Accessors ******* ***********
template<class Im3DValue>
Im3DValue Image3Dlinear<Im3DValue>::SafeValue(int pos) const 
{
    if ((pos < 0) || (pos >= this->GetNVoxels()))
        return zero;
    else return (*this)(pos);
}

template<class Im3DValue>
void
Image3Dlinear<Im3DValue>::ComplexImag(Image3Df &imag) const
{
    imag.Resize(this->Size());

    typename _Self::const_iteratorFast p;
    Image3Df::iteratorFast pImag;
    for (p = this->begin(),pImag=imag.begin(); p != this->end(); p++,pImag++)
    {
        *pImag=(*p).imag();
    }
}
template<class Im3DValue>
void
Image3Dlinear<Im3DValue>::ComplexReal(Image3Df &imag) const
{
    imag.Resize(this->Size());

    typename _Self::const_iteratorFast p;
    Image3Df::iteratorFast pImag;
    for (p = this->begin(),pImag=imag.begin(); p != this->end(); p++,pImag++)
    {
        *pImag=(*p).real();
    }
}
template<class Im3DValue>
void
Image3Dlinear<Im3DValue>::ComplexNorm2(Image3Df &imag) const
{
    imag.Resize(this->Size());

    typename _Self::const_iteratorFast p;
    Image3Df::iteratorFast pImag;
    for (p = this->begin(),pImag=imag.begin(); p != this->end(); p++,pImag++)
    {
        *pImag=norm(*p);
    }
}
template<class Im3DValue>
void
Image3Dlinear<Im3DValue>::ComplexModule(Image3Df &imag) const
{
    imag.Resize(this->Size());

    typename _Self::const_iteratorFast p;
    Image3Df::iteratorFast pImag;
    for (p = this->begin(),pImag=imag.begin(); p != this->end(); p++,pImag++)
    {
        *pImag=sqrt(norm(*p));
    }
}

template<class Im3DValue>
double
Image3Dlinear<Im3DValue>::L2Norm() const
{
    double energie = 0;
    typename _Self::const_iteratorFastMasked p;
    for (p=this->begin(); p!=this->end(); p++)
    {
        energie += (*p)*(*p);
    }
    return sqrt(energie);
}

template<class Im3DValue>
void
Image3Dlinear<Im3DValue>::GetComponent(int c,Image3Df &res) const 
{
    res.Resize(this->Size());
    typename _Self::const_iteratorFast p;
    Image3Df::iteratorFast rp;
    for(p=this->begin(),rp=res.begin();p!=this->end(); p++,rp++)
    {
        *rp=(*p)(c);
    }       
}
template<class Im3DValue>
void
Image3Dlinear<Im3DValue>::SetComponent(int c,const Image3Df &src)
{
    src.Resize(this->Size());
    typename _Self::iteratorFast p;
    Image3Df::const_iteratorFast rp;
    for(p=this->begin(),rp=src.begin();p!=this->end(); p++,rp++)
    {
        (*p)(c)=*rp;
    }       
}

#endif // _Image3Dlinear_hxx

---