Other RHESSI
Web Sites

VIEWGRAPHS

Outline .doc
Imager schematic .html
Tray photo .html
Grid photo ??
Grid parameters .xls
Efficiency vs E ??
Grid diagram .html + .ps
Modulation plots .ps
Modulation example .ps
Perspectives on Imaging Task .doc
Beam profiles .ps
Image Reconstruction Techniques

 .doc ***

Imaging flow .doc
Internal software features .doc
Aspect .doc
Calibration .doc
Normalization and photometry .doc
Data rate and volume .doc
Other Imaging issues .doc
Effect of grids on spectroscopy .doc
Effect of grids on light curves .doc
HESSI as an imager .doc
Factors for judging .doc
Tests for judging .doc
A Perspective on HESSI Imaging .doc

 

HESSI IMAGING – OUTLINE

 

bulletHow does HESSI imaging work ?
bulletThe image reconstruction task
bulletHESSI imaging topics
bulletEditorial comments on HESSI imaging
bulletEd Schmahl: - imaging in practice

 

 

PERSPECTIVES ON THE IMAGING

bulletResponse of an individual subcollimator
bulletAverage count rate is proportional to source intensity.
bulletModulates sources whose diameter < collimator FWHM.
bulletModulation frequency and phase depends on source location and pitch.
bulletLinear response to multiple sources.
bulletOver a small rotational interval, modulation measures one Fourier Component of source
bulletEquivalent to one baseline in a radio interferometer
bulletGrid rotation is equivalent to earth rotation synthesis
bulletCan use full set of image reconstruction tools developed for radio astronomy.
bulletOver a half-rotation, one subcollimator is a telescope with Bessel Function,

PSF ~ Jo(2pi * r / ang_PITCH) +…

bulletData set represents a set of counts measured in a large number of successive, short time bins.
bulletModulation pattern for each time bin is a map of the grid-pair transmission probability for a photon from the m’th pixel.

HESSI IMAGE RECONSTRUCTION TECHNIQUES

TECHNIQUE

BASIS

STRENGTHS

WEAKNESSES

BACKPROJECTION

Summed response of all timebins

Very robust
Fast
Inherently linear

Prominent sidelobes
Difficult photometry

CLEAN

Replace PSF of dominant sources with Gaussian

Removes sidelobes
Good heritage

Assumes point sources

MEM-Sato

Image contains minimum info consistent with data

Images look good.

More computer intensive.
Photometric uncertainties.

MEM-VIS

MEM, starting from visibility data

Images look good.
Efficient for long integrations

More computer intensive.
Photometric uncertainties.

PIXONS

MEM, with variable pixel size

Excellent image quality
Well-suited to extended sources

VERY computer intensive.

Forward Fitting

Optimizes parameters of assumed data model.

Excellent photometry and parameter determination.

Limited to relatively simple sources

 

 

TYPICAL IMAGING FLOW

User: Selects imaging parameters

Time
Energy
Subcollimators
Imaging algorithm
Imaging parameters
Etc

Hessi_image:

Locates data base
Reads event tags
Selects events based on user criteria
Bins events into short time bins
Locates or performs aspect solution
Associates (x,y,roll) aspect with each time bin
Associates live time with each time bin
Associates grid calibration with each time bin
Combines, aspect, grid calibration, etc into ‘phase’ for each time bin
Result is a ‘calibrated photon list’
Optionally, convert to visibilities
Uses selected algorithm
Displays and stores resulting image.


User: Manipulate display of resulting image

User: Iterate or repeat with other imaging parameters, etc.

 

SOME INTERNAL SOFTWARE FEATURES

bulletPolar coordinates
bulletUsed internally for core computations.
bulletExploits inherent symmetries.
bulletMajor saving in computation time.
bulletResults converted back to rectangular coordinates.
bulletRef. Ed Schmahl’s web site.

 

bulletUniversal Modulation Patterns
bulletProvides a very efficient way or representing and applying instantaneous subcollimator response.
bulletMajor savings in storage and computation time.

 

bulletGrid response is represented by a Fourier series,
bulletAnalytically, grid response is a complicated function, dependent on location and energy.
bulletFourier parameters vary slowly with energy, and offset location.
bulletParameters are average transmission, modulation amplitude and phase.
bulletFundamental, 2nd and 3rd harmonics are potentially relevant.

 

ASPECT

