Procédé d’une estimation de positions de sources irradiantes dans un champ d’observation (Ω, Ω2), à l’aide d’une gamma caméra (1), la gamma caméra comportant : un matériau de détection (11) ; des pixels (12), répartis selon une surface de détection du matériau de détection, chaque pixel étant configuré pour former un signal de détection sous l’effet d’une détection d’une interaction d’un photon ionisant dans le matériau de détection ; un collimateur disposé face au matériau de détection, et définissant un premier champ d’observation (Ω1), inclus dans le champ d’observation de la gamma caméra ; une unité de localisation (14), configurée pour assigner une position à chaque interaction dans le matériau de détection, à partir du signal de détection généré par ladite interaction; une unité de spectrométrie (15), configurée pour assigner une énergie à chaque interaction détectée, à partir du signal de détection généré par ladite interaction ; une unité de coïncidence (16), configurée pour identifier chaque interaction comme une interaction simple ou une interaction multiple, une interaction multiple étant formées de deux interactions détectées en coïncidence temporelle ; une mémoire (17), configurée pour mémoriser des paramètres d’interactions détectées; le procédé comportant les étapes suivantes: a) acquisition de signaux de détection par les pixels (12) durant une période d’acquisition, chaque signal de détection étant associé à une interaction détectée ; b) en fonction des signaux de détection, assignation d’une position et d’une énergie ( ) à chaque interaction détectée ; c) identification de chaque interaction détectée comme étant une interaction simple ou une interaction multiple ; d) mémorisation, par la mémoire, des paramètres de chaque interaction, les paramètres de chaque interaction comportant la position et l’énergie de ladite interaction ; le procédé étant caractérisé en ce qu’il comportant également e) sélection d’un isotope ( ) parmi une listes d’isotopes ; f) en fonction de l’énergie assignée à chaque interaction, simple ou multiple, calcul d’une probabilité d’émission par l’isotope ( ) ; g) à partir des paramètres de chaque interaction, calcul d’une probabilité d’émission du photon ayant généré l’interaction, depuis différents points du champ d’observation ( ; h) prise en compte d’une probabilité a priori d’émission d’un photon pour l’isotope sélectionné, depuis les différents points du champ d’observation ( ); i) en fonction des probabilités résultant de g) et h), calcul d’une distribution spatiale de l’intensité d’émission dans le premier champ d’observation pour l’isotope sélectionné lors de l’étape e). Procédé selon la revendication 1, dans lequel les étapes h) et i) sont effectuées de façon itérative, de telle sorte qu’à chaque itération : l’étape h) comporte une prise en compte d’une distribution spatiale de l’intensité d’émission de l’isotope initiale ou résultant d’une itération précédente ; l’étape i) comporte une mise à jour de la distribution spatiale de l’intensité d’émission pour l’isotope sélectionné en fonction : de la probabilité a priori d’émission résultant de l’étape h) de la même itération ; des probabilités résultant des étapes f) et g). Procédé selon l’une quelconque des revendications précédentes, dans lequel lors de l’étape e), l’isotope est sélectionné en fonction des énergies respectivement assignées à différentes interactions détectées ; ou sélectionné arbitrairement. Procédé selon l’une quelconque des revendications précédentes, dans lequel les étapes f) à i) sont effectuées successivement pour différents isotopes sélectionnés. Procédé selon l’une quelconque des revendications précédentes, dans lequel l’étape f) est mise en œuvre en calculant : une probabilité ( ) que le photon ayant généré l’interaction soit un photon direct, émis par l’isotope, selon une énergie d’émission de l’isotope et ayant atteint le matériau de détection en conservant l’énergie d’émission, sans traverser une partie absorbante du collimateur ; une probabilité ( ) que le photon ayant généré l’interaction soit un photon indirect, émis par l’isotope, selon une énergie d’émission de l’isotope, et ayant atteint le matériau de détection avec une énergie inférieure à l’énergie d’émission et/ou en traversant une partie absorbante du collimateur. Procédé selon la revendication 5, dans lequel l’étape f) comporte une prise en compte, pour chaque énergie d’émission de l’isotope : une distribution de probabilités de détection d’une interaction, dans le matériau détecteur, selon différentes énergies, en supposant que le photon ayant généré l’interaction est un photon direct ; une distribution de probabilités de détection d’une interaction, dans le matériau détecteur, selon différentes énergies, en supposant que le photon ayant généré l’interaction est un photon indirect ; les distributions de probabilités de détection étant préalablement calculées en considérant une interaction simple et une interaction multiple. Procédé selon l’une quelconque des revendications précédentes, dans lequel l’étape g) comporte : lorsque l’interaction est simple, sélection d’un premier modèle de réponse spatiale de la caméra, la probabilité de position d’émission depuis différents points du premier champ d’observation ( ) étant déterminée en fonction du premier modèle de réponse spatiale et de la position de l’interaction dans le matériau détecteur ; lorsque l’interaction est multiple, sélection d’un deuxième modèle de réponse spatiale de la caméra, la probabilité de position d’émission depuis différents points du champ d’observation ( ) étant déterminée en fonction du deuxième modèle de réponse spatiale et de l’énergie ainsi que la position des interactions détectées en coïncidence temporelle formant l’interaction multiple dans le matériau détecteur. Procédé selon la revendication 7, dans lequel le premier modèle de réponse spatiale détermine, pour différentes positions bidimensionnelles d’interactions simples ( ), établies parallèlement à la surface de détection, une probabilité de détection d’un photon émis à partir de différents points du premier champ d’observation (Ω1), chaque interaction simple étant une interaction non considérée comme multiple. Procédé selon la revendication 7, dans lequel le deuxième modèle de réponse spatiale détermine, pour différentes positions tridimensionnelles d’interactions détectées en coïncidence temporelle , formant l’interaction multiple, et pour différentes valeurs d’énergie desdites interactions ( , une probabilité de détection d’un photon émis à partir de différents points du champ d’observation (Ω). Gamma caméra (1), destinée à détecter une présence de sources irradiantes (5) dans un champ d’observation (Ω),la gamma caméra comportant : un matériau de détection (11) ; des pixels (12), répartis selon une surface de détection du matériau de détection, chaque pixel étant configuré pour former un signal de détection sous l’effet d’une détection d’une interaction d’un photon ionisant dans le matériau de détection ; un collimateur disposé face au matériau de détection, et définissant un premier champ d’observation (Ω1), inclus dans le champ d’observation de la gamma caméra ; une unité de localisation (14), configurée pour assigner une position à chaque interaction dans le matériau de détection, à partir du signal de détection généré par ladite interaction; une unité de spectrométrie (15), configurée pour assigner une énergie à chaque interaction détectée, à partir du signal de détection généré par ladite interaction ; une unité de coïncidence (16), configurée pour identifier chaque interaction comme une interaction simple ou une interaction multiple, une interaction multiple étant formées de deux interactions de diffusion détectées en coïncidence temporelle ; une mémoire (17), configurée pour mémoriser des paramètres d’interactions détectées; une unité de traitement (18), configurée pour traiter les interactions mémorisées dans la mémoire, l’unité de traitement étant configurée pour mettre en œuvre les étapes e) à i) d’un procédé selon l’une quelconque des revendications précédentes.