varmapToImage.hpp

Go to the documentation of this file.

/** \file varmapToImage.hpp
 * \author Jared R. Males (jaredmales@gmail.com)
 * \brief A utility to convert a wavefront variance map to an intensity image
 * \ingroup mxAO_files
 *
 */
 
#ifndef varmapToImage_hpp
#define varmapToImage_hpp
 
namespace mx
{
namespace AO
{
namespace analysis
{
 
/// Convert a wavefront variance map to an intensity image by convolving with the PSF.
/**
 * The PSF should have odd dimensions, and the peak pixel should be in the center pixel
 * defined by [0.5*psf.rows(), 0.5*psf.cols()].  The platescale (lambda/D per pixel) of the PSF should match that of
 * the variance map. The PSF should be normalized such that the center/peak pixel has value 1 (not so that the sum of
 * pixels is 1).
 *
 * \tparam imageT is an Eigen-array-like type.
 *
 * \ingroup mxAOAnalytic
 */
template <typename imageT>
void varmapToImage( imageT &im,     ///< [out] is the intensity image, resized to match varmap
                    imageT &varmap, ///< [in]  is the wavefront variance map
                    imageT &psf     ///< [in]  is the point spread function
)
{
    typedef typename imageT::Scalar floatT;
 
    im.resize( varmap.rows(), varmap.cols() );
    im.setZero();
 
    floatT psfNorm = 1.0; // psf.sum();
 
    int psf_i, psf_j;
    floatT psfVal;
 
#pragma omp parallel for
    for( int i = 0; i < im.rows(); ++i )
    {
        for( int j = 0; j < im.cols(); ++j )
        {
            for( int ii = 0; ii < im.rows(); ++ii )
            {
                psf_i = 0.5 * ( psf.rows() ) - ( ii - i );
 
                if( psf_i < 0 || psf_i >= psf.rows() )
                    continue;
 
                for( int jj = 0; jj < im.cols(); ++jj )
                {
                    psf_j = 0.5 * ( psf.cols() ) - ( jj - j );
 
                    if( psf_j >= 0 and psf_j < psf.cols() )
                    {
                        psfVal = psf( psf_i, psf_j ) / psfNorm;
                        im( i, j ) += varmap( ii, jj ) * psfVal;
                    }
                }
            }
        }
    }
}
 
} // namespace analysis
} // namespace AO
} // namespace mx
 
#endif // varmapToImage_hpp