bulletGoal is to make aspect solution transparent to user.
bulletAspect solution will either be calculated as needed or retrieved from a database
bulletBasic aspect solution software is working, but needs more bells and whistles
bulletAspect simulation software is coming soon…
bulletEnd-to-end test through hardware (including hardware) was successful
bulletCurrently, pointing behavior is built into photon simulations and taken out with a

‘virtual’ aspect solution

 

 

IMAGER CALIBRATION

bulletCalibration types
bulletDetector calibration
bulletDead-time calibration
bulletGrid calibration
bulletAspect system calibration
bulletAttenuator calibration
bulletImager Calibration Status
bulletGrids, aspect and attenuators were fully calibrated at the subsystem level
bulletGrid alignment was independently verified by ‘gridlet test’
bulletCalibration data will be installed in SSW data bases
bulletGoal is to understand instrument response at the ~1% level
bulletInflight Calibration
bulletAssumes grid, aspect and attenuator parameters are independent of time, energy and flare.
bulletRedundancies in data to be used to refine aspect, grid and attenuator parameters.
bulletAspect component positions
bulletRoll calibration using Crab pulsar
bulletGrid parameters using flare data redundancies
bulletAttenuator parameters using discontinuities at attenuator changes.
bulletImproved calibration data base can be applied retroactively
bulletImager performance will improve with time.

 

NORMALIZATION AND PHOTOMETRY

bulletHigh spectral resolution and imaging spectroscopy are drivers for photometric accuracy.
bulletGoal is ~1% photometry in optimum circumstances.
bulletMap units will correspond to photons / pixel / cm^2 (or equivalent).

è Integrate over a source component to extract photons/cm^2

bulletBack Projection
bulletFor a point source, peak in map indicates summed incident fluence
bulletUnderestimates flux by ~10% (to be fixed)
bulletFor multiple or extended sources, this is not a good choice since ‘hand corrections’ are required.
bulletUsed at present to support testing.

 

 

OTHER IMAGING ISSUES

bulletSnapshot imaging
bulletLargely untested up to now
bulletBackprojection works at subsecond intervals.
bulletVisibilities
bulletCalculated internally
bulletNot yet made conveniently accessible
bulletFormat will be table of Time, U, V, Amplitude, Phase, (Errors, Flags, tbd)
bulletConversion to AIPS-compatible format is under consideration
bulletSpectral calibration
bulletImaging uses only the diagonal elements of the spectral response matrix.
bulletImaging Spectroscopy
bulletBased on multi-image comparison
bulletBest option will be feature based rather than pixel-based
bulletWill be a significant reduction is spectral quality
bulletNeeds good photometry
bulletManual capability now
bulletGoal by launch is semi-automated, feature based.

 

DATA RATE AND VOLUME ISSUES

bulletBasic Capabilities
bulletOn board solid-state recorder has 4 Gbyte capacity (1 photon = 4 bytes)
bulletDownlink capability is ~1 Gbyte/day/station (6 x 10 minute passes)
bulletData volume dominated by background, except for very large flares
bulletMaximum data rate ~25000 counts/s/detector segment
bulletLarge flare could easily fill memory and have substantial dead time.
bulletUse of a second ground station (Wallops) during periods of activity increases downlink capability by ~75%
bulletAttenuators to control rates
bulletTwo aluminum attenuators (thin and thick).
bulletSimultaneous insertion for all detectors.
bulletControlled by on-board software, based on data rates and adjustable parameters.
bulletOne-time override mechanism to reduce effects of malfunction.
bulletDisabled in default position for first few weeks/months
bulletDecimation to control data volume in front segment.
bulletDiscards every fixed fraction of photons below preset energy.
bulletFour steps of decimation
bulletEffects compensated in analysis (except for statistics)
bulletFast rate mode to provide some imaging capability at high rates.
bulletAt high input rates, transmitted photon rates decrease due to deadtime.
bulletFast rate mode is triggered automatically to transmit binned rates in 4 energy channels.
bulletDetector-dependent time resolution is sufficient to permit imaging in broad spectral bands.
bulletDetector rear segment is largely unaffected by high rates.

 

 

EFFECT OF GRIDS ON SPECTROSCOPY

bulletDetectors view Sun through rotating grids.
bulletTransmission of grids depends on
bulletGrid
bulletTime
bulletEnergy
bulletSource location relative to pointing axis.
bulletDistorts spectra
bulletSpectral corrections depend on approximate location of source
bulletSpectral corrections are straightforward, for spectroscopy with integration times that are multiples of half-rotation time.
bulletImportant check
bulletWhen correctly calibrated, each detector should yield identical, independent spectra.

 

