Registro de imagens é uma operação importante para aplicações de sensoriamento remoto que basicamente envolve a identificação de muitos pontos de controle nas imagens. Como a identificação manual de pontos de controle pode ser demorada e tediosa foram desenvolvidas várias técnicas automáticas. O presente tutorial tem o intuito de demonstrar o processo de registro automático imagem a imagem no software ERDAS 9.1 através da extensão (modulo) AutoSync, que é bastante intuitivo. O processo de correção geométrica que será desenvolvido aqui necessita de imagens, tanto de entrada quanto de saída, com menos de 5% de nuvens por quadrante.

Para saber mais sobre registro automático acesse:

Sistema semi-automático de registro e mosaico de imagens. - Dmitri Fedorov

Automatic Image Registration and Mosaicking System

Hoje falaremos um pouco sobre o processo de Classificação Supervisionada no ERDAS 9.1, que é outra técnica comum no processamento digital de imagens. Diferentes métodos de classificação podem ser empregados. A técnica consiste basicamente na criação de assinaturas multi-espectrais, as quais constituem critérios para a classificação do uso do solo de acordo com modelos de distribuição estatística multivariada.

Para estudar mais:

Utilização de técnicas de classificação digital de imagens landsat ...

Comparação entre métodos de classificação supervisionada

Classificação automática de imagens de satélites ...

Um ótimo estudo a todos.

Por: Luis Sadeck

Depois de termos apresentado um tutorial sobre recorte regular no ERDAS 9.1, muitas perguntas chegaram sobre recorte irregular ou com limites de municípios ou áreas especiais. Apesar de já termos feito alguns tutoriais de outros softwares com essa temática, alguns técnicos só dispõem de licenças para outros softwares, por isso, vamos agora ao recorte irregular no ERDAS 9.1.

Bom trabalho!!!!!

Até a próxima.

Por: Luis Sadeck

Olá a todos, hoje vamos demonstrar como se faz o processo de recorte regular no Software ERDAS 9.1. Esse processo é bastante usado quando nossa área de trabalho não abrange a(s) imagem(s) toda(s) e precisamos retirar uma parte da informação que, nesse caso é desnecessária.

Optamos hoje pelo recorte regular, pensando na criação do Layout da carta imagem depois, ou seja, com o recorte regular você terá uma apresentação visual bem melhor.

Então vamos ao que importa.

Um abraço a todos e até a próxima.

Por: Luis Sadeck

Olá queridos, fazia tempo que não postava nada sobre o ERDAS, então para retomarmos os estudos nesse software que já está indo para a versão 9.2, vamos demonstrar um processo bem básico mais de extrema importância.

Quem trabalha com georreferenciamento através das imagens Geocover tem como produto uma imagem no sistema de projeção UTM e na maioria das vezes os dados em Shapefile estão em geográfica, isso gera alguns problemas, mesmo se o software de GIS conseguir sobrepor essas duas informações. Para que isso não aconteça e suas informações não corram o risco de erros básicos, vamos demonstrar o processo de reprojeção.

Espero que esse tutorial incentive o estudo das projeções.

Um Abraço.

Por: Luis Sadeck

Olá mais uma vez, falei agora a pouco sobre Layer Stack e não podia deixar de fora o Registro de Imagem, por isso, resolvi ficar mais um pouco e disponibilizar esse tutorial.

Sendo assim, para começar registro é uma transformação geométrica que relaciona coordenadas da imagem (linha e coluna) com coordenadas geográficas (latitude e longitude) de um mapa. Essa transformação elimina distorções existentes na imagem, causadas no processo de formação da imagem, pelo sistema sensor e por imprecisão dos dados de posicionamento da plataforma (aeronave ou satélite).

A necessidade de fazer o registro:
· Integração de imagens obtidas por sensores diferentes.
· Imagens obtidas em tempos diferentes. Análise temporal.
· Imagens tomadas em posições diferentes. Obter informação tridimensional
· Mosaico de imagens

Nesse tutorial usamos o Geocover 2000 para servir de referência para o nosso registro.

Espero poder ajudar

Um abraço

The values for pixels in one waveband typically ranges from 0 to 255, this being the range of values that can be represented by one byte of computer storage. Thus, the detectors on the sensing devices will be set such that the energy will be recorded digitally within that range. For any one image, the majority of the actual data values recorded may be concentrated within only one small part of the 0-255 range.

A consequence of this is that small variations in pixel values, which may have interest and which might distinguish between features, may be difficult to detect visually (Jones, 1998). Viewing remotely-sensed images is therefore concerned with manipulating the pixel values in such a way that features can be clearly discriminated between when the image is displayed. It is therefore clear that the amount of information that can be retrieved from an image lies not in the cell values per se, but in the way these values can be displayed.

