Como o empilhamento realmente funciona (e por que mais frames nem sempre é melhor)

Pergunte por que empilhamos astrofotografias e normalmente vão dizer “para reduzir o ruído”. É verdade e quase inútil, porque não diz quantos frames você precisa, por que a sua pilha de 200 parece pior que a pilha de 60 de outra pessoa, ou por que adicionar mais uma hora às vezes não muda absolutamente nada.

O argumento de verdade é curto, e uma vez que você o tem, a maioria das decisões de empilhamento se responde sozinha.

Sinal soma. Ruído não.

Aponte o telescópio para uma nebulosa e faça uma exposição. Cada pixel recebe um certo número de fótons do céu. Esse número não é fixo - a chegada de fótons é aleatória, e a contagem varia de exposição para exposição em aproximadamente a própria raiz quadrada. Colete 100 fótons e a variação entre frames é de uns 10. Colete 10.000 e é de uns 100.

Agora empilhe N frames. O sinal - a contagem real de fótons do objeto - soma linearmente: N vezes mais. O ruído é aleatório e se cancela em parte, crescendo apenas como a raiz quadrada de N.

Então a relação sinal/ruído melhora em √N.

Essa única relação explica a maior parte do que você observa:

• 4 frames é duas vezes melhor que 1.
• 16 frames é duas vezes melhor que 4.
• 64 frames é duas vezes melhor que 16.

Cada dobro de qualidade custa quatro vezes mais frames. É por isso que a diferença entre 10 e 40 frames é dramática e a diferença entre 200 e 260 é invisível. Você não está fazendo nada errado quando a segunda decepciona - você adicionou 30 % a mais de frames por 14 % de melhoria, e 14 % está abaixo do que o olho percebe com confiança.

Se quer um salto visível a partir de uma pilha de 200 frames, precisa de uns 800.

É por isso que o ruído de leitura decide a duração dos subs

Existe uma segunda fonte de ruído que não se comporta assim. Toda vez que o sensor é lido, a eletrônica adiciona uma penalidade fixa. Ela não depende da duração da exposição - um frame de 5 segundos e um de 500 pagam exatamente uma vez cada.

Então um dado tempo total de integração carrega ruído de leitura em proporção a em quantos frames você o divide. Noventa minutos em 540 subs de dez segundos pagam a penalidade 540 vezes. Os mesmos noventa minutos em 18 subs de cinco minutos pagam 18 vezes.

Esse é todo o argumento a favor de subexposições mais longas, e é por isso que um setup guiado se distancia de um telescópio inteligente em alvos fracos mesmo com tempo total idêntico. Não é que o equipamento guiado coleta mais luz. Ele lê o sensor trinta vezes menos.

A contrapressão é que subs longos precisam de bom rastreamento, e um frame de 300 segundos com rastros não vale nada. A duração certa é a mais longa que a sua montagem entrega de forma confiável - uma questão mecânica, não de processamento.

O alinhamento é onde as pilhas realmente se perdem

O argumento √N assume que cada frame contribui o mesmo sinal para o mesmo pixel. Isso exige que os frames estejam registrados, e é aqui que a maioria das pilhas decepcionantes falha de verdade.

Os frames se movem entre exposições. Deriva da montagem, refração perto do horizonte, flexão e - em qualquer montagem altazimutal, o que inclui todo telescópio inteligente - rotação de campo. Ao longo de uma hora, frames do início e do fim podem estar rotacionados dezenas de graus entre si.

Um alinhamento que só lida com translação não dá conta. O pior é como ele falha: não com uma mensagem de erro, mas com uma pilha em que tudo está levemente borrado. Parece foco ruim ou seeing ruim, e as pessoas passam uma temporada perseguindo o problema errado.

Um registro adequado trabalha a partir das próprias estrelas. Constrói triângulos com as estrelas mais brilhantes de cada frame e os descreve pelas razões dos lados - quantidades que sobrevivem a translação, rotação e escala. Triângulos correspondentes votam em qual estrela corresponde a qual, um ajuste robusto descarta os votos discordantes e uma passada de refinamento leva o resultado a precisão subpixel.

Essa última parte importa. Precisão subpixel não é luxo; um erro consistente de meio pixel na sua pilha custa o mesmo detalhe que teria sido perdido com um seeing um pouco pior.

Rejeição: a parte que remove o que a média não consegue

A média lida com ruído aleatório. Não lida com o que não é aleatório: satélites, aviões, raios cósmicos, faróis de carro varrendo o campo.

