Pular para o conteúdo principal

Fotometria do cometa Ison (Novembro de 2013)

Fig. 1 Imagem de 14 de Novembro com a presença do cometa Ison. Problema: com base nesta imagem, qual a magnitude do cometa? Imagem obtida em Barão Geraldo, Campinas, SP.
Finalmente, o dia 14 de Novembro permitiu o primeiro registro fotográfico do cometa Ison, como visto desde Campinas, SP (figura 1). A imagem mostra uma pequena mancha esverdeada em meio às estrelas da constelação da Virgem. Esse registro fotográfico, por mais simples que pareça, contém em si informação suficiente para determinar a magnitude (ou brilho) do cometa. 

Vamos mostrar aqui um método prático (fazendo jus ao nome do blog) de se estimar essa magnitude. A imagem da Figura 1 já identifica diversas estrelas na vizinhança do cometa. A magnitude dessas estrelas é conhecida. Como podemos fazer para estimar a do Ison? O método aqui descrito é semelhante à estimativas feitas pelo olho, quando se compara o brilho a ser inferido com duas estrelas de magnitude maior e menor que o do astro, obtendo-se um resultado por interpolação.

Sem muitos detalhes, passamos a descrever o método, que é baseado na extração de informação de brilho contida na imagem. Antes de tudo, convém listar as magnitudes (como dada pelo Stellarium) para as estrelas listadas na Fig. 1, em ordem decrescente de brilho:

m(X Vir) = 4.65
m( ψ Vir) = 4.75
m(g Vir) = 5.55
m(HIP 62743) = 6.45
m(HIP 63225) = 6.9
m(HIP 63240a)=7.25

Observe que colhemos uma amostra de estrelas com ampla gama de variação e sabemos que o cometa tem brilho no intervalo 4.65-7.25 claramente por inspeção visual da Fig. 1.

Extração da amostra de brilho de sensibilização da imagem.

Antigamente, a medida de brilho em placas contendo uma imagem era feita por meio de análise de microscópio. Cada estrela na imagem da Fig.1 era observada sob um microscópio e o analista media a área da "mancha" sensibilizada no filme. Como temos uma imagem digital, podemos fazer isso diretamente sobre a matriz da intensidades. Para tanto, convertemos a imagem da Fig.1 (que é colorida) em um padrão preto e branco (onde os tons de cinza vão de 0 (preto) a 255 (branco), ou formato 8-bits).

Extraímos um pedaço da imagem por meio de um algoritmo de análise. No nosso caso, usamos o software MathCad que permite a manipulação de imagens inteiras ou pedaços dela. Por exemplo, a Fig. 2 traz os mapas em tons de cinza do cometa na foto da Fig.1, de ψ Vir e de X Vir para um quadro com 20 X 20 pixels. Note que a imagem não aparece circular mas ovalada, por causa da exposição de 8 segundos e o uso de tripé sem acompanhamento.
Fig. 2 Matrizes de intensidade de amostras da imagem da Fig.1 convertida em preto e branco, centrado em cada elemento indicado e cobrindo uma área de 20 X 20 pixels. A estrutura 'ovalada' se deve ao movimento do objeto durante a exposição de 8 segundos.
Há uma outra maneira de representar o 'mapa' de intensidade da fig. 2. Através de um gráfico 2d de curvas de intensidade ou superfícies 3d de alturas. Nesse último caso, cada pixel representa um paralelepípedo com altura igual ao tom de cinza do pixel. Essa representação pode ser vista na fig. 3 à direita, junto com o mapa de contorno de intensidades em 2d para o cometa Ison (ou seja, a representação de tons de pixel do primeiro frame da Fig. 2).
Fig. 3 Esquerda: curvas de nível de intensidade para o cometa Ison (mapa 2d). Direita: representação 3d onde cada tom do pixel corresponde a uma altura. O brilho é dado, por exemplo, pela soma das alturas de todos os 'blocos' que tenham altura acima de um limiar escolhido.
O brilho fotográfico do objeto pode ser tomado como sendo igual a soma de todas as 'alturas' no mapa de intensidade que estão acima de certo limiar. No nosso caso, escolhemos esse limiar de forma que a função de brilho seja proporcional (ou inversamente proporcional) à magnitude das estrelas escolhidas.

O cálculo do brilho pode ser feito por meio da seguinte função, que chamamos de meas(image,limiar) (escrita em código MathCad, aqui image é a matriz de intensidade de entrada):


 Chamamos de B o brilho extraído a partir da contagem dos pixels submetidos à filtragem do limiar (ou seja, o valor resultado de meas(image, limiar)). Para a fig. 1, o resultado com limiar = 107 foi:

B(X Vir) = 5379
B( ψ Vir) = 5379
B(g Vir) = 3480
B(HIP 62743) = 1312
B(HIP 63225) = 915
B(HIP 63240a)=107

A curva de calibração da imagem e, finalmente, a magnitude do Ison

O próximo passo consistem em plotar em um gráfico os valores de B versus a magnitude (fornecida anteriormente). Cada ponto azul da Fig. 4 é uma estrela na lista usada para calibração. Assim, no eixo x podemos ler a magnitude visual usada e na ordenada o brilho (conforme calculado acima para a função de brilho). 

Fig. 4 Reta de interpolação dos dados (função de brilho versus magnitude). Com base na curva que melhor interpola os dados, a magnitude do Ison pode ser calculada a partir de seu brilho estimado da imagem.
A reta tem a equação B(mag) = a + b*mag, onde:

