Forum der Vereinigung der Sternfreunde

Hallo -
ich versuche derzeit mein Spektrum mit "Verlgleichs" Spektren aus dem ESO Sience Archive Portal zu kombinieren.
Zu diesem Zwecke versuche ich ein FITS-file aus dem Portal nach MIDAS zu laden...

mit INTAPE/FITS das File erstmal in interne Tabelle konvertiert.

INTAPE/FITS 1 eso ESO.fits
$mv eso0001.tbl eso.tbl

SHOW/TABLE liefert

SHOW/TABLE eso
Table : eso [Transposed format]
No.Columns : 6 No.Rows : 1
All.Columns: 909143 All.Rows: 8 Sel.Rows: 1
Sorted by :Sequence Reference:Sequence
Col.# 1:WAVE Unit:Angstrom Format:E15.5 R*8(129877)
Col.# 2:FLUX_REDUCED Unit:adu Format:E15.5 R*4(129877)
Col.# 3:ERR_REDUCED Unit:adu Format:E15.5 R*4(129877)
Col.# 4:BGFLUX_REDUCED Unit:adu Format:E15.5 R*4(129877)
Col.# 5:FLUX Unit:10^-16 erg/cm^2/ Format:E15.5 R*4(129877)
Col.# 6:ERR Unit:10^-16 erg/cm^2/ Format:E15.5 R*4(129877)


verstehen tu ich die Anzahl der Columns schon mal nicht...

READ/TAB eso
Table : eso

Sequence WAVE (array col)
-------- --------------------
(ERR) Problems in executing /usr/lib/eso-midas/17FEB/prim/exec/tdatatbl.exe


rumms....

war da schon mal jemand erfolgreicher ?

Liebe Grüße
Georg.

Kannst Du bitte einen Link zu der Fits-Datei posten.

Und:

echo/on

schaltet mehr Meldungen über die Arbeit eines Befehls ein.

echo/full

überschwemmt dich geradezu mit noch mehr Mitteilungen.

echo/off

schaltet den ganzen Wust wieder ab.

günter

_________________
Im längsten Frieden spricht der Mensch nicht so viel Unsinn und Unwahrheit als im kürzesten Kriege. (Jean Paul)

Hallo - ich habe mal versucht die Datei hier anzubinden.
Vielleicht kannst Du ja mal Dein Glück versuchen.

CS
Georg

Eine Fitsdatei besteht aus mindestens einer HDU, Header Data Unit. Mindestens! Es können auch mehrere sein, hier sind es zwei HDUs. Eine leere Image-HDU und eine Tabelle.

INDISK/MFITS Eso.fits

erzeugt die beiden Dateien toto0000.bdf und toto0001.tbl
SHOW/TABLE, CREATE/GRAPH

SET/GRAPH stype=0
PLOT/TABLE toto0001.tbl :wave :flux

viola

In einer Multiextension-Fits-Datei kann man die Originaldaten und die reduzierten speichern und auch Tabellen. Da lohnt es sich astropy zu lernen. Ceterum censeo, nix Carthago aber Python!

günter

_________________
Im längsten Frieden spricht der Mensch nicht so viel Unsinn und Unwahrheit als im kürzesten Kriege. (Jean Paul)

Hallo,
ich habe mal versucht, die Datei mit Iraf zu öffnen, was mir aber nur bzgl. des Headers gelang. Denb kann man allerdings auch mit dem Linux Textprogramm lesen. Danach enthält die Datei 15 Ordnungen. Vermutlich handelt es sich um ein ESO-typisches Multispecformat, was aber leider nicht durch Iraf geöffnet werden kann.
Viele Grüße
Christian

Okay Freunde - ich versinke in Demut.... (gebeut sich rückwärts zurückziehend . LOL)

das halt schon mal in der Tat geklappt. - Irgendwie musste das ja gehen. Die ESO legt bestimmt nicht Ihre Daten ab, so dass man es MIDAS nicht mehr lesen kann...
Ich dachte, ich wüsste einiges von MIDAS - falsch gedacht.... - naja, ist ja auch nur ca. 30 Jahre her...

Auf jeden Fall kann ich jetzt mal meine Spektren gegen dieses hier stellen und damit weiterbosseln... (more to come....) Vielen Dank.

Georg.

Hä,
Python delendam esse? Vere?
Christian

Wir haben noch ein bisschen mit der Datei ESO.fits herumgespielt, weil der Inhalt der eingebetteten Tabelle nicht so einfach zu lesen war.

Es stellte sich dann heraus, dass das eine dreidimensionale Tabelle war. Es gab die Spalten Wave, Flux etc. aber diese enthielten nicht einfach die Werte in aufeinander folgenden Zeilen, sondern jede Spalte hatte nur eine einzige Zeile. Die gesuchten Werte waren darin jeweils in einem Array, Vektor gespeichert. Dann waren die Werte aber auch mit den Mitteln von Midas in ein vernünftiges Format zu bringen.

Warum einfach, wenn es auch umständlich geht.

Mit AstroPy konnte man die Tabellen auch relativ leicht lesen. Vielleicht wurde die Tabelle auch so erzeugt.

Ja, Christian, ich bin im übrigen der Meinung, dass Python momentan das Mittel der Wahl ist. Ich sprach nicht von delendam esse, nur davon, dass Carthago nicht mehr so interessant ist wie Python.

günter

_________________
Im längsten Frieden spricht der Mensch nicht so viel Unsinn und Unwahrheit als im kürzesten Kriege. (Jean Paul)

Über diese hochkomplexe und aktuelle Thema solltest Du Dich vielleicht mit Marcus Porcius Cato Censorius, auch Cato der Ältere genannt, austauschen.
Viele Grüße
Christian

Ich habe mal ein kleines Python-Programm geschrieben, unter Verwendung von Astropy-Modulen, mit dem man die Datei ESO.fits lesen, ihre Struktur ergründen und auch plotten kann. Siehe anliegendes Bildschirmfoto.
Links das Programmskript, rechts im Textfeld die Ausdrucke, die das Programm bei der Ausführung wiedergibt und der Spektrumplot.
Die Tabelle enthält 6 Spalten WAVE bis ERR. Jede Spalte hat eine Zeile, die jeweils aus einer Liste von 129877 Werten inkl. einer Einheit besteht.

_________________
Herzliche Grüße / best regards

Lothar

https://lotharschanne.wordpress.com/