Как заполнить "дыры" в изображении?

У меня есть фотоизображения галактик. На этих изображениях есть некоторые нежелательные данные (например, звезды или полосы от самолетов), которые маскируются. Я хочу не просто заполнить замаскированные области неким средним значением, а интерполировать их в соответствии с окружающими данными. Как это сделать в python?

Мы пробовали различные функции из пакета SciPy.interpolate: RectBivariateSpline, interp2d, splrep/splev, map_coordinates, но все они, похоже, работают для поиска новых пикселей между существующими, мы не смогли заставить их заполнить произвольную "дыру" в данных.