References

Mather, P. (1999) Computer Processing of Remotely-Sensed Images, An Itroduction, 2nd ed.

Lillesand, T.M., and Keifer, R. W. (1994) Remote Sensing and Image Interpretation, 3rd ed

Jones, C. (1997) Geographical Information Systems and Computer Cartography

Queridos amigos, hoje estou passando rápido só para deixar mais um tutorial, dessa vez sobre o ERDAS, bem simples e de fácil compreensão. O tutorial é, mais especificamente sobre composição colorida ou como dizemos, Layer Stack.

Espero que tenha ficado bom! Aguardo a Contribuição de vocês.

Um abraço.

To distinguish different features in a TM image, the figure below left is of help, as it shows how different features generally display their reflectance values in the different bands of the TM sensor.

Comments on findings in “peak-tm.84”: Firstly, the peak in band 5 is not consistent with the generalised reflectance, which should be lower in band 5 than band 4. Secondly, the peat bog should, if it were wet, show a much lesser reflectance in bands 5 and 7. Even so, the lower overall reflectance value of the peat bog, compared to the other features.

Image Info

Image info yields valuable information about the image and the data contained in each band, such as:

Information about the raster itself, rows and columns, data type etc.
Statistical information, mean, max mode etc.
Image coordinates and pixel size

Image info as requested in practical 1:

Lond-pan94: #rows=6012, #columns=7000, meanDN=44.966, data type= Unsigned 8bit

Leic-tm94-geo:
#rows=836, #columns=872, meanDN in band 4 = 191.660, pixel size=35m

Leic-tm84-geo:
maxDN in band 4= 201, ULx=442045.0, ULy=32110.0, LRx=471025.0, Lry=290870.0

Flevoland:
Data type= Unsigned 16bit, ULx=0, Uly=0, LRx=487.0, LRy=-401.0

AVHRR-4: Data type= Unsigned 8bit

Spatial profile

Profile of band 4 in leic-tm92, across Cemetery, University campus and Victoria Park showing low level of reflectance (less vegetation) in built-up areas and campus, and high reflectance across Victoria Park.

Surface profile of band 4 in area near motorway junction. Assumptions that can be made is that the fields along the motorway are barren or dry (or with crops that have little reflectance), whereas the motorway is clearly flanked by trees and vegetation.

Blend, Swipe and Flicker

When two images have been georeferenced, they can be laid on top of each other. Blend allows blending the displays into each other, Swipe rolls like a curtain over the top layer, hiding or revealing the bottom layer. Flicker allows switching between the images. These functions can be operated manually using buttons/ scrollbars or automatically at a set rate.

References

Mather, P. (1999) Computer Processing of Remotely-Sensed Images, An Itroduction, 2nd ed.

Lillesand, T.M., and Keifer, R. W. (1994) Remote Sensing and Image Interpretation, 3rd ed

Jones, C. (1997) Geographical Information Systems and Computer Cartography

One of the most basic enhancement techniques is the contrast stretch. If the screen is set to display 256 grey scale colours, band3 in leic-tm92 will appear dark, as most values are clustered around the mean value of 29, with the overall values ranging from approximately 21 to 55, taking up only 34 different values or 14% of the full value range (0-255) that can be utilised.

Using stretch will distribute the pixel values according to the algorithm used so that the whole value range is used. In doing so the algorithm leaves the pixel value unchanged, but is assigned a display value in a lookup table (LUT). As an example, the pixel value at (276, -230) is 23 in the image; the LUT value is 23 without stretch and 47 with stretch applied.

Decorrelation stretch

The effect of decorrelation stretch is clearly visible in the following images. The effect of this stretch is greater colour differentiation, clearly seen in the histograms, where the pixel values are grouped and distinctively separated from each other.

Histogram - Without decorrelation stretch

Histogram - With decorrelation stretch

Image - Without decorrelation stretch

Image - With decorrelation stretch

Edge enhancement

The “Crisp” function passes an edge-detecting or “high-pass” filter across the image. A high-pass is a sharp change in the grey scale value at a particular pixel point, and may be interpreted as features of the built environment, roads or field boundaries. In the combined image below the right part shows more distinguishable features than the left part, which has not been enhanced. As can be seen in the upper right corner of each picture, differences in neighbouring pixel values make an impact on the crisp picture, whereas the area is more or less consistently blurry green in the original image.

(image to be posted)

References

Mather, P. (1999) Computer Processing of Remotely-Sensed Images, An Itroduction, 2nd ed.

