Lissage d'image dans R

Question

Comment est-ce que je peux lisser cette image dans R, de sorte qu'il ne reste que deux pics?

Si ce serait des données 1d, je ferais une moyenne en cours d'exécution ou d'ajuster une fonction de régression. Mais je n'ai pas trouvé d'informations très spécifiques sur l'application de ces méthodes sur une matrice 2D. Par exemple, j'ai essayé d'utiliser filter() à partir du package stats.

J'ai aussi pensé à kriging, mais c'est plutôt une question d'interpolation, n'est-ce pas?

Answer 1

Le paquet spatstat détient une fonction blur() que applicates un flou gaussien. Cela étouffe l'image d'une manière, que la plupart du bruit disparaît et que les deux pics principaux sont clairement distincts.

L'effet peut être vu dans l'image suivante et est assez remarquable, en particulier dans le tracé 3D.

Le code pour générer l'image était:

library(jpeg) 
library(fields) 
library(spatstat) 

picture <- readJPEG("~/Downloads/spectrogram.png.jpeg") 
picture2 <- as.matrix(blur(as.im(picture), sigma=6)) 

layout(matrix(c(1:4), nrow=2)) 
image.plot(picture, col=gray.colors(50), main="original image", asp=1) 
image.plot(picture2, col=gray.colors(50), main="blurred with sigma = 6", asp=1) 
drape.plot(1:nrow(picture), 1:ncol(picture), picture, border=NA, theta=0, phi=45, main="original spectrogram") 
drape.plot(1:nrow(picture), 1:ncol(picture), picture2, border=NA, theta=0, phi=45, main="blurred with sigma = 6") 

Answer 2

Je pense que vous devriez jeter un oeil à la fonction focal dans le package de trame. Par exemple (copié à partir de la documentation raster):

r <- raster(ncols=36, nrows=18, xmn=0) 
r[] <- runif(ncell(r)) 
# 3x3 mean filter 
r3 <- focal(r, w=matrix(1/9,nrow=3,ncol=3)) 

La documentation comprend plus de détails.

Answer 3

Vous voulez certainement jeter un coup d'œil au paquet EBImage. Il y a plusieurs fonctions pour lisser votre image.

Par exemple, un filtre médian:

# Load EBImage up 
require(EBImage) 
# Read in your image 
im = readImage('/path/to/your/image') 
# Apply a median filter with window size of 7 
im.med = medianFilter(im, size=7) 
# Display image 
display(im.med) 

ou vous pouvez essayer un flou gaussien:

# Apply a gaussian blur with sigma=4 
im.gaus = gblur(im, sigma=4) 
# Display image 
display(im.gaus) 

Hope this helps!