4

我正在研究关于 Sobel 边缘检测的 R 任务。不幸的是，我关注的视频教程将 R 用于其他所有任务，但切换到 python 进行图像处理 - 我猜他没有找到任何用于图像卷积类型工作的有用 R 包（本教程来自去年）。我已经尝试过EBImage和magick（这个似乎是新的），但没有找到太多。这个魔法小插图谈论image_convolve('Sobel')（大约在页面的一半处）但仅适用于垂直边缘，而不是水平边缘。有人可以推荐一些我可以使用的好材料吗？我对图像处理相当陌生。

更新：

我已经设法使用magick包（粘贴在下面的代码）分别检测垂直和水平边缘，但不知道如何将它们组合起来生成单个图像

library(magick)

# get image
img <- image_read("https://www.r-project.org/logo/Rlogo.png")
print(image_info(img))

# define horizontal and vertical Sobel kernel
Shoriz <- matrix(c(1, 2, 1, 0, 0, 0, -1, -2, -1), nrow = 3)
Svert <- t(Shoriz)

# get horizontal and vertical edges
imgH <- image_convolve(img, Shoriz)
imgV <- image_convolve(img, Svert)

print(plot(as.raster(img)))     # view original image
print(plot(as.raster(imgH)))    # view horizontal edges
print(plot(as.raster(imgV)))    # view vertical edges

在教程中，接下来我需要组合imgH并imgV计算这些边缘之间的欧几里得距离，但dist()不适用于图像对象本身。我需要从这些图像中获取数据，但不知道如何。类似于imageData()in package 的东西会有所帮助，但在.It 中EBImage找不到它。它有功能，但它的输出看起来很复杂。magickimage_data()

更新（2）：

我（希望）得到了我想要的EBImage包（下面的代码）。magick一旦我弄清楚如何从如上所述的边缘图像中获取像素数据，以及如何将最终的边缘数据转换回图像，我仍然希望它与包一起使用。

library(EBImage)

# get image
img <- readImage("https://www.r-project.org/logo/Rlogo.png")
print(img, short = T)

# define horizontal and vertical Sobel kernel
Shoriz <- matrix(c(1, 2, 1, 0, 0, 0, -1, -2, -1), nrow = 3)
Svert <- t(Shoriz)

# get horizontal and vertical edges
imgH <- filter2(img, Shoriz)
imgV <- filter2(img, Svert)

# combine edge pixel data to get overall edge data
hdata <- imageData(imgH)
vdata <- imageData(imgV)
edata <- sqrt(hdata^2 + vdata^2)

# transform edge data to image
imgE <- Image(edata, colormode = 2)
print(display(combine(img, imgH, imgV, imgE), method = "raster", all = T))

谢谢。

4

0 回答 0