SlideShare a Scribd company logo
A User’s Guide to Lucky
Imaging
(With a focus on mosaic cameras)
Tim Staley,
Craig Mackay & the Cambridge LI Group
RS Lucky Imaging Meeting
Chicheley Hall, April 2013
WWW: timstaley.co.uk
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
OUTLINE
WHY MOSAIC?
LUCKYCAM 2009 & THE NOT
REDUCTION
PERFORMANCE
TRADEOFFS, RECOMMENDATIONS
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
OUTLINE
WHY MOSAIC?
LUCKYCAM 2009 & THE NOT
REDUCTION
PERFORMANCE
TRADEOFFS, RECOMMENDATIONS
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
WHY MOSAIC?
Spatial sampling
High spatial resolution surveys require
many pixels on sky.
At the frame rates lucky imaging requires,
it becomes difficult to run a 1Kx1K CCD
without getting swamped by readout noise.
Unfortunately, current generation of
EMCCDs cannot be abutted (though this is
changing).
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
WHY MOSAIC?
Dynamic range
Independent gain control is a bonus.
Dynamic range on a single EMCCD is
limited to about 10 magnitudes by electron
well depth in multiplication register.
(Disclaimer: back of envelope calculation.)
Can’t guide well on a saturated star.
Ok if GS is Mag ≈ 16. Potentially
problematic if GS is Mag ≈ 10.
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
WHY MOSAIC?
But...
Significant hardware challenge – alignment,
synchronisation, multiple control systems...
High data rates.
More calibration.
(See below...)
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
OUTLINE
WHY MOSAIC?
LUCKYCAM 2009 & THE NOT
REDUCTION
PERFORMANCE
TRADEOFFS, RECOMMENDATIONS
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
A HARDWARE IMPLEMENTATION
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
A HARDWARE IMPLEMENTATION
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
A HARDWARE IMPLEMENTATION
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
A HARDWARE IMPLEMENTATION
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
FIELD TEST
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
FIELD TEST
Visitor Instrument @ 2.5m Nordic
Optical Telescope
Awarded 8 nights on sky, + technical time
for setup and take-down.
Off we went to La Palma.
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
FIELD TEST
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
FIELD TEST
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
RESULTS?
Minor software
glitches slowed
progress for first
couple of days.
Overall, system
performed very well.
Some good weather.
Now, we hads lots of
data to reduce
(∼ 8TB).
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
OUTLINE
WHY MOSAIC?
LUCKYCAM 2009 & THE NOT
REDUCTION
PERFORMANCE
TRADEOFFS, RECOMMENDATIONS
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
PHOTOMETRIC CALIBRATION
Specific to EMCCDs:
Accurate on-sky gain calibration via pixel
histograms is desirable (CCDs
independently adjustable).
Histogram based bias map estimation is the
most robust way to separate bias pedestal
and internally generated signal (i.e. dark
current + clock induced charge).
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
PHOTOMETRIC CALIBRATION
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
PIPELINE OVERVIEW
Custom pipeline written in C++.
A 2 stage process, comprising evaluation
and reduction. Almost always I/O limited
(∼60–80MB/sec).
Make use of thread-friendly STL based
containers and Intel Thread Building Blocks
library to implement a ‘pipeline of filters’
pattern.
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
PIPELINE OVERVIEW
First, frame evaluation:
Load frame, debias, gain normalize,
subtract internally generated signal (DC +
CIC).
Interpolate GS region, cross-correlate with
Airy core to estimate tip-tilt and (relative)
Strehl.
Save position and quality estimates to file.
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
PIPELINE OVERVIEW
Second, selection and recombination:
Sort frames by estimated quality, selecting
according to user-criteria.
Load and calibrate, optionally applying
thresholding to a copy.
Drizzle to create final image, writing to file
whenever a user-selection criteria is met.
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
PIPELINE OVERVIEW
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
PIPELINE OVERVIEW
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
ASTROMETRIC CALIBRATION
How do you calibrate a high-resolution mosaic
detector?
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
ASTROMETRIC CALIBRATION
How do you calibrate a high-resolution mosaic
detector?
Using standard binaries would take at least
(2n-1) pointings per focal lens.
So we observe a crowded field.
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
ASTROMETRIC CALIBRATION
How do you calibrate a high-resolution mosaic
detector?
Using standard binaries would take at least
(2n-1) pointings per focal lens.
So we observe a crowded field.
But: standard catalogues (2MASS, USNO)
do not have the resolution (2MASS pixels
are 2 arcseconds width).
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
ASTROMETRIC CALIBRATION
How do you calibrate a high-resolution mosaic
detector?
Using standard binaries would take at least
(2n-1) pointings per focal lens.
So we observe a crowded field.
But: standard catalogues (2MASS, USNO)
do not have the resolution (2MASS pixels
are 2 arcseconds width).
Solution: Create own catalogues from
HST/WFPC fields, cross match using
custom code.
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
ASTROMETRIC CALIBRATION
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
ASTROMETRIC CALIBRATION
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
OUTLINE
WHY MOSAIC?
LUCKYCAM 2009 & THE NOT
REDUCTION
PERFORMANCE
TRADEOFFS, RECOMMENDATIONS
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
What resolution improvement do we
obtain?
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
PSF PROFILES
1.0 0.5 0.0 0.5 1.0
Radius in arcseconds
0.000
0.005
0.010
0.015
0.020
0.025
0.030
0.035
0.040
Strehlratioatpeak
5% selection
100% selection
Seeing profile
How does this vary as we move away from GS?
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
FWHM ACROSS THE FIELD
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
FWHM ACROSS THE FIELD
0 10 20 30 40 50 60
Distance from guide star [arcsec]
100
200
300
400
500
FWHM[mas]
10% selection
100% selection
Seeing FWHM
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
How faint can we go?
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
FAINT LIMITS: GUIDE STARS
102 103
Source flux (photo-electrons per exposure)
0.00
0.05
0.10
0.15
0.20
0.25
Strehlratio
Normalised Airy reference CC
Standard Airy reference CC
Brightest pixel
Ideal case (upper horizontal)
Seeing disc (lower horizontal)
(Simulated)
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
FAINT LIMITS: GUIDE STARS
A real world example: the Einstein Cross.
(Left: HST. Right: LI. GS Mi ≈ 17)
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
FAINT LIMITS: SCIENCE TARGETS
3C 405 (Cygnus-A)
Left: HST. Right: LI.
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
FAINT LIMITS: SCIENCE TARGETS
3C 405 (Cygnus-A)
Left: HST. Right: LI.
One hour observation, 50% selection, ∼0.5”
seeing.
Galaxy approx 30” from guide star.
Weakly detect stars of estimated Mi ≥ 23.
Can do better with latest hardware.
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
OUTLINE
WHY MOSAIC?
LUCKYCAM 2009 & THE NOT
REDUCTION
PERFORMANCE
TRADEOFFS, RECOMMENDATIONS
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
OPTIMAL PIXEL SIZE AND
FRAMERATE TRADE-OFFS
Spatial sampling vs detector noise.
(Old issue with a new twist: Drizzle many
frames with subpixel offsets.)
Longer exposures for faint targets vs
atmospheric blurring above coherence time.
Data rates.
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
OBSERVATION PLANNING
Will there be a close / bright enough guide
star?
Will de-saturating my guide star ruin my
faint limit? (Can I position the mosaic to
avoid this?)
What resolution / magnitude depth do I
need to achieve, and what weather
conditions / selection criteria does this
require?
(For very bright sources:) Would I be better
off turning off the EM amplifier altogether?
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
ASTROMETRIC CALIBRATION
This is non-trivial!
Worth investigating:
http://www.astromatic.net/
software/scamp
(Cannot read custom catalogues when I last
checked, but this is probably an easy
feature to add.)
WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
OPEN QUESTIONS / RESEARCH
PROBLEMS
Use of multiple guide stars.
Crowded field photometry with a slowly
varying PSF.
Ad

