R中的Radon变换

Question

我正在寻找R中的Radon变换的实现。我所知道的唯一一个是在PET包中。但是，它没有给我我期望的结果。我正在使用wikipedia page

中的示例

require(png)
require(EBImage)
require(RCurl)
require(PET)

im <- readPNG( getBinaryURL( "http://upload.wikimedia.org/wikipedia/en/thumb/e/e5/Shepp_logan.png/170px-Shepp_logan.png" ) )[,,1]
rad = radon(im)$rData
# Normalize intensity values from 0-1
rad = normalize(rad)

display(t(im))
display(t(rad))

我得到这张图片：

当我期待：

我是否错误地使用了该功能的参数？我不太熟悉氡变换的细节

Answer 1

仔细观察数字：两个数字相同，只是偏移了90度。
回想一下Radon变换只是在不同角度上进行线积分，无论它是否转置都没关系。
另外，PET包装不见了，cran中似乎没有其他ra功能。
要安装最新的PET，您可以执行以下操作：

devtools::install_github('cran/PET')