Попробуйте этот програмный пакет (ESO MIDAS) http://www.eso.org/esomidas. Он включает в себя пакет DAOPHOТ ( фотометрия и точные координаты звезд на кадре), IMRES (программы восстановления изображений), INVENTORY( выделение и классификация различных объектов на изображениях),PEPSYS (планирование фотометрических наблюдений и полная их обработка), SURFPHOT( поверхностная фотометрия галактик, построение изофот...). Еще один пакет IRAF : http://iraf.noao.edu . И если ничего из этого не устроит, то поищите здесь: http://asds.stsci.edu.

Сам всем этим не пользовался, только ссылки знаю.

Спасибо, Ирис действительно рисует 

Изолиния - это множество точек, на которых некоторый параметр остается постоянным. Обычно изолинии изображаются семейством, то есть для параметра, по которому строится изолиния, выбирают множество дискретных значений, и для каждого значения строится изолиния.
Для изолиний якрости этот параметр - яркость.
Арина, если тебе надо построить изолинию для заданного уравнения, тебе надо последовательно приравнивать правую часть к некоторой константе и отыскивать множесто точек, являющихся решением данного уравнения.
То есть если задана функция
U(x, y), то, решая последовательно уравнения
U(x, y) = Ci, где Ci - некоторая константа, то множество точек {x, y}, являющихся решением данного уравнения, и будут формировать изолинию для уровня Ci.
Множество уровней (констант) Ci выбирается произвольно, но исходя из разумных критериев