Làm mịn hình ảnh trong R

Question

Làm thế nào tôi có thể làm mịn hình ảnh này trong R, sao cho chỉ có hai đỉnh vẫn còn?

Nếu đây là dữ liệu 1d, tôi sẽ thực hiện hoạt động trung bình hoặc phù hợp với chức năng hồi quy. Nhưng tôi không tìm thấy thông tin rất cụ thể về việc áp dụng các phương thức này trên ma trận 2d. Ví dụ: tôi đã cố gắng sử dụng filter() từ gói stats.

Tôi cũng nghĩ về kriging, nhưng đây là nhiều hơn về nội suy, phải không?

Answer 1

Gói spatstat giữ một chức năng blur() rằng applicates một gaussian blur. Điều này làm mờ hình ảnh theo một cách, mà hầu hết các tiếng ồn biến mất và hai đỉnh chính là rõ ràng rõ ràng.

Hiệu ứng có thể được nhìn thấy trong hình dưới đây và khá đáng chú ý, đặc biệt là trong ô 3D.

Mã để tạo ra hình ảnh là:

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

Tôi nghĩ bạn nên có một cái nhìn tại các focal chức năng trong gói raster. Ví dụ: (được sao chép từ tài liệu 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)) 

Tài liệu này bao gồm thêm chi tiết.

Answer 3

Bạn chắc chắn muốn xem gói EBImage. Có nhiều chức năng để làm mịn hình ảnh của bạn.

Ví dụ, một bộ lọc trung bình:

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

hoặc bạn có thể thử một gaussian blur:

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

Hope this helps!