Die eigentliche Mathematik ist ja kein Problem, es hapert an der Umsetzung der Projektionen. Bzw. Sollte man den Background bei jedem einzelnen Bild subtrahieren. Tja, und da geht das Problem schon los. Wie bringt man die Pioneer-Daten in die entsprechende Projektion/Ausschnitt/Rotation des aufgenommenen Bildes ?
Den ganzen Prozeß müßte man zusammen mit den Sky-subtraction oder statt dessen machen. Null Schimmer, wie ich das anstellen könnte.
musst du mir schon aufschreiben, ich habe keine Zeit, mir das nochmal zu überlegen, schon gar nicht wenn du es schon weißt.ganze Mathematik
meine ich je gerade: Du solltest ihn mir verständlich machen, wie ein Kochrezept.Prozess
Knickohr hat geschrieben:
ich glaube, ich (wir) verstehen es wirklich nicht. ich dachte eigentlich, das man so was ganz universell in Theli einbaut. So das es für alle nutzbar ist und angewendet werden kann, ganu gleich, welche Daten da rein geschmissen werden.
Aber OK, speziell das Knickohr-Verfahren
Pioneer-Daten hier :
http://www.stsci.edu/~kgordon/pioneer_i ... l_1_B.fits (blau)
http://www.stsci.edu/~kgordon/pioneer_i ... l_1_R.fits (rot)
Ein künstlicher Grünkanal müßte gerneriert werden (B+R/2 ???)
Ich habe das so gemacht :
# ic '%2 0. %1 1e-6 < ?' P_all_1_R_median.fits P_all_1_R_median_rotate.fits > P_all_1_R_black.fits
# swarp_theli P_all_1_R_median.fits P_all_1_R_black.fits -c resize.swarp
Siehe hier :
viewtopic.php?f=13&t=298&p=2320&hilit=pioneer#p2320
Die original Pioneer Daten liegen in der CAR-Projektion vor. Als ich diese In eine andere Projektion bringen wollte, oder "nur" das galaktische Zentrum in die Mitte brachte, gingen die eingentlichen Probleme los. Es entsteht ein Balken in der Mitte, bzw. es werden die "Ohren" bei z.B. MOL-Projektion abgeschnitten Das war dan auch der Punkt, wo ich gefrustet aufgegeben habe.
Das Ergebinis dieses Bildes (ohne den Streifen oder die Ohren) muß jetzt von meinem Bild subtrahiert werden (ebenfalls dann CAR oder MOL Projektion). Ja richtig, mein Bild ist bis einschließlich Kalibration durch die Pipeline durch. Es liegt also entweder als RGB und L vor oder aber als R, B, G, L einzeln.
Das was da dann raus kommt (verbleibender Lichtsiff & Sterne) muß kräfitg weichgezeichnet (Median-Filter ???) werden und das Ergebnis dann nochmal von meinem Bild subtrahiert werden.
Das Ganze natürlich für alle Farbkanäle. R und B liegen als Rohdaten der Pioneer-Sonden vor, der Grünkanal muß synthetisiert werden (wie schon oben angesprochen).
Da ja bei swarp die eigenltiche Auflösung der Pioneer-Daten keine Rolle spielt, müßte sich das ja mit jeder beliebigen Größe meines Bildes rechnen lassen. Oder mache ich da einen Denkfehler ?
ich dachte eigentlich, das man so was ganz universell in Theli einbaut.
ist ja . Aber trivial ist es nicht.stressen
h so. Dann müsste es ja mit einem Pioneer Bild für alle Bilder an allen Positionen gehen.
Ideal wäre es, wenn man das quasi schon als Auswahlbox bei der Skysubstraktion auswählen kann. Also quasi die Pioneer-Daten so in Theli einbindet, das sie zusätzlich zu der üblichen Skysubstraktion auch noch den Lichtsiff durch die Pioneer-Daten eliminiert
Vielleicht sollte man die ganze Sache nicht am Endbild machen, sondern in der Pipeline, also an jedem einzelnen Bild. Da fehlt mir aber der komplette Workflow
# ########################################################################
# Convert Galactic coordinates to Equatorial coordinates in the epoch J2000
# ########################################################################
def gal2eq(l, b):
"""
Convert Galatic coordinates to Equatorial Coordinates in the epch J2000.
könne wir uns sparen?gal2eq(l, b):
Knickohr hat geschrieben:Die original Pioneer Daten um jeweils 180° rotieren und dann alles, was im Original schwarz ist mit dem rotierten ersetzen. Somit füllen wir den gesamten Himmel aus, also auch das fehlende Stückchen, wo die Sonden wegen der Sonne keine Aufnahmen machen konnte. Ein leichtes Median-Filter noch drüber laufen lassen, damit die Übergänge "glatt" sind und die Störungen im Bild draußen.
#---------------------------------- Output ----------------------------------------------------
IMAGEOUT_NAME test.fits
WEIGHTOUT_NAME test.weight.fits
#------------------------------- Input Weights ------------------------------------------------
WEIGHT_TYPE NONE # Not used here
WEIGHT_SUFFIX .weight.fits
WEIGHT_IMAGE
#------------------------------- Co-addition --------------------------------------------------
COMBINE_TYPE SUM # This coaddition is for illustration
# only: the weight-map will contain a sum of field footprints
#-------------------------------- Astrometry --------------------------------------------------
CELESTIAL_TYPE EQUATORIAL # Coordinate system forced to equatorial
PROJECTION_TYPE CAR # Code for Plate Caree
CENTER_TYPE MANUAL # Imposed to alpha = delta = 0.0
CENTER 00:00:00.0, +00:00:00.0 #
PIXELSCALE_TYPE MANUAL # The full sky area will exceed the fraction that contains the fields
PIXEL_SCALE 900 # in arcsec at image center, on both axes
IMAGE_SIZE 1460,730 # plus a margin
#-------------------------------- Resampling --------------------------------------------------
RESAMPLING_TYPE BILINEAR # For illustration purposes, no need for a sophisticated interpolation
OVERSAMPLING 3 # A small oversampling only to have pretty, antialiased field limits
INTERPOLATE N
GAIN_KEYWORD GAIN
GAIN_DEFAULT 0.0
#--------------------------- Background subtraction -------------------------------------------
SUBTRACT_BACK N # No background subtraction
BACK_TYPE AUTO
BACK_DEFAULT 0.0
BACK_SIZE 128
BACK_FILTERSIZE 3
#------------------------- Virtual memory management ------------------------------------------
VMEM_DIR .
VMEM_MAX 32000
MEM_MAX 32000 # 32000 MB should be enough to avoid swapping
#------------------------------ Miscellaneous -------------------------------------------------
DELETE_TMPFILES Y # Delete temporary resampled FITS files
VERBOSE_TYPE NORMAL
Zurück zu Updates und Wunschzettel
Mitglieder in diesem Forum: 0 Mitglieder und 3 Gäste