调用matlab引擎通过txt绘图

sunminmin

在matlab的commond窗口中键入下面的代码能得到map.txt对应的图形：(注:map.txt中每一行为一个点的坐标（x,y,z)。）

1. M = dlmread('map.txt'); % 读取map.txt文件中的数据
   
2. x = M(:,1); % 取第一列
   
3. y = M(:,2); % 取第二列
   
4. z = M(:,3); % 取第三列
   
5. x = reshape(x,65,65); % 重构为size=65*65的方格
   
6. y = reshape(y,65,65); % 重构为size=65*65的方格
   
7. z = reshape(z,65,65); % 重构为size=65*65的方格
   
8. surf(x,y,z) % 绘图
   
9. shading flat %各小曲面之间不要网格
   
10. xlabel('X(单位：千米)')
    
11. ylabel('Y(单位：千米)')
    
12. zlabel('Z(单位：千米)')

复制代码

但是在VC++中建立MFC AppWizard基于基本对话框的工程TxtToMap,调用matlab引擎绘不出map.txt对应的图形：

   

1. void CTxtToMapDlg::OnButtonBrowsing()
   
2. {
   
3. // TODO: Add your control notification handler code here
   
4. char szFilters[] = "MyType Files (*.txt)|*.txt|All Files (*.*)|*.*||";
   
5. 
   
6. CFileDialog fileDlg (TRUE, "txt", "*.txt",OFN_FILEMUSTEXIST| OFN_HIDEREADONLY, szFilters, this);
   
7. 
   
8. // Display the file dialog. When user clicks OK, fileDlg.DoModal()
   
9. // returns IDOK.
   
10. if( fileDlg.DoModal()==IDOK )
    
11. {
    
12. CString pathName = fileDlg.GetPathName();
    
13. 
    
14. m_FilePath = pathName; //m_FilePath里保存txt文件的路径
    
15. 
    
16. UpdateData(FALSE);
    
17. }
    
18. }
    

复制代码

1. void CTxtToMapDlg::OnButtonTxtToMap()
   
2. {
   
3. // TODO: Add your control notification handler code here
   
4. 
   
5. 
   
6. Engine *ep;
   
7. 
   
8. if (!(ep = engOpen(NULL))) //打开Matlab引擎
   
9. 
   
10. {
    
11. 
    
12. MessageBox ("Can't start MATLAB engine", MB_OK);
    
13. 
    
14. exit(-1);
    
15. }
    
16. 
    
17. int N = 4225;
    
18. mxArray *M = mxCreateDoubleMatrix(N,3,mxREAL);
    
19. mxArray *x = mxCreateDoubleMatrix(1,N,mxREAL);
    
20. mxArray *y = mxCreateDoubleMatrix(1,N,mxREAL);
    
21. mxArray *z = mxCreateDoubleMatrix(1,N,mxREAL);
    
22. engPutVariable(ep,"M",M);
    
23. engPutVariable(ep,"x",x);
    
24. engPutVariable(ep,"y",y);
    
25. engPutVariable(ep,"z",z);
    
26. engEvalString(ep,"M = dlmread('m_FilePath')"); //m_FilePath保存map.txt的路径，读取map.txt文件中的数据
    
27. engEvalString(ep,"x = M(:,1)"); //取第一列
    
28. engEvalString(ep,"y = M(:,2)"); // 取第二列
    
29. engEvalString(ep,"z = M(:,3)"); // 取第三列
    
30. engEvalString(ep,"x = reshape(x,65,65)"); //重构为size=65*65的方格
    
31. engEvalString(ep,"y = reshape(y,65,65)"); //重构为size=65*65的方格
    
32. engEvalString(ep,"z = reshape(z,65,65)"); //重构为size=65*65的方格
    
33. engEvalString(ep,"surf(x,y,z)"); //绘图
    
34. 
    
35. }
    

复制代码

我刚开始学vc++与matlab混合编程调用matlab引擎，我自己也知道上面void CTxtToMapDlg::OnButtonTxtToMap() 中的代码写的有错，请问如何修改才能得到map.txt对应的图形？

map.rar (29.55 KB, 下载次数: 2)

2011-5-9 20:55 上传

点击文件名下载附件
map.txt

sunminmin

