Kanzelhöhe Observatory
 for solar and environmental research
logo uni graz
 
 
 Home   Presentation   Instruments   Observations   Sun & Climate   Gallery    Documents    Links  
 
 
 

KPDC - Continuum Photoheliograms

2k x 2k CCD Full Disk Continuum (Broadband) Photoheliograms at 546 nm (10 nm FWHM)

»go direct to the KPDC synoptic data archives directory...

» Instrument Operations Timeline

» Synoptic Full Disk Continuum Photoheliograms on film - PhoKa (before July 2007)

 

Data Description - Contents

» Observing Programme

» Coordinate System and Image Orientation

» Data Levels

» Data Processing Pipeline

» Data Formats and File Names

» Daily Synoptic Data Archives Directory Structure"

»Kanzelhöhe Specific FITS Keywords and FITS Sample Headers

 

more information

» KPDC - Instrument and Data Processing Procedures Description...

 

Observing Programme

A.Patrol of the solar photosphere. Obtain a continuous time series of full disk continuum images whenever weather conditions permit. Instrumental cadence standard is 1 image/min. The cadence in the long term archives may be reduced to 1 image/hour on periods with low or no visible solar activity by selecting the "best image" of that time interval.
B.Publish a representative daily synoptic image of the solar photosphere out of the time series.

Coordinate System and Image Orientation

Coordinate System
Basically the images are mapped as the solar disk would be seen in the sky with the eye, i. e. N is up, E is to the left. The correction angle Θ specifies the alignment between the x-axis of the image array and the E-W direction. For details about Θ see the » KPDC - Instrument and Data Processing Procedures Description
Coordinate Origin
The first pixel in the data array is the left bottom edge of the image.
Heliographic Coordinates
P is the position angle of the solar rotation axis (solar N), reckoned positive from (celestial) N towards E.

Data Levels

A. Raw LevelData as written by the image acquisition software. Basically all data are in FITS, specific keywords denote properties which have to be recorded at observation time and cannot be restored or recalculated at a later time. A dark current frames (DC) is recorded each day.
B. Archive LevelBasically raw data, i. e. the image itself or the pixel data remain unchanged but prepared for archiving and better usability. The processing comprises a quality check of the individual images and a cadence reduction. We avoid any rebinning, averaging or other pixel value manipulation except a) a shifting of the solar disk towards the center coordinates of the field of view (pixel array) by shifting an integer number of pixels and b) setting pixels far outside the solar disk to 0 for a better compression of the data. More keywords are inserted into the FITS headers to make the data compliant with common accepted data standards.
C. Synoptic LevelA gain table correction (flatfield) which assumes circular symmetry on large scales on the solar disk ("Burlov"-method) is applied to archive level data. This method delivers also a center to limb variation (CLV) profile and therefore a "quiet Sun map". Dividing the flatfielded image by this "quiet Sun map" yields a normalized ("flat") solar disk image where the quiet sun level corresponds to 1 in the data array. In particular for the JPEG copies the contrast is enhanced by unsharp masking and a linear rescaling of the actual intensity range (e.g. 0...120% of the int. at disk center) to full byte scale (0...255), for individual properties of the data products see the formats and the data processing pipeline section.
The daily synoptic image is processed according to this level and published in various formats.

Data Processing Pipeline

The data processing pipeline describes the work flow from the data acqusition to the final data repository in the archives and substantial processing steps.

At Observation Time

-- expose to constant "density"vary exposure time to adapt for equal average intensities in AOI in each frame
-- frame selection on pixel contrastcalculate standard deviation of intensities in AOI and select frame with max. std. dev. within a certain number of frames
--» A. Raw Level DataFITS, normally not available

After Observation Processing

-- discard corrupted imagescheck for reasonable disk radii and intensities in all quadrants
-- cadence reductionselect "best of hour" by the Optimum Window Method, or keep all during high activity periods
-- determination of actual disk center coord. and radii
-- shift disk to center of FoVonly by integer number of pixel to avoid any data averaging
-- calculate properties and parameters for achive level datae. g. angle Θ, image scale,...
-- mask outside disk = 0set pixels far outside the disk to zero for a better data compression
--» B. Archive Level DataFITS, optional JPEG available
-- application of a gain table (flatfield)"Burlov"-method, which assumes and reduces the image to a circular symmetry on large scales and delivers also a CLV profile
-- calculate properties and parameters for synoptic level data e. g. P, B0, L0 and SOHO standard FITS keys for synoptic observations
-- contrast enhancement and derotationunsharp masking, intensity range rescaling, rotated by P and Θ to have solar N up
--» C1 Synoptic Level Datadenominated as "low contrast", FITS and JPEG available.
-- normalizationdivide by CLV (quiet Sun map)
-- contrast enhancement and derotationunsharp masking, linear intensity range rescaling, rotated by P and Θ to have solar N up
--» C2 Synoptic Level Data - normalizeddenominated as "high contrast", FITS and JPEG available.

