Revista Brasileira de Geofı́sica (2006) 24(4): 559-572

© 2006 Sociedade Brasileira de Geofı́sica

ISSN 0102-261X
[Link]/rbg

INTEGRAÇÃO E ANÁLISE DE DADOS AEROGEOFÍSICOS POR MEIO DA APLICAÇÃO DE TÉCNICAS

DE PROCESSAMENTO DIGITAL DE IMAGENS E CLASSIFICAÇ ÃO NÃO SUPERVISIONADA:
O EXEMPLO DO GREENSTONE BELT RIO DAS VELHAS, QUADRILÁTERO FERRÍFERO, MG

Alexandre de Amorim Teixeira1 , Adalene Moreira Silva2 , Augusto César Bittencourt Pires3 ,
Roberto Alexandre Vitória de Moraes4 e Carlos Roberto de Souza Filho5
Recebido em 25 agosto, 2006 / Aceito em 17 novembro, 2006
Received on August 25, 2006 / Accepted on November 17, 2006

ABSTRACT. Poor outcrop and deep weathering characterize the Rio das Velhas Greenstone Belt in the southeastern from Brazil. This paper summarizes the use of
high-density airborne survey for mineral exploration studies based on interpretation enhancements of magnetic, radiometric and frequency domain electromagnetic data
using image-processed methods and an unsupervised classification. The generated products provide new insights and an excellent tool for mapping and trace individual
lithological units, improving the information content of the single geophysical channels. The geophysical images were processed using different combinations. The best
product was the analytical signal amplitude and inclination integrated by IHS transformation. Information extracted from this image maps the geology and lineament
patterns at both regional and local scales. The K-means technique using ten classes was also applied to the geophysical data. These results enhance the lithologies
mapped by the field geologists at the 1:100.000. Also shows important host rocks and different gold mineralized geological domains. Such domains host the known
gold mineralization, illustrating the utility of these techniques to improve the geological knowledge in the study area.

Keywords: Digital Imaging Processing, gold mineralization, Rio das Velhas Greenstone Belt.

RESUMO. Afloramentos escassos e um intemperismo acentuado caracterizam o Greenstone Belt Rio das Velhas no sudeste do Brasil. Este artigo sumariza a utilização
de dados aerogeofı́sicos de alta densidade de amostragem aplicados à exploração mineral baseado no realce e na interpretação de dados magnéticos, radiométricos e
eletromagnéticos no domı́nio da freqüência através de métodos de processamento digital de imagens e classificação não supervionada. Os produtos gerados forneceram
novos insights e uma excelente ferramenta para mapeamentos dos diferentes litotipos, melhorando o conteúdo da informação presente nos canais individuais. As
imagens geofı́sicas foram processadas utilizando diferentes combinações. O melhor produto obtido foi a integração via IHS da amplitude e da inclinação do sinal
analı́tico. Informações extraı́das desta imagem mapeiam tanto a geologia quanto lineamentos na escalas regionais e locais. Aplicou-se a técnica de classificação não
supervisionada conhecida com média K aos dados aerogeofı́sicos. O resultado realça litologias mapeadas por geólogos de campo na escala 1:100.000. Ela mapeia
também rochas hospedeiras e diferentes domı́nios associados com a mineralização aurı́fera. Tais domı́nios são conhecidos por hospedar a mineralização aurı́fera,
ilustrando a utilidade destas técnicas para enriquecer o conhecimento geológico da área estudada.

Palavras-chave: Processamento Digital de Imagens, mineralização aurı́fera, Greenstone Belt Rio das Velhas.

1 Agência Nacional de Águas – ANA, Setor Policial Sul, Área 05 Quadra 03, Bloco “B” – Sala 105, 70610-200 Brası́lia, DF, Brasil. Tel: +55 (61) 2109-5378

– E-mail: [Link]@[Link]
2 Laboratório de Geofı́sica Aplicada, Instituto de Geociências, Universidade de Brası́lia, 70910-900 Brası́lia, DF, Brasil. Tel: +55 (61) 3307-2877; Tel./Fax: +55 (61)

3340-4759 – E-mail: adalene@[Link]

3 Laboratório de Geofı́sica Aplicada, Instituto de Geociências, Universidade de Brası́lia, 70910-900 Brası́lia, DF, Brasil. Tel: +55 (61) 3307-2877; Tel./Fax: +55 (61)

3340-4759 – E-mail: acbpires@[Link]

4 HGeo/InterGeo, SCLN 111, Bloco C, Salas 219/221, 70754-530 Brası́lia, DF, Brasil. Tel./Fax: +55 (61) 3274 4533 – E-mail: rmoraes@[Link]
5 Instituto de Geociências (IG), Universidade Estadual de Campinas (UNICAMP), Caixa Postal 6152, 13083-970 Campinas, SP, Brasil. Tel: +55 (19) 3788-4535;

Fax: +55 (19) 3289-1097 – E-mail: beto@[Link]

560 INTEGRAÇÃO E ANÁLISE DE DADOS AEROGEOFÍSICOS

INTRODUÇÃO Os terrenos granito-gnáissicos do QF e de suas imediações