Lillesand, T.M., and Keifer, R. W. (1994) Remote Sensing and Image Interpretation, 3rd ed

Jones, C. (1997) Geographical Information Systems and Computer Cartography

From the general spectral reflectance the following observation can be made:

Vegetation - NIR/R >>> 1, Water - NIR/R < 1, Soil – NIR/R > 1.

NIR/R–images can serve as a crude classifier of images, and indicate vegetated areas in particular. Therefore this ratio has been developed into a range of different vegetation indices.

Original image

NIR/R (TM4/TM3), with values ranging from 0 to 8, soil to vigorous vegetation, vegetated areas shown in white.

NIR/R, stretched, with values ranging from 16 to 254, giving a smoother transition between surface types.

References

Mather, P. (1999) Computer Processing of Remotely-Sensed Images, An Itroduction, 2nd ed.

Lillesand, T.M., and Keifer, R. W. (1994) Remote Sensing and Image Interpretation, 3rd ed

Jones, C. (1997) Geographical Information Systems and Computer Cartography

Part of image with missing scan line

Regarding the position of the missing scan line, to find the correct row number, it must considered that the image peak-tm84 has 512 rows and 512 columns according to it’s image info, with coordinates upper left 1/1(y/x) and lower right 512/-510 (y/x).

Since the missing scan line is at x = –200, this is row number 201.

Sensor de-striping

Destriped image (left) and original image (right).

The destriped image shows a small improvement compared to the original image.

Mather (1999) describes two ways of destriping, using a linear method or histogram matching. In the first, the mean and standard deviation of the histograms of every 16 bands are similarly forced to be equal to the mean and the standard deviation of all the pixels in the image. In the latter a target histogram is constructed, which is the cumulative histogram of the entire image. Using the cumulative histogram of one of the detectors, meaning the cumulative histogram of every 16th line, a comparison is made, and an output pixel value is calculated. ERDAS uses a convolution method, which passes an averaging matrix across the image, which may account for stripes still being apparent in the image.

Atmospheric path luminance correction

Original image with default 50/50% --- And 100/70% brightness/contrast

The pixel values in the bands used for display (1-3) take up only about the first 40 of the possible 256 values, which is why the image seems dark. The brightness tool moves the colours into white, while contrast stretches the histogram for better colour differentiation. Thus, the dark bluish pixel values are replaced by bright tones in the right image, speckled with red tones.

Corrected image with default 50/50% ---- and 100/70% brightness/contrast

After correcting for luminance, all pixel values start more or less from 0, which is why the default image is almost totally black compared to the original image.

Due to the effect of path radiance the histogram initially does not have its lowest value at zero. This is why the 100/70% contrast brightness is bluish bright, since the path radiance is highest in band 1, which is chosen as the blue colour component. After subtracting the minimum DN, the corrected image appears with a normalised true colour approximation. Bands 5 and 7 normally need not to be corrected (Mather, 1999), which is also apparent in the table below.

The table shows the minimum and maximum values before (B) and after (A) the dark pixel subtraction.

Geometric correction

-- Using Ground Control Point (GCP)

The GCPs should be spread evenly over the image, covering the whole image, and be placed as much as possible into the corners of the image, to give best coverage for calculating the transformation.

This is not he case in the provided data, as no corners are used, and there is a slight over-representation of GCP in the lower half of the image. Unfortunately, adding GCP from the given map is also limited, as the map covers much less area than the actual image. Clustering GCPs in the centre of the image will not give accuracy to the transformation of coordinates in the extremities of the image.

Adding GCP from a map is easy and straightforward. However, care should be taken when choosing the correct position of a GCP, as a scanned map might add to the RMS error, because roads tend to be drawn wider than they actually are, or even slightly displaced, i.e. University Road is said to be approximately 30m wide. In addition, the pixel size of the map is 3m; the image itself operates with 30m, which is another source of error in correct placement of GCP.

-- Image-to-image transform
Image to image transform is especially useful where geo-referencing is not necessary for the analysis of multiple images of the same area. Again, the same rule applies as above, namely spreading the GCPs as much as possible. For the reasons already mentioned, areas that are not covered by both images and that do not contain GCPs should not be considered as “matching” and should therefore not be used in comparison studies.

References

Mather, P. (1999) Computer Processing of Remotely-Sensed Images, An Itroduction, 2nd ed.

Lillesand, T.M., and Keifer, R. W. (1994) Remote Sensing and Image Interpretation, 3rd ed

Jones, C. (1997) Geographical Information Systems and Computer Cartography

