<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.0 Transitional//EN"> <html><head><meta name="robots" content="noindex"> <meta http-equiv="Content-Type" content="text/html;charset=iso-8859-1"> <title>ITK: itkExpectationMaximizationMixtureModelEstimator.h Source File</title> <link href="DoxygenStyle.css" rel="stylesheet" type="text/css"> </head><body bgcolor="#ffffff"> <!-- Section customized for INSIGHT : Tue Jul 17 01:02:45 2001 --> <center> <a href="index.html" class="qindex">Main Page</a> <a href="modules.html" class="qindex">Groups</a> <a href="namespaces.html" class="qindex">Namespace List</a> <a href="hierarchy.html" class="qindex">Class Hierarchy</a> <a href="classes.html" class="qindex">Alphabetical List</a> <a href="annotated.html" class="qindex">Compound List</a> <a href="files.html" class="qindex">File List</a> <a href="namespacemembers.html" class="qindex">Namespace Members</a> <a href="functions.html" class="qindex">Compound Members</a> <a href="globals.html" class="qindex">File Members</a> <a href="pages.html" class="qindex">Concepts</a></center> <!-- Generated by Doxygen 1.5.9 --> <div class="navpath"><a class="el" href="dir_02e23d841c198393a6e961e9a0678fff.html">Code</a> » <a class="el" href="dir_2c33022d19e8077162952ea9ed57656d.html">Numerics</a> » <a class="el" href="dir_e82fe55e9fab4124d6625d716e091a5f.html">Statistics</a> </div> <div class="contents"> <h1>itkExpectationMaximizationMixtureModelEstimator.h</h1><a href="itkExpectationMaximizationMixtureModelEstimator_8h.html">Go to the documentation of this file.</a><div class="fragment"><pre class="fragment"><a name="l00001"></a>00001 <span class="comment">/*=========================================================================</span> <a name="l00002"></a>00002 <span class="comment"></span> <a name="l00003"></a>00003 <span class="comment"> Program: Insight Segmentation & Registration Toolkit</span> <a name="l00004"></a>00004 <span class="comment"> Module: $RCSfile: itkExpectationMaximizationMixtureModelEstimator.h,v $</span> <a name="l00005"></a>00005 <span class="comment"> Language: C++</span> <a name="l00006"></a>00006 <span class="comment"> Date: $Date: 2005-07-26 15:54:55 $</span> <a name="l00007"></a>00007 <span class="comment"> Version: $Revision: 1.12 $</span> <a name="l00008"></a>00008 <span class="comment"></span> <a name="l00009"></a>00009 <span class="comment"> Copyright (c) Insight Software Consortium. All rights reserved.</span> <a name="l00010"></a>00010 <span class="comment"> See ITKCopyright.txt or http://www.itk.org/HTML/Copyright.htm for details.</span> <a name="l00011"></a>00011 <span class="comment"></span> <a name="l00012"></a>00012 <span class="comment"> This software is distributed WITHOUT ANY WARRANTY; without even </span> <a name="l00013"></a>00013 <span class="comment"> the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR </span> <a name="l00014"></a>00014 <span class="comment"> PURPOSE. See the above copyright notices for more information.</span> <a name="l00015"></a>00015 <span class="comment"></span> <a name="l00016"></a>00016 <span class="comment">=========================================================================*/</span> <a name="l00017"></a>00017 <span class="preprocessor">#ifndef __itkExpectationMaximizationMixtureModelEstimator_h</span> <a name="l00018"></a>00018 <span class="preprocessor"></span><span class="preprocessor">#define __itkExpectationMaximizationMixtureModelEstimator_h</span> <a name="l00019"></a>00019 <span class="preprocessor"></span> <a name="l00020"></a>00020 <span class="preprocessor">#include "<a class="code" href="itkMembershipFunctionBase_8h.html">itkMembershipFunctionBase.h</a>"</span> <a name="l00021"></a>00021 <span class="preprocessor">#include "<a class="code" href="itkMixtureModelComponentBase_8h.html">itkMixtureModelComponentBase.h</a>"</span> <a name="l00022"></a>00022 <a name="l00023"></a>00023 <span class="keyword">namespace </span>itk{ <a name="l00024"></a>00024 <span class="keyword">namespace </span>Statistics{ <a name="l00025"></a>00025 <a name="l00052"></a>00052 <span class="keyword">template</span>< <span class="keyword">class</span> TSample > <a name="l00053"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html">00053</a> <span class="keyword">class </span><a class="code" href="itkWin32Header_8h.html#7d4384322557205d75f0719c012dded6">ITK_EXPORT</a> <a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html" title="This class generates the parameter estimates for a mixture model using expectation...">ExpectationMaximizationMixtureModelEstimator</a> : <span class="keyword">public</span> <a class="code" href="classitk_1_1Object.html" title="Base class for most itk classes.">Object</a> <a name="l00054"></a>00054 { <a name="l00055"></a>00055 <span class="keyword">public</span>: <a name="l00057"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#3fbabdc1094fb40aef2b3ae6797b2794">00057</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html" title="This class generates the parameter estimates for a mixture model using expectation...">ExpectationMaximizationMixtureModelEstimator</a> <a class="code" href="classitk_1_1LightObject.html" title="Light weight base class for most itk classes.">Self</a>; <a name="l00058"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#bfe6957469d79941d7d0f1f94e4f53f7">00058</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1Object.html" title="Base class for most itk classes.">Object</a> <a class="code" href="classitk_1_1LightObject.html" title="Light weight base class for most itk classes.">Superclass</a>; <a name="l00059"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#3ca2ad314e0d7572288b8ffa3009828f">00059</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1SmartPointer.html">SmartPointer< Self ></a> <a class="code" href="classitk_1_1SmartPointer.html">Pointer</a>; <a name="l00060"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#b93d0366cd110726c688ab588c4a0149">00060</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1SmartPointer.html">SmartPointer<const Self></a> <a class="code" href="classitk_1_1SmartPointer.html">ConstPointer</a>; <a name="l00061"></a>00061 <a name="l00063"></a>00063 <a class="code" href="itkMacro_8h.html#574c1f86661d409dd32c742307a0e96a">itkTypeMacro</a>(<a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html" title="This class generates the parameter estimates for a mixture model using expectation...">ExpectationMaximizationMixtureModelEstimator</a>, <a name="l00064"></a>00064 <a class="code" href="classitk_1_1Object.html" title="Base class for most itk classes.">Object</a>); <a name="l00065"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#383e2498afe3ea7ee4315b9ea4b13076">00065</a> <a class="code" href="itkMacro_8h.html#0fdc15685afa271857c57e7cd482f30c">itkNewMacro</a>(<a class="code" href="classitk_1_1LightObject.html" title="Light weight base class for most itk classes.">Self</a>) ; <a name="l00067"></a>00067 <a name="l00069"></a>00069 <span class="keyword">typedef</span> <span class="keyword">typename</span> TSample::MeasurementType MeasurementType ; <a name="l00070"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#20fe0b7a2a07992cebbe66c41c5c8230">00070</a> <span class="keyword">typedef</span> <span class="keyword">typename</span> TSample::MeasurementVectorType MeasurementVectorType ; <a name="l00071"></a>00071 <a name="l00073"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#887c82f8bf88331628ced16faddc047f">00073</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1Statistics_1_1MixtureModelComponentBase.html" title="base class for distribution modules that supports analytical way to update the distribution...">MixtureModelComponentBase< TSample ></a> <a class="code" href="classitk_1_1Statistics_1_1MixtureModelComponentBase.html" title="base class for distribution modules that supports analytical way to update the distribution...">ComponentType</a> ; <a name="l00074"></a>00074 <a name="l00076"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#d3004a8761d64ed85abc5e033333bc68">00076</a> <span class="keyword">typedef</span> std::vector< ComponentType* > ComponentVectorType ; <a name="l00077"></a>00077 <a name="l00079"></a>00079 <span class="keyword">typedef</span> <a class="code" href="classitk_1_1Statistics_1_1MembershipFunctionBase.html" title="MembershipFunctionBase class declares common interfaces for membership functions...">MembershipFunctionBase< MeasurementVectorType ></a> <a name="l00080"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#51135253e8e6473d8c2be2473ca02614">00080</a> <a class="code" href="classitk_1_1Statistics_1_1MembershipFunctionBase.html" title="MembershipFunctionBase class declares common interfaces for membership functions...">ComponentMembershipFunctionType</a> ; <a name="l00081"></a>00081 <a name="l00083"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#2cf4267a6004982f615c9fe17e78a26d">00083</a> <span class="keyword">typedef</span> <a class="code" href="classitk_1_1Array.html">Array< double ></a> <a class="code" href="classitk_1_1Array.html">ProportionVectorType</a> ; <a name="l00084"></a>00084 <a name="l00086"></a>00086 <span class="keywordtype">void</span> SetSample(<span class="keyword">const</span> TSample* sample) ; <a name="l00087"></a>00087 <a name="l00089"></a>00089 <span class="keyword">const</span> TSample* GetSample() <span class="keyword">const</span>; <a name="l00090"></a>00090 <a name="l00093"></a>00093 <span class="keywordtype">void</span> SetInitialProportions(<a class="code" href="classitk_1_1Array.html">ProportionVectorType</a> &propotion) ; <a name="l00094"></a>00094 <a class="code" href="classitk_1_1Array.html">ProportionVectorType</a>* GetInitialProportions() ; <a name="l00096"></a>00096 <a name="l00098"></a>00098 <a class="code" href="classitk_1_1Array.html">ProportionVectorType</a>* GetProportions() ; <a name="l00099"></a>00099 <a name="l00104"></a>00104 <span class="keywordtype">void</span> SetMaximumIteration(<span class="keywordtype">int</span> numberOfIterations) ; <a name="l00105"></a>00105 <span class="keywordtype">int</span> GetMaximumIteration() ; <a name="l00107"></a>00107 <a name="l00109"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#96a8dbddc1c8ff9aece12690a8b88cf6">00109</a> <span class="keywordtype">int</span> GetCurrentIteration() <a name="l00110"></a>00110 { <span class="keywordflow">return</span> m_CurrentIteration ; } <a name="l00111"></a>00111 <a name="l00113"></a>00113 <span class="keywordtype">int</span> AddComponent(ComponentType* component) ; <a name="l00114"></a>00114 <a name="l00116"></a>00116 <span class="keywordtype">unsigned</span> <span class="keywordtype">int</span> GetNumberOfComponents() ; <a name="l00117"></a>00117 <a name="l00119"></a>00119 <span class="keywordtype">void</span> Update() ; <a name="l00120"></a>00120 <a name="l00122"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#01d72eb1f529a40a889b37205f591ba8">00122</a> <span class="keyword">enum</span> <a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#01d72eb1f529a40a889b37205f591ba8">TERMINATION_CODE</a> { CONVERGED = 0, NOT_CONVERGED = 1 } ; <a name="l00123"></a>00123 <a name="l00125"></a>00125 TERMINATION_CODE GetTerminationCode() ; <a name="l00126"></a>00126 <a name="l00129"></a>00129 ComponentMembershipFunctionType* GetComponentMembershipFunction(<span class="keywordtype">int</span> componentIndex) ; <a name="l00130"></a>00130 <a name="l00131"></a>00131 <span class="keyword">protected</span>: <a name="l00132"></a>00132 <a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html" title="This class generates the parameter estimates for a mixture model using expectation...">ExpectationMaximizationMixtureModelEstimator</a>() ; <a name="l00133"></a><a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#2b05829b0f858ae9ba2843be18f5be6f">00133</a> <span class="keyword">virtual</span> <a class="code" href="classitk_1_1Statistics_1_1ExpectationMaximizationMixtureModelEstimator.html#2b05829b0f858ae9ba2843be18f5be6f">~ExpectationMaximizationMixtureModelEstimator</a>() {} <a name="l00134"></a>00134 <span class="keywordtype">void</span> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#920ae03048a5ad715d803ca2e92b52de">PrintSelf</a>(std::ostream& os, <a class="code" href="classitk_1_1Indent.html" title="Control indentation during Print() invocation.">Indent</a> indent) <span class="keyword">const</span> ; <a name="l00135"></a>00135 <a name="l00136"></a>00136 <span class="keywordtype">bool</span> CalculateDensities() ; <a name="l00137"></a>00137 <span class="keywordtype">double</span> CalculateExpectation() ; <a name="l00138"></a>00138 <span class="keywordtype">bool</span> UpdateComponentParameters() ; <a name="l00139"></a>00139 <span class="keywordtype">bool</span> UpdateProportions() ; <a name="l00140"></a>00140 <a name="l00142"></a>00142 <span class="keywordtype">void</span> <a class="code" href="namespaceHardConnectedComponentImageFilter.html#f75a5945cdd8315676941024a0b8c1e9">GenerateData</a>() ; <a name="l00143"></a>00143 <a name="l00144"></a>00144 <span class="keyword">private</span>: <a name="l00146"></a>00146 <span class="keyword">const</span> TSample* m_Sample ; <a name="l00147"></a>00147 <a name="l00148"></a>00148 <span class="keywordtype">int</span> m_MaxIteration ; <a name="l00149"></a>00149 <span class="keywordtype">int</span> m_CurrentIteration ; <a name="l00150"></a>00150 TERMINATION_CODE m_TerminationCode ; <a name="l00151"></a>00151 ComponentVectorType m_ComponentVector ; <a name="l00152"></a>00152 ProportionVectorType m_InitialProportions ; <a name="l00153"></a>00153 ProportionVectorType m_Proportions ; <a name="l00154"></a>00154 } ; <span class="comment">// end of class</span> <a name="l00155"></a>00155 <a name="l00156"></a>00156 <a name="l00157"></a>00157 } <span class="comment">// end of namespace Statistics </span> <a name="l00158"></a>00158 } <span class="comment">// end of namespace itk</span> <a name="l00159"></a>00159 <a name="l00160"></a>00160 <a name="l00161"></a>00161 <span class="preprocessor">#ifndef ITK_MANUAL_INSTANTIATION</span> <a name="l00162"></a>00162 <span class="preprocessor"></span><span class="preprocessor">#include "itkExpectationMaximizationMixtureModelEstimator.txx"</span> <a name="l00163"></a>00163 <span class="preprocessor">#endif</span> <a name="l00164"></a>00164 <span class="preprocessor"></span> <a name="l00165"></a>00165 <span class="preprocessor">#endif</span> <a name="l00166"></a>00166 <span class="preprocessor"></span> <a name="l00167"></a>00167 <a name="l00168"></a>00168 <a name="l00169"></a>00169 <a name="l00170"></a>00170 <a name="l00171"></a>00171 <a name="l00172"></a>00172 <a name="l00173"></a>00173 </pre></div></div> <hr><address><small> Generated at Thu May 7 22:38:56 2009 for ITK by <a href="http://www.stack.nl/~dimitri/doxygen/index.html"> <img src="http://www.stack.nl/~dimitri/doxygen/doxygen.png" alt="doxygen" align="middle" border=0 width=110 height=53> </a> 1.5.9 written by <a href="mailto:dimitri@stack.nl">Dimitri van Heesch</a>, © 1997-2000</small></address> </body> </html>