arrays - 2Dnumpy配列内に長方形を描画する

Question

センサーの各ピクセルからの個々のデータを含む2Dnumpy配列があります。画像は、カメラからのライブフィードを使用してGUIに表示されます。画面の領域を区別するために、画像上に長方形を描画できるようにしたいと思います。画像の側面に平行な長方形を描くのはとても簡単なようですが、最終的には長方形を回転できるようにしたいと思います。長方形を回転させたときに、長方形がどのピクセルをカバーしているかをどのように知ることができますか？

Answer 1

依存関係を気にしない場合は、PythonImagingLibraryを使用できます。2D numpy配列data、およびpolyポリゴン座標の配列（形状（n、2））が与えられた場合、これにより、配列に値0で満たされたポリゴンが描画されます。

img = Image.fromarray(data)
draw = ImageDraw.Draw(img)
draw.polygon([tuple(p) for p in poly], fill=0)
new_data = np.asarray(img)

これが自己完結型のデモです。

import numpy as np
import matplotlib.pyplot as plt

# Python Imaging Library imports
import Image
import ImageDraw


def get_rect(x, y, width, height, angle):
    rect = np.array([(0, 0), (width, 0), (width, height), (0, height), (0, 0)])
    theta = (np.pi / 180.0) * angle
    R = np.array([[np.cos(theta), -np.sin(theta)],
                  [np.sin(theta), np.cos(theta)]])
    offset = np.array([x, y])
    transformed_rect = np.dot(rect, R) + offset
    return transformed_rect


def get_data():
    """Make an array for the demonstration."""
    X, Y = np.meshgrid(np.linspace(0, np.pi, 512), np.linspace(0, 2, 512))
    z = (np.sin(X) + np.cos(Y)) ** 2 + 0.25
    data = (255 * (z / z.max())).astype(int)
    return data


if __name__ == "__main__":
    data = get_data()

    # Convert the numpy array to an Image object.
    img = Image.fromarray(data)

    # Draw a rotated rectangle on the image.
    draw = ImageDraw.Draw(img)
    rect = get_rect(x=120, y=80, width=100, height=40, angle=30.0)
    draw.polygon([tuple(p) for p in rect], fill=0)
    # Convert the Image data to a numpy array.
    new_data = np.asarray(img)

    # Display the result using matplotlib.  (`img.show()` could also be used.)
    plt.imshow(new_data, cmap=plt.cm.gray)
    plt.show()

このスクリプトは、次のプロットを生成します。