No que tange à integração de dados multifontes para mapeamento apresentam, às vezes, uma forma dômica e compreendem di-
geológico, não existe uma técnica que contemple a complexidade versos complexos: Caeté, Bonfim, Belo Horizonte e Santa Rita.
dos diversos tipos de rochas e suas alterações intempéricas ou São polideformados e os contatos com as unidades supracrus-
hidrotermais. Assim, esforços têm sido despendidos no intuito tais adjacentes são de natureza tectônica. Esses complexos me-
de gerar novas metodologias e procedimentos que auxiliem na tamórficos são constituı́dos, majoritariamente, por ortognaisses
definição de domı́nios definidos nestes dados e que se correlaci- do tipo TTG, migmatizados ou não. Apresentam intrusões de
onem espacialmente e espectralmente com a unidade litológica, granitóides, anfibolitos e pegmatitos. As rochas dos complexos
esteja ela alterada ou não. foram geradas preferencialmente nos ciclos Arqueano e Transa-
A região do Quadrilátero Ferrı́fero permite esse tipo de pes- mazônico (Cordani et al., 1980 e 1989; Teixeira, 1985; Machado
quisa, devido ao grande acervo de dados de geologia básica, et al., 1989; Romano, 1989; Machado & Carneiro 1992; Carneiro
aerogeofı́sica e de prospecção geoquı́mica, gerados por órgãos et al., 1995; Noce, 1995; Noce et al., 1992; Teixeira et al., 1996).
governamentais e empresas de mineração. O greenstone belt As rochas do Greenstone Belt Arqueano (Schorscher 1978)
Rio das Velhas hospeda alguns dos mais importantes depósitos são agrupadas no Supergrupo Rio das Velhas (SGRV), dividido
aurı́feros do paı́s, bem como apresenta uma grande diversidade nos Grupos Nova Lima e Maquiné. Ladeira (1980) propôs dividir
litológica caracterizada por rochas vulcânicas e metassedimenta- o Grupo Nova Lima em três unidades, da base para o topo: uni-
res, com expressivas zonas de alteração hidrotermal. dade metavulcânica, unidade metassedimentar quı́mica, unidade
Produtos integrados por meio de Processamento Digital de metassedimentar clástica. Segundo Schrank & Silva (1993), parte
Imagens e Classificadores Hierárquicos permitem a geração de importante do pacote sedimentar é representado por seqüências
temas que possibilitam uma nova compreensão do caráter estru- turbidı́ticas de composição básica a félsica. Sobreposto a esse
tural e litológico da região. No entanto, ainda é necessário re- conjunto encontra-se o Grupo Maquiné que foi dividido em
finar, desenvolver, aplicar e comparar novas técnicas neste sen- duas formações, da base para o topo: Palmital e Casa Forte
tido. O foco deste trabalho é auxiliar a interpretação geológica e (Dorr, 1969).
a compreensão das mineralizações aurı́feras no greenstone belt Não existe uma subdivisão válida para o SGRV como um todo.
Rio das Velhas. No que tange a uma continuação no processa- Pinto & Silva (1996), Zuchetti et al. (1996) propuseram uma nova
mento de dados aerogeofı́sicos já desenvolvido (Silva, 1999), a subdivisão estratigráfica para o Greenstone Belt Rio das Velhas.
deconvolução de Euler representa uma tentativa de se compreen- Pinto & Silva (1996) propõem elevar as unidades a categoria de
der o comportamento espacial das fontes magnéticas relevantes e formações. Segundo Zuchetti et al. (1996) estas unidades fo-
sua relação com arcabouço estrutural e mineralizações associa- ram definidas por observações de campo, com o apoio de estu-
das. Dentro desta ótica será apresentada uma discussão sobre os dos petrológicos e geoquı́micos. Um empilhamento estratigráfico
métodos utilizados, visando apresentar as vantagens e desvan- clássico das unidades não foi efetuado devido à alta magnitude de
tagens de cada produto integrado no processo de interpretação deformação, o elevado grau de intemperismo e extensão da ba-
estrutural e litológica para o ambiente estudado. cia de deposição. Ainda, de acordo com Pinto & Silva (1996),
as unidades foram agrupadas em blocos tectônicos, denomina-
dos de: Nova Lima, Caeté, Santa Bárbara e São Bartolomeu, que
CONTEXTO GEOLÓGICO segundo os autores, retratam ambientes petrogenéticos distintos
O Quadrilátero Ferrı́fero (QF) situa-se na borda sul do Cráton São (Figura 2).
Francisco e representa uma parte preservada de um antigo núcleo Dentro dessas modificações foi mantida a designação formal
maior denominado de Cráton do Paramirim (Almeida, 1977). De- de SGRV, com os Grupos Nova Lima (base) e Maquiné (topo). No
vido ao seu potencial mineral o QF é alvo de estudos geológicos entanto, o Complexo Córrego dos Boiadeiros e o Grupo Quebra
sistemáticos desde o final do século XIX. Osso foram redefinidos como formações e englobados dentro do
Os litotipos aflorantes abrangem uma área de aproximada- Grupo Nova Lima. O Complexo Córrego dos Boiadeiros foi sepa-
mente 7190 km2 e são constituı́dos de terrenos granitos-gnáissi- rado como um sill ultrabásico intrusivo na base do Grupo Nova
cos; Greenstone Belt Arqueano; seqüências supracrustais do Pa- Lima (Zucchetti et al., 1996; Pinto & Silva, 1996).
leoproterozóico e do Mesoproterozóico (Chemale Jr et al., 1991 O Supergrupo Minas (SGM) sobrepõe-se aos dos TTG’s e
e 1994) (Figura 1). SGRV, geralmente em contato tectônico. Este divide-se em três

Revista Brasileira de Geofı́sica, Vol. 24(4), 2006

 A. AMORIM TEIXEIRA, A. MOREIRA SILVA, A.C. BITTENCOURT PIRES, R.A. VITORIA DE MORAES e C.R. DE SOUZA FILHO 561

Figura 1 – Mapa Geológico do Quadrilátero Ferrı́fero, MG (Modificado de Dorr, 1969 in Silva, 1999).

grupos, da base para o topo: Caraça, Itabira e Piracicaba (Dorr, Gandarela, no topo do grupo, exibe mármores dolomı́ticos, su-
1969; Cordani et al., 1980 e 1989; Ladeira & Viveiros, 1984; Alk- bordinadamente, itabiritos e filitos.
mim & Marshak, 1998; Carneiro et al., 1995; Noce, 1995). A unidade basal do Grupo Piracicaba, Formação Cercadinho,
O Grupo Caraça apresenta, na base, a Formação Moeda, caracteriza-se pela alternância de quartzitos e filitos, freqüente-
constituı́da por quartzitos com intercalações de filitos e nı́veis mente ferruginosos. A Formação Fecho do Funil é constituı́da
conglomeráticos aurı́feros-uranı́feros na porção basal. Esta tran- por filitos quartzosos, filitos dolomı́ticos e lentes de dolomito.
siciona para a Formação Batatal, onde predominam filitos se- As Formações Taboões (ortoquartzitos) e Barreiras (filitos grafi-
ricı́ticos, por vezes carbonosos ou ferruginosos. tosos) são de ocorrência restrita.
O Grupo Itabira inicia pela Formação Cauê, composta por A unidade superior do Grupo Piracicaba, Formação Sabará,
formações ferrı́feras bandadas do tipo Lago Superior e, subor- é constituı́da por clorita xistos, filitos, metagrauvacas, metatufos,
dinadamente, por filitos ferruginosos e dolomitos. A Formação metaconglomerados, quartzitos e formação ferrı́fera subordinada.

Brazilian Journal of Geophysics, Vol. 24(4), 2006

562 INTEGRAÇÃO E ANÁLISE DE DADOS AEROGEOFÍSICOS

Figura 2 – Nova subdivisão estratigráfica para o Greenstone Belt Rio das Velhas. As unidades foram agrupadas em
blocos tectônicos, denominados de: Nova Lima, Caeté, Santa Bárbara e São Bartolomeu, que segundo Pinto & Silva
(1996), Zuchetti et al. (1996) retratam ambientes petrogenéticos distintos.

