How I can have the volume of a DICOM image

dose this line of this code is correct to heve the sum of all voxels in image?

if(exist('filename','var')==0) [filename, pathname] = uigetfile('*.nii', 'Read nii-file'); filename = [pathname filename]; end s1=sum(filename(:))

[Rik]

I would advise a Matlab tutorial (e.g. this one). This code will not implicitly read the file, so the sum is calculated over the name itself.

[Vania krm]

Thanks Rik can you explain a bit more, because I have no access to the Matlab tutorial.

[Rik]

There are many free online courses to learn Matlab.

That filename is a char array (a vector made up of letters), each letter has a value, so you can calculate the sum of a name: 'Hi!' is [72 105 33], so sum('Hi!') is 210.

[Vania krm]

Aha understand! Thanks with this code is correct? https://it.mathworks.com/matlabcentral/fileexchange/29344-read-medical-data-3d

function V = nii_read_volume(info) % function for reading volume of NifTi ( .nii ) volume file % nii_read_header(file-info) % % volume = nii_read_volume(file-header) % % examples: % 1: info = nii_read_header() % V = nii_read_volume(info); % imshow(squeeze(V(:,:,round(end/2))),[]); % % 2: V = nii_read_volume('test.nii');

if(~isstruct(info)) info=nii_read_header(info); end

% Open v3d file fid=fopen(info.Filename,'rb');

% Seek volume data start
datasize=prod(info.Dimensions)*(info.BitVoxel/8);
fsize=info.Filesize;
fseek(fid,fsize-datasize,'bof');

% Read Volume data
switch(info.DataTypeStr)
    case 'INT8'
      V = int8(fread(fid,datasize,'int8'));
    case 'UINT8'
      V = uint8(fread(fid,datasize,'uint8'));
    case 'INT16'
      V = int16(fread(fid,datasize,'int16'));
    case 'UINT16'
      V = uint16(fread(fid,datasize,'uint16'));
    case 'INT32'
      V = int32(fread(fid,datasize,'int32'));
    case 'UINT32'
      V = uint32(fread(fid,datasize,'uint32'));
    case 'INT64'
      V = int64(fread(fid,datasize,'int64'));
    case 'UINT64'
      V = uint64(fread(fid,datasize,'uint64'));    
    otherwise
      V = uint8(fread(fid,datasize,'uint8'));
end
fclose(fid);

% Reshape the volume data to the right dimensions V = reshape(V,info.Dimensions);

if I do : S=sume (V(:))

[Rik]

If you use that function to read your data, and that function has no bugs, then S=sum(V(:)); will give you the sum.

To correctly format code on this forum, select all the code you pasted and the click the {}Code button.

[Vania krm]

Rik you really help me. very very thanks and wish you best.