More Related Content

Similar to A user's guide to lucky imaging (6)

Energy Efficient Communications Solutions for Nomadic Nodes within a Wireless...
Energy Efficient Communications Solutions for Nomadic Nodes within a Wireless...Energy Efficient Communications Solutions for Nomadic Nodes within a Wireless...
Energy Efficient Communications Solutions for Nomadic Nodes within a Wireless...
luca-bencini
 
Multi-channel Detector Readout Integrated Circuits with ADCs for X-ray and Ga...
Multi-channel Detector Readout Integrated Circuits with ADCs for X-ray and Ga...Multi-channel Detector Readout Integrated Circuits with ADCs for X-ray and Ga...
Multi-channel Detector Readout Integrated Circuits with ADCs for X-ray and Ga...
Gunnar Maehlum
 
Realtime pothole detection system using improved CNN Models
Realtime pothole detection system using improved CNN ModelsRealtime pothole detection system using improved CNN Models
Realtime pothole detection system using improved CNN Models
nithinsai2992
 
Cmc Laboratories Client Presentation 2009
Cmc Laboratories Client Presentation 2009Cmc Laboratories Client Presentation 2009
Cmc Laboratories Client Presentation 2009
CMC Laboratories, Inc.
 
Zuxin2009
Zuxin2009Zuxin2009
Zuxin2009
zuxin
 