a = 15100
b = -2082

Esses valores são extraídos a partir de uma interpolação linear feita sobre os pontos do gráfico da Fig. 4. Pela aplicação do superfície de intensidade do cometa (fig. 3) usando o mesmo limiar obtemos:

B(Ison)= 3233

Portanto, a magnitude do cometa será dada por:

mag(Ison)=[B(Ison)-a] / b

Substituindo os valores na equação obtemos:

mag(Ison) = 5.7.

Essa é a magnitude estimada com base na imagem obtida da Fig. 1.

Premissas e observações 

Alguns comentários são importantes sobre o que aqui descrevemos:
  • Usamos valores de magnitude visual, enquanto que a CCD da câmera fornece brilho proporcional à magnitude fotográfica;
  • Há diferenças, portanto, entre as magnitudes e as correspondentes funções de brilho por causa das cores dos objetos fotografados (para o olho humano um objeto com certa cor pode ser mais brilhante do que para a CCD);
  • Se mais de uma imagem fosse novamente analisada, é provável que o brilho estimado não seria exatamente igual ao obtido. Isso ocorre por causa de erros no registro, flutuações de intensidade incontroláveis no CCD;
  • O processo de conversão para preto e branco altera as intensidades de entrada;
  • Podemos usar outras funções para estimativa do brilho, mas a relação entre a magnitude e o brilho pode não ser linear como mostrado na Fig. 4;
  • Se usarmos outro valor de Limiar, p. ex, 106, a magnitude obtida é igual a 5,68 (ao invés de 5,70). O método é sensível ao valor de intensidade usado, mas a variação está no intervalo +/-0.03;
  • O correto seria tomar várias imagens e aplicar, de forma automática, o método para cada uma delas. Tomar a média e estimar o erro. 

Nosso objetivo aqui foi apenas mostrar como é feito - de forma grosseira - as estimativas de magnitude sem dependência com o olho humano. Esse método pode ser aplicado para outros objetos, por exemplo, estrelas variáveis, com vantagens porque é uma estimativa totalmente numérica e não depende de julgamento psicológico do observador.

Nossa estimativa para a magnitude do cometa Ison em 14/11/13 as 07:15 UT foi, portanto, 5,7+/-0,03.

Referências






Comentários

Postagens mais visitadas deste blog

As Três Marias

Aspecto do Cinturão de Órion com as "Três Marias" conforme visto em um pequeno binóculo. Por Roberto Mura via Wikipedia. Ver também: A constelação de Órion (segundo gregos e egípcios) N ão é incomum a curiosidade em torno das famosas "Três Marias". Que agrupamento de estrelas é esse? Não podem ser encontradas com essa designação em nenhum mapa celeste, então elas se referem a quê? Na verdade, as constelações são designadas sob diversos nomes, mas apenas alguns são conhecidos como "oficiais". As "Três Marias" é um caso não oficial. Trata-se de, como é fácil de se pesquisar na internet, de um nome dado ao " Cinturão de Órion ", na constelação de mesmo nome. O "Dicionário Enciclopédico de Astronomia e Astronáutica" (1) de R. R. de Freitas Mourão, assim define esse nome: Três Marias . Asterismo na constelação de Órion, formado por três estrelas brilhantes, em linha reta, e igualmente espaçadas; Três Irmãs, Três Rei...

A questão dos nomes astronomia, astrofilia e o significado da astronomia amadora.

D. Pedro II (1825-1891) O dia 2 de Dezembro é conhecido nacionalmente como 'dia do astrônomo'. Esse é também o dia de nascimento de Pedro de Alcântara ou D. Pedro II que é reconhecido (merecidamente) como patrono da astronomia brasileira. Sempre nesse dia surge a dúvida se a comemoração deve incluir não somente os astrônomos profissionais mas os amadores também. Para tentar resolver a dúvida, devemos antes discutir um pouco sobre o significado próprio dos nomes usados para designar a ciência da astronomia e a prática da observação do céus "sem compromisso", como realizada por amadores. Por que as ciências que estudam fenômenos, ocorrências relacionados à vida se chama "biologia", e aquela que estuda fenômenos e ocorrências relacionados aos astros se chama "astronomia"? Além disso, não se fala em uma 'física amadora' ou uma 'biologia amadora'; não se ouve falar (em clubes) de 'físicos amadores' ou 'biólogos a...

A constelação de Órion (segundo gregos e egípcios)

Uma representação de Órion como visto no hemisfério sul, pelo sotware Stellarium . Alguns dizem que a constelaçao de Canis Major é de Órion, porque esse era caçador e o cão foi colocado junto a ele no céu. (Pseudo-Hyginus, Astronomica 2. 35, [1]) Órion na mitologia grega Segundo a mitologia grega, Órion (Ωριων), não se sabe se filho de Ireus ou do deus Poseidon com Euriale, natural da Beócia, era um caçador gigante muito bonito conhecido dos habitantes daquela região como Caldaon. Tendo peregrinado até a ilha de Quios, infestada de animais perigosos, sob influência da ninfa Hélice apaixonou-se por Mérope, filha do rei Enopion. Depois de caçar e matar todas as bestas de Quios, trouxe os espolios da caça aos pés da princesa como presente. Entretanto, o rei não queria o casamento. Depois de várias tentativas de aproximação da princesa, Órion conseguiu entrar no quarto de Mérope. Enopion pediu ajuda ao deus Dionísio que fez com que Órion caisse, sob influência dos sátiros...