图像分割 高斯混合模型 马尔科夫 image segmentation MRF GMM ICM

image segmentation MRF GMM ICM 图像分割 马尔科夫随机场 高斯混合模型 ICM算法 另附详细的说明pdf 讲清GMM，MRF做分割的详细原理，可作为参考。

/**
 * Copyright （C） 2012 Quan Wang 
 * Signal Analysis and Machine Perception Laboratory
 * Department of Electrical Computer and Systems Engineering
 * Rensselaer Polytechnic Institute Troy NY 12180 USA
 */

/**
 * This function expands the matrix using mirror boundary condition. 
 * For example 
 *
 * A = [
 *     1  2  3  11
 *     4  5  6  12
 *     7  8  9  13
 *     ]
 *
 * B = BoundMirrorExpand（A） will yield
 *
 *     5  4  5  6  12  6
 *     2  1  2  3  11  3
 *     5  4  5  6  12  6 
 *     8  7  8  9  13  9 
 *     5  4  5  6  12  6
 */

#include “mex.h“
#include 
#include 
#include 
#include 

double *boundMirrorExpand（double *A int m int n）
{
    double *B=new double[（m+2）*（n+2）];
    B[0+0*（m+2）]=A[1+1*m];
    B[（m+1）+0*（m+2）]=A[（m-2）+1*m];
    B[0+（n+1）*（m+2）]=A[1+（n-2）*m];
    B[（m+1）+（n+1）*（m+2）]=A[（m-2）+（n-2）*m];
    for（int i=1;i    {
        B[0+i*（m+2）]=A[1+（i-1）*m];
        B[m+1+i*（m+2）]=A[m-2+（i-1）*m];
    }
    for（int i=1;i    {
        B[i+0*（m+2）]=A[（i-1）+1*m];
        B[i+（n+1）*（m+2）]=A[i-1+（n-2）*m];
    }
    for（int i=1;i    {
        for（int j=1;j        {
            B[i+j*（m+2）]=A[i-1+（j-1）*m];
        }
    }
    return B;
}

void copyMatrix（double *A double *B int m int n）
{
    for（int i=0;i    {
        B[i]=A[i];
    }
}

void mexFunction（ int nlhs mxArray *plhs[]
        int nrhs const mxArray *prhs[]）
{
    if（nrhs!=1）
    {
        mexErrMsgIdAndTxt（ “MATLAB:BoundMirrorExpand:invalidNumInputs“
                “One input required.“）;
    }
    if（nlhs!=1）
    {
        mexErrMsgIdAndTxt（ “MATLAB:BoundMirrorExpand:invalidNumOutputs“
                “One output required.“）;
    }
    
    double *A=mxGetPr（prhs[0]）;
    int m=mxGetM（prhs[0]）;
    int n=mxGetN（prhs[0]）;
    
    if（m<3 || n<3）
    {
        mexErrMsgIdAndTxt（ “MATLAB:BoundMirrorExpand:matrixTooSmall“
                “Input matrix is too small.“）;
    }
    
    plhs[0] = mxCreateDoubleMatrix（m+2 n+2 mxREAL）;
    
    double *B=mxGetPr（plhs[0]）;
    
    double *C=boundMirrorExpand（Amn）;
    
    copyMatrix（C B m+2 n+2）;
    
    delete[] C;

}

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件     3064399  2012-12-20 23:09  GMM-HMRF_v1.1\1212.4527v1.pdf
     文件       84296  2007-03-25 01:38  GMM-HMRF_v1.1\code\color-image\385028.jpg
     文件        2211  2014-03-16 03:48  GMM-HMRF_v1.1\code\color-image\BoundMirrorExpand.cpp
     文件        1852  2014-03-16 03:49  GMM-HMRF_v1.1\code\color-image\BoundMirrorShrink.cpp
     文件        1251  2012-12-20 23:22  GMM-HMRF_v1.1\code\color-image\HMRF_EM.m
     文件        2273  2012-12-20 23:22  GMM-HMRF_v1.1\code\color-image\MRF_MAP.m
     文件         888  2014-03-16 05:26  GMM-HMRF_v1.1\code\color-image\demo.m
     文件         441  2014-03-16 04:27  GMM-HMRF_v1.1\code\color-image\gaussianBlur.m
     文件         481  2012-12-20 23:22  GMM-HMRF_v1.1\code\color-image\get_GMM.m
     文件         452  2012-12-17 10:16  GMM-HMRF_v1.1\code\color-image\image_kmeans.m
     文件         174  2012-04-26 09:06  GMM-HMRF_v1.1\code\color-image\ind2ij.m
     文件        1229  2012-12-20 23:27  GMM-HMRF_v1.1\code\three-dimensional\HMRF_EM.m
     文件      125000  2012-12-20 23:27  GMM-HMRF_v1.1\code\three-dimensional\Image.raw
     文件        2384  2012-12-20 23:27  GMM-HMRF_v1.1\code\three-dimensional\MRF_MAP.m
     文件         921  2012-12-20 23:27  GMM-HMRF_v1.1\code\three-dimensional\demo.m
     文件         618  2012-12-20 23:27  GMM-HMRF_v1.1\code\three-dimensional\generate_3D_image.m
     文件         392  2012-12-20 23:27  GMM-HMRF_v1.1\code\three-dimensional\get_GMM.m
     文件         437  2012-12-17 10:12  GMM-HMRF_v1.1\code\three-dimensional\image_kmeans.m
     文件         280  2012-12-17 10:11  GMM-HMRF_v1.1\code\three-dimensional\ind2ijq.m
     文件      125000  2012-12-20 23:28  GMM-HMRF_v1.1\code\three-dimensional\kmeans.raw
     文件      125000  2012-12-20 23:28  GMM-HMRF_v1.1\code\three-dimensional\segmentation.raw
     文件        1309  2014-03-15 17:36  license.txt