How to decommission a data centre
How to decommission a data centreHow to decommission a data centre
How to decommission a data centre
Databarracks
 
Energy Efficient Communications Solutions for Nomadic Nodes within a Wireless...
Energy Efficient Communications Solutions for Nomadic Nodes within a Wireless...Energy Efficient Communications Solutions for Nomadic Nodes within a Wireless...
Energy Efficient Communications Solutions for Nomadic Nodes within a Wireless...
luca-bencini
 
Multi-channel Detector Readout Integrated Circuits with ADCs for X-ray and Ga...
Multi-channel Detector Readout Integrated Circuits with ADCs for X-ray and Ga...Multi-channel Detector Readout Integrated Circuits with ADCs for X-ray and Ga...
Multi-channel Detector Readout Integrated Circuits with ADCs for X-ray and Ga...
Gunnar Maehlum
 
Realtime pothole detection system using improved CNN Models
Realtime pothole detection system using improved CNN ModelsRealtime pothole detection system using improved CNN Models
Realtime pothole detection system using improved CNN Models
nithinsai2992
 
Cmc Laboratories Client Presentation 2009
Cmc Laboratories Client Presentation 2009Cmc Laboratories Client Presentation 2009
Cmc Laboratories Client Presentation 2009
CMC Laboratories, Inc.
 
Zuxin2009
Zuxin2009Zuxin2009
Zuxin2009
zuxin
 
How to decommission a data centre
How to decommission a data centreHow to decommission a data centre
How to decommission a data centre
Databarracks
 

More from Tim Staley (8)

Which transient when? - A utility function for transient follow-up scheduling
Which transient when? - A utility function for transient follow-up schedulingWhich transient when? - A utility function for transient follow-up scheduling
Which transient when? - A utility function for transient follow-up scheduling
Tim Staley
 
How to build a TraP: An image-plane transient-discovery tool
How to build a TraP: An image-plane transient-discovery toolHow to build a TraP: An image-plane transient-discovery tool
How to build a TraP: An image-plane transient-discovery tool
Tim Staley
 
From gamma-ray to radio: Multi-wavelength follow-up in the first five minutes
From gamma-ray to radio: Multi-wavelength follow-up in the first five minutesFrom gamma-ray to radio: Multi-wavelength follow-up in the first five minutes
From gamma-ray to radio: Multi-wavelength follow-up in the first five minutes
Tim Staley
 
Tunable algorithms for transient follow-up
Tunable algorithms for transient follow-upTunable algorithms for transient follow-up
Tunable algorithms for transient follow-up
Tim Staley
 
A brief introduction to version control systems
A brief introduction to version control systemsA brief introduction to version control systems
A brief introduction to version control systems
Tim Staley
 
Training your astronomy robots to work as a team
Training your astronomy robots to work as a teamTraining your astronomy robots to work as a team
Training your astronomy robots to work as a team
Tim Staley
 
Fast radio follow-up
Fast radio follow-upFast radio follow-up
Fast radio follow-up
Tim Staley
 
Fast radio follow-up of GRBs
Fast radio follow-up of GRBsFast radio follow-up of GRBs
Fast radio follow-up of GRBs
Tim Staley
 
