3
3
Program: Insight Segmentation & Registration Toolkit
4
4
Module: $RCSfile: itkRGBGibbsPriorFilter.txx,v $
6
Date: $Date: 2005/08/10 18:13:43 $
7
Version: $Revision: 1.46 $
6
Date: $Date: 2008-01-27 18:29:24 $
7
Version: $Revision: 1.47 $
9
9
Copyright (c) Insight Software Consortium. All rights reserved.
10
10
See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.
309
309
bool changeflag = (k > 3);
311
for(unsigned int j = 0; j < 2; j++)
311
for(unsigned int jj = 0; jj < 2; jj++)
314
energy[j] += GibbsEnergy(i, 0, j);
315
energy[j] += GibbsEnergy(i+rowsize+1, 1, j);
316
energy[j] += GibbsEnergy(i+rowsize, 2, j);
317
energy[j] += GibbsEnergy(i+rowsize-1, 3, j);
318
energy[j] += GibbsEnergy(i-1, 4, j);
319
energy[j] += GibbsEnergy(i-rowsize-1, 5, j);
320
energy[j] += GibbsEnergy(i-rowsize, 6, j);
321
energy[j] += GibbsEnergy(i-rowsize+1, 7, j);
322
energy[j] += GibbsEnergy(i+1, 8, j);
323
if ( m_LabelStatus[i] == j )
314
energy[jj] += GibbsEnergy(i, 0, jj);
315
energy[jj] += GibbsEnergy(i+rowsize+1, 1, jj);
316
energy[jj] += GibbsEnergy(i+rowsize, 2, jj);
317
energy[jj] += GibbsEnergy(i+rowsize-1, 3, jj);
318
energy[jj] += GibbsEnergy(i-1, 4, jj);
319
energy[jj] += GibbsEnergy(i-rowsize-1, 5, jj);
320
energy[jj] += GibbsEnergy(i-rowsize, 6, jj);
321
energy[jj] += GibbsEnergy(i-rowsize+1, 7, jj);
322
energy[jj] += GibbsEnergy(i+1, 8, jj);
323
if ( m_LabelStatus[i] == jj )
642
642
ChangedPixelVec = inputImageIt.Get();
645
const std::vector<double> & dist =
645
const std::vector<double> & distValue =
646
646
m_ClassifierPtr->GetPixelMembershipValue( ChangedPixelVec );
648
648
LabelType pixLabel;
649
if (dist[1] > m_ObjectThreshold) pixLabel = 0;
649
if (distValue[1] > m_ObjectThreshold) pixLabel = 0;
650
650
else pixLabel = 1;
651
651
labelledImageIt.Set( pixLabel );