Optimum Window Method (see Giammanco, 2000, Astron. Astrophys. Suppl. Ser. 147, 1)


 

Data Formats and File Names

FileExplanation
A. Raw Level Data

wlYYYYMMDD_hhmmss.fts.gz
raw level image frame, 2 Byte integer, compressed FITS

dcYYYYMMDD_hhmmss.fts.gz
raw level dark current frame, 2 Byte integer, compressed FITS

YYYYMMDD.NAM
list of available raw level images which passed the quality checking and cadence reduction procedure of the corresponding date, ASCII text
B. Archive Level Data

YYYYMMDD_hhmmss.fts.gz
archive level image frame, 2 Byte integer, compressed FITS

YYYYMMDD_hhmmss.jpg
archive level image frame, half scale, JPEG
C. Synoptic Level Data

kanz_bband_fd_YYYYMMDD_hhmm.ar.fts.gz
archive level representation of the daily synoptic image, 2 Byte integer, compressed FITS

kanz_bband_fd_YYYYMMDD_hhmm.clv.fts.gz
center to limb variation (CLV) profile, derived by the flatfielding procedure, compressed FITS ASCII table

kanz_bband_fd_YYYYMMDD_hhmm.hi.fts.gz
normalized image (CLV contrast compensated), rescaled to by use of BSCALE=0.01 to 1 Byte integer, compressed FITS

kanz_bband_fd_YYYYMMDD_hhmm.hi.jpg
normalized image (CLV contrast compensated), unsharp masking, intensity rescaling 40...120% -> 0...255, derotated for P and angle Θ, heliograph. coordinate grid, half scale, JPEG

kanz_bband_fd_YYYYMMDD_hhmm.lo.fts.gz
gain table corrected image, 2 Byte integer, compressed FITS

kanz_bband_fd_YYYYMMDD_hhmm.lo.jpg
gain table corrected image, unsharp masking, intensity rescaling 0...120% -> 0...255, derotated for P and angle Θ, half scale, JPEG

YYYYMMDD_hhmm[ss]...date and time of the observation, daily synoptic files may appear without seconds [ss] in filename

The file names comply with the SOHO file name convention for synoptic data.

 

Daily Synoptic Data Archives Directory Structure (rel. to web-root)

»go to the KPDC synoptic data archives directory...

DirectoryExplanation
http://cesar.kso.ac.at/phokaD/...
...KPDC synoptic root
                            ./FITS/...
...FITS format subdirectory
                                 ./arch/...
...synoptic images, archive level format
                                 ./clv/...
...CLV tables
                                 ./high/...
...synoptic level, normalized image (high contrast)
                                 ./low/...
...synoptic level, low (normal) contrast
                                     ./YYYY/...
...each level has annual subdirectories
                            ./JPEG/...
...JPEG format subdirectory
                                 ./high/...
...synoptic level, high contrast
                                 ./low/...
...synoptic level, low (normal) contrast
                                     ./YYYY/...
...each level has annual subdirectories
                                 ./recent/...
...recent (last week) synoptic level, low and high contrast

YYYY... subdirectory according to the year of observation

 

Kanzelhöhe Specific FITS Keywords and FITS Sample Headers

In the list below there are explanations of Kanzelhöhe specific FITS keywords and essential hints for the interpretation of FITS standard keywords, some other keywords might be self explaining. The FITS data files comply fully with the FITS standard and make use of the proposed keywords for SOHO. Please refer to these documentations for a comprehensive explanation of standard or SOHO specific keywords.
» FITS standard
» Proposed keywords for SOHO
 

Explanations