Which transient when? - A utility function for transient follow-up scheduling
Which transient when? - A utility function for transient follow-up schedulingWhich transient when? - A utility function for transient follow-up scheduling
Which transient when? - A utility function for transient follow-up scheduling
Tim Staley
 
How to build a TraP: An image-plane transient-discovery tool
How to build a TraP: An image-plane transient-discovery toolHow to build a TraP: An image-plane transient-discovery tool
How to build a TraP: An image-plane transient-discovery tool
Tim Staley
 
From gamma-ray to radio: Multi-wavelength follow-up in the first five minutes
From gamma-ray to radio: Multi-wavelength follow-up in the first five minutesFrom gamma-ray to radio: Multi-wavelength follow-up in the first five minutes
From gamma-ray to radio: Multi-wavelength follow-up in the first five minutes
Tim Staley
 
Tunable algorithms for transient follow-up
Tunable algorithms for transient follow-upTunable algorithms for transient follow-up
Tunable algorithms for transient follow-up
Tim Staley
 
A brief introduction to version control systems
A brief introduction to version control systemsA brief introduction to version control systems
A brief introduction to version control systems
Tim Staley
 
Training your astronomy robots to work as a team
Training your astronomy robots to work as a teamTraining your astronomy robots to work as a team
Training your astronomy robots to work as a team
Tim Staley
 
Fast radio follow-up
Fast radio follow-upFast radio follow-up
Fast radio follow-up
Tim Staley
 
Fast radio follow-up of GRBs
Fast radio follow-up of GRBsFast radio follow-up of GRBs
Fast radio follow-up of GRBs
Tim Staley
 
Ad

Recently uploaded (20)

Skin_Glands_Structure_Secretion _Control
Skin_Glands_Structure_Secretion _ControlSkin_Glands_Structure_Secretion _Control
Skin_Glands_Structure_Secretion _Control
muralinath2
 
Infrastructure for Tracking Information Flow from Social Media to U.S. TV New...
Infrastructure for Tracking Information Flow from Social Media to U.S. TV New...Infrastructure for Tracking Information Flow from Social Media to U.S. TV New...
Infrastructure for Tracking Information Flow from Social Media to U.S. TV New...
Himarsha Jayanetti
 
biochemistry amino acid from chemistry to life machinery
biochemistry amino acid from chemistry to life machinerybiochemistry amino acid from chemistry to life machinery
biochemistry amino acid from chemistry to life machinery
chaitanyaa4444
 
Preparation of Permanent mounts of Parasitic Protozoans.pptx
Preparation of Permanent mounts of Parasitic Protozoans.pptxPreparation of Permanent mounts of Parasitic Protozoans.pptx
Preparation of Permanent mounts of Parasitic Protozoans.pptx
Dr Showkat Ahmad Wani
 
whole ANATOMY OF EYE with eye ball .pptx
whole ANATOMY OF EYE with eye ball .pptxwhole ANATOMY OF EYE with eye ball .pptx
whole ANATOMY OF EYE with eye ball .pptx
simranjangra13
 
Body temperature_chemical thermogenesis_hypothermia_hypothermiaMetabolic acti...
Body temperature_chemical thermogenesis_hypothermia_hypothermiaMetabolic acti...Body temperature_chemical thermogenesis_hypothermia_hypothermiaMetabolic acti...
Body temperature_chemical thermogenesis_hypothermia_hypothermiaMetabolic acti...
muralinath2
 
Quiz 3 Basic Nutrition 1ST Yearcmcmc.pptx
Quiz 3 Basic Nutrition 1ST Yearcmcmc.pptxQuiz 3 Basic Nutrition 1ST Yearcmcmc.pptx
Quiz 3 Basic Nutrition 1ST Yearcmcmc.pptx
NutriGen
 
RAPID DIAGNOSTIC TEST (RDT) overviewppt.pptx
RAPID DIAGNOSTIC TEST (RDT)  overviewppt.pptxRAPID DIAGNOSTIC TEST (RDT)  overviewppt.pptx
RAPID DIAGNOSTIC TEST (RDT) overviewppt.pptx
nietakam
 
