0

実行に時間がかかる一連のコード セクションがあります。Mathworks サイトのベクトル化のページを読みました。一部でまだ少し混乱しています。plane_intersect を実行する部分をベクトル化することはできますか?

ベクトル化されていない

for N = 1:sizeDimages
    imPos = anaInfoSat(N).ImagePositionPatient;
    A2Z = imPos(3);
    A2Y = imPos(2);
    A2X = imPos(1);

    [upP1(N,:)]=plane_intersect([cosSx(1),cosSy(1),cosSz(1)],[uppersatX1,uppersatY1,uppersatZ1],crossS,[A2X,A2Y,A2Z]);
    [loP1(N,:)]=plane_intersect([cosSx(1),cosSy(1),cosSz(1)],[lowersatX1,lowersatY1,lowersatZ1],crossS,[A2X,A2Y,A2Z]);
end

私のベクトル化の試みは、upP1 が Nx3 行列であることです。upP1 マトリックスを事前に割り当てます。以下のコードは、次元の不一致に関するエラーを返します。ImagePosition は 1x3 行列です。

N = 1:sizeDimages;
imPos = anaInfoSat(N).ImagePositionPatient;
A2Z = imPos(3);
A2Y = imPos(2);
A2X = imPos(1);

[upP1(N,:)]=plane_intersect([cosSx(1),cosSy(1),cosSz(1)],[uppersatX1,uppersatY1,uppersatZ1],crossS,[A2X,A2Y,A2Z]);
[loP1(N,:)]=plane_intersect([cosSx(1),cosSy(1),cosSz(1)],[lowersatX1,lowersatY1,lowersatZ1],crossS,[A2X,A2Y,A2Z]);

これは plane_intersect コードの一部であり、それが何をするかを知るのに十分なはずです。

function [P,N,check]=plane_intersect(N1,A1,N2,A2)
%plane_intersect computes the intersection of two planes(if any) 
% Inputs: 
%       N1: normal vector to Plane 1
%       A1: any point that belongs to Plane 1
%       N2: normal vector to Plane 2
%       A2: any point that belongs to Plane 2
%
%Outputs:
%   P    is a point that lies on the interection straight line.
%   N    is the direction vector of the straight line
% check is an integer (0:Plane 1 and Plane 2 are parallel' 
%                              1:Plane 1 and Plane 2 coincide
%                              2:Plane 1 and Plane 2 intersect)
%
% Example:
% Determine the intersection of these two planes:
% 2x - 5y + 3z = 12 and 3x + 4y - 3z = 6
% The first plane is represented by the normal vector N1=[2 -5 3]
% and any arbitrary point that lies on the plane, ex: A1=[0 0 4]
% The second plane is represented by the normal vector N2=[3 4 -3]
% and any arbitrary point that lies on the plane, ex: A2=[0 0 -2]
%[P,N,check]=plane_intersect([2 -5 3],[0 0 4],[3 4 -3],[0 0 -2]);
4

1 に答える 1

3

ベクトル化されたコードanaInfoSat(N).ImagePositionPatient;では、単一の回答は返されませんが、複数の回答が返されます。結果を単一の変数に割り当てると、最初の回答のみを受け取ります。これが、寸法の不一致エラーが発生する理由です。

マトリックスに結合できるデータクラスに応じて

imPos = [anaInfoSat(N).ImagePositionPatient];

またはセル配列に

imPos = {anaInfoSat(N).ImagePositionPatient};

同時に複数の変数に割り当てることもできます。

[A2X, A2Y, A2Z] = anaInfoSat(1:3).ImagePositionPatient;
于 2013-02-27T18:52:56.940 に答える