Segundo Dorr (1969), a Formação Sabará possui uma espes- O Supergrupo Espinhaço é representado por um pacote ro-
sura máxima de 3000 m, aproximando-se da espessura máxima chas quartzı́ticas a nordeste do Quadrilátero Ferrı́fero, na Serra
atingida pelo conjunto das demais unidades do SGM, que é de das Cambotas. A posição estratigráfica desta unidade foi sem-
3500 m. Barbosa (1968) e Ladeira (1980) propõem a elevação da pre controvertida por incluir pacotes distintos de rochas, tecto-
Formação Sabará à categoria de grupo, devido a sua espessura e nicamente justapostos (Freitas et al., 1991; Crocco-Rodrigues et
presença de litotipos diferenciados. Segundo Reis et al. (2002), al., 1991).
o Grupo Sabará representa uma bacia de antepaı́s compartimen- Dentro do quadro geológico regional do Quadrilátero Fer-
tada, gerada durante a deformação, soerguimento e erosão dos rı́fero, deve-se destacar também a existência de uma densidade
depósitos Arqueanos e Paleoproterozóicos do Cinturão Mineiro. razoável de corpos ı́gneos máficos intrudindo as unidades Arque-

Revista Brasileira de Geofı́sica, Vol. 24(4), 2006

 A. AMORIM TEIXEIRA, A. MOREIRA SILVA, A.C. BITTENCOURT PIRES, R.A. VITORIA DE MORAES e C.R. DE SOUZA FILHO 563

anas a Fanerozóicas. Duas gerações de diques máficos de direção grafia, para o gamaespectrômetro de 45 metros para o mag-
predominantemente NS cortam os terrenos granito-gnáissicos netômetro/VLF e de 30 metros para o eletromagnetômetro. A
(TTG) do QF e regiões adjacentes (Carneiro, 1992). São encontra- taxa média de amostragem do eletromagnetômetro é de ± 5 me-
dos também sills e diques de rochas básicas de direção NS, ora tros (5 leituras/segundo), a do magnetômetro é de ± 5 metros
deformados, ora metamorfizados, com idades em torno de 1,7 G.a (5 leituras/segundo), a do VLF – EM é de ± 5 metros (5 leitu-
(Carneiro et al., 1995; Silva et al., 1995). A intrusão de diques ras/segundo) e, finalmente, o gamaespectrômetro possui taxa de
máficos na região do QF foi mais expressiva durante a abertura do amostragem de ± 25 metros (1 leitura/segundo). A velocidade
protooceano Brasiliano/Panafricano (Silva et al., 1995). O último de operação do helicóptero utilizado como plataforma esteve en-
evento, de idade Mesocenozóica aflora na forma de diques de tre 70 a 110 km/h. O eletromagnetômetro empregou três bobinas
direções variáveis, não apresentando deformação nem metamor- coplanares com transmissores centrados nas freqüências de 500,
fismo. A origem destes diques está relacionada à fragmentação 4.175 e 33.000 Hz e duas bobinas coaxiais com 935 e 4600 Hz.
do Supercontinente Gondwana (Silva et al., 1995). O afastamento (transmissor/receptor) foi de 7 metros.
As rochas dos terrenos granito-gnáissicos foram metamorfi- Os dados aerogeofı́sicos utilizados nesse trabalho foram pro-
zadas em condições de fácies anfibolito a granulito (Hertz, 1978) cessados por Silva (1999) e reprocessados por Fuck (2001).
e retrometamorfizadas em xisto verde, enquanto que as supra- O processamento destes dados envolveu a edição e junção dos
crustais do SGRV e do SGM evidenciam paragênese metamór- quatro blocos representativos da amostragem feita numa base de
fica de fácies xisto verde a anfibolito (Hertz, 1978; Hoefs et al., dados única.
1982; Ladeira, 1980; Marshak & Alkmim, 1989). Localmente, as Após esta etapa, os dados foram interpolados. Para os
condições metamórficas do SGM atingem apenas o anquimeta- dados magnetométricos e eletromagnetométricos (resistividade
morfismo. Hertz (1978) e Jordt-Evangelista et al. (1991) descre- aparente), o método de interpolação mais eficiente foi o bi-
vem auréolas de metamorfismo de contato em torno Complexo direcional (implementado no Oasis Montaj como bi-grid ), com
Bação e na porção norte da Serra do Curral, respectivamente. célula quadrada de 50 metros de lado. No caso dos dados ga-
maespectrométricos utilizou-se o método da curvatura mı́nima
(Fuck, 2001). Os métodos de interpolação foram testados exaus-
DADOS AEROGEOFÍSICOS tivamente e, posteriormente, analisados não só através da análise
O projeto Rio das Velhas consiste de aerolevantamento geofı́sico visual, mas também através do espectro de potência. A seguir,
sistemático utilizando os métodos magnetométrico, gamaespec- as diferentes malhas interpoladas foram micronivelados com o
trométrico e eletromagnetométrico (domı́nio da freqüência) (in- emprego de algoritmo otimizado (Fuck, 2001), o qual tratou es-
cluindo VLF) de alta densidade, com o objetivo de identificar pecificamente os problemas de desnivelamentos encontrados nas
novos alvos aurı́feros em zonas de Greenstone belt (Hildenbrand diferentes áreas abrangidas pelo aerolevantamento.
& Perez da Gama, 1993).
A escolha desses métodos teve por base a boa resposta
geofı́sica esperada para as rochas do Greenstone belt Rio das Magnetometria
Velhas que apresentam contraste de susceptibilidade magnética A análise do relevo referente ao campo magnético anômalo foi
com as rochas máfico-ultramáficas/formações ferrı́feras e zonas auxiliada por suas transformações lineares. A amplitude do sinal
de cisalhamento. Já a gamaespectrometria, contribui com a diver- analı́tico foi usada na demarcação da posição horizontal das fon-
sidade litológica causada pela variação de concentração dos ra- tes magnéticas. A fase do sinal analı́tico possui um padrão textu-
dioelementos K, U e Th, e sua relação com eventos hidrotermais. ral que auxilia a caracterização dos domı́nios magnéticos e feições
O método eletromagnético indica a presença de corpos conduto- lineares. A utilização das derivadas, principalmente das verticais,
res como sulfetos maciços e filitos carbonosos. ajudou na compreensão das posições espaciais relativas destas
O aerolevantamento foi realizado em 1992 e as linhas de vôo fontes e as horizontais, ajudaram na demarcação em mapa.
foram orientadas perpendicularmente à direção preferencial das A Figura 3 apresenta-se um dos mapas com as soluções da
rochas do Supergrupo Rio das Velhas no interior da área (Fi- equação de Euler gerados neste trabalho que apresentou a me-
gura 1) com espaçamento das linhas de vôo de 250, amarra- nor dispersão (ı́ndice estrutural N = 1, Linha de pólos – tipo di-
dos por linhas de controle de 5000 metros. Os sensores fo- que semi-infinito), feito na tentativa de se compreender o com-
ram mantidos a uma altura de 60 metros, por sobre a topo- portamento dos falhamentos. Seguindo a proposta de Reid et

