5
Давайте предположим, что я хочу реализовать следующий код питонаМожно ли векторизовать функцию, которая обращается к различным элементам в массиве numpy?
Эта функция принимает изображение как 1 одномерный массив, и перебирает отдельных элементов в массиве (пикселей в исходном изображении), который влияет на выходной массив, который является также изображение, представленное в виде 1 одномерного массива
пример: один пиксель на входном изображении (красный) влияет 8 окружающих пикселей в (оранжевый) 
базовая реализация в C является
/* C version
* Given an input image create an
* output image that is shaped by individual pixels
* of the input image
*/
int image[width * height]; //image retrieved elsewhere
int output [width * height]; //output image
int y = 0, x = 0;
for(y = 1; y < height-1 ; ++ y) {
for(x = 1; x < width-1; ++ x) {
if (image[y * width + x] > condition) {
/* pixel affects the surrounding 8 pixels in the output image */
output[(y-1) * width + x - 1]++; /* upper left */
output[(y-1) * width + x ]++; /* above */
output[(y-1) * width + x + 1]++; /* upper right */
output[y * width + x + 1 ]++; /* right */
output[y * width + x - 1 ]++; /* left */
output[(y+1) * width + x - 1]++; /* lower left */
output[(y+1) * width + x ]++; /* below */
output[(y+1) * width + x + 1]++; /* lower right */
}
}
}
наивный подход питона будет использовать точно такой же поэлементно доступ, как показано ниже
#Python version
input = blah # formed elsewhere
output = np.zeros(width * height)
for y in xrange(1, height-1):
for x in xrange(1, width-1):
if input[y * width + x] > condition:
output[(y-1) * width + x - 1]+= 1; # upper left
output[(y-1) * width + x ]+= 1; # above
output[(y-1) * width + x + 1]+= 1; # upper right
output[y * width + x + 1 ]+= 1; # right
output[y * width + x - 1 ]+= 1; # left
output[(y+1) * width + x - 1]+= 1; # lower left
output[(y+1) * width + x ]+= 1; # below
output[(y+1) * width + x + 1]+= 1; # lower right
есть лучший способ осуществить это? Можно ли векторизовать эту функцию?