In the true colour image, with bands 1, 2 and 3, there seem to be little difference between principal component analysis (PCA) and tasselled cap transform (TCT). Even though the colour differs, the same pattern can be discerned.

The resulting images from PCA and TCT are similar, but not directly comparable.

TCT defines three functions or axes a priori: “wetness”, “greenness” and “brightness”.

TCT is particularly useful in defining soils (wetness) and vegetation (greenness) from Landsat MSS or TM data.

In PCA, the axes are dependent on statistical relationships among the spectral bands of the particular image being analysed.

Thus, inter-image comparison is more difficult when using PCA compared to TCT, because the rotation of the principal axes is dependent on the image that is used.

PCA will identify the number of dimensions in a data set: If all bands were highly correlated, the data would practically speaking have only one dimension.

In PCA, the axes with the most variability (or the least correlation) are found in PC1 (or principal component 1). Thus, this image contains most of the information that can be inferred from the image.

Principal component 1

Axes with most variability in the data are located in the first principal component (PC1). Thus, this component contains most of the information that can be inferred from the data and shows many distinguishable features.

PC1 and PC2

PC2 shows other hidden information, i.e. a clear separation of water bodies from surroundings, something that is not so obvious from PC1.

TCT band 5 and band 3

Band 5 in TCT emphasises wetness, so water bodies will stand out from the background. From the TCT coefficients in Mather (1999), green areas in band 3 will be dark, whereas wet areas will appear bright.

References

Mather, P. (1999) Computer Processing of Remotely-Sensed Images, An Itroduction, 2nd ed.

Lillesand, T.M., and Keifer, R. W. (1994) Remote Sensing and Image Interpretation, 3rd ed

Jones, C. (1997) Geographical Information Systems and Computer Cartography

When performing an unsupervised classification it is necessary to find the right number of classes that are to be found. Too many, and the image will not differ noticeable from the original, too few and the selection will be too coarse.

Original image

Unsupervised classification, 10 classes

Unsupervised classification, 6 classes

The difference between 6 and 10 unsupervised classes is the merger of urban and residential as well as agricultural fields. The table below summarises the convergence for every iteration, depending on the number of classes. The convergence describes how many of the pixels stay in the same cluster between one iteration and the next.

The classes were determined by referring to the Ordnance Survey Landranger map of Leicester. Still, not all classes are consistent, i.e. water tends to appear many places. Forest, which one would expect to stand out in summer, when this picture is taken, does not have its own class.

Dry grassland was inferred from the area around the racecourse. Agricultural areas were named so because of their typical pattern. Parkland areas were consistent with park areas on the map.

Looking at the separability table, some class are clearly separable, others not so clearly. Class 1 (city centre) and class 5 (residential) show the lowest average separability. The same goes for class 2 (roads) and 4 (railways). What is a bit puzzling is the fact that small water bodies like rivers have no distinct class, as they seem to mingle with other classes, especially roads and railways. This is due to their low reflectance, making it difficult to distinguish between them.

Supervised classification

There are different approaches to supervised classification.

Supervised classification requires a priori knowledge of the number of classes, as well as knowledge concerning statistical aspects of the classes. All methods start with establishing training samples, which are areas that are assumed or verified to be of a particular type. The classification algorithms will sent “sort” the pixels in the image accordingly.

Minimum distance

Maximum Likelihood

Maximum Likelihood/ Parallelepiped

The Minimum Distance algorithm allocates each cell by its minimum Euclidian distance to the respective centroid for that group of pixels, which is similar to Thiessen polygons. The Maximum Likelihood classifier applies the rule that the geometrical shape of a set of pixels belonging to a class often can be described by an ellipsoid. Consequently, pixels are grouped according to their position in the influence zone of the class ellipsoid. The Parallelepiped method is a rectangle, where the lowest and highest pixel values for the class in each band make up the boundary.

The Maximum Likelihood classification above has determined more water bodies and identified dense urban areas better than the Minimum Distance. Still, both methods display apparent errors. Choosing better training areas can improve this, but not necessarily so. The incorrect classification of water bodies in the Parallelepiped is due to the method’s inherent inaccuracy as the minimum and maximum value in a class seldom are representative values for a particular class.

To compare the unsupervised and supervised classification above is difficult, because their signature files do not show the same classes. In an image with high separability unsupervised classification may be used , whereas low separability will need the aid of supervision.

References

Mather, P. (1999) Computer Processing of Remotely-Sensed Images, An Itroduction, 2nd ed.

Lillesand, T.M., and Keifer, R. W. (1994) Remote Sensing and Image Interpretation, 3rd ed

Jones, C. (1997) Geographical Information Systems and Computer Cartography