Brazilian Journal of Geophysics, Vol. 24(4), 2006

564 INTEGRAÇÃO E ANÁLISE DE DADOS AEROGEOFÍSICOS

Figura 3 – Mapa de estimativa de profundidade de fontes magnéticas por deconvolução de Euler. Legenda: Muito Raso, Raso
(ciano claro), Pouco Raso (ciano escuro), Pouco Profundo (azul), Profundo (vermelho), Médio Profundo (amarelo claro), Muito
Profundo (amarelo escuro).

al. (1990), utilizaram-se janela de 10 × 10 células da malha do te classes (µ ± nσ ; n = −3, −2, −1, 0, 1, 2, 3) usando-
campo magnético anômalo com dimensões espaciais de 50, 125, se distribuição log-normal (Tabela 1), visto que visualmente o
250 e 500 metros. Essa ação visa a determinação de corpos em gráfico de freqüência das estimativas de profundidade apresenta
maiores profundidades à medida que se aumenta o tamanho da ja- esse tipo de distribuição.
nela. Para tanto, também se calculou a profundidade dos corpos
com tolerâncias de 5, 10 e 15% para todas as malhas. Quanto Tabela 1 – Reclassificação das estimativas de profundidade.
maior esta e menor é o tamanho da célula da malha, maior é o Intervalo da
número de ocorrências calculadas. Desse modo, a escolha da to- profundidade estimada Classificação
lerância visou limitar a grande quantidade de ocorrências geradas (metros)
pelas células menores e ampliar o número das ocorrências calcu- 0–62 Muito Raso
ladas pelas células maiores. Foram escolhidas as tolerâncias de 62–119 Raso
5% para malhas com células de 50 e 125 metros; de 10% para 119–230 Pouco Raso
aquelas com células de 250 metros e, por fim, 15% para as de 230–436 Pouco Profundo
500 metros. Todas as estimativas feitas foram corrigidas para a 436–852 Profundo
superfı́cie do terreno. 852–1643 Médio Profundo
As estimativas de profundidades foram classificadas em se- 1643–3706 Muito Profundo

Revista Brasileira de Geofı́sica, Vol. 24(4), 2006

 A. AMORIM TEIXEIRA, A. MOREIRA SILVA, A.C. BITTENCOURT PIRES, R.A. VITORIA DE MORAES e C.R. DE SOUZA FILHO 565

Gamaespectrometria a freqüência utilizada e a resistividade do meio (Silva, 1999).

O procedimento utilizado na interpretação dos dados gamaespec- Cálculos baseados em Zonge & Hughes (1991) mostram uma
trométricos envolveu: a) a comparação da imagem correspon- profundidade estimada de 20 metros para o canal de 33.000 Hz,
dente a cada canal discriminado (K, U, e Th) com a topografia 50 metros para o canal médio (4175 Hz) e de 100 metros para o
para estudar a influência do relevo sobre estes dados; b) a análise canal de 935 Hz.
do canal de contagem total para a definição de grandes domı́nios O canal de 4175 Hz revelou ser o melhor produto, pois
gamaespectrométricos, bem como a análise da contribuição rela- apresenta alta razão sinal/ruı́do (Silva, 1999) e foi, portanto, o
tiva de cada um dos outros canais; c) o uso de composições em único a ser utilizado na integração de dados. A profundidade de
falsa cor (imagem ternária) do tipo RBG e CMY, com os canais observação foi calculada e estimada em torno de 50 metros.
de K, Th e U, usando ou não o modelo digital de terreno como
uma quarta componente e seu estudo espacial para definir uni- PROCESSAMENTO DIGITAL DE IMAGENS
dades e/ou domı́nios com assinaturas gamaespectrométricas se- AEROGEOFÍSICAS (PDI)
melhantes; d) uso das razões U/Th, Th/K e U/K na análise da sua Neste trabalho utilizaram-se as seguintes técnicas de PDI:
correlação com os diversos domı́nios gamaespectrométricos e Composição Colorida RGB, IHS e Principais Componentes.
com as regiões onde estão mapeadas as mineralizações na região.
Composição colorida RGB
Eletromagnetometria As composições coloridas falsa-cor no espaço RGB utilizadas
para realçar e combinar os dados aerogeofı́sicos são: RGB (K, Th
As imagens com a resistividade aparente foram geradas a partir
e U) e RGB (Amplitude do Sinal Analı́tico, Condutividade-canal
dos valores interpretados existentes. Vale ressaltar que estas ima-
de 4.175 Hz e Canal de Contagem Total) (Figura 4).
gens podem não só apontar corpos condutores importantes para
Na Figura 4, as áreas vermelhas representam regiões com alta
exploração mineral, mas podem ser uma boa ferramenta para ma-
susceptibilidade magnética, baixa condutividade e baixos valores
peamento geológico. Os dados com os resultados das resistivi-
de radiação gama. Já as áreas amarelas representam áreas com
dades calculadas para as freqüências de 935, 4175 e 33000 Hz
alta susceptibilidade magnética e alta condutividade, mas com
foram interpolados utilizando o método bi-direcional com célula
baixa radiação.
de 50 metros e a seguir, micronivelados semelhantemente à
magnetometria. Depois, foram convertidos para condutividade
Sistema IHS
(Silva, 1999).
Na verdade, as assinaturas identificáveis nas imagens com Utilizou-se, também, a transformação IHS para integrar as infor-
os resultados da resistividade, podem ser detalhadas pelas da mações aerogeofı́sicas desse trabalho. Assim, fundiu-se a fase do
condutividade e vice-versa. As imagens de condutividade mos- sinal analı́tico e a primeira derivada vertical do campo magnético
traram pequenas variações nesta propriedade nos litotipos do anômalo com as imagens de amplitude do sinal analı́tico, condu-
Greenstone Belt Rio das Velhas, o que parece ser tı́pico de ter- tividade (4.175 Hz) e composição colorida falsa-cor RGB (K, Th
renos Pré-cambrianos (Palacky, 1991). As imagens com os loga- e U) (ex. Figura 5). Os melhores resultados são aqueles proveni-
ritmos destes valores detalharam melhor estas assinaturas, visto entes da integração dos dados eletromagnéticos e amplitude do
que em rochas pré-cambrianas as diferenças entre os valores de sinal analı́tico, pois realçam domı́nios magnéticos e condutivos
condutividade são muito pequenas (Palacky, 1991; Silva, 1999). que são sabidamente conhecidos por hospedar mineralizações
De acordo com Palacky (1991), as freqüências eletromagnéticas aurı́feras em BIFs fácies óxido na região.
baixas e médias constituem canais melhores para mapear con-
dutores múltiplos e pequenos alvos em profundidade, devido a Principais componentes
problemas associados a ruı́dos. Em terrenos tropicais, como Aplicou-se a análise por principais componentes aos dados ga-
aqueles do presente trabalho, o canal médio mostra uma ima- maespectrométricos (K, Th, U), juntamente com a amplitude do
gem com menor quantidade de ruı́do. O canal de 33.000 Hz mos- sinal analı́tico e o canal de condutividade de 4.175Hz, utilizando
tra forte contribuição antrópica, apesar da alta razão sinal/ruı́do a matriz de covariância. A Tabela 2 apresenta os autovetores ge-
(Silva, 1999). Uma das questões chaves versa sobre a profundi- rados por principais componentes e a Tabela 3 apresenta a matriz
dade das fontes. A quantificação depende de muitos fatores como de correlação entre as variáveis geofı́sicas.

