[romh]

Am tot cautat ceva librarii sau proiecte folosite de comunitatea de astronomi. Am gasit pt python proiectul asta numit astropy. E dezvoltat de diverse personaje de la diferite universitati. In mare parte sunt wrapere care depind de alte librarii ca numpy sau scipy, dar sunt si niste lucruri specifice domeniului.

Este format din mai multe submodule, pachete mentinute de echipe diferite.

- o colectie de pachete este dedicata constantelor (viteza luminii, raza pamantului samd), librarii diverse pentru conversie si formatare timp la diferite zone, unitati de masura, coordonate(ra/dec sau altaz+rotatie plus conversie). Pare comod de apelat o constanta dupa nume in loc sa o definesti de fiecare data.

- niste librarii pentru modelare sau wrappere de la librarii dependente… Mai mult regresii liniare din ce am vazut.

- pentru structuri de date avem ceva numit Nddata - un fel de sir n dimensional care poate sa tina cam orice tip de date. Mai avem ceva numit DataTable - util de folosit cand e vorba de incarcat date dintr-un fisier fits si de prelucrat/sortat/filtrat sub forma de tabel.

- exista niste librarii pt citire/scriere nu foarte diferite de libraria implicita io dar standardizeaza ceva legat de formatul fisierelor, iti permite sa citesti un fisier fits direct intr-un datatable fara sa iti bati capul cu parsarea care imi imaginez ca este anevoioasa.

- o librarie cu statistici, niste intervale de incredere, distributii poisson si alte lucruri pe care nu le inteleg momentan. Exista o librarie de convolutii cu transformate fourier, filtre gausiene si alte lucruri. Exista o librarie care automatizeaza downloadul de date de la nu stiu ce website IERS cu date despre rotatia pamantului…

Si cam asta este un sumar foarte succint. Nu cred ca o sa folosesc prea curand dar mi s-a parut interesant ca multe din librariile mentionate sunt o abstractizare in jurul acestui format numit fits pe care nu l-am prea intalnit pana acum.

Referitor la prelucrarea imaginilor nu am vazut nimic special - imaginile inca sunt tratate ca simple matrici pe care le aduni daca vrei sa le stivuiesti. Dar sunt niste librarii care te ajuta de exemplu sa adaugi un grid cu coordonatele galactice peste poza, reiesind doar din metadatele fits.

Mai mult un proiect care sa standardizeze lucrul astronomilor profesionisti care prelucreaza date de la jwst sau alte observatoare....

____________
https://www.astropy.org/

[romh]

Colectia asta de librarii nu e chiar asa de inutila. O sa incerc sa pun niste exemple cu tot felul de experimente pe care le fac si care sper ca pot fi utile pentru alti amatori. Iar notebookurile sunt simple tutoriale gasite online si adaptate un pic la nevoile unui amator.

Fotometrie

Este procesul de masurare a intensitatii si luminozitatii printre altele si pentru calculul magnitudinii aparente
Vom folosi libraria Photutils din colectia Astropy pentru a identifica stelele dintr-o poza jpeg si apoi vom plota si coordonatele acestei stele.

Utilitate:

- Extragerea de informatii din pozele astronomilor amatori
- Detectarea de duble
- Calculul magnitudinii aparente

Principala provocare o reprezinta faptul ca nu lucrez .. inca cu poze fips. E un format pe care trebuie sa il inteleg mai bine. Dar am vrut sa iau o poza jpeg facuta cu telefonul acum ceva timp si sa vad daca reusesc sa extrag stelele. As putea sa fac ceva simplu dar libraria de mai sus deja foloseste niste modele citate de tipi de la nasa... imi imaginez ca e mai buna decat orice pot sa inventez eu.

Notebookul Jupyter se gaseste aici, puteti sa rulati toate celulele si sa le analizati dupa
https://colab.research.google.com/drive/1a...sp=sharing

Modelul de detectie se numeste daostarfinder. pare sa fie cel putin 3 alte modele.

[nobody]

