查找3D橢球matlab

Question

我在Matlab中有以下問題：我有一個二維三維體積（NxNxN）對應於3D ROI，我想要找到最適合此ROI的橢球。我找到了中心的三維座標，三維半球的長度以及擬合橢球的三維角度。結果矢量類似於：[195.1126 169.3114 62.4193 28.2725 23.0191 9.3104 52.9536 -5.7639 1.7606]其中三個第一個元素代表中心座標，元素編號4,5,6代表半徑半徑，元素編號爲7,8,9分別代表x，y，z中以度數表示的旋轉。現在我想要找到一個二進制3D體（如3D矩陣：NxNxN，原始ROI的大小相同），其中0不在橢圓體中，1是橢圓體。 任何人都可以幫我解決這個問題嗎？ 非常感謝。

Answer 1

試試這個：的解釋是在評論

 V = [195.1126 169.3114 62.4193,28.2725 23.0191 9.3104,52.9536 -5.7639 1.7606]; 
c = V(1:3); % The centroid 
l = V(4:6); % Axis Length 
phi = -V(7:9); % Angles 
E = zeros(512,512,85); % initialize Elipse matrix 
R = rotx(phi(1))*roty(phi(2))*rotz(phi(3)); % Create rotation matrix around phi 
%% 
[X,Y,Z] = meshgrid(1:size(E,2),1:size(E,1),1:size(E,3)); % Create a Grid 
XYZr = ([X(:)-c(1),Y(:)-c(2),Z(:)-c(3)])*R; % rotate Grid and centralize around centroid 
%% Calculate the elipse equation 
inx = (((XYZr(:,1)))./l(1)); % 
iny = (((XYZr(:,2)))./l(2)); 
inz = (((XYZr(:,3)))./l(3)); 
ind = [inx,iny,inz]; 

E(:) = sum(ind.^2,2)<1; 
%% Display 
figure; imagesc(E(:,:,62)); 
figure; imagesc(squeeze(E(:,195,:))); 
figure; imagesc(squeeze(E(170,:,:)));