Brazilian Journal of Geophysics, Vol. 24(4), 2006

566 INTEGRAÇÃO E ANÁLISE DE DADOS AEROGEOFÍSICOS

Figura 4 – Imagem da fusão da Amplitude do Sinal Analı́tico e Fase do sinal analı́tico (intensidade).

Tabela 2 – Tabela de autovetores da transformada por principais componentes.

K Th U ASA C(4.175Hz) Autovalores (%)

PC1 –0,64 –0,37 –0,17 +0,47 +0,427 37,70
PC2 +0,03 +0,36 +0,19 +0,79 –0,432 23,90
PC3 –0,56 –0,21 –0,04 –0,27 –0,752 18,65
PC4 +0,51 –0,74 –0,23 +0,24 –0,245 14,40
PC5 –0,02 –0,34 +0,93 –0,03 +0,07 5,3

Tabela 3 – Matriz de correlação das variáveis geofı́sicas.

K Th U ASA C(4.175Hz)
K 1 0,28 0,21 –0,24 –0,19
Th 0,28 1 0,44 –0,05 –0,21
U 0,21 0,44 1 –0,006 –0,22
ASA –0,24 –0,05 –0,006 1 0,09
C(4.175 Hz) –0,19 –0,21 –0,22 0,09 1

Revista Brasileira de Geofı́sica, Vol. 24(4), 2006

 A. AMORIM TEIXEIRA, A. MOREIRA SILVA, A.C. BITTENCOURT PIRES, R.A. VITORIA DE MORAES e C.R. DE SOUZA FILHO 567

Figura 5 – Composição colorida falsa-cor RGB (Amplitude do sinal Analı́tico, Condutividade no canal 4.175Hz e Canal de Contagem Total).

Observa-se pela tabela de autovetores (Tabela 2) que, de uma dos os canais com exceção do canal condutivo. A escolha dessa
forma geral, os canais gamaespectrométricos apresentam peso componente no canal azul é para homogeneizar a composição
semelhante ao longo de todas as principais componentes. Isso dos planos de informação na imagem colorida. Essa composição
ocorre principalmente por causa da afinidade geofı́sica desses colorida, com apenas três bandas, sintetiza 76% da variância
elementos em comparação com as propriedades fı́sicas dos ou- total das informações gamaespectrométricas, magnetométrica
tros planos de informação, como susceptibilidade magnética e e condutiva.
condutividade e pode ser comprovado pelos valores observados
na matriz de correlação (Tabela 3). CLASSIFICAÇÃO NÃO SUPERVISIONADA VIA MÉDIA K
A melhor composição de imagem utilizando as principais A maioria das técnicas hierárquicas pode necessitar da criação
componentes compreende a composição colorida falsa-cor RGB, e manipulação matrizes de similaridade extremamente grandes
para realce litológico, com a PC1 no vermelho, a PC4 no Verde se o número de amostras for elevado, o que é comum nos da-
e a PC2 no Azul (Figura 6). Observa-se pela tabela de autoveto- dos produzidos pelas técnicas de geofı́sica aérea. Por isso, foram
res (Tabela 2) que a PC1 possui maior peso da amplitude do sinal concebidos procedimentos de agrupamento usando um número
analı́tico e da condutividade (+0,47 e +0,42, respectivamente), em limitado de centros arbitrários de grupos visando reduzir a difi-
contrapartida, a contribuição dos canais gamaespectrométricos é culdade computacional. Sem dúvida, um dos mais usados dentre
negativa. Esses domı́nios são estratégicos para exploração mine- estes é o procedimento da média-K (K-means). Nesta técnica,
ral na região do Grupo Nova Lima. Já a PC4 possui peso positivo k pontos caracterizados por m variáveis são designados como
para o canal de potássio e para a amplitude do sinal analı́tico, centróides iniciais dos grupos. Uma matriz de similaridades en-
realçando assim, domı́nios magnéticos com possı́vel alteração tre os k centróides e as n amostras é calculada e as amostras
hidrotermal. Por fim, a PC2 possui contribuição positiva de to- mais similares são agrupadas ao centróide mais próximo. No-

Brazilian Journal of Geophysics, Vol. 24(4), 2006

568 INTEGRAÇÃO E ANÁLISE DE DADOS AEROGEOFÍSICOS

Figura 6 – Composição colorida falsa-cor RGB (PC1, PC4 e PC2).

vos centróides são então calculados e o processo segue intera- A Figura 8 mostra a assinatura geofı́sica das classes resul-
tivamente até todas as amostras haverem sido classificadas. Em tantes da classificação não supervisionada por média-k. Todos
princı́pio, o centróide se deslocará rapidamente em direção ao os planos de informação, bem como todas as classes pertencen-
centro verdadeiro de um grupo em crescimento, na medida em tes aos mesmos, foram reclassificados em relação à sua média
que informações verdadeiras sobre este se sobrepuserem àquela em baixo, baixo-médio, médio, médio-alto e alto.
inicialmente arbitrária.
Neste artigo foram testadas diversas combinações de infor- DISCUSSÕES E CONCLUSÕES
mações visando definir o que melhor caracterizasse os diferen- O melhor produto resultante da transformação IHS para inter-
tes litotipos da região estudada, e conseqüentemente, as rochas pretação estrutural é a amplitude do sinal analı́tico com a fase
hospedeiras da mineralização aurı́fera. Os dados foram exausti- do sinal analı́tico no canal de intensidade (Figura 5). Assim, ao
vamente testados para diferentes classes e diferentes interações. mesmo tempo em que se identifica a posição da fonte magnética,
As classificações estatı́sticas efetuadas mostraram que os dados pode-se identificar uma assinatura do mergulho do gradiente da
possuem em média dez classes. Portanto, foi efetuado um grupa- anomalia. Observa-se que a Zona de Cisalhamento São Vicente
mento em dez classes, após dez interações utilizando as seguintes varia o mergulho de falha ao longo da sua extensão, o que é vali-
informações: canais de potássio, tório, urânio, bem como a am- dado pelos trabalhos anteriores na área (Araújo, 2001). As áreas
plitude do sinal analı́tico e o canal de condutividade (4.175 Hz). mineralizadas ao longo da falha apresentam assinatura padrão
Todos os dados foram reescalonados para 8 bits, ou 256 tons de com valores de fase do sinal analı́tico próximo a zero. As es-
cinza (Figura 7). timativas de profundidade das fontes magnéticas por meio da