EFFECT OF GRIDS ON LIGHT CURVES

Grids introduce significant artifacts into light curves on three time scales.

1.  Modulation time scales (~1 to ~500 milliseconds)
bulletTimescale depends on grid pitch and source radial offset from pointing axis
2.  Periodic at twice rotation rate
bulletDue to internal shadowing in grids
bulletTime of maximum depends on grid orientation direction from source to pointing axis
bulletAmplitude of slow modulation depends on energy, source offset, and grid pitch:thickness ratio

3.  Periodic at rotation rate

bulletDue to shadowing in grids, combined with grid tilt
bulletCorrection of slow variations is straightforward, but requires knowledge of approximate position of source
bulletCorrection for modulation is possible only for the summed response over several grids.
bullet‘Demodulation’ will be an analysis option.
bulletTimescales for modulation for various grids do not overlap
bulletResult will be a single lightcurve, (as a function of energy) valid on all timescales.

 

FACTORS TO CONSIDER FOR JUDGING WHICH FEATURES TO TRUST ?

bulletCounting statistics
bulletStatistical s/n in the peak of a BPmap is s ~ 2 * SQRT(total counts)
bulletComplexity of image (simpler is better)
bulletFraction of total counts in feature of interest (a smaller fraction is more susceptible to systematic calibration errors) bulletSpecial circumstances
bulletImportance of dead time
bulletFast-rate mode
bulletAttenuator changes
bulletSpecial source locations (relative to pointing axis or imaging field of view)
bulletSuspicious symmetries (arcs relative to pointing axis or strong image sources)
bulletCharacteristics of aspect solution
bulletKnown hardware anomalies
bulletKnown software anomalies
bulletUnusual symmetries
bulletUnusual source variability
bulletBelievability of feature

 

TESTS FOR JUDGING WHICH FEATURES TO TRUST ?

bulletReimage with:
bulletdifferent algorithms
bulletdifferent FOV
bulletpixel size
bulletimage center
bulletsubcollimator combinations
bulletenergy range
bulletenergy bin sizes
bullettime ranges
bulletcalibration parameters
bulletodd/even half-rotations

 

bulletApply redundancy tests
bulletCompare spectra/lightcurves from different detectors
bulletCheck rotational symmetry of light curve

 

bulletTrace back the feature to:
bulletFourier components
bulletObserved modulation curves

 

bulletSimulate and analyze similar sources

 

HESSI AS AN IMAGER

bulletStrengths
bulletNominal performance requirements (angular- and energy-resolution, imaging spectroscopy…)
bulletImage location
bulletColocation at different energies
bulletEnergy calibration
bulletPhotometry (usually)
bulletDynamically adaptable to a wide range of image scales
bulletDynamically adaptable to large range of source strengths
bulletLimitations
bulletLimited image complexity (Will NOT provide TRACE-like images!)
bulletLimited dynamic range within an individual image (goal ~100:1 in favorable cases)
bulletExtended sources can be invisible
bulletAll sources contribute to noise of each source feature.
bulletLimiting Factors
bulletStatistics
bulletMay dominate for weak flares, short integrations, narrow energy windows
bulletImage complexity
bulletMay be important for spatially complex cases
bulletSystematic errors – In other cases, limitation may be set by
bulletKnowledge of instrument response
bulletAlgorithm limitations

 

A PERSPECTIVE ON HESSI IMAGING

Consider HESSI as an imager which you can configure.

bullet Variables you control:
bullet Integration time
bullet Energy range
bullet Choice/weighting of subcollimators
bullet Software algorithms and parameters
bullet Factors to consider:
bullet Science objective
bullet Imaging, spectra or light curves
bullet Spatial scales
bullet Energy
bullet Characteristics of flare

Performance of imager will be sensitive to your choices.

We will all be on a steep learning curve.

We should plan on doing easy science first.

 

Link to NASA Home Page

Responsible NASA Official:  Gordon D. Holman

Web Design:  Merrick Berg, Brian Dennis, Gordon Holman, & Gilbert Prevost

Heliophysics Science Division
NASA/Goddard Space Flight Center
Laboratory for Solar Physics/ Code 671
Greenbelt, MD, 20771, USA
Gordon.D.Holman@nasa.gov

Link to NASA/Goddard Home Page

+ NASA Privacy Statement, Disclaimer, and Accessibility Certification

This site last updated November 10, 2008.