ShapeAnalysis.cpp

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.
 */

#include <ImLib3D/ShapeAnalysis.hpp>
#include <ImLib3D/Draw.hpp>
#include <ImLib3D/Display.hpp>
#include <ImLib3D/TestPatterns.hpp>


void
IP3D::BoundingBox(const Mask3D& mask, RectZone3Di &bbox)
{
    //      RectZone3Di box(_Base::Height(), _Base::Width(), _Base::Depth(), 0, 0, 0);
    bbox.SetP0(Vect3Di(mask.Height(), mask.Width(), mask.Depth()));
    bbox.SetP1(Vect3Di(0, 0, 0));
    Mask3D::const_iteratorXYZ p;
    bool noValues = true;
    for (p=mask.begin(); p!=mask.end(); p++)
    {
        if ((*p)>0)
        {
            if (p.x < bbox.x0)
            bbox.x0 = p.x;
            if (p.y < bbox.y0)
            bbox.y0 = p.y;
            if (p.z < bbox.z0)
            bbox.z0 = p.z;
            noValues = false;
        }
    }
    int n = mask.GetNVoxels();
    p=mask.begin(); p+=n-1;//No p=end() assignment possible at the moment
    for (; p!=mask.begin(); p--)
    {
        if ((*p)>0)
        {
            if (p.x > bbox.x1)
            bbox.x1 = p.x;
            if (p.y > bbox.y1)
            bbox.y1 = p.y;
            if (p.z > bbox.z1)
            bbox.z1 = p.z;
            noValues = false;
        }
    }
    if (noValues)
    {
        bbox.x0 = 0; bbox.y0 = 0; bbox.z0 = 0;
    }
}

---