Revista Brasileira de Geofı́sica, Vol. 24(4), 2006

 A. AMORIM TEIXEIRA, A. MOREIRA SILVA, A.C. BITTENCOURT PIRES, R.A. VITORIA DE MORAES e C.R. DE SOUZA FILHO 569

Figura 7 – Classes resultantes da aplicação da técnica de classificação não supervisionada por média-k. Legenda:
Classe 1: vermelho, Classe 2: verde claro, Classe 3: azul marinho, Classe 4: amarelo, Classe 5: ciano, Classe 6:
magenta, Classe 7: marrom escuro, Classe 8: verde escuro, Classe 9: roxo e Classe 10: marrom claro.

deconvolução de Euler auxiliam o entendimento do comporta- ultramáficos com textura spinifex e sills cumuláticos, correla-
mento geométrico dessa falha em profundidade. No geral, a fa- ciona com a classe 1 e apresenta baixos valores radiométricos e
lha apresenta pouca profundidade (200 m a 400 m) com algumas alta condutividade. Localmente, correlaciona-se a classe 2 com
áreas profundas (400 m a 800 m). A análise qualitativa da fusão essa formação onde ela é mais magnética, porém, não condutiva.
Fase/Amplitude do sinal analı́tico (Figura 5) e da deconvolução A Formação Morro Vermelho caracteriza a classe 8 e diferencia-se
de Euler (Figura 3) permite a determinação de zonas homólogas da classe 1 por possuir médios valores de radiação, porém, muito
a partir da textura da fase do sinal analı́tico, dos valores de am- condutiva. Essa formação caracteriza-se por pequenas lentes de
plitude do sinal analı́tico e dos lineamentos das estimativas de sedimentos quı́micos exalativos como BIFs do tipo algoma, me-
profundidade da deconvolução de Euler. tacherts e xistos carbonáticos. A classe 7 correlaciona-se com
A partir das classes obtidas pela classificação por média-k a Formação Ribeirão Vermelho composta por rochas metassedi-
dos dados aerogeofı́sicos, pôde-se associar a assinatura de al- mentares vulcanoclásticas e constituem-se de tufos dacı́ticos pi-
gumas classes (Figuras 7 e 8) com unidades litoestatigráficas do roclásticos e aglomerados de baixo fluxo. Essa classe caracteriza-
Supergrupo Rio das Velhas mapeadas no projeto Rio das Velhas se por altos valores de radiação de K, Th e U e por médio a bai-
(Pinto, 1996) na escala 1:100.000. Assim, os BIFs, de uma forma xos valores magnéticos e condutivos. A Formação Córrego do
geral, se correlacionam com a classe 10, caracterizada princi- Sı́tio representa as rochas metassedimentares de fontes mistas do
palmente por baixos valores de radiação e alta amplitude do si- Grupo Nova Lima (Pinto & Silva, 1996) e caracteriza-se por finas
nal analı́tico. A Formação Quebra-Osso, composta por magmas intercalações de quartzo-sericita-clorita xisto e feldspato-quartzo

Brazilian Journal of Geophysics, Vol. 24(4), 2006

570 INTEGRAÇÃO E ANÁLISE DE DADOS AEROGEOFÍSICOS

Figura 8 – Assinatura Geofı́sica das classes resultantes da classificação não supervisionada por média-k.

xisto interpretados como metapelitos e metagrauvacas com fi- CARNEIRO MA. 1992. O Complexo Metamórfico Bonfim Setentrio-
nas intercalações de xistos carbonáticos e BIFs. Essa unidade nal (Quadrilátero Ferrı́fero, Minas Gerais): Litoestratigrafia e Evolução
correlaciona-se com as classes 4 e 5. Na parte norte há uma Geológica de um Segmento de Crosta Continental do Arqueano. São
predominância da classe 4 que possui assinatura com médios a Paulo, 233p. (Tese de Doutoramento, Instituto de Geociências da Uni-
versidade de São Paulo).
médio-baixos valores dos dados geofı́sicos. Porém, na região
ao longo da Zona de cisalhamento São Vicente que apresenta CARNEIRO MA, NOCE CM & TEIXEIRA W. 1995. Evolução policı́clica
mineralizações do tipo load , há a predominância da classe 5, que do Quadrilátero Ferrı́fero: uma análise fundamentada no conhecimento
reflete, provavelmente, a alteração hidrotermal da classe 4 carac- atual da geocronologia U-Pb e geoquı́mica isotópica Sm-Nd. Revista da
terizada pelo alto valor relativo de K e baixo valor de Tório. Pode- Escola de Minas, 48 (4): 264-273.

se correlacionar a classe 5 com a PC4 da Análise por Principais CHEMALE JR F, ROSIÈRE CA & ENDO I. 1991. Evolução tectônica do
Componentes dos dados aerogeofı́sicos descrita anteriormente. Quadrilátero Ferrı́fero, Minas Gerais – Um modelo (Tectonic Evolution
of the Quadrilátero Ferrı́fero, Minas Gerais. A Model). Pesquisas, 18:
104–127.
REFERÊNCIAS
CHEMALE JR F, ROSIÈRE CA & ENDO I. 1994. Tectonic Evolution of the
ALKMIM FF & MARSHAK S. 1998. Transamazonian Orogeny in the Quadrilátero Ferrı́fero, MG, Brazil. Precambrian Research, 65: 25–54.
Southern São Francisco Craton Region, Minas Gerais, Brazil: evidence CORDANI UG, KAWASHITA K, MUELLER G, QUADE H, REIMER V &
for Paleoproterozoic collision and collapse in the Quadrilátero Ferrı́fero. ROESER HP. 1980. Interpretação tectônica e petrológica de dados geo-
Precambrian Research, 90: 29–58. cronológicos do embasamento do bordo sudeste do Quadrilátero Fer-
rı́fero, Minas Gerais. Anais da Academia Brasileira de Ciências, Rio de
ALMEIDA FFM. 1977. O Cráton São Francisco. RBG, 7(4): 349–364.
Janeiro, 52: 785–799.
ARAÚJO JGM. 2001. Influência das zonas de cisalhamento de São Vi- CORDANI UG, TEIXEIRA W, SIGA JR O. 1989. Geocronologia do Qua-
cente e Tapera na mineralização aurı́fera do Quadrilátero Ferrı́fero, Mi- drilátero Ferrı́fero. In: XXI Semana de Estudos de Ouro Preto, SICEG,
nas Gerais. Instituto de Geociências. Universidade de Brası́lia. Brası́lia. Ouro Preto, Bol. 21: 27–44.
Dissertação de mestrado. 100 p
CROCCO-RODRIGUES FA, RIBEIRO-RODRIGUES LC & ROSIERE CA.
BARBOSA ALM. 1968. Contribuições recentes à geologia do Qua- 1991. Análise crı́tica sobre a controvérsia do Grupo Tamanduá, Qua-
drilátero Ferrı́fero. IGC-UFMG (1985), 44p. drilátero Ferrı́fero. In: 6◦ Simpósio de Geologia de Minas Gerais/ XXX

