Документ взят из кэша поисковой машины. Адрес оригинального документа : http://xmm.vilspa.esa.es/sas/7.0.0/doc/imgrad/node2.html
Дата изменения: Thu Jun 29 02:11:07 2006
Дата индексирования: Sat Dec 22 07:53:27 2007
Кодировка:
Поисковые слова: южная атлантическая аномалия

XMM-Newton Science Analysis System


imgrad (tools-1.63) [xmmsas_20060628_1801-7.0.0]

Meta Index / Home Page

Description

This task approximates the gradient at each pixel of an image by doing least-squares fitting of a plane to the 3x3 array of values centred on each pixel. Since gradient is a vector quantity, two output images are required. The user can choose between a cartesian representation (in which the output images record the x and y components of the gradient) or a polar representation (in which the output images are the magnitude and azimuth of the gradient.).

`Least-squares fitting of a plane' sounds very grand but in fact the algebra boils down to the following:

$$(\nabla I)_{x,y} = \frac{1}{6} \left( \begin{array}{c} \sum_... ...i,y+1} - \sum_{i=x-1}^{x+1} I_{i,y-1} \\ \end{array} \right).$$

Here $I_{x,y}$ represents the image value at the ($x,y$)th pixel.

No gradient value is calculated at the edge: nulls are stored in the output at these pixels. If the input image contains null-valued pixels, all 8 nearest neighbours of such pixels are set to null in the output.

The azimuth, where this is to be calculated, is $\arctan(\nabla_y/\nabla_x)$.