Antonie van Leeuwenhoek- Father of Microbiology
Antonie van Leeuwenhoek- Father of MicrobiologyAntonie van Leeuwenhoek- Father of Microbiology
Antonie van Leeuwenhoek- Father of Microbiology
Anoja Kurian
 
Chapter 4_Part 2_Infection and Immunity.ppt
Chapter 4_Part 2_Infection and Immunity.pptChapter 4_Part 2_Infection and Immunity.ppt
Chapter 4_Part 2_Infection and Immunity.ppt
JessaBalanggoyPagula
 
Polytene chromosomes. A Practical Lecture.pptx
Polytene chromosomes. A Practical Lecture.pptxPolytene chromosomes. A Practical Lecture.pptx
Polytene chromosomes. A Practical Lecture.pptx
Dr Showkat Ahmad Wani
 
SuperconductingMagneticEnergyStorage.pptx
SuperconductingMagneticEnergyStorage.pptxSuperconductingMagneticEnergyStorage.pptx
SuperconductingMagneticEnergyStorage.pptx
BurkanAlpKale
 
Turkey Diseases and Disorders Volume 2 Infectious and Nutritional Diseases, D...
Turkey Diseases and Disorders Volume 2 Infectious and Nutritional Diseases, D...Turkey Diseases and Disorders Volume 2 Infectious and Nutritional Diseases, D...
Turkey Diseases and Disorders Volume 2 Infectious and Nutritional Diseases, D...
Ali Raei
 
VERMICOMPOSTING A STEP TOWARDS SUSTAINABILITY.pptx
VERMICOMPOSTING A STEP TOWARDS SUSTAINABILITY.pptxVERMICOMPOSTING A STEP TOWARDS SUSTAINABILITY.pptx
VERMICOMPOSTING A STEP TOWARDS SUSTAINABILITY.pptx
hipachi8
 
Polymerase Chain Reaction (PCR).Poer Pint
Polymerase Chain Reaction (PCR).Poer PintPolymerase Chain Reaction (PCR).Poer Pint
Polymerase Chain Reaction (PCR).Poer Pint
Dr Showkat Ahmad Wani
 
Gel Electrophorosis, A Practical Lecture.pptx
Gel Electrophorosis, A Practical Lecture.pptxGel Electrophorosis, A Practical Lecture.pptx
Gel Electrophorosis, A Practical Lecture.pptx
Dr Showkat Ahmad Wani
 
Causes of mortalities of eggs and spawn and remedies.pptx
Causes of mortalities of eggs and spawn and remedies.pptxCauses of mortalities of eggs and spawn and remedies.pptx
Causes of mortalities of eggs and spawn and remedies.pptx
anshumanmohanty9090
 
Structure formation with primordial black holes: collisional dynamics, binari...
Structure formation with primordial black holes: collisional dynamics, binari...Structure formation with primordial black holes: collisional dynamics, binari...
Structure formation with primordial black holes: collisional dynamics, binari...
Sérgio Sacani
 
Skin function_protective_absorptive_Presentatation.pptx
Skin function_protective_absorptive_Presentatation.pptxSkin function_protective_absorptive_Presentatation.pptx
Skin function_protective_absorptive_Presentatation.pptx
muralinath2
 
Introduction to Mobile Forensics Part 1.pptx
Introduction to Mobile Forensics Part 1.pptxIntroduction to Mobile Forensics Part 1.pptx
Introduction to Mobile Forensics Part 1.pptx
Nivya George
 
Skin_Glands_Structure_Secretion _Control
Skin_Glands_Structure_Secretion _ControlSkin_Glands_Structure_Secretion _Control
Skin_Glands_Structure_Secretion _Control
muralinath2
 
Infrastructure for Tracking Information Flow from Social Media to U.S. TV New...
Infrastructure for Tracking Information Flow from Social Media to U.S. TV New...Infrastructure for Tracking Information Flow from Social Media to U.S. TV New...
Infrastructure for Tracking Information Flow from Social Media to U.S. TV New...
Himarsha Jayanetti
 
biochemistry amino acid from chemistry to life machinery
biochemistry amino acid from chemistry to life machinerybiochemistry amino acid from chemistry to life machinery
biochemistry amino acid from chemistry to life machinery
chaitanyaa4444
 
