In digital image processing according to the invention, a weight function F where the distance between pixels is a variable and a weight function H where the difference between the pixel values of neighboring pixels is a variable are set in filtering that uses a bilateral filter (step S4). Here, the PET image to be processed is not used； instead, a CT image is used as another digital image to set the weight function H, which gives a weight dependent on the edge intensity (the difference between the pixel values of the pixel of interest and an adjacent pixel) of the CT image. The CT image is also used in this manner as another digital image to determine filter coefficients (step S6). The determined filter coefficients are used in filtering of the PET image, which is the digital image to be processed (step S7). Therefore, filtering can be performed without being affected by the noise level of the PET image to be processed. As a result, spatial resolution is maintained while noise is reduced.Dans un traitement d'image numérique selon l'invention, une fonction de poids F dans laquelle la distance entre des pixels est une variable et une fonction de poids H dans laquelle la différence entre les valeurs de pixel de pixels voisins est une variable sont réglées dans un filtrage qui utilise un filtre bilatéral (étape S4). Dans la présente invention, l'image PET à traiter n'est pas utilisée ； à la place, une image CT est utilisée en tant qu'autre image numérique pour régler la fonction de poids H, ce qui donne un poids qui dépend de l'intensité de contour (la différence entre les valeurs de pixel du pixel d'intérêt et d'un pixel adjacent) de l'image CT. L'image CT est également utilisée de cette manière en tant qu'autre image numérique pour déterminer des coefficients de filtre (étape S6). Les coefficients de filtre déterminés sont utilisés dans le filtrage de l'image PET, qui est l'image numérique à traiter (étape S7). Par conséquent, un filtrage peut être réalisé sans être
