0

Mathematicaのコードを理解しようとしていますが、その構文に精通していません。ここの論理に頭を悩ませているようには見えません。

anamorph3[img_, angle_: 270 Degree, imgWidth_: 512] :=
Module[{data, f, matrix, dim, rOuter, rInner = 1.},
dim = ImageDimensions[img];
rOuter = rInner (1 + angle #2/#1 & @@ dim);
data = Table[
ListInterpolation[#[[All, All, i]], 
{{rOuter, rInner}, {-angle/2, angle/2}}], {i, 3}] &@ImageData[img];
f[i_, j_] := If[Abs[j] <= angle/2 && rInner <= i <= rOuter, 
Through[data[i, j]], {1., 1., 1.}];
Image@Table[f[Sqrt[i^2 + j^2], ArcTan[i, -j]], 
{i, -rOuter, rOuter, 2 rOuter/(imgWidth - 1)},
{j, -rOuter, rOuter, 2 rOuter/(imgWidth - 1)}]]

このコードの詳細については、こちらをご覧ください。

ImageTransformationでアナモフィックバージョンの画像を生成する方法。

誰かがこのための擬似コード/Pythonコードを書くことができれば素晴らしいでしょう。

4

1 に答える 1

1

私はそれが次のようなものになると思います：

def anamorph3(img, angle = math.radians(270), imgWidth = 512):
    rInner = 1.
    dim = ImageDimensions(img)
    rOuter = rInner * (1 + angle * dim[2]/dim[1])
    data = list()
    for i in range(3):
        channelData = [[x[i] for x in v] for v in ImageData(img)]
        data.append(ListInterpolate(channelData, [[rOuter, rInner],[-angle/2,angle/2]]))
    def f(i, j):
        if abs(j) <= angle/2 && rInner <= i <= rOuter:
            l = list()
            for fun in data:
                l.append(fun(i, j))
            return l
        else:
            return [1, 1, 1]
    newImageData = list()
    for i in range(-rOuter, rOuter, 2*rOuter/(imgWidth-1)):
        l = list()
        for j in range(-rOuter, rOuter, 2*rOuter/(imgWidth-1)):
            l.append(f(math.sqrt(i**2 + j**2), math.atan(i, -j)))
        newImageData.append(l)
    return Image(newImageData)

ここでImageDimensionsは、画像の寸法をImageData返し、0-1の範囲に正規化された各ピクセルの各カラーチャネルのカラーデータを含む3Dマトリックスを返し、ListInterpolate指定されたxおよびy座標範囲内の数値のセット間を補間する関数を作成します。 。Image3Dマトリックスを取得し、RGB画像に変換します。

于 2013-01-20T23:56:53.990 に答える