将笛卡尔图像转换为极坐标、外观差异

Question

我试图对下面的第一幅图像进行极坐标变换，最后得到第二幅。然而，我的结果是第三张图片。我有一种感觉，这与我选择的“原产地”有关，但我不确定。

  radius = sqrt(width**2 + height**2)
  nheight = int(ceil(radius)/2)
  nwidth = int(ceil(radius/2))
 for y in range(0, height):
   for x in range(0, width):
     t =  int(atan(y/x))
     r = int(sqrt(x**2+y**2)/2)
     color = getColor(getPixel(pic, x, y))
     setColor( getPixel(radial,r,t), color)

Answer 1

有一些差异/错误：

1. 他们使用图像的中心作为原点
2. 它们适当地缩放轴。在您的示例中，您绘制的是角度（在0和π之间），而不是使用图像的整个高度。
3. 您使用了错误的atan函数（atan2在这种情况下工作得更好：））
4. 不太重要，但你不必要地舍入了很多，这会稍微降低准确性，会减慢速度。

这是结合我建议的改进的代码。它不是非常高效，但它应该有望工作：）

  maxradius = sqrt(width**2 + height**2)/2
  rscale = width / maxradius
  tscale = height / (2*math.pi)
  for y in range(0, height):
   dy = y - height/2
   for x in range(0, width):
     dx = x - width/2
     t =  atan2(dy,dx)%(2*math.pi)
     r = sqrt(dx**2+dy**2)
     color = getColor(getPixel(pic, x, y))
     setColor( getPixel(radial,int(r*rscale),int(t*tscale)), color)

具体而言，它通过以下方式解决了上述问题：

1. 我们使用dx=x-width/2作为距离中心的度量，与dy类似。然后，在整个计算过程中，我们使用它们来代替x和y

有问题尽管问！