coordinate system
pixel 0,0...left, bottom
NAXIS1
CCD x-axis, roughly parallel to E-W
NAXIS2
CCD y-axis, roughly parallel to N-S
CTYPE1, CTYPE2
description of the physical meaning of the coordinate axes
CRPIX1, CRPIX2
position of physical coordinate origin within the data array, here position of the solar disk center
(standard FITS keyword)
CRVAL1, CRVAL2
optional offsets of the origin defined by CRPIX1, CRPIX2
CDELT1, CDELT2
scale of the physical units per pixel [arcsec/px]
CENTER_X, CENTER_Y
position of solar disk center in coordinates of the data array
SOLAR_R
actual solar radius in pixels
ANGLE
derotation angle P + Θ: rotate the image by ANGLE to have the solar rotational axis parallel to NAXIS2 (solar N up)
SOLAR_P0, SOLAR_B0
actual physical ephemeris P and B0 at the time of observation
pixel values
ADU = BZERO + BSCALE * D,
ADU...intensity level at the CCD in data units (for a 10 bit Cam from 0...1023),
D...value (as a 2 byte integer) of the corresponding array element in the FITS data array
DATAMIN, DATAMAX
min and max of of the CCD pixel values in data units (ADUs)
descriptive
observing conditions and data processing
EXP_TIME
actual exposure (integration) time used for the frame either in msec or sec
EXP_MODE
exposure mode:
0...exposure time is controlled by the image acquisition application to reach a preselected level (PRE_INT) of pixel values in the area of interest (a rectangle defined by A_O_INT),
1...the actual expsoure time is twice the value which is needed to reach the PRE_INT level,
2...the actual exposure time is fixed to the value given by EXP_TIME,
3...combination of 1 and 2
PRE_INT
preselected level (mean of pixel values in the area of interest) which should be reached by adapting the integration time
A_O_INT
the vector defines the area of interest (AOI) in pixel coordinates. The AOI is used for the calculation of a mean value and a standard deviation needed of exposure time control and frame selection
TYPE-DP
data processing level:
raw...raw data level
archive...archive data level
low...synoptic data level, normal contrast
high...synoptic data level, normalized (CLV compensated)
TELESCOP
telescope used for observation, original string 'KANZELHOEHE PHOTOSPHERE TELESCOPE' was changed to 'KPT'
ORIGIN
observing programme which created the observations

 

Reading und Understanding the CLV Profile from the FITS ASCII Data Tables

The CLV profile is saved as a ACSII data table which is defined in an extension of the FITS header (see section D.). The meaning of the standard keywords is a little bit different as for image data. The 42 characters of each line form 3 columns of E14.7 format numbers written digit behind digit, TBCOLi gives the starting column of the i-th number of the line. There are 2045 lines of data representing the steps for the profile calculation from solar disk center towards the limb and beyond. The first number of each line gives the distance r from disk center in units of solar radii, the second the same distance measured in μ (μ2 = 1-r2), and the third the calculated intensity level of a quiet sun in data units (ADUs) at that position.

 

A. Raw Level Data: FITS Sample Header

SIMPLE  =                    T / file does conform to FITS standard             
BITPIX = 16 / number of bits per data pixel
NAXIS = 2 / number of data axes
NAXIS1 = 2048 / length of data axis 1
NAXIS2 = 2048 / length of data axis 2
BZERO = 32768 / offset data range to that of unsigned short
BSCALE = 1 / default scaling factor
DATE_OBS= '2007-11-12T10:30:24Z' / Date of observation
EXP_TIME= 5.886 / Exposure Time [ms]
EXP_MODE= 0 / Exp. Mode (0=auto,1=dbl,2=fix,3=both)
PRE_INT = 850 / Preselected PixInt in AOI
A_O_INT = '900,1000,1100,1009' / Rect. for PixInt [X0,Y0,X1,Y1]
OBJECT = 'SUN ' / NOAA # or SUN for full disk
OBSTYPE = 'BBAND '
TYPE-DP = 'RAW ' / Data Processing Type
TELESCOP= 'KANZELHOEHE PHOTOSPHERE TELESCOPE'
INSTRUME= 'TM4100 ' / Camera Type
ORIGIN = 'KANZELHOEHE PHOTOSPHERE PATROL'
HISTORY Raw data from aquisition software, no processing applied
COMMENT coord-origin: left bottom
COMMENT Kanzelhoehe Observatory, A-9521 Treffen
COMMENT Univ. of Graz, Austria
END

B. Archive Level Data: FITS Sample Header

