**問題描述：**輸入nifti圖像，截取圖像的一部分并輸出，比如截取圖像的101010這一塊，并導出為nii文件

inputFile = 'D:\aa\dcm\input.nii'; % 輸入文件路徑subsetSize = [10 10 10]; % 截取的圖像塊大小
subsetStart = [1 1 1]; % 截取的起始位置
% 讀取NIfTI圖像
img = niftiread(inputFile);% 獲取圖像尺寸
imgSize = size(img);% 檢查截取區(qū)域是否超出圖像邊界
if any(subsetSize + subsetStart > imgSize)error('Subset exceeds image boundaries.');
end% 截取圖像
subset = img(subsetStart(1):subsetStart(1)+subsetSize(1)-1, ...subsetStart(2):subsetStart(2)+subsetSize(2)-1, ...subsetStart(3):subsetStart(3)+subsetSize(3)-1);% 保存導出,導出為nifti
niiFile = 'output.nii';  
niftiwrite(subset,niiFile);

3D展示：

itk-snap – edit – 3D Panel – Toggle Volume Rendering V