Roberts Kreuz ist etwas knifflig, weil es keine ungerade Größe hat (2x2 statt 3x3 oder 5x5). Ich habe es mit numpy + scipy mit einer gepolsterten 3x3-Faltungsmaske gemacht.
import sys
import numpy as np
from scipy import ndimage
import Image
roberts_cross_v = np.array( [[ 0, 0, 0 ],
[ 0, 1, 0 ],
[ 0, 0,-1 ]] )
roberts_cross_h = np.array( [[ 0, 0, 0 ],
[ 0, 0, 1 ],
[ 0,-1, 0 ]] )
def load_image( infilename ) :
img = Image.open( infilename )
img.load()
# note signed integer
return np.asarray( img, dtype="int32" )
def save_image( data, outfilename ) :
img = Image.fromarray( np.asarray( np.clip(data,0,255), dtype="uint8"), "L" )
img.save( outfilename )
def roberts_cross( infilename, outfilename ) :
image = load_image( infilename )
vertical = ndimage.convolve( image, roberts_cross_v )
horizontal = ndimage.convolve( image, roberts_cross_h )
output_image = np.sqrt( np.square(horizontal) + np.square(vertical))
save_image( output_image, outfilename )
infilename = sys.argv[1]
outfilename = sys.argv[2]
roberts_cross( infilename, outfilename )
Aus dem Wikipedia-Eintrag zu Roberts Kreuz. http://en.wikipedia.org/wiki/Roberts_Cross
Die Ausgabe meines Skripts.