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 poids F dans laquelle la distance entre pixels est une variable et une fonction poids H dans laquelle la différence entre les valeurs de pixel de pixels voisins est une variable sont établies dans un filtrage qui utilise un filtre bilatéral (étape S4). Ici, l'image PET à traiter n'est pas utilisée ； à la place, une image CT est utilisée à titre d'autre image numérique pour établir la fonction poids H, qui donne un poids dépendant de l'intensité des contours (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 à titre d'autre image numérique pour déterminer des coefficients de filtre (étape S6). Les coefficients de filtre déterminés sont utilisés dans un filtrage de l'image PET, qui est l'image numérique à traiter (étape S7). En conséquence, un filtrage peut être effectué sans subir l'influence du niveau de bruit de l'
