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