SIMPLE  =                    T                                                  
BITPIX = 16
NAXIS = 2
NAXIS1 = 2048
NAXIS2 = 2048
BUNIT = 'CCD COUNTS'
BSCALE = 1.00000
BZERO = 32768.0
DATAMIN = 0
DATAMAX = 957
CTYPE1 = 'SOLAR X ' / [arcsec/px]
CTYPE2 = 'SOLAR_Y ' / [arcsec/px]
CRPIX1 = 1025.22
CRPIX2 = 1025.31
CRVAL1 = 0.00000
CRVAL2 = 0.00000
CDELT1 = 1.02376
CDELT2 = 1.02376
DATE_OBS= '2007-11-21T08:09:23Z'
CENTER_X= 1025.22 / [pix]
CENTER_Y= 1025.31 / [pix]
SOLAR_R = 948.802 / [pix]
ANGLE = 18.8725 / [deg]
SOLAR_P0= 19.6450 / [deg]
SOLAR_B0= 2.14354 / [deg]
OBJECT = 'SUN '
TYPE-DP = 'ARCHIVE ' / Data Processing Type
OBS_TYPE= 'BBAND '
TELESCOP= 'KPT '
INSTRUME= 'TM4100 '
WAVELNTH= 5460 / [ANG], FWHM=100 [ANG]
EXPTIME = 0.0143880 / [sec]
EXP_MODE= 0 / Exp. Mode (0=auto,1=dbl,2=fix,3=both)
PRE_INT = 850 / Preselected PixInt in AOI
A_O_INT = '900,1000,1100,1009' / Rect. for PixInt [X0,Y0,X1,Y1]
ORIGIN = 'KANZELHOEHE PHOTOSPHERE PATROL'
DATE = '2007-11-21' / Date of FITS Creation
COMMENT Orientation: N up, W right, first pix is left bottom
COMMENT Pixels far outside disk are set to 0
COMMENT KANZELHOEHE SOLAR OBSERVATORY
COMMENT A-9521 TREFFEN, AUSTRIA
END

C. Synoptic Level Data: FITS Sample Header

SIMPLE  =                    T                                                  
BITPIX = 16
NAXIS = 2
NAXIS1 = 2048
NAXIS2 = 2048
BUNIT = 'CCD COUNTS'
BSCALE = 1.00000
BZERO = 32768.0
DATAMIN = 0
DATAMAX = 902
CTYPE1 = 'SOLAR X ' / [arcsec/px]
CTYPE2 = 'SOLAR_Y ' / [arcsec/px]
CRPIX1 = 1025.76
CRPIX2 = 1025.26
CRVAL1 = 0.00000
CRVAL2 = 0.00000
CDELT1 = 1.02324
CDELT2 = 1.02324
DATE_OBS= '2007-11-21T07:58:16Z'
CENTER_X= 1025.76 / [pix]
CENTER_Y= 1025.26 / [pix]
SOLAR_R = 949.285 / [pix]
ANGLE = 22.8817 / [deg]
SOLAR_P0= 19.6474 / [deg]
SOLAR_B0= 2.14447 / [deg]
OBJECT = 'SUN '
TYPE-DP = 'LOW ' / Data Processing Type
OBS_TYPE= 'BBAND '
TELESCOP= 'KPT '
INSTRUME= 'TM4100 '
WAVELNTH= 5460 / [ANG], FWHM=100 [ANG]
EXPTIME = 0.0113370 / [sec]
ORIGIN = 'KANZELHOEHE PHOTOSPHERE PATROL'
DATE = '2007-11-21' / Date of FITS Creation
HISTORY Large scale asymmetries corrected by fitting
HISTORY polynomials into rings of the Solar disk
COMMENT Orientation: N up, W right, first pix is left bottom
COMMENT KANZELHOEHE SOLAR OBSERVATORY
COMMENT A-9521 TREFFEN, AUSTRIA
END

D. Synoptic Level Data: FITS Sample Extension Header for CLV Data Tables

XTENSION= 'TABLE   '                                                            
BITPIX = 8
NAXIS = 2
NAXIS1 = 42
NAXIS2 = 2045
PCOUNT = 0
GCOUNT = 1
TFIELDS = 3
TBCOL1 = 1
TFORM1 = 'E14.7 '
TBCOL2 = 15
TFORM2 = 'E14.7 '
TBCOL3 = 29
TFORM3 = 'E14.7 '
BUNIT1 = 'FRACS OF DISK'
BUNIT2 = '1-R*R '
BUNIT3 = 'CCD COUNTS'
TTYPE1 = 'RADIUS '
TTYPE2 = 'MIU '
TTYPE3 = 'CLVSA_R '
END

 


 


footer bild
 
 
UNIGRAZonline   NEWS    SITEMAP   SEARCH   GERMAN   UNI GRAZ         EDITOR