Preparation of Permanent mounts of Parasitic Protozoans.pptx
Preparation of Permanent mounts of Parasitic Protozoans.pptxPreparation of Permanent mounts of Parasitic Protozoans.pptx
Preparation of Permanent mounts of Parasitic Protozoans.pptx
Dr Showkat Ahmad Wani
 
whole ANATOMY OF EYE with eye ball .pptx
whole ANATOMY OF EYE with eye ball .pptxwhole ANATOMY OF EYE with eye ball .pptx
whole ANATOMY OF EYE with eye ball .pptx
simranjangra13
 
Body temperature_chemical thermogenesis_hypothermia_hypothermiaMetabolic acti...
Body temperature_chemical thermogenesis_hypothermia_hypothermiaMetabolic acti...Body temperature_chemical thermogenesis_hypothermia_hypothermiaMetabolic acti...
Body temperature_chemical thermogenesis_hypothermia_hypothermiaMetabolic acti...
muralinath2
 
Quiz 3 Basic Nutrition 1ST Yearcmcmc.pptx
Quiz 3 Basic Nutrition 1ST Yearcmcmc.pptxQuiz 3 Basic Nutrition 1ST Yearcmcmc.pptx
Quiz 3 Basic Nutrition 1ST Yearcmcmc.pptx
NutriGen
 
RAPID DIAGNOSTIC TEST (RDT) overviewppt.pptx
RAPID DIAGNOSTIC TEST (RDT)  overviewppt.pptxRAPID DIAGNOSTIC TEST (RDT)  overviewppt.pptx
RAPID DIAGNOSTIC TEST (RDT) overviewppt.pptx
nietakam
 
Antonie van Leeuwenhoek- Father of Microbiology
Antonie van Leeuwenhoek- Father of MicrobiologyAntonie van Leeuwenhoek- Father of Microbiology
Antonie van Leeuwenhoek- Father of Microbiology
Anoja Kurian
 
Chapter 4_Part 2_Infection and Immunity.ppt
Chapter 4_Part 2_Infection and Immunity.pptChapter 4_Part 2_Infection and Immunity.ppt
Chapter 4_Part 2_Infection and Immunity.ppt
JessaBalanggoyPagula
 
Polytene chromosomes. A Practical Lecture.pptx
Polytene chromosomes. A Practical Lecture.pptxPolytene chromosomes. A Practical Lecture.pptx
Polytene chromosomes. A Practical Lecture.pptx
Dr Showkat Ahmad Wani
 
SuperconductingMagneticEnergyStorage.pptx
SuperconductingMagneticEnergyStorage.pptxSuperconductingMagneticEnergyStorage.pptx
SuperconductingMagneticEnergyStorage.pptx
BurkanAlpKale
 
Turkey Diseases and Disorders Volume 2 Infectious and Nutritional Diseases, D...
Turkey Diseases and Disorders Volume 2 Infectious and Nutritional Diseases, D...Turkey Diseases and Disorders Volume 2 Infectious and Nutritional Diseases, D...
Turkey Diseases and Disorders Volume 2 Infectious and Nutritional Diseases, D...
Ali Raei
 
VERMICOMPOSTING A STEP TOWARDS SUSTAINABILITY.pptx
VERMICOMPOSTING A STEP TOWARDS SUSTAINABILITY.pptxVERMICOMPOSTING A STEP TOWARDS SUSTAINABILITY.pptx
VERMICOMPOSTING A STEP TOWARDS SUSTAINABILITY.pptx
hipachi8
 
Polymerase Chain Reaction (PCR).Poer Pint
Polymerase Chain Reaction (PCR).Poer PintPolymerase Chain Reaction (PCR).Poer Pint
Polymerase Chain Reaction (PCR).Poer Pint
Dr Showkat Ahmad Wani
 
Gel Electrophorosis, A Practical Lecture.pptx
Gel Electrophorosis, A Practical Lecture.pptxGel Electrophorosis, A Practical Lecture.pptx
Gel Electrophorosis, A Practical Lecture.pptx
Dr Showkat Ahmad Wani
 
