7

R の画像 (行列として表される) を、原点が 0,0 (左上隅) の極座標空間に変換しようとしています。次のような215x215行列が与えられます。x

ここに画像の説明を入力

x0 = as.vector(col(x))
y0 = as.vector(行(x))

r = sqrt( (x0^2) + (y0^2) )#x
a = atan(y0/x0)#y
m = as.matrix(data.frame(y=a, x=r))

m = ラウンド (m)
m[m>215] = NA
m[m==0] = NA

xn = x[m]
xn = 行列(xn, 215, 215)

ただし、 xn は次のようになります。

ここに画像の説明を入力

私がこれを期待するとき:

ここに画像の説明を入力

私が間違っていることは何ですか?

4

1 に答える 1

10

角度に問題があります:atan角度をラジアンで返します。丸めてしまうと、あまり情報が残らない…

試してみてください:

a = atan(y0/x0) * 215 / (pi/2)

変形画像

これは明らかに逆変換であり、中心が画像の中央にあります。

# Load the image
library(png)
library(RCurl)
d <- readPNG( getBinaryURL( "http://i.stack.imgur.com/rMR3C.png" ) )
image(d, col=gray(0:255/255))

# Origin for the polar coordinates
x0 <- ncol(d)/2
y0 <- nrow(d)/2

# The value of pixel (i,j) in the final image 
# comes from the pixel, in the original image, 
# with polar coordinates (r[i],theta[i]).
# We need a grid for the values of r and theta.
r <- 1:ceiling(sqrt( max(x0,nrow(d))^2 + max(y0,ncol(d))^2))
theta <- -pi/2 + seq(0,2*pi, length = 200)
r_theta <- expand.grid(r=r, theta=theta)

# Convert those polar coordinates to cartesian coordinates:
x <- with( r_theta, x0 + r * cos(theta) )
y <- with( r_theta, y0 + r * sin(theta) )
# Truncate them
x <- pmin( pmax( x, 1 ), ncol(d) )
y <- pmin( pmax( y, 1 ), nrow(d) )

# Build an empty matrix with the desired size and populate it
r <- matrix(NA, nrow=length(r), ncol=length(theta))
r[] <- d[cbind(x,y)]
image(r, col=gray(0:255/255))
于 2013-05-27T17:20:31.797 に答える