24小时热门版块排行榜

>论坛更新日志 (4056)
>虫友互识 (422)
>文献求助 (396)
>导师招生 (373)
>硕博家园 (123)
>考博 (111)
>论文投稿 (109)
>休闲灌水 (104)
>博后之家 (99)
>招聘信息布告栏 (97)
>催化 (59)
>公派出国 (55)
>基金申请 (51)
>考研 (51)
>教师之家 (45)
>论文道贺祈福 (44)

返回列表

eeds136

金虫 (著名写手)

应助: 2 (幼儿园)
金币: 397
散金: 970
红花: 6
沙发: 6
帖子: 1281
在线: 60.7小时
虫号: 2569615
注册: 2013-07-28
性别: GG
专业: 模式识别

[求助] Matlab可以把生成的数据导入excel,VC调用Matlab引擎却不可以,为什么

程序如下：
void CMyDlg::OnHuiDuGongSheng()
{
// TODO: Add your command handler code here

engEvalString(ep,"Gray=imread('121.bmp');");
engEvalString(ep,"[M,N,O] = size(Gray);");
engEvalString(ep,"M = 128;");
engEvalString(ep,"N = 128;");

//为了减少计算量，对原始图像灰度级压缩，将Gray量化成16级
engEvalString(ep,"for i = 1:M;\
for j = 1:N;\
for n = 1:256/16;\
if (n-1)*16<=Gray(i,j)&&Gray(i,j)<=(n-1)*16+15;\
Gray(i,j) = n-1;\
end;\
end;\
end;\
end;");

//计算四个共生矩阵P,取距离为1，角度分别为0,45,90,135
engEvalString(ep,"P = zeros(16,16,4);");
engEvalString(ep,"for m = 1:16;\
   for n = 1:16;\
   for i = 1:M;\
   for j = 1:N;\
   if j<N&&Gray(i,j)==m-1&&Gray(i,j+1)==n-1;\
   P(m,n,1) = P(m,n,1)+1;\
   P(n,m,1) = P(m,n,1);\
   end;\
   if i>1&&j<N&&Gray(i,j)==m-1&&Gray(i-1,j+1)==n-1;\
   P(m,n,2) = P(m,n,2)+1;\
   P(n,m,2) = P(m,n,2);\
   end;\
   if i<M&&Gray(i,j)==m-1&&Gray(i+1,j)==n-1;\
   P(m,n,3) = P(m,n,3)+1;\
   P(n,m,3) = P(m,n,3);\
   end;\
   if i<M&&j<N&&Gray(i,j)==m-1&&Gray(i+1,j+1)==n-1;\
   P(m,n,4) = P(m,n,4)+1;\
   P(n,m,4) = P(m,n,4);\
   end;\
end;\
end;\
if m==n;\
P(m,n,:) = P(m,n,:)*2;\
end;\
end;\
end;");

//对共生矩阵归一化
engEvalString(ep,"for n = 1:4;\
P(:,:,n) = P(:,:,n)/sum(sum(P(:,:,n)));\
end;");
//对共生矩阵计算能量、熵、惯性矩、相关4个纹理参数
engEvalString(ep,"H = zeros(1,4);");
engEvalString(ep,"I = H;");
engEvalString(ep,"Ux = H;");
engEvalString(ep,"Uy = H;");
engEvalString(ep,"deltaX= H;");
engEvalString(ep,"deltaY = H;");
engEvalString(ep,"C =H;");
engEvalString(ep,"for n = 1:4;\
E(n) = sum(sum(P(:,:,n).^2));\
for i = 1:16;\
for j = 1:16;\
if P(i,j,n)~=0;\
H(n) = -P(i,j,n)*log(P(i,j,n))+H(n);\
end;\
I(n) = (i-j)^2*P(i,j,n)+I(n);\
Ux(n) = i*P(i,j,n)+Ux(n);\
Uy(n) = j*P(i,j,n)+Uy(n);\
end;\
end;\
end;");
   engEvalString(ep,"for n = 1:4;\
                                 for i = 1:16;\
   for j = 1:16;\
   deltaX(n) = (i-Ux(n))^2*P(i,j,n)+deltaX(n);\
   deltaY(n) = (j-Uy(n))^2*P(i,j,n)+deltaY(n);\
   C(n) = i*j*P(i,j,n)+C(n);\
end;\
end;\
C(n) = (C(n)-Ux(n)*Uy(n))/deltaX(n)/deltaY(n);\
   end;");

   //求能量、熵、惯性矩、相关的均值和标准差作为最终8维纹理特征
   engEvalString(ep,"T(1) = mean(E);");
   engEvalString(ep,"T(2) = mean(H);");
   engEvalString(ep,"T(3) = mean(I);");
   engEvalString(ep,"T(4) = mean(C);");
//    engEvalString(ep,"fprintf('%f,%f,%f,%f',T(1),T(2),T(3),T(4));");
   engEvalString(ep,"xlswrite('testdata.xls', [T(1) T(2) T(3) T(4)], 'B2:H2');");

回复此楼