Wygładzanie obrazu w R

Question

Jak mogę wygładzić to zdjęcie w R, aby pozostały tylko dwa piki?

Gdyby to były dane 1d, zrobiłbym średnią średnią lub dopasowałem do niej funkcję regresji. Ale nie znalazłem bardzo konkretnych informacji na temat stosowania tych metod na matrycy 2d. Na przykład próbowałem użyć filter() z pakietu stats.

Też myślałem o kriging, ale to jest bardziej o interpolacji, prawda?

Answer 1

Pakiet spatstat posiada funkcję blur() że applicates rozmycie Gaussa. To sprawia, że ​​obraz w pewnym sensie znika, a dwa główne szczyty są wyraźnie widoczne.

Efekt można zobaczyć na poniższym obrazku i jest dość niezwykły, zwłaszcza w działce 3d.

Kod do generowania obrazu było:

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

myślę, że należy przyjrzeć się funkcji w pakiecie rastrowej focal. Na przykład (skopiowany z dokumentacją 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)) 

Dokumentacja zawierać więcej szczegółów.

Answer 3

Na pewno chcesz rzucić okiem na pakiet EBImage. Istnieje wiele funkcji do wygładzania obrazu.

Przykładowo, mediana filtr:

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

lub można spróbować rozmycie Gaussa:

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

nadzieję, że to pomaga!