Faça a média de 100 frames onde um tem um rastro de satélite e o rastro ainda está lá com 1 % de intensidade - fraco, mas visível contra um céu escuro, e sobrevive ao esticamento.

A rejeição o remove olhando ao longo da pilha em cada pixel e descartando valores que não pertencem ali. Os métodos diferem no que “não pertencem” significa, e escolher entre eles não é arbitrário:

Sigma clipping pega a média e o desvio padrão em cada pixel e descarta valores além de um limiar. Funciona bem e é o padrão por um bom motivo - mas assume que cada frame amostra a mesma distribuição.

Sigma clipping winsorizado puxa os outliers em direção ao grosso em vez de removê-los. Mais seguro em pilhas rasas, onde estimar um desvio padrão a partir de oito amostras e depois descartar valores tende a jogar fora sinal real.

Linear-fit clipping responde a uma pergunta diferente. Ordena os valores de cada pixel ao longo da pilha e ajusta uma reta. Uma dispersão suave - do tipo produzido pela transparência mudando ao longo da noite - é esperada por esse modelo e não custa nada, enquanto um rastro de satélite continua sendo um valor longe da reta.

Essa distinção importa mais do que parece. Quando nuvens altas chegam, ou seu alvo desce em direção à cúpula de luz, os frames genuinamente diferem. A dispersão é dado real, não ruído. O sigma clipping vê os frames nas pontas da faixa como outliers e começa a descartá-los - e os da ponta escura costumam ser os frames mais limpos da sessão. Você perde seus melhores dados para um método que assumiu que estavam errados.

O sinal a observar é quanto o nível de céu por frame deriva. Se é estável, sigma clipping é o certo. Se se move mais que alguns pontos percentuais, um ajuste linear vai servir melhor. O Akastroid mede isso e troca automaticamente acima de cerca de 8 % de deriva, uma decisão que não exige conhecimento de astronomia para acertar e exige bastante para acertar à mão.

Pilhas rasas precisam de tratamento diferente

Com oito frames você não consegue estimar uma distribuição. Qualquer método de rejeição está trabalhando com poucas amostras, e um clipping agressivo vai remover sinal real enquanto afirma remover outliers.

Abaixo de uns oito frames, média simples ou uma mediana costuma ser a escolha honesta. Entre oito e quinze, winsorizar. Acima disso, o clipping de verdade fica seguro.

O modo de falha aqui é sutil: rejeição agressiva numa pilha rasa produz uma imagem que parece limpa e à qual falta, silenciosamente, detalhe fraco. Não parece quebrada. Parece levemente vazia, e não há nada na tela que diga por quê.

Quando mais frames realmente para de ajudar

Três casos em que a próxima hora não compra nada:

Você está limitado pelo ruído de leitura. Com subs muito curtos, o piso de ruído é definido pela leitura, não pelos fótons do céu. Empilhar ainda ajuda, mas você está lutando contra um custo fixo que mais frames não removem. Subs mais longos removem.

Você está limitado pelo céu num alvo fraco. Sob céus claros o fundo domina, e as regiões externas fracas do seu alvo podem simplesmente nunca subir acima dele. A relação sinal/ruído ainda melhora como √N - mas √N de “mal detectável” ainda é mal detectável. Céus mais escuros ou filtros de banda estreita mudam isso; mais tempo não.

Seu alinhamento é o limite. Se o registro está deixando erros de meio pixel, adicionar frames adiciona mais dados levemente desregistrados. A pilha fica mais suave e não mais nítida. Vale conferir isso antes de culpar o céu.

Uma resposta prática

• Duração do sub: a mais longa que a sua montagem rastreia com confiança. Isso domina o ruído de leitura, e o ruído de leitura domina os alvos fracos.
• Quantidade de frames: pense em fatores de quatro. Se quer uma melhoria visível, planeje quadruplicar, não adicionar um pouco.
• Registro: precisa lidar com rotação. Inegociável em qualquer montagem altazimutal.
• Rejeição: adapte à profundidade da sua pilha e à estabilidade da sua noite.
• Conheça seu limite: ruído de leitura, brilho do céu ou alinhamento. Adicionar frames só ajuda nos dois primeiros, e devagar.

Nada disso exige entender a matemática além de √N. Mas explica por que a pilha com a qual você está insatisfeito está insatisfatória - e isso costuma ser mais útil de saber do que mais um controle deslizante.