A system for change detection in medical images, comprising: - an image data interface configured to access a first medical image (022, 200) and a second medical image (024, 210); - a processor configured to: - generate a difference image representing intensity differences between the first medical image (022, 200) and the second medical image (024, 210); - determine an intensity distribution of the difference image; - fit a mixture model to the intensity distribution to identify a plurality of probability distributions (315, 320) which collectively model the intensity distribution; - determine a plurality of intensity ranges in the intensity distribution, wherein each one of the plurality of intensity ranges is determined as a function of a respective one of the plurality of probability distributions (315, 320); characterized in that each of the plurality of probability distributions represents a different type of change; and each one of the plurality of intensity ranges represents the different type of change; wherein the processor is further configured to - label image data of the difference image by determining into which of the plurality of intensity ranges said labeled image data falls. - an image data interface configured to access a first medical image (022, 200) and a second medical image (024, 210); - a processor configured to: - generate a difference image representing intensity differences between the first medical image (022, 200) and the second medical image (024, 210); - determine an intensity distribution of the difference image; - fit a mixture model to the intensity distribution to identify a plurality of probability distributions (315, 320) which collectively model the intensity distribution; - determine a plurality of intensity ranges in the intensity distribution, wherein each one of the plurality of intensity ranges is determined as a function of a respective one of the plurality of probability distributions (315, 320); characterized in that each of the plurality of probability distributions represents a different type of change; and each one of the plurality of intensity ranges represents the different type of change; wherein the processor is further configured to - label image data of the difference image by determining into which of the plurality of intensity ranges said labeled image data falls. - generate a difference image representing intensity differences between the first medical image (022, 200) and the second medical image (024, 210); - determine an intensity distribution of the difference image; - fit a mixture model to the intensity distribution to identify a plurality of probability distributions (315, 320) which collectively model the intensity distribution; - determine a plurality of intensity ranges in the intensity distribution, wherein each one of the plurality of intensity ranges is determined as a function of a respective one of the plurality of probability distributions (315, 320); characterized in that each of the plurality of probability distributions represents a different type of change; and each one of the plurality of intensity ranges represents the different type of change; wherein the processor is further configured to - label image data of the difference image by determining into which of the plurality of intensity ranges said labeled image data falls.