[转]Visual C++与Matlab混合编程方法在图像处理中的应用
1   引言
Matlab是美国Mathworks公司于20世纪80年代中期推出的一套数值计算软件，可以实现数值
分析、矩阵运算、自动控制、信号处理和图像处理等若干领域的计算和图形显示功能。它
不仅包含大量高度集成的函数可供直接调用以解决各种复杂的计算，而且提供了简洁的人
机界面、丰富的图形用户界面GUI(Graphical User Interfaces)开发功能以及求解特定学
科问题的求解工具箱ToolBox。与此同时编写程序Matlab所需时间往往能比使用其他编程语
言缩短许多倍。自面世以来，Matlab在教学和科研等领域受到了广泛的重视和应用, 在解
决复杂的数学计算和新算法的研究中具有极大的优势。
1.1   软件的特点归纳[4]
(1) 易学易用性  
Matlab是一门编程语言，其语法规则和C语言大同小异，因而具有一般编程语言基础的用户
很快就可以掌握。  
(2) 代码短小高效  
Matlab将不同数学分支的算法以函数的形式分类成库，使用时直接调用这些函数并赋予实
际参数就可以解决问题，快速而准确。
(3) 计算功能强大  
该软件具有强大的矩阵计算功能，除了一般的加、减、乘、除、转置和求逆运算外，还非
常适合处理稀疏矩阵等特殊矩阵运算以及有限元等大型数值算法的编程。
(4) 可扩展性能  
Matlab最重要的特点是易于扩展，它允许用户自行编写指定功能的m文件，组成自己的工具
箱。当前支持Matlab的商用ToolBox(工具箱)有数百种之多。Matlab支持DDE和ActiveX自动
化等机制，可以与同样支持该技术的应用程序接口。利用m文件或MEX文件可以实现也VB、
VC等程序的进程内无缝集成;利用web服务器，可实现Matlab与网络程序的接口;利用API函
数，还可以实现Matlab与硬件的接口。
1.2   Matlab自身所存在的某些不足也限制了它的应用范围
(1) Matlab是一种解释性语言，所以它的代码执行速度慢，这对于实时性要求较高的领域
，如自动控制、信号处理等，其实时效率是较差的。
(2) Matlab程序不能脱离其环境运行，其功能只能在它本身所提供的平台上使用，因而不
能被用于开发商用软件。
(3) Matlab程序可被直接看到源代码，因而算法和数据的保密性欠佳。
Visual C++(以下简称VC)是Microsoft公司推出的强大的可视化集成编程环境, 从底层软件
到上层直接面向用户的软件都可以用VC开发, 而且强大的调试功能也为大型复杂软件的开
发提供了高效的排错手段。使用VC开发的系统具有界面友好、代码效率高和执行速度快等
一系列优点. 同时C++语言支持面向对象的程序设计, 而利用面向对象方法设计的软件具有
良好的可重用性、可维护性和可扩充性等特点。它是目前综合性{zg}、最强大、也是最复
杂的软件开发工具之一, 应用极为广泛。
因此实现VC与Matlab混合编程，使两者结合起来，协同工作，将大大减少编程的工作量、
提高程序执行效率。利用VC这一优秀编程工具可以弥补Matlab在处理具体问题时的缺陷,利
用Matlab的强大数值运算等功能则可以增强VC对信号处理、三维处理、自动控制等方面的
能力。两者结合将会充分发挥各自的优势，必将提高软件开发效率，使开发的软件具有更
高的性能，更广阔的应用领域。

2   VC++与Matlab的接口技术简介
2.1   在VC环境中调用Matlab程序主要方法
为将Matlab的强大功能融入到各种应用程序中，通过高级语言编译器将Matlab的m文件编译
为二进制代码已成为研究热点。
(1) 在VC中启用MATLAB EN GINE(引擎)，采用客户机/服务器(Client/Server)的计算模式
。这种方式需要Matlab在后台运行，离不开Matlab环境。但是它可以调用Matlab工具箱函
数和图形函数。
(2) 利用MATLAB COMPIL ER(编译器)，将Matlab 的函数编译成脱离Matlab环境的可执行程
序(.exe文件)， 在VC中调用。
(3) 利用Matlab的数学函数和图形库。6.1版本的Matlab软件包中提供了C/C++的数学函数
和图形库，通过其编译器可以将Matlab中编写的m文件转换成以C/C++代码的文件，而且可
以将m文件生成dll库，甚至可以直接调用其库函数，生成不依赖Matlab的可执行文件。

