有一个问题进行向量化

1视图(30天)

显示旧的评论

若谢泼德 2016年7月12日

0
链接

这个问题直接联系

//www.tatmou.com/matlabcentral/answers/295166-how-to-vectorize-a-problem

大家好,我想要一些建议关于vectorizing问题。

我计算图像配准和我使用的是两个for循环以浏览整个图像和需要多少时间。所以我想知道我可以vectorize它。

这是我的部分代码:

                         [a, b] =大小(我);
                        
                         为我= 1:2 * a - 1
                        
                         %翻译向量
                        
                         tx = a: 1:——+ 1;
                        
                         tz = b: 1: - b + 1;
                        
                         %设置我为每一个循环
                        
                         I_trans =我;
                        
                         % x翻译
                        
                         I_trans = imtranslate (I_trans [0 tx (i) 1]);
                        
                         为j = 1:2 * b - 1
                        
                         %设置I_trans_x每个循环
                        
                         I_trans_x = I_trans;
                        
                         % z翻译
                        
                         I_trans_z = imtranslate (I_trans_x [tz (j) 1 0]);
                        
                         %删除行和列的0
                        
                         I_trans_z (~ (I_trans_z 2):) = [];
                        
                         I_trans_z (:, ~ (I_trans_z, 1)) = [];
                        
                         [crop_x, crop_z] =大小(I_trans_z);
                        
                         如果(crop_x > =我)
                        
                         c = b-crop_z + 1;
                        
                         d = a-crop_x + 1;
                        
                         rect_J = [c d crop_z crop_x];
                        
                         J_crop = imcrop (J, rect_J);
                        
                         其他的
                        
                         rect_J = [0 0 crop_z crop_x];
                        
                         J_crop = imcrop (J, rect_J);
                        
                         结束
                        
                         % GI计算
                        
                         胃肠道= GI_1 (I_trans_z J_crop);
                        
                         GI_tab (i, j) =胃肠道;
                        
                         %图
                        
                         % imshowpair (I_trans_z J_crop,“蒙太奇”);
                        
                         %标题(int2str (GI));
                        
                         结束
                        
                         结束