A system for change detection in medical images, comprising: - an image data interface configured to access a first medical image and a second medical image ; - a processor configured to: - generate a difference image representing intensity differences between the first medical image and the second medical image ; - determine an intensity distribution of the difference image ; - fit a mixture model to the intensity distribution to identify a plurality of probability distributions 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 ; 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 and a second medical image ; - a processor configured to: - generate a difference image representing intensity differences between the first medical image and the second medical image ; - determine an intensity distribution of the difference image ; - fit a mixture model to the intensity distribution to identify a plurality of probability distributions 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 ; 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 and the second medical image ; - determine an intensity distribution of the difference image ; - fit a mixture model to the intensity distribution to identify a plurality of probability distributions 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 ; 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. <EOS>