(4) MIDEVA(Matcom)是{dy}个由Matlab到C++的编译开发软件平台.使用Matcom编译器可以
将Matlab源代码译成同等功能的C++代码，并编译为EXE或DLL文件。既保持了Matlab的优良
算法，又提高了执行速度。它还支持一定的图形显示，生成代码的可读性很好。因其简单
便捷、功能强大、应用灵活，本文采用这种实现方案。

3   Matcom在图像处理程序中的具体使用方法和技巧
下面以实现图像处理中绘制直方图和快速傅立叶变换的编程过程为例，介绍通过Matcom实
现VC对Matlab程序的调用过程，运行环境为:WindowsXP、Matlab6.5、Matcom4.5、Visual
C++6.0。    
3.1   编译环境设置[2]    
首先安装MIDEVA，然后运行Visual C++6.0，从菜单条中选择Tools/Customize/Add-ins a
nd Macro Files，在出现的对话框中选中Visual Matcom Add-in，关闭对话框，这时在Vi
sual C++[1]的开发环境中看到一个Visual Matcom工具栏，表明安装成功。

3.2   代码编写
在MATCOM或者Matlab环境中编写实现绘制图像直方图的文件myhist.m代码如下:
function h=myhist(x)   ％绘制灰度图像的直方图
h=zeros(1,256);[m,n]=size(x)
for i=1:m
for j=1:n
h(1,x(i,j)+1)=h(1,x(i,j)+1)+1;
end
end
plot(0:255,h)
建立一个新的VC++工程或者使用现有的图像处理的MFC工程Dip.dsw，用C++实现DIB图像的
读取，获取图像相关信息，图像显示，保存等基本功能[3]。再建立一个用于图像频域处理
的类CFreqPro，将来在此类中添加有关频域处理的成员函数。点击Visual Matcom工具栏上
的m++图标，选择保存过的Matlab文件myhist.m进行转化。这时会在VC中出现一个转换完毕
的说明文件，文件中如果报告有错误，可以双击C++files文件夹下的myhist.m直接进行修
改，不必回到Matlab环境，再重新转化直到没有错误报告为止。在需要使用myhist()函数
的文件中添加头文件matlib.h和myhist.h，加入代码#include”matlib.h”，#include”
myhist.h”。
在调用matcom的Matrix<LIB> C++库函数之前用initM(MATCOM_VERSION)初始化类库调用，
并用exitM()结束对类库的调用。在MFC工程文件中，一般是在一个类(例如CDipDoc)的构造
函数中添加initM(MATCOM_VERSION)以完成类库的初始化工作，在该类的析构函数中调用e
xitM()结束对类库的调用。
3.3   创建图像数据矩阵和将数据矩阵赋给图像数据区的功能核心代码
在matcom中，提供了一个双精度Matrix类型－Mm，因为所有的操作均为矩阵运算。为利用
矩阵运算完成图像变换，首先应将其图像数据赋给矩阵变量，而且能将运算后的矩阵变量
再回赋到图像数据区[5]。为此，在CFreqPro中添加成员函数GetMatData()和SetMatData(
)，分别实现创建图像数据矩阵和将数据矩阵赋给图像数据区的功能。这两个函数对常用的
8位灰度图像处理具有较好的通用性和参考价值。核心代码如下(有删节):
Mm CFreqPro::GetMatData()
{ DWORD SizeImage = nWidth * nHeight;
//创建图像数据矩阵，并将其元素初始值设为0
m_matBits = zeros(1, SizeImage);
//默认的矩阵数据类型是double，首先将其转换为unsigned char型
//以便用memcpy内存拷贝命令将图像数据赋给矩阵。
m_matBits = muint8(m_matBits);
//通过Matrix<LIB>C++库的.addr()函数返回矩阵变量的内存指针，完成对矩阵单元的访问

//用memcpy内存拷贝命令将图像数据赋给矩阵。
memcpy(m_matBits.addr(), pBits, SizeImage);
//由于Mm类型的矩阵是按先列后行的顺序排列，
//在此用reshape()函数将创建的一维矩阵m_matBits变为nWidth×nHeight的二维矩阵。

//再用rot90()函数将二维矩阵逆时针旋转90度，将矩阵变为nHeight×nWidth的二维矩阵
，
//并使矩阵的第nHeight行对应图像数据的{dy}行数据。
m_matBits = rot90(reshape(m_matBits, nWidthBytes, nHeight));
//将矩阵由unsigned char型转换为double型，以便进行运算
m_matBits = mdouble(m_matBits);
}
///////////////////////////////////
BOOL SetMatData()
{ //将矩阵数据范围限定于(0—255)
m_matBits = mabs(m_matBits);
Mm L = m_matBits > 255.0;
Mm NotL = !L;
m_matBits = times(m_matBits, NotL) + L * 255.0;
//将矩阵由double型转换为unsigned char型，以便将图像数据赋给矩阵
m_matBits = muint8(m_matBits);
//对矩阵进行转置操作
m_matBits = rot90(m_matBits, -1);
//将图像数据赋给矩阵
memcpy(pBits, m_matBits.addr(), (nWidth * nHeight)*sizeof(unsigned char));
return( TRUE );}
3.4   图像的傅立叶变换和绘制直方图函数
图像数据矩阵建立后，就可以利用Matrix<Lib>C++库函数通过各种矩阵运算完成图像变换
的工作。下面的函数仅实现图像的傅立叶变换和绘制直方图，用同样方法，离散余弦变换
、离散沃尔什－哈达玛变换也都能实现。在图像频域处理类CFreqPro中添加成员函数FFT2
()
Mm CFreqPro::FFT2(CDib *pDibObject)
{if( pDib!= NULL ) m_pDib = pDib;
//创建图像数据矩阵
GetMatData();
//获得矩阵的行数和列数
Mm mSize = size(m_matBits);
//调用Matrix<Lib>C++库函数fft2()完成二维离散傅立叶变换
Mm ff = fft2(m_matBits);
Mm matTransed = ff;
//调用Matrix<Lib>C++库函数fftshift()将频域中心移到矩阵中心
ff = fftshift(ff);
//调用Matrix<Lib>C++库函数mabs()计算频谱
m_matBits = mabs(ff) / sqrt(mSize.r(1,1)*mSize.r(1,2));
//绘制直方图
myhist(m_matBits);  
//将矩阵数据赋给图像数据区
SetMatData();
return matTransed;}
运用本方法对一8位灰度图像(如图2)进行了离散傅立叶变换，结果如图3所示，利用myhis
t()函数实现直方图绘制。

