24小时热门版块排行榜    

查看: 511  |  回复: 0

eeds136

金虫 (著名写手)

[求助] 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');");
回复此楼

» 猜你喜欢

一切都似那么神奇。。。。
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

找到一些相关的精华帖子,希望有用哦~

科研从小木虫开始,人人为我,我为人人
相关版块跳转 我要订阅楼主 eeds136 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 333求调剂 +5 87639 2026-03-21 7/350 2026-03-21 19:31 by ColorlessPI
[考研] 求调剂 +3 13341 2026-03-20 3/150 2026-03-21 18:28 by 学员8dgXkO
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 0703化学297求调剂 +3 Daisy☆ 2026-03-20 3/150 2026-03-21 17:45 by ColorlessPI
[考研] 生物学一志愿985,分数349求调剂 +3 zxts12 2026-03-21 3/150 2026-03-21 16:34 by 33来了真来了
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 0856材料专硕353求调剂 +3 NIFFFfff 2026-03-20 3/150 2026-03-21 10:23 by luoyongfeng
[考研] 二本跨考郑大材料306英一数二 +3 z1z2z3879 2026-03-17 3/150 2026-03-21 02:29 by JourneyLucky
[考研] 304求调剂 +7 司空. 2026-03-18 7/350 2026-03-20 23:08 by JourneyLucky
[考研] 求调剂,一志愿:南京航空航天大学大学 ,080500材料科学与工程学硕,总分289分 +4 @taotao 2026-03-19 4/200 2026-03-20 22:14 by JourneyLucky
[考研] 一志愿中南化学(0703)总分337求调剂 +8 niko- 2026-03-19 9/450 2026-03-20 21:57 by luoyongfeng
[考研] 材料与化工专硕调剂 +7 heming3743 2026-03-16 7/350 2026-03-20 19:31 by zhukairuo
[考研] 工科材料085601 279求调剂 +7 困于星晨 2026-03-17 9/450 2026-03-20 17:38 by 无懈可击111
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +5 1孙悟空 2026-03-17 6/300 2026-03-19 23:46 by zcl123
[考研] 081700化工学硕调剂 +3 【1】 2026-03-16 3/150 2026-03-19 23:40 by edmund7
[考研] 一志愿福大288有机化学,求调剂 +3 小木虫200408204 2026-03-18 3/150 2026-03-19 13:31 by houyaoxu
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
信息提示
请填处理意见