回复 1 # sunminmin 的帖子

如何调用matlab引擎把map.txt load进matlab，让matlab绘出对应的三维图？关键是这个数据量太大，怎样拟合才能通过matlab引擎把图绘出来？
上面的void CTxtToMapDlg::OnButtonTxtToMap()怎么修改，那样写错在哪里？

sunminmin

能这样用吗？   engEvalString(ep,"M = dlmread('m_FilePath')"); //m_FilePath保存map.txt的路径，读取map.txt文件中的数据

sunminmin

还是靠自己解决了。

1. 
   
2.   Engine *ep;
   
3. 
   
4. if (!(ep = engOpen(NULL))) //打开Matlab引擎
   
5.   
   
6. {              
   
7.   
   
8.   MessageBox ("Can't start MATLAB engine", MB_OK);
   
9.   
   
10.   exit(-1);
    
11.   
    
12. }
    
13. 
    
14. 
    
15. int N = 4225; //txt文件中有4225行数据
    
16. mxArray *M = mxCreateDoubleMatrix(N,3,mxREAL);
    
17. mxArray *x = mxCreateDoubleMatrix(1,N,mxREAL);
    
18. mxArray *y = mxCreateDoubleMatrix(1,N,mxREAL);
    
19. mxArray *z = mxCreateDoubleMatrix(1,N,mxREAL);
    
20. mxArray *xx = mxCreateDoubleMatrix(65,65,mxREAL);
    
21. mxArray *yy = mxCreateDoubleMatrix(65,65,mxREAL);
    
22. mxArray *zz = mxCreateDoubleMatrix(65,65,mxREAL);
    
23. engPutVariable(ep,"M",M);
    
24. engPutVariable(ep,"x",x);
    
25. engPutVariable(ep,"y",y);
    
26. engPutVariable(ep,"z",z);
    
27. engPutVariable(ep,"xx",x);
    
28. engPutVariable(ep,"yy",y);
    
29. engPutVariable(ep,"zz",z);
    
30. 
    
31. CString s1("M = load('");
    
32. CString s2("')");
    
33. CString s = s1+m_FilePath+s2;
    
34. char * p = (LPSTR)(LPCTSTR)s;
    
35. engEvalString(ep,p);
    
36. engEvalString(ep,"x = M(:,1)"); //取第一列
    
37. engEvalString(ep,"y = M(:,2)"); // 取第二列
    
38. engEvalString(ep,"z = M(:,3)"); // 取第三列
    
39. engEvalString(ep,"xx = reshape(x,65,65)"); //重构为size=65*65的方格
    
40. engEvalString(ep,"yy = reshape(y,65,65)"); //重构为size=65*65的方格
    
41. engEvalString(ep,"zz = reshape(z,65,65)"); //重构为size=65*65的方格
    
42. engEvalString(ep,"mesh(xx,yy,zz)"); //绘图
    
43. engEvalString(ep,"xlabel('X(单位：千米)')");
    
44. engEvalString(ep,"ylabel('Y(单位：千米)')");
    
45. engEvalString(ep,"zlabel('Z(单位：百米)')");

复制代码

sunminmin

解决问题之后回来分享经验

在那个调用引擎画图的那个消息响应函数里有错，即engEvalString(ep,"M = load('m_FilePath')");有错。
就是调用m_filename变量的时候，传的是“M = load(m_filename.txt)"是字符串啊，没办法调出m_filename的值来
M = load('m_FilePath')"这句话再matlab里是不能够知道m_FilePath是什么的，得把几个字符串拼接起来
M = load('+ m_FilePath+ ')" 把这3段字符串拼接然后传入matlab 。


现在的问题就是
比如
m_FilePath = "C:\New folder\my.txt"
engEvalString(ep,"M = load('m_FilePath')");这句话运行后相当于在matlab里运行这句话M = load('m_FilePath')，
但是实际上应该传入这句话M = load('C:\New folder\my.txt')


用下面的语句
CString s1("M = load('");
CString s2("')");
Ctring s = s1+m_FilePath+s2;
char * p = (LPSTR)(LPCTSTR)s;
engEvalString(ep,p);
替换engEvalString(ep,"M = load('+m_FilePath+')");
这句


问题就解决了！