4   结束语
在图像处理程序中，既可以直接调用调用Matrix<Lib>C++库函数，也可以在VC环境中转换
.m函数文件。如果需要转换的.m文件不是一个函数，而是脚本文件(Script)，则要在工程
目录下找到转换文件的.cpp文件，将其中的C代码拷贝到需要调用它的函数里面,也可以按
照MATCOM的语法(类似Matlab的语法)直接进行编写。Matcom还可以实现函数的嵌套。当所
编译的.m文件依赖于其他.m文件时，只要把被调用的.m文件与要编译的.m文件放在同一目
录下，把生成的被调用文件的.h和.cpp文件插入到VC开发的工程中就可以了。作者在图像
识别程序中使用这种方法来提高程序编写的效率，取得了非常好的效果。   

参考文献
[1] 张志涌等. 精通MATLAB6.5版. 北京:北京航空航天大学出版社,2003
[2] MathTools Ltd,MIDEVA.MATCOM&Visual MATCOM Installtion Guide 1999
[3] 郎锐. 数字图像处理学Visual c++实现. 北京:北京希望电子出版社,2003
[4] 阮沈勇等. MATLAB程序设计. 北京:电子工业出版社,2004
[5] Duane Hanselman,B ruce littlefield. 精通MATLAB综合辅导与指南. 西安:西安交通
大学出版社,1998作者简介
刘馨月(1978-)   女   硕士研究生   研究方向:图像处理与模式识别。