Levigatura immagine in R

Question

Come posso lisciare questa immagine in R, in modo che rimangano solo due picchi?

Se si trattasse di dati 1d, eseguirei una media di esecuzione o adattare una funzione di regressione ad esso. Ma non ho trovato informazioni molto specifiche sull'applicazione di questi metodi su una matrice 2d. Ad esempio, ho provato a utilizzare filter() dal pacchetto stats.

Ho pensato anche al kriging, ma questo riguarda più l'interpolazione, vero?

Answer 1

Il pacchetto spatstat detiene una funzione blur() che applicates una sfocatura gaussiana. In questo modo l'immagine viene attenuata, la maggior parte del rumore scompare e i due picchi principali sono chiaramente distinguibili.

L'effetto può essere visto nella figura seguente ed è piuttosto notevole, specialmente nella trama 3D.

Il codice per generare l'immagine era:

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

penso che si dovrebbe avere uno sguardo alla funzione di focal nel pacchetto raster. Per esempio (copiato dalla documentazione 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 documentazione comprende ulteriori dettagli.

Answer 3

Si consiglia di dare un'occhiata al pacchetto EBImage. Ci sono più funzioni per lisciare l'immagine.

Ad esempio, un filtro mediano:

# 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) 

o si potrebbe provare una sfocatura gaussiana:

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

Spero che questo aiuti!