kphilipp

L = 400;
D = 80;
speckleSize = L/D
img = zeros(L, L);
pad = L/2; % padding
for k=pad-D:D+pad
  for l=pad-D:D+pad
    if abs((pad-k)^2+(pad-l)^2) < D^2/4
      img(k,l) = exp(unifrnd(-pi,pi));
    end

kphilipp

# further information:
# * http://stackoverflow.com/questions/11114349/how-to-visualize-descriptor-matching-using-opencv-module-in-python
# * http://docs.opencv.org/doc/tutorials/features2d/feature_homography/feature_homography.html#feature-homography
# * http://stackoverflow.com/questions/9539473/opencv-orb-not-finding-matches-once-rotation-scale-invariances-are-introduced
# * OpenCV 2 Computer Vision Application Programming Cookbook, Chapter 9
import cv2
import scipy as sp
import numpy as np