Changed pixel neighborhood array to image region iterator.
Updated image border handling.
Used variable length vectors to represent joint pixel.