Revista Brasileira de Geofı́sica, Vol. 24(4), 2006

 A. AMORIM TEIXEIRA, A. MOREIRA SILVA, A.C. BITTENCOURT PIRES, R.A. VITORIA DE MORAES e C.R. DE SOUZA FILHO 571

Semana de Estudos de Ouro Preto, SICEG, Ouro Preto. Anais ..., Revista NOCE CM. 1995. Geocronologia dos Eventos Magmáticos, Sedimenta-
da Escola de Minas, 45(1/2): 47–48. res e Metamórficos do Quadrilátero Ferrı́fero, Minas Gerais. São Paulo.
128p. Instituto de Geociências da Universidade de São Paulo. Tese de
DORR II JVN. 1969. Physiographic, Stratigraphic, and Structural De-
Doutorado.
velopment of the Quadrilátero Ferrı́fero, Minas Gerais. U.S.G.S. Paper,
641-A: 1–110. NOCE CM, PINHEIRO SO, LADEIRA EA, FRANCA CR & KATTAH S.
1992. A seqüência vulcanossedimentar do Grupo Nova Lima na região de
FREITAS ME, CROCCO-RODRIGUES FA & ROSIÈRE CA. 1991. A cor-
Piedade do Paraopeba, oeste do Quadrilátero Ferrı́fero, Minas Gerais.
relação do Quartzito Cambotas com a Seqüência Espinhaço apoiada por
Revista Brasileira de Geociências, 22(2): 175–183.
mapeamento geológico em escala 1:10.000. In: 6◦ Simpósio de Geo-
logia de Minas Gerais / XXX Semana de Estudos de Ouro Preto, SICEG. PALACKY GJ. 1991. Resistivity characteristics of geologic targets. In:
Anais..., Revista da Escola de Minas, 45(1/2): 51–53. NABIGHIAN MN (Ed.). Electromagnetic Methods in Applied Geophysics.
Society of Exploration Geophysicists, Tulsa, Oklahoma, 1: 53–129.
FUCK RF. 2001. Reprocessamento dos Dados Aerogeofı́sicos do Projeto
Rio das Velhas. Relatório bolsa RHAE. CNPq. PINTO PC. 1996. Explanatory notes / Rio das Velhas project. Rel. téc.
DNPM – SUREG/BH.
HERTZ N. 1978. Metamorphic rocks of the Quadrilátero Ferrı́fero, Minas
Gerais, Brazil. U.S.G.S. Professional Paper, 641-C: 1–81. PINTO PC & SILVA LC. 1996. Contrasting Tectono-Stratigraphic Do-
mains in the Rio das Velhas Greenstone Belt (RVGB), MG, Brazil. In: SBG,
HILDENBRAND JD & PEREZ DA GAMA MF. 1993. Levantamento Aéreo
Symposium Archaean Terranes of the South American Platform. Brası́lia,
Magnetométrico, Gamaespectrométrico e Eletromagnético. Aquisição,
Anais, 23–25.
Processamento e Interpretação dos Dados. Projeto Rio das Velhas. Re-
latório Final. Prospec S/A (inédito), 1: 47. REID AB, ALLSOP JM, GRANSER H, MILLETT AJ & SOMERTON IW.
1990. Magnetic interpretation in three dimensions using Euler decon-
HOEFS J, MULLER G & SCHUSTER AK. 1982. Polymetamorphic re-
volution. 58th Annual International Meeting. Society of Exploration Geo-
lations in iron ores from the Iron Quadrangle, Brazil: the correlation
physicists.
of oxygen isotope variations with deformation history. Contributions to
Mineralogy and Petrology, 79: 241–251. REIS LA, MARTINS-NETO MA, GOMES NS, ENDO I & JORDT-EVAN-
GELISTA H. 2002. A bacia de antepaı́s paleoproterozóica Sabará, Qua-
JORDT-EVANGELISTA H, ALKMIM FF & MARSHAK S. 1991. Metamor-
drilátero Ferrı́fero, Minas Gerais. Revista Brasileira de Geologia, 32(1):
fismo progressivo e a ocorrência dos três polimorfos de Al2 O3 SiO5 (cia-
27–42.
nita, andaluzita e silimanita) na Formação Sabará em Ibirité, Quadrilátero
Ferrı́fero, MG. In: 6◦ Simpósio de Geologia de Minas Gerais/ XXX ROMANO AW. 1989. Evolution tectonique de la région NW du Quadri-
Semana de Estudos de Ouro Preto, SICEG. Anais ..., Revista da Escola latère Ferrifère – Minas Gerais, Brésil. Nancy, 259p. (Thèse Docteur,
de Minas, 45(1 e 2): 157–160. Université de Nancy I).

LADEIRA EA. 1980. Metallogenesis of Gold at the Morro Velho Mine, SCHRANK A & SILVA MG. 1993. Os greenstone belts do Cráton do São
and in Nova Lima District, Quadrilátero Ferrı́fero, Minas Gerais, Brazil. Francisco, Brasil. In: DOMINGUEZ JML & MISI A (Ed.). O Cráton do
Unpubl. PhD thesis, University of Western Ontario, London, 272 pp. São Francisco. Salvador, SBG/SGM/CNPq, pp. 85–118.

LADEIRA EA & VIVEIROS JFM. 1984. Hipótese sobre a estruturação do SCHORSCHER HD. 1978. Komatiitos na estrutura greenstone belt Série
Quadrilátero Ferrı́fero com base em dados disponı́veis. Boletim da So- Rio das Velhas. Quadrilátero Ferrı́fero, Minas Gerais, Brasil. XXX Congr.
ciedade Brasileira de Geologia-Núcleo Minas Gerais, 4: 1–14. Bras. Geol., Anais... Recife, Bol. Esp. 1: 292–293.