Causes of mortalities of eggs and spawn and remedies.pptx
Causes of mortalities of eggs and spawn and remedies.pptxCauses of mortalities of eggs and spawn and remedies.pptx
Causes of mortalities of eggs and spawn and remedies.pptx
anshumanmohanty9090
 
Structure formation with primordial black holes: collisional dynamics, binari...
Structure formation with primordial black holes: collisional dynamics, binari...Structure formation with primordial black holes: collisional dynamics, binari...
Structure formation with primordial black holes: collisional dynamics, binari...
Sérgio Sacani
 
Skin function_protective_absorptive_Presentatation.pptx
Skin function_protective_absorptive_Presentatation.pptxSkin function_protective_absorptive_Presentatation.pptx
Skin function_protective_absorptive_Presentatation.pptx
muralinath2
 
Introduction to Mobile Forensics Part 1.pptx
Introduction to Mobile Forensics Part 1.pptxIntroduction to Mobile Forensics Part 1.pptx
Introduction to Mobile Forensics Part 1.pptx
Nivya George
 
Ad

A user's guide to lucky imaging

  • 1. A User’s Guide to Lucky Imaging (With a focus on mosaic cameras) Tim Staley, Craig Mackay & the Cambridge LI Group RS Lucky Imaging Meeting Chicheley Hall, April 2013 WWW: timstaley.co.uk
  • 2. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS OUTLINE WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
  • 3. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS OUTLINE WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
  • 4. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS WHY MOSAIC? Spatial sampling High spatial resolution surveys require many pixels on sky. At the frame rates lucky imaging requires, it becomes difficult to run a 1Kx1K CCD without getting swamped by readout noise. Unfortunately, current generation of EMCCDs cannot be abutted (though this is changing).
  • 5. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS WHY MOSAIC? Dynamic range Independent gain control is a bonus. Dynamic range on a single EMCCD is limited to about 10 magnitudes by electron well depth in multiplication register. (Disclaimer: back of envelope calculation.) Can’t guide well on a saturated star. Ok if GS is Mag ≈ 16. Potentially problematic if GS is Mag ≈ 10.
  • 6. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS WHY MOSAIC? But... Significant hardware challenge – alignment, synchronisation, multiple control systems... High data rates. More calibration. (See below...)
  • 7. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS OUTLINE WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
  • 8. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS A HARDWARE IMPLEMENTATION
  • 9. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS A HARDWARE IMPLEMENTATION
  • 10. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS A HARDWARE IMPLEMENTATION
  • 11. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS A HARDWARE IMPLEMENTATION
  • 12. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS FIELD TEST
  • 13. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS FIELD TEST Visitor Instrument @ 2.5m Nordic Optical Telescope Awarded 8 nights on sky, + technical time for setup and take-down. Off we went to La Palma.
  • 14. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS FIELD TEST
  • 15. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS FIELD TEST
  • 16. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS RESULTS? Minor software glitches slowed progress for first couple of days. Overall, system performed very well. Some good weather. Now, we hads lots of data to reduce (∼ 8TB).
  • 17. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS OUTLINE WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
  • 18. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS PHOTOMETRIC CALIBRATION Specific to EMCCDs: Accurate on-sky gain calibration via pixel histograms is desirable (CCDs independently adjustable). Histogram based bias map estimation is the most robust way to separate bias pedestal and internally generated signal (i.e. dark current + clock induced charge).
  • 19. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS PHOTOMETRIC CALIBRATION
  • 20. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS PIPELINE OVERVIEW Custom pipeline written in C++. A 2 stage process, comprising evaluation and reduction. Almost always I/O limited (∼60–80MB/sec). Make use of thread-friendly STL based containers and Intel Thread Building Blocks library to implement a ‘pipeline of filters’ pattern.
  • 21. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS PIPELINE OVERVIEW First, frame evaluation: Load frame, debias, gain normalize, subtract internally generated signal (DC + CIC). Interpolate GS region, cross-correlate with Airy core to estimate tip-tilt and (relative) Strehl. Save position and quality estimates to file.
  • 22. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS PIPELINE OVERVIEW Second, selection and recombination: Sort frames by estimated quality, selecting according to user-criteria. Load and calibrate, optionally applying thresholding to a copy. Drizzle to create final image, writing to file whenever a user-selection criteria is met.
  • 23. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS PIPELINE OVERVIEW
  • 24. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS PIPELINE OVERVIEW
  • 25. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS ASTROMETRIC CALIBRATION How do you calibrate a high-resolution mosaic detector?
  • 26. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS ASTROMETRIC CALIBRATION How do you calibrate a high-resolution mosaic detector? Using standard binaries would take at least (2n-1) pointings per focal lens. So we observe a crowded field.
  • 27. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS ASTROMETRIC CALIBRATION How do you calibrate a high-resolution mosaic detector? Using standard binaries would take at least (2n-1) pointings per focal lens. So we observe a crowded field. But: standard catalogues (2MASS, USNO) do not have the resolution (2MASS pixels are 2 arcseconds width).
  • 28. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS ASTROMETRIC CALIBRATION How do you calibrate a high-resolution mosaic detector? Using standard binaries would take at least (2n-1) pointings per focal lens. So we observe a crowded field. But: standard catalogues (2MASS, USNO) do not have the resolution (2MASS pixels are 2 arcseconds width). Solution: Create own catalogues from HST/WFPC fields, cross match using custom code.
  • 29. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS ASTROMETRIC CALIBRATION
  • 30. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS ASTROMETRIC CALIBRATION
  • 31. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS OUTLINE WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
  • 32. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS What resolution improvement do we obtain?
  • 33. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS PSF PROFILES 1.0 0.5 0.0 0.5 1.0 Radius in arcseconds 0.000 0.005 0.010 0.015 0.020 0.025 0.030 0.035 0.040 Strehlratioatpeak 5% selection 100% selection Seeing profile How does this vary as we move away from GS?
  • 34. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS FWHM ACROSS THE FIELD
  • 35. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS FWHM ACROSS THE FIELD 0 10 20 30 40 50 60 Distance from guide star [arcsec] 100 200 300 400 500 FWHM[mas] 10% selection 100% selection Seeing FWHM
  • 36. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS How faint can we go?
  • 37. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS FAINT LIMITS: GUIDE STARS 102 103 Source flux (photo-electrons per exposure) 0.00 0.05 0.10 0.15 0.20 0.25 Strehlratio Normalised Airy reference CC Standard Airy reference CC Brightest pixel Ideal case (upper horizontal) Seeing disc (lower horizontal) (Simulated)
  • 38. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS FAINT LIMITS: GUIDE STARS A real world example: the Einstein Cross. (Left: HST. Right: LI. GS Mi ≈ 17)
  • 39. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS FAINT LIMITS: SCIENCE TARGETS 3C 405 (Cygnus-A) Left: HST. Right: LI.
  • 40. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS FAINT LIMITS: SCIENCE TARGETS 3C 405 (Cygnus-A) Left: HST. Right: LI. One hour observation, 50% selection, ∼0.5” seeing. Galaxy approx 30” from guide star. Weakly detect stars of estimated Mi ≥ 23. Can do better with latest hardware.
  • 41. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS OUTLINE WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS
  • 42. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS OPTIMAL PIXEL SIZE AND FRAMERATE TRADE-OFFS Spatial sampling vs detector noise. (Old issue with a new twist: Drizzle many frames with subpixel offsets.) Longer exposures for faint targets vs atmospheric blurring above coherence time. Data rates.
  • 43. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS OBSERVATION PLANNING Will there be a close / bright enough guide star? Will de-saturating my guide star ruin my faint limit? (Can I position the mosaic to avoid this?) What resolution / magnitude depth do I need to achieve, and what weather conditions / selection criteria does this require? (For very bright sources:) Would I be better off turning off the EM amplifier altogether?
  • 44. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS ASTROMETRIC CALIBRATION This is non-trivial! Worth investigating: http://www.astromatic.net/ software/scamp (Cannot read custom catalogues when I last checked, but this is probably an easy feature to add.)
  • 45. WHY MOSAIC? LUCKYCAM 2009 & THE NOT REDUCTION PERFORMANCE TRADEOFFS, RECOMMENDATIONS OPEN QUESTIONS / RESEARCH PROBLEMS Use of multiple guide stars. Crowded field photometry with a slowly varying PSF.