$d$次元Fourier変換を \begin{equation} \mathcal{F}u(\xi) = \frac{1}{(2\pi)^{d/2}}\int_{\R^d}\exp(- i x \cdot \xi) u(x) dx \end{equation}で定義する。 ここでは、2次元Fourier変換を考え、$u(x) = \chi_{[-1,1]^2}(x)$ つまり、 \begin{equation} \ch…

pyplot.imread でグレイスケールで読ませる方法を調べるのが面倒だったので、cv2 を使う。グレースケールだと、FFTとIFFT - らんだむな記憶で触れたような色相の問題も一応なさそう。 #! /usr/bin/env python # -*- coding: utf-8 -*- import sys, os import…

高速Fourier変換と逆高速Fourier変換を画像にかましてみる。 import sys import numpy as np from matplotlib import pyplot as plt img = plt.imread(sys.argv[1], 1) # FFT f = np.fft.fft2(img) # IFFT img_back = np.fft.ifft2(f) # complex -> uint8 im…

cv2.imread の動きを見てみたい。 画像とビデオの読み込みと書き込み — opencv v2.1 documentationによると、第2引数はグレースケール or 3チャンネルかを制御できる。I = [0,255]とする。 import sys import cv2 img = cv2.imread(sys.argv[1], 0) print im…

numpyの親戚でしょ？うふふふって $ pip install scipy したら終わらねーorzPythonで数値計算を始める（with Homebrew and pip） - Qiitaによると、結構時間がかかるらしい。 裏でコンパイラが爆走。10分以上経過した気がする...。 $ pip install scipy Coll…

OpenGLやWebGLも少し気になるなぁ、という気もしないでもないが、だるぅ～～～いのでOpenCVを入れてみる。 よく分かってないので、Install OpenCV 3 and Python 2.7+ on Ubuntu - PyimageSearchのまんまでビルドする。頭を使ったら負けだ！ $ sudo apt-get u…