A method for determining metastases within an anatomical body part, wherein the metastases exhibits an enhanced signal in an image of the anatomical body part generated by a medical imaging method using a contrast agent, the method being designed to be performed by a computer and comprising the following steps: acquiring an enhanced image (T1w+c) of the anatomical body part obtained by a medical imaging method using the contrast agent; acquiring a non-enhanced image (T1w-c) of the anatomical body part obtained by a medical imaging method without using the contrast agent, wherein both the enhanced image and the non-enhanced image include two-dimensional layers of voxels positioned one above the other; spatially correlating the enhanced image (T1w+c) and the non-enhanced image (T1w-c) with respect to each other; acquiring non-enhanced sub-regions in each of the enhanced image (T1w+c) and the non-enhanced image (T1w-c), wherein the non-enhanced sub-regions are assumed or known to not exhibit an enhanced signal; determining a normalized difference image by normalizing the spatially correlated enhanced (T1w+c) and non-enhanced images (T1w-c) with respect to each other on the basis of the intensities in the respective non-enhanced sub-regions and by performing a difference operation on the normalized spatially correlated enhanced (T1w+c) and non-enhanced images (T1w-c); acquiring an intensity condition for differentiating between the structure of interest and other structures which do not exhibit an enhanced signal; determining a three-dimensional binary image on the basis of the intensity condition and the normalized difference image, wherein the binary image comprises a plurality of image layers comprising voxels representing either an enhanced signal intensity or a non-enhanced signal intensity; wherein the method further comprises the steps of: acquiring metastasis geometry data which describe geometrical properties of a metastasis in the image layers; determining a set of adjacent voxels which represent the enhanced signal intensity and are integrally closed and apart from other sets and/or voxels which represent the enhanced signal intensity; analyzing the 2D properties of the set; and determining that the set represents a metastasis if the set exhibits two-dimensional geometrical properties which comply with the geometrical properties described by the metastasis geometry data, characterized in that the geometrical properties include relative geometrical properties which describe the relative variation in the geometry of the metastasis of one image layer with respect to another image layer, and wherein the two-dimensional properties of the set are determined to comply with the geometrical properties described by the metastatis geometry data if the center of a circle representing a borderline of the set deviates from the one image layer to the other image layer within a predetermined range of variation. acquiring an enhanced image (T1w+c) of the anatomical body part obtained by a medical imaging method using the contrast agent; acquiring a non-enhanced image (T1w-c) of the anatomical body part obtained by a medical imaging method without using the contrast agent, wherein both the enhanced image and the non-enhanced image include two-dimensional layers of voxels positioned one above the other; spatially correlating the enhanced image (T1w+c) and the non-enhanced image (T1w-c) with respect to each other; acquiring non-enhanced sub-regions in each of the enhanced image (T1w+c) and the non-enhanced image (T1w-c), wherein the non-enhanced sub-regions are assumed or known to not exhibit an enhanced signal; determining a normalized difference image by normalizing the spatially correlated enhanced (T1w+c) and non-enhanced images (T1w-c) with respect to each other on the basis of the intensities in the respective non-enhanced sub-regions and by performing a difference operation on the normalized spatially correlated enhanced (T1w+c) and non-enhanced images (T1w-c); acquiring an intensity condition for differentiating between the structure of interest and other structures which do not exhibit an enhanced signal; determining a three-dimensional binary image on the basis of the intensity condition and the normalized difference image, wherein the binary image comprises a plurality of image layers comprising voxels representing either an enhanced signal intensity or a non-enhanced signal intensity; acquiring metastasis geometry data which describe geometrical properties of a metastasis in the image layers; determining a set of adjacent voxels which represent the enhanced signal intensity and are integrally closed and apart from other sets and/or voxels which represent the enhanced signal intensity; analyzing the 2D properties of the set; and determining that the set represents a metastasis if the set exhibits two-dimensional geometrical properties which comply with the geometrical properties described by the metastasis geometry data,