Probabil imaginea trebuie curatata si calibrata.
La ce limitari CMOS te referi ?

[romh]

Salut

De cred ca e o problema cu magnitudinile pt cmos - cand am citit recent cate ceva despre diferentele ccd/cmos am vazut ca cmos are cate un adc nu foarte performant pe fiecare linie orizontala de celule fotosensitive. Diferentele de calitate intre atatea adc-uri pot rezulta in luminozitati diferite de la o linie la alta.

Dar... poate m-am grabit cu concluzia. Uite aici e poza originala

https://drive.google.com/uc?id=17wZCqWSwADyTLr9UuKDq7O62pZ5zmRtx

Ar trebui sa fie o poza cu Arcturus pe care am luat-o cu tlefonul la primele observatii - e steaua mare din mijloc. Am actualizat notebookul ca sa imi ploteze stelele cu cele mai mari 5 magnitudini sau ma rog... ce este aperture_sum ... trebuie sa mai citesc.

Cea mai stralucitoare este intr-adevar Arcturus, apoi imi ia niste stele din coltul dreapta... Si aici cred ca ai dreptate - aduna poluarea luminoasa la acele stele, sau cel putin asa mi se pare... Am nevoie sa filtrez poluarea luminoasa din soft.

Pana acum singurul filtru necesar pt poze este alb-negru. Poza este salvata in alb-negru apoi culorile sunt inversate apoi este extrasa media din matricea originala, apoi este aplicat acel algoritm dao. Deci exista deja un pic de filtrare si pregatire, mai mult nu este necesar din ce am citit sumar in documentatie.

Sau poate fi vorba de parametrii pe care ii dau acestui model dao - fwhm si threshold... Mai experimentez...

[nobody]

Poza originala din link are doar 54 KB ?

[romh]

Da, are tocmai 54 kb ...

Ca sa scap de poluarea luminoasa am aplicat un filtru de blur gausian apoi am scazut matricea rezultata din matricea originala.

matricea originala
https://drive.google.com/file/d/1e-o9nQ99F...share_link

matricea blurata
https://drive.google.com/file/d/1YdV3stUeZ...share_link

matricea oringala - matrice blurata
https://drive.google.com/file/d/15ipQHW_rJ...share_link

Atasez stelele identificate mai jos. Cu cat mai mare este dimensiunea matricei de blur cu atat mai putine stele vad... Si da, poluarea a disparut, dar impreuna cu ea si stelele care au fost mascate de poluare...

Stelele astea 5... tinand cont ca imaginea e cu susul in jos.
Arcturus e cea din mijloc, apoi cu el in linie dreapta Epsilon Bootis si Delta Bootis. Apoi sus se vede cred ca Alphecca din Corona Borealis. Iar langa Arcturus in dreapta nu stiu ce stea e aia. Nu corespunde la nimic din ce vad pe stelarium.

[Erwin]

Într-adevăr, astropy este destinată în primul rând uzului profesional. Există o mulțime de aplicații ale formatului FIT(S) în știință în general, pot fi stocate date numerice de la diverse instrumente, nu numai imagini de la camere. În HEADER sunt date despre: autor, instituție, instrumentul cu care au fost obținute, timpul, diferite caracteristici etc. Headerul poate fi extins după nevoi și există unele care nu respectă întru-totul standardele. În zona a doua numită DATA este un array (matrice) uni sau multi-dimensională, imaginile au de obicei 2 dimensiuni dar se pot folosi și 3D sau chiar mai mult. Am folosit pyFITS pentru manevrarea imaginilor FITS cu care am lucrat. Mai sunt și alte colecții de biblioteci utile, mai folosim de exemplu sunpy pentru Soare. Multe dintre ele se folosesc pentru prelucrarea imaginilor de la observatoarele aflate pe sateliți așa că au o destinație precisă.

Avantajul python este că e gratuit și înlocuiește cu succes IDL-ul care costă o căruță de bani pe an o singură licență. Treptat s-au construit biblioteci de funcții care emulează funcțiile disponibile în IDL și care sunt accesibile oricui.