質問: 半径r1, 中心座標 (x1, y1, z1) の球体の内側のピクセルの値をゲットしたい。 3d画像 (xyz)内の半径r1, 中心座標 (x1, y1, z1) の球体の内側のピクセルの値をゲットしたいです。 個々のピクセルの値には興味はなく、球体内のピクセルの色強度の平均値・sdといった統計量をゲットしたいです。 球体は1つではなく数百以上あるため、効率的に定量する方法を求めています。 教えて下さい。 半径 (r1, r2, r3 ...)と中心座標 ((x1, y1, z1), (x2, y2, z2), (x3, y3, z3) ...)の値は持っています。 私の説明が不十分であればお尋ねください。 よろしくお願いします。
回答: clear;clc;close all load mri d = squeeze(d); volshow(d) [x, y, z]=meshgrid([1:siz(1)],[1:siz(2)],[1:siz(3)]); r1=10; cx=100; cy=50; cz=13; % obtain the index of the voxel, where (x-cx)^2+(y-cy)^2+(z-cz)^2