MACHADO N & CARNEIRO MA. 1992. U-Pb evidence of late Archaean SILVA AM. 1999. Geophysical and Geological Data Integration Using a
tectono-thermal activity in the southern São Francisco shield, Brazil. Ca- New Statistical Approach for Mineral Target Selection Applied to the Rio
nadian Journal of Earth Sciences, 29(11): 2341–2346. das Velhas Greenstone Belt, Quadrilátero Ferrı́fero. Inst. de Geociências,
Universidade de Brası́lia, Brası́lia, Tese de Doutoramento.
MACHADO N, NOCE CM, BELO DE OLIVEIRA OA & LADEIRA EA. 1989.
Evolução geológica do Quadrilátero Ferrı́fero Arqueano e Proterozóico SILVA AM, CHEMALE JR F & KUYUMJIAN RM. 1995. O Ambiente
Inferior, com base em geocronologia U/Pb. In: 5◦ Simpósio de Geologia tectônico dos diques máficos do Quadrilátero Ferrı́fero e Espinhaço Me-
do Núcleo Minas Gerais – 1◦ Simpósio de Geologia do Núcleo. ridional e suas implicações para a evolução do Cráton São Francisco.
In: Simpósio Nacional de Estudos Tectônicos, 5, Gramado. SBG-CPGQ-
MARSHAK S & ALKMIM FF. 1989. Proterozoic Extension/Contraction
UFRGS, 1996. Boletim de Resumos Expandidos..., p. 209–210.
Tectonics of the São Francisco Craton and Adjacent Regions, Minas Ge-
rais, Brazil: A Kinematic Model Relating Quadrilátero Ferrı́fero, São Fran- TEIXEIRA W. 1985. A Evolução Geotectônica da Porção Meridional do
cisco Basin and Cordilheira do Espinhaço. Tectonics, 8 (3): 555–571. Cráton São Francisco, com Base em Interpretações Geocronológicas.

Brazilian Journal of Geophysics, Vol. 24(4), 2006

572 INTEGRAÇÃO E ANÁLISE DE DADOS AEROGEOFÍSICOS

São Paulo. 207p. (Tese de Doutorado, Instituto de Geociências da Uni- in Applied Geophysics. Society of Exploration Geophysicists, Tulsa,
versidade de São Paulo). Oklahoma, 2: 713–810.
TEIXEIRA W, CARNEIRO MA, NOCE CM, MACHADO N, SATO K &
TAYLOR PN. 1996. Pb, Sr and Nd isotope constrains on the Archaean ZUCCHETTI M, BALTAZAR OF & RAPOSO FO. 1996. Estratigrafia. Pro-
evolution of gneissic-granitoid complexes in São Francisco Craton, Bra- jeto Rio das Velhas. Mapa Geológico Integrado. Escala: 1:100000. Texto
zil. Precambrian Research, 78: 151–164. Explicativo. Programa de Estudos Mineiros. Convênio DNPM-CPRM.
ZONGE KL & HUGHES LJ. 1991. Controlled Source Audio-Frequency DNPM-Secretaria de Minas e Energia-Ministério de Minas e Energia.
Magnetotellurics. In: NABIGHIAN MN (Ed.). Electromagnetic Methods Belo Horizonte, 1996, p. 13–42.

NOTAS SOBRE OS AUTORES

Alexandre de Amorim Teixeira. Geólogo pela Universidade de Brası́lia – UnB (2000) e Mestre (2003) em processamento de dados em geologia e análise ambien-
tal pela Universidade de Brası́lia. Atualmente, é especialista em geoprocessamento pela Agência Nacional de Águas – ANA.

Adalene Moreira Silva. Engenheira Geóloga (1989) graduada pela Universidade Federal de Ouro Preto, Mestre (1992) e Doutora em Geologia (1999) pela Uni-
versidade de Brası́lia. Profissionalmente, atuou como pesquisadora do Serviço Geológico Americano (1997-1999), Professora Visitante da Universidade de Brası́lia
(2000-2002), Pesquisadora (05/2002-04/2003) e Professora do Instituto de Geociências da UNICAMP (04/2003-10/2005). Atualmente é Professora do Instituto de
Geociências da Universidade de Brası́lia e vem pesquisando sobre novas técnicas de processamento, interpretação e integração de dados aplicados a exploração mine-
ral, hidrogeologia e ambiente.

Augusto César Bittencourt Pires. Geólogo e geofı́sico, formado em 1968 pela UFRJ. Ph.D. em Geofı́sica pela Colorado School of Mines em 1975 e com Pós-
Doutorado em Geofı́sica pela mesma instituição. Foi responsável pelo departamento de Geofı́sica da ENCAL; Chefe do Departamento de Geofı́sica do Observatório
Nacional; Chefe do Centro de Geofı́sica Aplicada do DNPM. No CNPq exerceu funções de Supervisor a Diretor de Programas. Foi Diretor de Programas da CAPES
e Coordenador de Polı́ticas e Planos do MCT. Foi Consultor Legislativo do Senado. É Professor Titular, Coordenador do Laboratório de Geofı́sica Aplicada da UnB e
consultor sênior da HGeo/ InterGeo.

Roberto Alexandre Vitória de Moraes. Geólogo (1968) pela UFPe; Ph.D. em Geofı́sica (1997) Colorado School of Mines, Golden, CO, EUA. Geofı́sico do DNPM,
Rio de Janeiro-RJ; geólogo/geofı́sico da DOCEGEO, Belo Horizonte-MG; geofı́sico-chefe da PROSPEC S/A, Rio de Janeiro-RJ; geofı́sico-chefe do SEMEPO da Paulipetro,
São Paulo-SP; geofı́sico-chefe da Petromisa, Rio de Janeiro-RJ; geofı́sico sênior da Petrobras, Rio de Janeiro-RJ, pesquisador sênior do IPT, São Paulo-SP e Professor
Adjunto do IG da UnB. Especialista nos métodos magnético, gravimétrico, elétrico, eletromagnético (aéreos e terrestres), gamaespectrométrico e em métodos numéricos
para modelagem e inversão em Geofı́sica. Diretor Cientı́fico da HGeo e da InterGeo.

Carlos Roberto de Souza Filho. Engenheiro Geólogo (Universidade Federal de Ouro Preto, 1988), Mestre em Metalogênese (UNICAMP, 1991), PhD (Open Uni-
versity, Inglaterra, 1995, Jovem-Pesquisador (FAPESP-UNICAMP, 1995-1997), Professor-Doutor (UNICAMP, 1997-2002). Desde 2002, é Professor Livre Docente do
Departamento de Geologia e Recursos Naturais do Instituto de Geociências da UNICAMP. Atualmente é Coordenador do Programa de Pós-graduação de Geociências e
lı́der do Grupo de Geotecnologias do IG. Pesquisador Nı́vel 1 do CNPq, é responsável pelos laboratórios de pesquisa em Fs-LA-ICP-MS (IsoGeo), Espectroscopia de
Reflectância (LER) e Processamento de Informações Georreferenciadas (LAPIG). É pesquisador da NASA nos programas ASTER e Mars-Earth Analogs, Editor-Associado
das Revistas Computers & Geosciences e Mathematical Geology.

Revista Brasileira de Geofı́sica, Vol. 24(4), 2006