A processor-implemented method of contouring three-dimensional medical image data, the method comprising:
- receiving, via a first interface, source data in the three-dimensional medical image data, which source data represent an anatomic structure of interest that is to be contoured, and tissues adjoining the anatomic structure of interest;
- selecting a first subset of the source data defining a first two-dimensional graphic representation of the source data, the first subset being arranged in a first image slice through the anatomic structure of interest at a first axial position in the source data;
- outputting, via a second interface, the first subset of the source data for presentation on a graphical display;
- receiving a first set of instructions in a third interface in the form of user commands;
- identifying, in response to the first set of instructions received via the third interface, a first edge of the anatomic structure of interest in the first image slice, and after having received the first set of instructions:
- selecting a second subset of the source data defining a second two-dimensional graphic representation of the source data, the second subset being arranged in a second image slice through the anatomic structure of interest at a second axial position in the source data, which second axial position is different from the first axial position;
- outputting, via the second interface, the second subset of the source data for presentation on the graphical display;
- identifying, in response to a second set of instructions received via the third interface, a second edge of the anatomic structure of interest in the second image slice;
- calculating a three-dimensional shell based on the first and second edges and the source data, the three-dimensional shell representing an approximation of a delimiting surface of the anatomic structure of interest;
- calculating an initial estimate of the three-dimensional shell between the first and second axial positions in the source data by interpolating between first and second sets of data points, the first set of data points describing a first outline of the first edge and the second set of data points describing a second outline of the second edge;
- analyzing a first data portion of the first subset of the source data in proximity to the first edge;
- analyzing a second data portion of the second subset of the source data in proximity to the second edge, and based thereon
- deriving a set of image features defining an estimated perimeter of the anatomic structure of interest, and
- calculating the subsequent estimate of the three-dimensional shell on the further basis of the set of image features;
- and
- interpolating the set of image features defining the estimated perimeter over the initial estimate by using harmonic functions, and where image features derived from the first and second data portions represent boundary conditions.