A computer-implemented method for lesion marking and quantitative analysis of nuclear medicine images of a human subject, the method comprising: (a) accessing, by a processor of a computing device, a bone scan image set for the human subject, said bone scan image set obtained following administration of an agent to the human subject; (b) automatically segmenting, by the processor, each image in the bone scan image set to identify one or more skeletal regions of interest, each skeletal region of interest corresponding to a particular anatomical region of a skeleton of the human subject, thereby obtaining an annotated set of images; (c) automatically detecting, by the processor, an initial set of one or more hotspots, each hotspot corresponding to an area of elevated intensity in the annotated set of images, said automatically detecting comprising: using ( i) intensities of pixels in the annotated set of images and ( ii) a plurality of preliminary threshold values to detect a set of potential hotspots; computing a global threshold scaling factor using the set of potential hotspots; adjusting the plurality of preliminary threshold values using the global threshold scaling factor, thereby obtaining a plurality of adjusted threshold values; and using ( i) intensities of pixels in the annotated set of images and ( ii) the plurality of adjusted threshold values to identify the initial set of hotspots; (d) for each hotspot in the initial set of hotspots, extracting, by the processor, a set of hotspot features associated with the hotspot; (e) for each hotspot in the initial set of hotspots, calculating, by the processor, a metastasis likelihood value corresponding to a likelihood of the hotspot representing a metastasis, based on the set of hotspot features associated with the hotspot; and (f) causing, by the processor, rendering of a graphical representation of at least a portion of the initial set of hotspots for display within a graphical user interface. (a) accessing, by a processor of a computing device, a bone scan image set for the human subject, said bone scan image set obtained following administration of an agent to the human subject; (b) automatically segmenting, by the processor, each image in the bone scan image set to identify one or more skeletal regions of interest, each skeletal region of interest corresponding to a particular anatomical region of a skeleton of the human subject, thereby obtaining an annotated set of images; (c) automatically detecting, by the processor, an initial set of one or more hotspots, each hotspot corresponding to an area of elevated intensity in the annotated set of images, said automatically detecting comprising: using ( i) intensities of pixels in the annotated set of images and ( ii) a plurality of preliminary threshold values to detect a set of potential hotspots; computing a global threshold scaling factor using the set of potential hotspots; adjusting the plurality of preliminary threshold values using the global threshold scaling factor, thereby obtaining a plurality of adjusted threshold values; and using ( i) intensities of pixels in the annotated set of images and ( ii) the plurality of adjusted threshold values to identify the initial set of hotspots; using ( i) intensities of pixels in the annotated set of images and ( ii) a plurality of preliminary threshold values to detect a set of potential hotspots; computing a global threshold scaling factor using the set of potential hotspots; adjusting the plurality of preliminary threshold values using the global threshold scaling factor, thereby obtaining a plurality of adjusted threshold values; and using ( i) intensities of pixels in the annotated set of images and ( ii) the plurality of adjusted threshold values to identify the initial set of hotspots; (d) for each hotspot in the initial set of hotspots, extracting, by the processor, a set of hotspot features associated with the hotspot; (e) for each hotspot in the initial set of hotspots, calculating, by the processor, a metastasis likelihood value corresponding to a likelihood of the hotspot representing a metastasis, based on the set of hotspot features associated with the hotspot; and (f) causing, by the processor, rendering of a graphical representation of at least a portion of the initial set of hotspots for display within a graphical user interface.