projects:micado:astrom [OmegaCEN]

This page is read only. You can view the source, but not change it. Ask your administrator if you think this is wrong.
{{NEWPAGE}}
/********************** DON'T EDIT ABOVE THIS LINE **********************/

= Astrometry with MICADO

== Action items

* fully project statistical model into analysis: tip-tilt jitter deliveres indepedent measurements every ~3sec, AO keeps PSF in shape at KHz timescales. So exptimes<3sec do not help you to beat tilt-jitter. This needs to be confirmed with 1sec spacing simulations. If confirmed: can convert N exposures, to beat down tilt-jitter with sqrt(N), into brightness / exposure times (use phase A numbers) 
* Need distribution of stellar colors in SIMs, otherwise no effect of Atmospheric Dispersion on relative positions between stars, only broadening of their PSF.
* Centroiding accuracy constraint could come from 1sec spaced exposures as they are sampling more similar centroid. 


== Authors 

* Remko Stuik (**RS**), Gijs Verdoes Kleijn (**GVK**), Ramon Navarro (**RN**) and Eline Tolstoy (**ET**), John McFarland (**JMF**)

== Executive summary 

We have performed an astrometric analsyis of simulated MICADO images to assess the astrometric error budget 

* Set 1: 2.2 micron mono-chromatic images (so in K-band) to test relative astrometry without atmospheric dispersion effects.

Conclusions are:


== Simulated data

We have constructed simulated images by combining AO PSFs, provided by LESIA, with our ADC modelling. 

The LESIA simulation was made in-order, i.e. it was run at 1 KHz and every 1000-th frame
was saved. This simulation was run for 1000 seconds, leading to 1000 frames. The
wind blows the turbulence over the telescope in 2.5 seconds, so every 2.5 seconds
you have a totally different piece of atmosphere, while within 2.5 seconds, you
still have parts of similar turbulence on the telescope and parts of the turbulence
will be correlated, leading to similar features in the PSF. So in principle image 0,
1 and 2 should have common components, while image 0 and 3 should have totally
different features, etc.


^Simulation spec^Description^
| M1 diameter, windspeed (temp,humid perhaps handy as well) | 38.5m, f/17.48 |
| epoch | single epoch with total duration over 1 hour|
| wavelength | monochromatic at 2.2micron |
| PSF | FWHM at 2.2 microns: 1.22*2.2e-6/39 = 14 mas|
| R_0 | 0.129m, typical scale air bubbles|
| zenith angle | 0 deg |
| field angle | 0deg? angle between AO star and centre of FoV|
| Guidestar | mag<5 , 1 NGS on axis|
| noise | added bias level of 1000 with sigma=10ADU |
| Instrumental geom. distorsion | Not taken into account: assumed zero geometric distorsion |
|Telescope instabilities | Not taken into account |
|Achromatic differential refraction |applied, Trippe expects to be ~15mas over FoV in linear term and ~2micro-arcsec in second order term.|
|Chromatic differential refraction \\ aka Atmospheric Dispersion | applied, ADC study #1 computes ~10-30mas  dispersion in absence of ADC (Trippe expects induces ~1mas astrometric error in absence of ADC.) |
|AO Guidestar Measurement errors | Applied, but are negligible given brightness of AO star. NB, Trippe expects order N-1 spatial distorsions for N GSs. So N-1 order coordinate transformation needed on each individual exposure. |
|Atmos. Differential tilt jitter | yes|
|Anisoplanatism SCAO | yes|
| Galaxies as calibrators | N/A|
| Sky projected pixelscale | no |
|Detector| Central chip Hawaii-4RG: 4k x 4k, 15 micron/pixels. \\ Pixel Scale 3 mas/pix (1.5 mas zoom?) \\ --> Half chip = 2048 pix -> maximum 6.144"* SQRT(2) = 8.6" (corner)|
|pixelsampling | **RS:** give caveats in modelling this effect.| 
| gain | 1 ADU / e-  **TBC**|
|Stars | **RS**: describe stellar spectra used|
|Star positions| Two stars at 5arcsec and 10arcsec from AO star.|
|Exposure times| 1 KHz snapshots, separated by 3sec ("sampling A" ) and 1sec ("sampling B")|
|Stacking procedure| Straight sum |
|Shack Hartmann | 1KHz of 77x77 subapertures each 4x4pixels in size|
|windspeed | 14.5 m/s|

**RS**: add description of ADCs and their modelling:
   * baseline ADC (counterrotating, pupil plane)
   * linear ADC

=== Data sets

We have 1 set of simulated MICADO images that take into account:  
* set 1: atmosphere+telecope+SCAO 

For each set we have 100 snapshots of the wavefront that are converted to 100 images. Each image is formed from a 1 KHz (i.e., 0.001s)  snapshot and the time interval between snapshots is 3s (i.e., of order the wind crossing time of the E-ELT). Additionally we have constructed series of 50 coadds of 2 time-adjacent images ("50 x 2-coadds"), and in similar fashion 25 x 4-coadds, 20 x 5-coadds, etcetera.  


^Set^Simulated images^ Sourcelists ^ Summary catalog ^
|no ADC, 3sec sampling| [[http://dbview.astro-wise.org/DbView?project=MICADO%40E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=500&fcontext_project_only=on&show_expanded=no&Exportselect=NoDownload&class_str=ReducedScienceFrame&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AsS%2527expanded%2527%250Ap3%250A%2528dp4%250AS%2527ReducedScienceFrame%2527%250Ap5%250A%2528lp6%250AS%2527open%2527%250Ap7%250AaS%2527root%2527%250Ap8%250Aag5%250AasS%2527ReducedScienceFrame.chip%2527%250Ap9%250A%2528lp10%250AS%2527prefilled%2527%250Ap11%250AaS%2527link_property%2527%250Ap12%250AaS%2527Chip%2527%250Ap13%250AasS%2527ReducedScienceFrame.instrument%2527%250Ap14%250A%2528lp15%250Ag11%250AaS%2527link_property%2527%250Ap16%250AaS%2527Instrument%2527%250Ap17%250AasS%2527ReducedScienceFrame.filter%2527%250Ap18%250A%2528lp19%250Ag11%250AaS%2527link_property%2527%250Ap20%250AaS%2527Filter%2527%250Ap21%250AassS%2527derived%2527%250Ap22%250A%2528dp23%250As.&expanded_leaf=&expanded_leaf_type=&expanded_leaf_class=&derived_leaf=&node_state=&ReducedScienceFrame.expanded=open&ReducedScienceFrame.chip.expanded=prefilled&ReducedScienceFrame.instrument.expanded=prefilled&ReducedScienceFrame.filter.expanded=prefilled&ReducedScienceFrame.creation_date.op=%3C&ReducedScienceFrame.creation_date=2014-02-18&ReducedScienceFrame.EXPTIME.op=%3D&ReducedScienceFrame.EXPTIME=0.001&ReducedScienceFrame.is_valid.op=%3E&ReducedScienceFrame.is_valid=0&ReducedScienceFrame.filter.name=MICADO_SIM_MONO&ReducedScienceFrame.filter.name.op=%3D&ReducedScienceFrame.instrument.name=MICADO&ReducedScienceFrame.instrument.name.op=%3D | 100 snapshots]] \\ [[http://dbview.astro-wise.org/DbView?project=MICADO%40E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=500&fcontext_project_only=on&show_expanded=no&Exportselect=NoDownload&class_str=ReducedScienceFrame&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AsS%2527expanded%2527%250Ap3%250A%2528dp4%250AS%2527ReducedScienceFrame%2527%250Ap5%250A%2528lp6%250AS%2527open%2527%250Ap7%250AaS%2527root%2527%250Ap8%250Aag5%250AasS%2527ReducedScienceFrame.chip%2527%250Ap9%250A%2528lp10%250AS%2527prefilled%2527%250Ap11%250AaS%2527link_property%2527%250Ap12%250AaS%2527Chip%2527%250Ap13%250AasS%2527ReducedScienceFrame.instrument%2527%250Ap14%250A%2528lp15%250Ag11%250AaS%2527link_property%2527%250Ap16%250AaS%2527Instrument%2527%250Ap17%250AasS%2527ReducedScienceFrame.filter%2527%250Ap18%250A%2528lp19%250Ag11%250AaS%2527link_property%2527%250Ap20%250AaS%2527Filter%2527%250Ap21%250AassS%2527derived%2527%250Ap22%250A%2528dp23%250As.&expanded_leaf=&expanded_leaf_type=&expanded_leaf_class=&derived_leaf=&node_state=&ReducedScienceFrame.expanded=open&ReducedScienceFrame.chip.expanded=prefilled&ReducedScienceFrame.instrument.expanded=prefilled&ReducedScienceFrame.filter.expanded=prefilled&ReducedScienceFrame.creation_date.op=%3C&ReducedScienceFrame.creation_date=2014-02-19&ReducedScienceFrame.EXPTIME.op=%3D&ReducedScienceFrame.EXPTIME=0.002&ReducedScienceFrame.is_valid.op=%3E&ReducedScienceFrame.is_valid=0&ReducedScienceFrame.filter.name=MICADO_SIM_MONO&ReducedScienceFrame.filter.name.op=%3D&ReducedScienceFrame.instrument.name=MICADO&ReducedScienceFrame.instrument.name.op=%3D | 3 x50 2-coadds]] \\ 3x25 4-coadds \\ 3x20 5-coadds \\  3x10 10-coadds \\ 3x5 20-coadds \\ [[http://dbview.astro-wise.org/DbView?project=MICADO%40E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=500&fcontext_project_only=on&show_expanded=no&Exportselect=NoDownload&class_str=ReducedScienceFrame&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AsS%2527expanded%2527%250Ap3%250A%2528dp4%250AS%2527ReducedScienceFrame%2527%250Ap5%250A%2528lp6%250AS%2527open%2527%250Ap7%250AaS%2527root%2527%250Ap8%250Aag5%250AasS%2527ReducedScienceFrame.chip%2527%250Ap9%250A%2528lp10%250AS%2527prefilled%2527%250Ap11%250AaS%2527link_property%2527%250Ap12%250AaS%2527Chip%2527%250Ap13%250AasS%2527ReducedScienceFrame.instrument%2527%250Ap14%250A%2528lp15%250Ag11%250AaS%2527link_property%2527%250Ap16%250AaS%2527Instrument%2527%250Ap17%250AasS%2527ReducedScienceFrame.filter%2527%250Ap18%250A%2528lp19%250Ag11%250AaS%2527link_property%2527%250Ap20%250AaS%2527Filter%2527%250Ap21%250AassS%2527derived%2527%250Ap22%250A%2528dp23%250As.&expanded_leaf=&expanded_leaf_type=&expanded_leaf_class=&derived_leaf=&node_state=&ReducedScienceFrame.expanded=open&ReducedScienceFrame.chip.expanded=prefilled&ReducedScienceFrame.instrument.expanded=prefilled&ReducedScienceFrame.filter.expanded=prefilled&ReducedScienceFrame.creation_date.op=%3C&ReducedScienceFrame.creation_date=2014-02-19&ReducedScienceFrame.EXPTIME.op=%3D&ReducedScienceFrame.EXPTIME=0.025&ReducedScienceFrame.is_valid.op=%3E&ReducedScienceFrame.is_valid=0&ReducedScienceFrame.filter.name=MICADO_SIM_MONO&ReducedScienceFrame.filter.name.op=%3D&ReducedScienceFrame.instrument.name=MICADO&ReducedScienceFrame.instrument.name.op=%3D | 3x4 25-coadds]] \\ 3x2 50-coadds \\ 3x1 100-coadds  | 100 snapshots \\ 50 2-coadds \\ 25 4-coadds \\ 20 5-coadds \\ [[ http://dbview.astro-wise.org/DbView?project=MICADO%40E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=100&fcontext_project_only=on&show_expanded=no&Exportselect=NoDownload&class_str=SourceList&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AS%2527SourceList.instrument.name%2527%250Ap3%250A%2528S%2527%253D%2527%250Ap4%250AS%2527MICADO%2527%250Ap5%250Atp6%250AssS%2527expanded%2527%250Ap7%250A%2528dp8%250AS%2527SourceList.frame.instrument%2527%250Ap9%250A%2528lp10%250AS%2527prefilled%2527%250Ap11%250AaS%2527link_property%2527%250Ap12%250AaS%2527Instrument%2527%250Ap13%250AasS%2527SourceList.frame.filter%2527%250Ap14%250A%2528lp15%250Ag11%250AaS%2527link_property%2527%250Ap16%250AaS%2527Filter%2527%250Ap17%250AasS%2527SourceList.chip%2527%250Ap18%250A%2528lp19%250AS%2527prefilled%2527%250Ap20%250AaS%2527link_property%2527%250Ap21%250AaS%2527Chip%2527%250Ap22%250AasS%2527SourceList.filter%2527%250Ap23%250A%2528lp24%250AS%2527prefilled%2527%250Ap25%250AaS%2527link_property%2527%250Ap26%250AaS%2527Filter%2527%250Ap27%250AasS%2527SourceList.frame%2527%250Ap28%250A%2528lp29%250AS%2527open%2527%250Ap30%250AaS%2527link_property%2527%250Ap31%250AaS%2527ReducedScienceFrame%2527%250Ap32%250AasS%2527SourceList.instrument%2527%250Ap33%250A%2528lp34%250AS%2527prefilled%2527%250Ap35%250AaS%2527link_property%2527%250Ap36%250AaS%2527Instrument%2527%250Ap37%250AasS%2527SourceList%2527%250Ap38%250A%2528lp39%250Ag30%250AaS%2527root%2527%250Ap40%250Aag38%250AasS%2527SourceList.frame.chip%2527%250Ap41%250A%2528lp42%250Ag11%250AaS%2527link_property%2527%250Ap43%250AaS%2527Chip%2527%250Ap44%250AassS%2527derived%2527%250Ap45%250A%2528dp46%250AS%2527SourceList.frame%2527%250Ap47%250Ag32%250Ass.&expanded_leaf=&expanded_leaf_type=&expanded_leaf_class=&derived_leaf=&node_state=&SourceList.frame.instrument.expanded=prefilled&SourceList.frame.filter.expanded=prefilled&SourceList.chip.expanded=prefilled&SourceList.filter.expanded=prefilled&SourceList.frame.expanded=open&SourceList.instrument.expanded=prefilled&SourceList.expanded=open&SourceList.frame.chip.expanded=prefilled&SourceList.is_valid.op=%3E&SourceList.is_valid=0&SourceList.detection_frame.base=BaseFrame&SourceList.frame.base=ReducedScienceFrame&SourceList.frame.creation_date.op=%3C&SourceList.frame.creation_date=2014-02-19&SourceList.frame.EXPTIME.op=%3D&SourceList.frame.EXPTIME=0.020&SourceList.frame.is_valid.op=%3E&SourceList.frame.is_valid=0&SourceList.frame.instrument.name=MICADO&SourceList.frame.instrument.name.op=%3D&SourceList.instrument.name=MICADO&SourceList.instrument.name.op=%3D | 3x5 20-coadds]] \\ [[http://dbview.astro-wise.org/DbView?project=MICADO%40E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=100&fcontext_project_only=on&show_expanded=no&Exportselect=NoDownload&class_str=SourceList&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AS%2527SourceList.instrument.name%2527%250Ap3%250A%2528S%2527%253D%2527%250Ap4%250AS%2527MICADO%2527%250Ap5%250Atp6%250AssS%2527expanded%2527%250Ap7%250A%2528dp8%250AS%2527SourceList.frame.instrument%2527%250Ap9%250A%2528lp10%250AS%2527prefilled%2527%250Ap11%250AaS%2527link_property%2527%250Ap12%250AaS%2527Instrument%2527%250Ap13%250AasS%2527SourceList.frame.filter%2527%250Ap14%250A%2528lp15%250Ag11%250AaS%2527link_property%2527%250Ap16%250AaS%2527Filter%2527%250Ap17%250AasS%2527SourceList.chip%2527%250Ap18%250A%2528lp19%250AS%2527prefilled%2527%250Ap20%250AaS%2527link_property%2527%250Ap21%250AaS%2527Chip%2527%250Ap22%250AasS%2527SourceList.filter%2527%250Ap23%250A%2528lp24%250AS%2527prefilled%2527%250Ap25%250AaS%2527link_property%2527%250Ap26%250AaS%2527Filter%2527%250Ap27%250AasS%2527SourceList.frame%2527%250Ap28%250A%2528lp29%250AS%2527open%2527%250Ap30%250AaS%2527link_property%2527%250Ap31%250AaS%2527ReducedScienceFrame%2527%250Ap32%250AasS%2527SourceList.instrument%2527%250Ap33%250A%2528lp34%250AS%2527prefilled%2527%250Ap35%250AaS%2527link_property%2527%250Ap36%250AaS%2527Instrument%2527%250Ap37%250AasS%2527SourceList%2527%250Ap38%250A%2528lp39%250Ag30%250AaS%2527root%2527%250Ap40%250Aag38%250AasS%2527SourceList.frame.chip%2527%250Ap41%250A%2528lp42%250Ag11%250AaS%2527link_property%2527%250Ap43%250AaS%2527Chip%2527%250Ap44%250AassS%2527derived%2527%250Ap45%250A%2528dp46%250AS%2527SourceList.frame%2527%250Ap47%250Ag32%250Ass.&expanded_leaf=&expanded_leaf_type=&expanded_leaf_class=&derived_leaf=&node_state=&SourceList.frame.instrument.expanded=prefilled&SourceList.frame.filter.expanded=prefilled&SourceList.chip.expanded=prefilled&SourceList.filter.expanded=prefilled&SourceList.frame.expanded=open&SourceList.instrument.expanded=prefilled&SourceList.expanded=open&SourceList.frame.chip.expanded=prefilled&SourceList.is_valid.op=%3E&SourceList.is_valid=0&SourceList.detection_frame.base=BaseFrame&SourceList.frame.base=ReducedScienceFrame&SourceList.frame.creation_date.op=%3C&SourceList.frame.creation_date=2014-02-19&SourceList.frame.EXPTIME.op=%3D&SourceList.frame.EXPTIME=0.025&SourceList.frame.is_valid.op=%3E&SourceList.frame.is_valid=0&SourceList.frame.instrument.name=MICADO&SourceList.frame.instrument.name.op=%3D&SourceList.instrument.name=MICADO&SourceList.instrument.name.op=%3D | 3x4 25-coadds]] \\ etc.. | csv files | 
|no ADC, 1sec sampling| [[http://dbview.astro-wise.org/DbView?project=MICADO%40E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=500&fcontext_project_only=on&show_expanded=no&Exportselect=NoDownload&class_str=ReducedScienceFrame&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AsS%2527expanded%2527%250Ap3%250A%2528dp4%250AS%2527ReducedScienceFrame%2527%250Ap5%250A%2528lp6%250AS%2527open%2527%250Ap7%250AaS%2527root%2527%250Ap8%250Aag5%250AasS%2527ReducedScienceFrame.chip%2527%250Ap9%250A%2528lp10%250AS%2527prefilled%2527%250Ap11%250AaS%2527link_property%2527%250Ap12%250AaS%2527Chip%2527%250Ap13%250AasS%2527ReducedScienceFrame.instrument%2527%250Ap14%250A%2528lp15%250Ag11%250AaS%2527link_property%2527%250Ap16%250AaS%2527Instrument%2527%250Ap17%250AasS%2527ReducedScienceFrame.filter%2527%250Ap18%250A%2528lp19%250Ag11%250AaS%2527link_property%2527%250Ap20%250AaS%2527Filter%2527%250Ap21%250AassS%2527derived%2527%250Ap22%250A%2528dp23%250As.&expanded_leaf=&expanded_leaf_type=&expanded_leaf_class=&derived_leaf=&node_state=&ReducedScienceFrame.expanded=open&ReducedScienceFrame.chip.expanded=prefilled&ReducedScienceFrame.instrument.expanded=prefilled&ReducedScienceFrame.filter.expanded=prefilled&ReducedScienceFrame.creation_date.op=between&ReducedScienceFrame.creation_date=2014-02-18%2C+2014-02-20&ReducedScienceFrame.EXPTIME.op=%3D&ReducedScienceFrame.EXPTIME=0.001&ReducedScienceFrame.is_valid.op=%3E&ReducedScienceFrame.is_valid=0&ReducedScienceFrame.filter.name=MICADO_SIM_MONO&ReducedScienceFrame.filter.name.op=%3D&ReducedScienceFrame.instrument.name=MICADO&ReducedScienceFrame.instrument.name.op=%3D | 100 snapshots]] | etc..| csv files |
|For perfect pupil plane ADC, z=0 deg, 2 wavelengths | [[http://dbview.astro-wise.org/DbView?project=MICADO@E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=500&fcontext_project_only=on&show_expanded=no&Exportselect=NoDownload&class_str=ReducedScienceFrame&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AsS%2527expanded%2527%250Ap3%250A%2528dp4%250AS%2527ReducedScienceFrame%2527%250Ap5%250A%2528lp6%250AS%2527open%2527%250Ap7%250AaS%2527root%2527%250Ap8%250Aag5%250AasS%2527ReducedScienceFrame.chip%2527%250Ap9%250A%2528lp10%250AS%2527prefilled%2527%250Ap11%250AaS%2527link_property%2527%250Ap12%250AaS%2527Chip%2527%250Ap13%250AasS%2527ReducedScienceFrame.instrument%2527%250Ap14%250A%2528lp15%250Ag11%250AaS%2527link_property%2527%250Ap16%250AaS%2527Instrument%2527%250Ap17%250AasS%2527ReducedScienceFrame.filter%2527%250Ap18%250A%2528lp19%250Ag11%250AaS%2527link_property%2527%250Ap20%250AaS%2527Filter%2527%250Ap21%250AassS%2527derived%2527%250Ap22%250A%2528dp23%250As.&expanded_leaf=&expanded_leaf_type=&expanded_leaf_class=&derived_leaf=&node_state=&ReducedScienceFrame.expanded=open&ReducedScienceFrame.chip.expanded=prefilled&ReducedScienceFrame.instrument.expanded=prefilled&ReducedScienceFrame.filter.expanded=prefilled&ReducedScienceFrame.AIRMEND.op=%3D&ReducedScienceFrame.AIRMEND=1&ReducedScienceFrame.creation_date.op=%3C&ReducedScienceFrame.creation_date=2014-02-24&ReducedScienceFrame.is_valid.op=%3E&ReducedScienceFrame.is_valid=0&ReducedScienceFrame.filter.name=MICADO_SIM_MULT&ReducedScienceFrame.filter.name.op=%3D&ReducedScienceFrame.instrument.name=MICADO&ReducedScienceFrame.instrument.name.op=%3D | 217 stacks]] | [[http://dbview.astro-wise.org/DbView?project=MICADO@E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=500&fcontext_project_only=on&show_expanded=yes&Exportselect=NoDownload&class_str=SourceList&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AS%2527SourceList.instrument.name%2527%250Ap3%250A%2528S%2527%253D%2527%250Ap4%250AS%2527MICADO%2527%250Ap5%250Atp6%250AsS%2527SourceList.is_valid%2527%250Ap7%250A%2528S%2527%253E%2527%250Ap8%250AI0%250Atp9%250AsS%2527SourceList.creation_date%2527%250Ap10%250A%2528S%2527%253C%2527%250Ap11%250Acdatetime%250Adatetime%250Ap12%250A%2528S%2527%255Cx07%255Cxde%255Cx02%255Cx18%255Cx00%255Cx00%255Cx00%255Cx00%255Cx00%255Cx00%2527%250Ap13%250Atp14%250ARp15%250Atp16%250AssS%2527expanded%2527%250Ap17%250A%2528dp18%250AS%2527SourceList.frame.instrument%2527%250Ap19%250A%2528lp20%250AS%2527prefilled%2527%250Ap21%250AaS%2527link_property%2527%250Ap22%250AaS%2527Instrument%2527%250Ap23%250AasS%2527SourceList.frame.filter%2527%250Ap24%250A%2528lp25%250Ag21%250AaS%2527link_property%2527%250Ap26%250AaS%2527Filter%2527%250Ap27%250AasS%2527SourceList.chip%2527%250Ap28%250A%2528lp29%250AS%2527prefilled%2527%250Ap30%250AaS%2527link_property%2527%250Ap31%250AaS%2527Chip%2527%250Ap32%250AasS%2527SourceList.filter%2527%250Ap33%250A%2528lp34%250AS%2527prefilled%2527%250Ap35%250AaS%2527link_property%2527%250Ap36%250AaS%2527Filter%2527%250Ap37%250AasS%2527SourceList.frame%2527%250Ap38%250A%2528lp39%250AS%2527open%2527%250Ap40%250AaS%2527link_property%2527%250Ap41%250AaS%2527ReducedScienceFrame%2527%250Ap42%250AasS%2527SourceList.instrument%2527%250Ap43%250A%2528lp44%250AS%2527prefilled%2527%250Ap45%250AaS%2527link_property%2527%250Ap46%250AaS%2527Instrument%2527%250Ap47%250AasS%2527SourceList%2527%250Ap48%250A%2528lp49%250Ag40%250AaS%2527root%2527%250Ap50%250Aag48%250AasS%2527SourceList.frame.chip%2527%250Ap51%250A%2528lp52%250Ag21%250AaS%2527link_property%2527%250Ap53%250AaS%2527Chip%2527%250Ap54%250AassS%2527derived%2527%250Ap55%250A%2528dp56 | z=0, 217 sourcelists]] |
| For perfect pupil plane ADC z=30 deg | [[http://dbview.astro-wise.org/DbView?project=MICADO@E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=500&fcontext_project_only=on&show_expanded=no&Exportselect=NoDownload&class_str=ReducedScienceFrame&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AsS%2527expanded%2527%250Ap3%250A%2528dp4%250AS%2527ReducedScienceFrame%2527%250Ap5%250A%2528lp6%250AS%2527open%2527%250Ap7%250AaS%2527root%2527%250Ap8%250Aag5%250AasS%2527ReducedScienceFrame.chip%2527%250Ap9%250A%2528lp10%250AS%2527prefilled%2527%250Ap11%250AaS%2527link_property%2527%250Ap12%250AaS%2527Chip%2527%250Ap13%250AasS%2527ReducedScienceFrame.instrument%2527%250Ap14%250A%2528lp15%250Ag11%250AaS%2527link_property%2527%250Ap16%250AaS%2527Instrument%2527%250Ap17%250AasS%2527ReducedScienceFrame.filter%2527%250Ap18%250A%2528lp19%250Ag11%250AaS%2527link_property%2527%250Ap20%250AaS%2527Filter%2527%250Ap21%250AassS%2527derived%2527%250Ap22%250A%2528dp23%250As.&expanded_leaf=&expanded_leaf_type=&expanded_leaf_class=&derived_leaf=&node_state=&ReducedScienceFrame.expanded=open&ReducedScienceFrame.chip.expanded=prefilled&ReducedScienceFrame.instrument.expanded=prefilled&ReducedScienceFrame.filter.expanded=prefilled&ReducedScienceFrame.AIRMEND.op=%3D&ReducedScienceFrame.AIRMEND=1.1547&ReducedScienceFrame.creation_date.op=%3C&ReducedScienceFrame.creation_date=2014-02-24&ReducedScienceFrame.is_valid.op=%3E&ReducedScienceFrame.is_valid=0&ReducedScienceFrame.instrument.name=MICADO&ReducedScienceFrame.instrument.name.op=%3D | 217 stacks]] | [[http://dbview.astro-wise.org/DbView?project=MICADO@E-ELT&set_cookie=1&SQL=WITH+SQ30+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22frame.instrument.object_id%22+FROM+AWOPER.%22Instrument%2B%22++%0D%0A+WHERE+%22name%22+%3D+%27MICADO%27%29%0D%0A+%2CSQ31+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22frame.filter.object_id%22+FROM+AWOPER.%22Filter%2B%22++%0D%0A+WHERE+%22name%22+%3D+%27MICADO_SIM_MULT%27%29%0D%0A+%2CSQ20+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22frame.object_id%22%2C%22AIRMEND%22+AS+%22frame.AIRMEND%22%2C%22AIRMSTRT%22+AS+%22frame.AIRMSTRT%22%2C%22creation_date%22+AS+%22frame.creation_date%22%2C%22DATE%23%22+AS+%22frame.DATE%22%2C%22DATE_OBS%22+AS+%22frame.DATE_OBS%22%2C%22EXPTIME%22+AS+%22frame.EXPTIME%22%2C%22filename%22+AS+%22frame.filename%22%2C%22is_valid%22+AS+%22frame.is_valid%22%2C%22NAXIS1%22+AS+%22frame.NAXIS1%22%2C%22NAXIS2%22+AS+%22frame.NAXIS2%22%2C%22OBJECT%23%22+AS+%22frame.OBJECT%22%2C%22psf_radius%22+AS+%22frame.psf_radius%22%2C%22quality_flags%22+AS+%22frame.quality_flags%22%2C%22scale_factor%22+AS+%22frame.scale_factor%22%2CSQ30.*%2CSQ31.*+FROM+AWOPER.%22ReducedScienceFrame%2B%22++%0D%0A+JOIN+SQ30+ON+%22instrument@%22+%3D+SQ30.%22frame.instrument.object_id%22%0D%0A+JOIN+SQ31+ON+%22filter@%22+%3D+SQ31.%22frame.filter.object_id%22+WHERE+%22AIRMEND%22+%3D+1.1547%29%0D%0A+%2CSQ21+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22filter.object_id%22%2C%22name%22+AS+%22filter.name%22+FROM+AWOPER.%22Filter%2B%22++%0D%0A%29%0D%0A+%2CSQ22+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22chip.object_id%22%2C%22name%22+AS+%22chip.name%22+FROM+AWOPER.%22Chip%2B%22++%0D%0A%29%0D%0A+%2CSQ23+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22instrument.object_id%22%2C%22name%22+AS+%22instrument.name%22+FROM+AWOPER.%22Instrument%2B%22++%0D%0A%29%0D%0A+%2CSQ10+AS+%28%0D%0A+SELECT+%22%2BPROJECT%22+AS+%22project_id%22%2C%27quality+view%27+%22quality+view%22%2C%22%2BPRIVILEGES%22%2C%22object_id%22%2C%27re-process%27+%22re-process%22%2CSQ23.%22instrument.name%22%2CSQ22.%22chip.name%22%2CSQ21.%22filter.name%22%2C%22associatelist%22%2C%22COMBINE_METHOD%22%2C%22creation_da | 30deg 217 sourcelists]] | |
|For perfect pupil plane ADC z=35 deg | [[http://dbview.astro-wise.org/DbView?project=MICADO@E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=500&fcontext_project_only=on&show_expanded=no&Exportselect=NoDownload&class_str=ReducedScienceFrame&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AsS%2527expanded%2527%250Ap3%250A%2528dp4%250AS%2527ReducedScienceFrame%2527%250Ap5%250A%2528lp6%250AS%2527open%2527%250Ap7%250AaS%2527root%2527%250Ap8%250Aag5%250AasS%2527ReducedScienceFrame.chip%2527%250Ap9%250A%2528lp10%250AS%2527prefilled%2527%250Ap11%250AaS%2527link_property%2527%250Ap12%250AaS%2527Chip%2527%250Ap13%250AasS%2527ReducedScienceFrame.instrument%2527%250Ap14%250A%2528lp15%250Ag11%250AaS%2527link_property%2527%250Ap16%250AaS%2527Instrument%2527%250Ap17%250AasS%2527ReducedScienceFrame.filter%2527%250Ap18%250A%2528lp19%250Ag11%250AaS%2527link_property%2527%250Ap20%250AaS%2527Filter%2527%250Ap21%250AassS%2527derived%2527%250Ap22%250A%2528dp23%250As.&expanded_leaf=&expanded_leaf_type=&expanded_leaf_class=&derived_leaf=&node_state=&ReducedScienceFrame.expanded=open&ReducedScienceFrame.chip.expanded=prefilled&ReducedScienceFrame.instrument.expanded=prefilled&ReducedScienceFrame.filter.expanded=prefilled&ReducedScienceFrame.AIRMEND.op=%3D&ReducedScienceFrame.AIRMEND=1.22077&ReducedScienceFrame.creation_date.op=%3C&ReducedScienceFrame.creation_date=2014-02-24&ReducedScienceFrame.is_valid.op=%3E&ReducedScienceFrame.is_valid=0&ReducedScienceFrame.instrument.name=MICADO&ReducedScienceFrame.instrument.name.op=%3D | 217 stacks = 100x1, 50x2, 25x4, 20x5, 10x10, 4x25, 2x50, 1x100 ]] | [[http://dbview.astro-wise.org/DbView?project=MICADO@E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=500&fcontext_project_only=on&show_expanded=yes&Exportselect=NoDownload&class_str=SourceList&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AS%2527SourceList.instrument.name%2527%250Ap3%250A%2528S%2527%253D%2527%250Ap4%250AS%2527MICADO%2527%250Ap5%250Atp6%250AsS%2527SourceList.is_valid%2527%250Ap7%250A%2528S%2527%253E%2527%250Ap8%250AI0%250Atp9%250AsS%2527SourceList.creation_date%2527%250Ap10%250A%2528S%2527%253C%2527%250Ap11%250Acdatetime%250Adatetime%250Ap12%250A%2528S%2527%255Cx07%255Cxde%255Cx02%255Cx18%255Cx00%255Cx00%255Cx00%255Cx00%255Cx00%255Cx00%2527%250Ap13%250Atp14%250ARp15%250Atp16%250AssS%2527expanded%2527%250Ap17%250A%2528dp18%250AS%2527SourceList.frame.instrument%2527%250Ap19%250A%2528lp20%250AS%2527prefilled%2527%250Ap21%250AaS%2527link_property%2527%250Ap22%250AaS%2527Instrument%2527%250Ap23%250AasS%2527SourceList.frame.filter%2527%250Ap24%250A%2528lp25%250Ag21%250AaS%2527link_property%2527%250Ap26%250AaS%2527Filter%2527%250Ap27%250AasS%2527SourceList.chip%2527%250Ap28%250A%2528lp29%250AS%2527prefilled%2527%250Ap30%250AaS%2527link_property%2527%250Ap31%250AaS%2527Chip%2527%250Ap32%250AasS%2527SourceList.filter%2527%250Ap33%250A%2528lp34%250AS%2527prefilled%2527%250Ap35%250AaS%2527link_property%2527%250Ap36%250AaS%2527Filter%2527%250Ap37%250AasS%2527SourceList.frame%2527%250Ap38%250A%2528lp39%250AS%2527open%2527%250Ap40%250AaS%2527link_property%2527%250Ap41%250AaS%2527ReducedScienceFrame%2527%250Ap42%250AasS%2527SourceList.instrument%2527%250Ap43%250A%2528lp44%250AS%2527prefilled%2527%250Ap45%250AaS%2527link_property%2527%250Ap46%250AaS%2527Instrument%2527%250Ap47%250AasS%2527SourceList%2527%250Ap48%250A%2528lp49%250Ag40%250AaS%2527root%2527%250Ap50%250Aag48%250AasS%2527SourceList.frame.chip%2527%250Ap51%250A%2528lp52%250Ag21%250AaS%2527link_property%2527%250Ap53%250AaS%2527Chip%2527%250Ap54%250AassS%2527derived%2527%250Ap55%250A%2528dp56 | 35deg, 217 sourcelists]] | |
|For perfect pupil plane ADC z=60 deg, 2 wavelengths | [[http://dbview.astro-wise.org/DbView?project=MICADO@E-ELT&set_cookie=&QSort=creation_date&QDesc=descending&2ndQSort=creation_date&2ndQDesc=descending&numrows=500&fcontext_project_only=on&show_expanded=no&Exportselect=NoDownload&class_str=ReducedScienceFrame&mode=table_res&depth=1&python_dict=%2528dp0%250AS%2527query%2527%250Ap1%250A%2528dp2%250AsS%2527expanded%2527%250Ap3%250A%2528dp4%250AS%2527ReducedScienceFrame%2527%250Ap5%250A%2528lp6%250AS%2527open%2527%250Ap7%250AaS%2527root%2527%250Ap8%250Aag5%250AasS%2527ReducedScienceFrame.chip%2527%250Ap9%250A%2528lp10%250AS%2527prefilled%2527%250Ap11%250AaS%2527link_property%2527%250Ap12%250AaS%2527Chip%2527%250Ap13%250AasS%2527ReducedScienceFrame.instrument%2527%250Ap14%250A%2528lp15%250Ag11%250AaS%2527link_property%2527%250Ap16%250AaS%2527Instrument%2527%250Ap17%250AasS%2527ReducedScienceFrame.filter%2527%250Ap18%250A%2528lp19%250Ag11%250AaS%2527link_property%2527%250Ap20%250AaS%2527Filter%2527%250Ap21%250AassS%2527derived%2527%250Ap22%250A%2528dp23%250As.&expanded_leaf=&expanded_leaf_type=&expanded_leaf_class=&derived_leaf=&node_state=&ReducedScienceFrame.expanded=open&ReducedScienceFrame.chip.expanded=prefilled&ReducedScienceFrame.instrument.expanded=prefilled&ReducedScienceFrame.filter.expanded=prefilled&ReducedScienceFrame.AIRMEND.op=%3D&ReducedScienceFrame.AIRMEND=2&ReducedScienceFrame.creation_date.op=%3C&ReducedScienceFrame.creation_date=2014-02-24&ReducedScienceFrame.is_valid.op=%3E&ReducedScienceFrame.is_valid=0&ReducedScienceFrame.instrument.name=MICADO&ReducedScienceFrame.instrument.name.op=%3D | 217 stacks]] | [[http://dbview.astro-wise.org/DbView?project=MICADO@E-ELT&set_cookie=1&SQL=WITH+SQ30+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22frame.instrument.object_id%22+FROM+AWOPER.%22Instrument%2B%22++%0D%0A+WHERE+%22name%22+%3D+%27MICADO%27%29%0D%0A+%2CSQ31+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22frame.filter.object_id%22+FROM+AWOPER.%22Filter%2B%22++%0D%0A+WHERE+%22name%22+%3D+%27MICADO_SIM_MULT%27%29%0D%0A+%2CSQ20+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22frame.object_id%22%2C%22AIRMEND%22+AS+%22frame.AIRMEND%22%2C%22AIRMSTRT%22+AS+%22frame.AIRMSTRT%22%2C%22creation_date%22+AS+%22frame.creation_date%22%2C%22DATE%23%22+AS+%22frame.DATE%22%2C%22DATE_OBS%22+AS+%22frame.DATE_OBS%22%2C%22EXPTIME%22+AS+%22frame.EXPTIME%22%2C%22filename%22+AS+%22frame.filename%22%2C%22is_valid%22+AS+%22frame.is_valid%22%2C%22NAXIS1%22+AS+%22frame.NAXIS1%22%2C%22NAXIS2%22+AS+%22frame.NAXIS2%22%2C%22OBJECT%23%22+AS+%22frame.OBJECT%22%2C%22psf_radius%22+AS+%22frame.psf_radius%22%2C%22quality_flags%22+AS+%22frame.quality_flags%22%2C%22scale_factor%22+AS+%22frame.scale_factor%22%2CSQ30.*%2CSQ31.*+FROM+AWOPER.%22ReducedScienceFrame%2B%22++%0D%0A+JOIN+SQ30+ON+%22instrument@%22+%3D+SQ30.%22frame.instrument.object_id%22%0D%0A+JOIN+SQ31+ON+%22filter@%22+%3D+SQ31.%22frame.filter.object_id%22+WHERE+%22AIRMEND%22+%3D+2%29%0D%0A+%2CSQ21+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22filter.object_id%22%2C%22name%22+AS+%22filter.name%22+FROM+AWOPER.%22Filter%2B%22++%0D%0A%29%0D%0A+%2CSQ22+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22chip.object_id%22%2C%22name%22+AS+%22chip.name%22+FROM+AWOPER.%22Chip%2B%22++%0D%0A%29%0D%0A+%2CSQ23+AS+%28%0D%0A+SELECT+%22object_id%22+AS+%22instrument.object_id%22%2C%22name%22+AS+%22instrument.name%22+FROM+AWOPER.%22Instrument%2B%22++%0D%0A%29%0D%0A+%2CSQ10+AS+%28%0D%0A+SELECT+%22%2BPROJECT%22+AS+%22project_id%22%2C%27quality+view%27+%22quality+view%22%2C%22%2BPRIVILEGES%22%2C%22object_id%22%2C%27re-process%27+%22re-process%22%2CSQ23.%22instrument.name%22%2CSQ22.%22chip.name%22%2CSQ21.%22filter.name%22%2C%22associatelist%22%2C%22COMBINE_METHOD%22%2C%22creation_date%22 | 60deg, 217 sourcelists]]|

**team internal notes**

* The data as provided by LESIA+Remko Stuik can be accessed in this Dropbox [[ https://www.dropbox.com/sh/hxyiem2dkbbm1z4/8Ja_YWAxeU | simulations repository ]] **JMF**: add other locations, if they exist.
* Data format
  * Fits data with (subset?) of following header item: {{:projects:micado:header_simulations.txt|}}
  * A more trimmed down version with indication of a preferred minimum set of header keywords {{:projects:micado:minimum_header_for_simulations.txt|}}


== Analysis



=== Source extraction and centroiding

The simulated exposures contain a background with a gaussian noise of mean of 1000 ADU and sigma of 10 ADU. The two stars have peak pixel values approaching 60000 and 45000 ADU, have total fluxes of ~1E6 ADU covering 200-300 pixels and do not contain poisson errors. Thus their formal S/N>5500. 

For a well-sampled PSF the limit to centroid measurement accuracy sigma_c_m is (e.g., Kuijken et al 2012):
 * sigma_c_m ~0.7 FWHM / (S/N)
which for our simulated images translates to sigma_c_m 1 micro-arcsec for star at 2.2 micron.

We expect an additional systematic error due to: 
  * the speckled, not well-sampled PSF in individual exposures (see e.g., [[http://imageview.astro-wise.org/?Queryselect=WebImage&Formselect=ReducedScienceFrame&filename=Sci-JMCFARLAND-MICADO-------MICADO_SIM_MONO-MICADO_SIM_CHIP1-Red---Sci-56699.5999502-c46ce11f1a20e610a11bb857261f399dfd201dd3.fits&subdir=DB_img&resizewindow=1 | 1 exposure]] which smooth out when summed: [[http://imageview.astro-wise.org/?Queryselect=WebImage&Formselect=ReducedScienceFrame&filename=Sci-JMCFARLAND-MICADO-------MICADO_SIM_MONO-MICADO_SIM_CHIP1-Red---Sci-56699.6001573-627799c2a63b817455d03fd6109c55eaa28a512b.fits&subdir=DB_img&resizewindow=1 | 100 exposure]]. 
* simplistic centroiding: we extract sources using SExtractor measuring the position in two ways: the barycenter of the isophotal footprint and Gaussian fit to the flux distribution (i.e., X_IMAGE and XWIN_IMAGE). 
 
The single exposures have "speckled" PSFs being 0.001s snapshots. Due to this Sextractor detects a varying number of spurious sources in the wings of the PSFs. However, the central detection is dominant in terms of flux (contains at least ~5 times more flux than any other source). 

We use the following information from the source extractions:

* Position measures:
   * X_IMAGE/Y_IMAGE, ERRX2_IMAGE/ERRY2_IMAGE = barycenter (flux-weighted first moment) and their error (1sigma)
   * XWIN_IMAGE/YWIN_IMAGE = Gaussian windowed measurement of these parameters

* Profile measures:
   * FWHM_IMAGE (assuming gaussian profile)

* Flux measures:
   * FLUX_APER: flux inside a 25 pixel diameter aperture.


{{:projects:micado:centroiding5arcsecsource.jpg?400|}}
{{:projects:micado:centroiding10arcsecsource.jpg?400|}}

Fig1. Examples of centroid determinations in single exposure images for source 1 (left) and 2 (right). Shown are the central few pixels. The tiny green square indicates the measured centroid.


First we compared the distance measured from the barycentric method and Gaussian fitting. Figure 1 shows they typically deliver the same centroid to within 0.05 pixel with a few outliers for the single exposure images. We verified the standard deviation of the measured distances between source 1 and 2 in single exposure images is the same for both centroiding methods within 1%. This standard deviation is also the same within 1% if one takes only distances of exposures in which both centroiding methods for both sources agree to within 0.1 pixel. For this reason we restrict further analysis to barycentric centroiding measurements only. 



{{:projects:micado:barycentercentroidingvsgaussianfit.jpg?800|}}

Fig 1. Difference in X and Y coordinate of centroid determined from a barycentric and Gaussian fit. 


**Exposures sampled every 3sec**:
---------------
^ ID ^ source ^ exposures per stack ^ number of stacks ^ mean number of spurious detections per stack ^ mean X-distance (sigma) ^ mean Y-distance (sigma) ^ mean X_IMAGE  (sigma) ^ mean Y_IMAGE (sigma) ^ mean FWHM_IMAGE (sigma) ^ mean FLUX_APER (sigma) ^ NPIX ^
^    ^        ^                     ^                  ^                                              ^ [pixel]                 ^ [pixel]                 ^ [pixel]               ^ [pixel]              ^  [pixel]                ^ [ADU]                  ^      ^
|  1  | 5arcsec off-axis  | 1 | 100 |      |             |             | 1459.947 (0.550) | 1459.689 (0.354) | 3.871 (0.093) | 1259574 (105752) | 162 |
|  2  | 10arcsec off-axis | 1 | 100 |      |             |             | 2638.259 (0.705) | 2638.231 (0.427) | 3.902 (0.163) | 949341 (101678) | 150 |
| 2-1 |                   |    |    | 5.0 | 1178.312 (0.858) | 1178.541 (0.465) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 2 | 150 |      |             |             | 1459.905 (0.452) | 1459.733 (0.246) | 3.981 (0.188) | 1371170 (120839) | 200 |
|  2  | 10arcsec off-axis | 2 | 150 |      |             |             | 2638.297 (0.572) | 2638.260 (0.291) | 4.117 (0.312) | 1090340 (94302) | 216 |
| 2-1 |                   |    |    | 2.5 | 1178.391 (0.702) | 1178.527 (0.344) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 4 | 75 |      |             |             | 1459.889 (0.326) | 1459.734 (0.175) | 3.992 (0.121) | 1453868 (79965) | 238 |
|  2  | 10arcsec off-axis | 4 | 75 |      |             |             | 2638.312 (0.465) | 2638.255 (0.202) | 4.134 (0.206) | 1156714 (63395) | 253 |
| 2-1 |                   |    |    | 1.0 | 1178.422 (0.557) | 1178.521 (0.242) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 5 | 60 |      |             |             | 1459.875 (0.310) | 1459.717 (0.147) | 4.000 (0.123) | 1454038 (88757) | 239 |
|  2  | 10arcsec off-axis | 5 | 60 |      |             |             | 2638.314 (0.351) | 2638.236 (0.160) | 4.139 (0.197) | 1173078 (61769) | 270 |
| 2-1 |                   |    |    | 0.7 | 1178.439 (0.472) | 1178.519 (0.188) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 10 | 30 |      |             |             | 1459.880 (0.232) | 1459.749 (0.079) | 3.965 (0.053) | 1471525 (65031) | 253 |
|  2  | 10arcsec off-axis | 10 | 30 |      |             |             | 2638.313 (0.257) | 2638.252 (0.106) | 4.102 (0.090) | 1181422 (60024) | 274 |
| 2-1 |                   |    |    | 0.4 | 1178.433 (0.344) | 1178.502 (0.124) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 20 | 15 |      |             |             | 1459.889 (0.146) | 1459.743 (0.058) | 3.955 (0.028) | 1491539 (21963) | 275 |
|  2  | 10arcsec off-axis | 20 | 15 |      |             |             | 2638.308 (0.151) | 2638.255 (0.068) | 4.065 (0.056) | 1192187 (19476) | 284 |
| 2-1 |                   |    |    | 0.0 | 1178.419 (0.205) | 1178.512 (0.084) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 25 | 12 |      |             |             | 1459.888 (0.156) | 1459.742 (0.057) | 3.948 (0.025) | 1491653 (20958) | 274 |
|  2  | 10arcsec off-axis | 25 | 12 |      |             |             | 2638.309 (0.129) | 2638.259 (0.056) | 4.063 (0.043) | 1192162 (19058) | 284 |
| 2-1 |                   |    |    | 0.0 | 1178.421 (0.195) | 1178.517 (0.083) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 50 | 6 |      |             |             | 1459.890 (0.061) | 1459.743 (0.045) | 3.953 (0.023) | 1491299 (16119) | 276 |
|  2  | 10arcsec off-axis | 50 | 6 |      |             |             | 2638.305 (0.096) | 2638.256 (0.023) | 4.052 (0.026) | 1191872 (13888) | 284 |
| 2-1 |                   |    |    | 0.0 | 1178.415 (0.130) | 1178.513 (0.044) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 100 | 3 |      |             |             | 1459.890 (0.053) | 1459.740 (0.028) | 3.949 (0.004) | 1490996 (4364) | 275 |
|  2  | 10arcsec off-axis | 100 | 3 |      |             |             | 2638.308 (0.086) | 2638.257 (0.022) | 4.050 (0.014) | 1192059 (3823) | 285 |
| 2-1 |                   |    |    | 0.0 | 1178.417 (0.120) | 1178.517 (0.033) |             |             |             |             |      |

Table 1. Statistics of extracted source parameters for set 1 (no ADC) for the sets of single exposures and coadds of 2, 4, 5, 10, 20, 25, 50 and 100 single exposures sampled every 3 seconds. Averages and standard deviations are listed for most parameters. Description of parameters is given in the text.  The position measurements use SExtractor's default barycentric method.

**Exposures sampled every 1sec**:
---------------
^ ID ^ source ^ exposures per stack ^ number of stacks ^ mean number of spurious detections per stack ^ mean X-distance (sigma) ^ mean Y-distance (sigma) ^ mean X_IMAGE  (sigma) ^ mean Y_IMAGE (sigma) ^ mean FWHM_IMAGE (sigma) ^ mean FLUX_APER (sigma) ^ NPIX ^
^    ^        ^                     ^                  ^                                              ^ [pixel]                 ^ [pixel]                 ^ [pixel]               ^ [pixel]              ^  [pixel]                ^ [ADU]                  ^      ^
|  1  | 5arcsec off-axis  | 1 | 100 |      |             |             | 1459.868 (0.524) | 1459.732 (0.406) | 3.884 (0.141) | 1254071 (119759) | 163 |
|  2  | 10arcsec off-axis | 1 | 100 |      |             |             | 2638.264 (0.757) | 2638.255 (0.420) | 3.930 (0.214) | 966125 (110074) | 158 |
| 2-1 |                   |    |    | 5.1 | 1178.396 (0.893) | 1178.524 (0.475) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 2 | 50 |      |             |             | 1459.878 (0.417) | 1459.768 (0.275) | 3.987 (0.177) | 1366792 (128284) | 198 |
|  2  | 10arcsec off-axis | 2 | 50 |      |             |             | 2638.258 (0.656) | 2638.276 (0.363) | 4.086 (0.267) | 1095880 (88204) | 224 |
| 2-1 |                   |    |    | 2.5 | 1178.380 (0.805) | 1178.507 (0.435) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 4 | 25 |      |             |             | 1459.863 (0.357) | 1459.769 (0.194) | 3.969 (0.097) | 1439098 (94023) | 231 |
|  2  | 10arcsec off-axis | 4 | 25 |      |             |             | 2638.240 (0.454) | 2638.251 (0.306) | 4.188 (0.239) | 1150882 (87792) | 250 |
| 2-1 |                   |    |    | 1.3 | 1178.377 (0.567) | 1178.483 (0.273) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 5 | 20 |      |             |             | 1459.873 (0.362) | 1459.772 (0.193) | 4.023 (0.111) | 1435484 (109567) | 232 |
|  2  | 10arcsec off-axis | 5 | 20 |      |             |             | 2638.252 (0.529) | 2638.231 (0.150) | 4.152 (0.241) | 1163133 (65725) | 263 |
| 2-1 |                   |    |    | 1.1 | 1178.379 (0.626) | 1178.459 (0.245) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 10 | 10 |      |             |             | 1459.876 (0.194) | 1459.752 (0.102) | 3.980 (0.073) | 1485786 (62257) | 259 |
|  2  | 10arcsec off-axis | 10 | 10 |      |             |             | 2638.218 (0.361) | 2638.290 (0.133) | 4.147 (0.085) | 1191883 (42554) | 278 |
| 2-1 |                   |    |    | 0.4 | 1178.343 (0.385) | 1178.539 (0.200) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 20 | 5 |      |             |             | 1459.840 (0.168) | 1459.779 (0.060) | 3.951 (0.045) | 1446750 (113709) | 241 |
|  2  | 10arcsec off-axis | 20 | 5 |      |             |             | 2638.249 (0.130) | 2638.253 (0.051) | 4.083 (0.061) | 1196224 (22430) | 288 |
| 2-1 |                   |    |    | 0.2 | 1178.409 (0.237) | 1178.474 (0.094) |             |             |             |             |      |
|  1  | 5arcsec off-axis  | 25 | 4 |      |             |             | 1459.834 (0.138) | 1459.766 (0.046) | 3.948 (0.013) | 1491999 (19009) | 274 |
|  2  | 10arcsec off-axis | 25 | 4 |      |             |             | 2638.250 (0.146) | 2638.251 (0.026) | 4.085 (0.029) | 1195797 (16665) | 288 |
| 2-1 |                   |    |    | 0.0 | 1178.417 (0.162) | 1178.485 (0.054) |             |             |             |             |      |

Table 2. Equivalent of Table 1. but now for exposures sampled every 1 second.  Also, coadds of 50 and 100 single exposures have been omitted for statistical reasons (sample size too low).


In all images the sources have a FWHM ~ 4 pix (12 mas), i.e., close to the diffraction limited of 13 mas at 2.2 micron. The fact that it is slightly smaller could be due to the fact that the extraction does not identify the full wings of the profile as part of a single source. The second order moments X2 and Y2 indicate that the sources are elongated in the X direction compared to the Y direction. **Any logic why this might be the case?**. 

The difference in peak values (60000 ADU for the source at 5arcsec from the AO star and 45000 ADU for the source at 10 arcsec) is expected due to anisoplanatism (e.g. Trippe et al.) 

Figure 2 show that the scatter in the barycenters is anisotropic: it is a factor 2 higher in the X direction compared to the Y direction. **Any logic why scatter is higher in the direction parallel to the wind?**
 

{{:projects:micado:xpossource1.jpg?400|}}
{{:projects:micado:ypossource1.jpg?400|}}
{{:projects:micado:xpossource2.jpg?400|}}
{{:projects:micado:ypossource2.jpg?400|}}

Fig 2. Extracted barycenters in X and Y for source 1 and 2 in the single exposure images. 



{{:projects:micado:dxdy.jpg?400|}}
{{:projects:micado:dist.jpg?400|}}

Fig3. **Left**: Pixel position difference delta(X_IMAGE) (red) and delta(Y_IMAGE) (blue) for source 1 and 2 in single exposures. **Right**: distance between source 1 and 2 in milliarcsec. 




  
We want to determine if there is a systematic floor to the accuracy with which we can determine the distance between the two stars. 
In Fig 2. the red circles show that the standard deviation of the distance measurements across series of different coadds (consisting of 1, 2, 4, 5, 10 and 20 exposures) scales as 1/sqrt(#exposures). All coadd combinations of the 100 single exposures give a similarly accurate estimate of the distance with a standard deviation of ~200-250 microarcsec (shown as green diamonds). The exposure are statistically independent samplings of the distance and summing more and more is expected to decrease the standard deviation. In this fashion a standard deviation of 10 microarcsec is expected to be reached with (200/10)^2 * 100 = 40000 exposures. 
 

{{:projects:micado:twopointastrometricaccuracy.jpg?800|}}

Figure 2: Standard deviation in microarcsec of distance measurements measured in single exposure images and in coadds of 2,4,5, 10 and 20 exposures. Both the standard deviation of the ensemble is shown (red circles) and the standard deviation of the mean of the ensemble (green diamonds).      

 
Contributions to this error budget are coming from systematics in centroiding and differential tilt jitter

^Source of systematic error ^ size (micro-arcsec)^comment^
|Centroiding |  unknown| |
|Differential tilt jitter |  54| predicted |
| - |  200| observed|

The prediction for the differential tilt jitter is done by scaling the description in Trippe et al. to the specifications used in our simulations. 

The difference between the two error budgets is that differential tilt jitter is expected to broaden the PSF on timescales longer than the wind crossing time (~3sec) while systematic centroiding errors do not cause an intrinsic broadening. For a coadd of 100 exposures which are sampled over a period of 100 or 300seconds the FWHM = 3.904 pix for source 1 and FWHM = 4.05 pix for source 2. For the single exposures it is FWHM = 3.880 and FWHM = 3.905 pix. The <0.15 pix difference is much smaller than the observed scatter in the single exposures of standard deviation of ~0.6-0.7 pix in X and ~0.3-0.4 pix in Y. Especially since FWHM ~ 2 sigma for typical PSFs. The most likely explanation is that systematic errors in the centroiding dominate.

Given that the single exposures are 0.001s snapshots the 40000 exposures would correspond to 4s exposure times.

If the 3s samplings of 0.001s snapshots are representative of 3sec exposures this would amount to 3hrs of exposure time. 

----
{{:projects:micado:singleradial.png?300|}}
{{:projects:micado:singlex.png?300|}}
{{:projects:micado:singley.png?300|}}

Figures ?-?: Radial and profile plots (IRAF imexamine) of a single simulated exposure (0.001 second image)

{{:projects:micado:coadd100radial.png?300|}}
{{:projects:micado:coadd100x.png?300|}}
{{:projects:micado:coadd100y.png?300|}}

Figures ?-?: Radial and profile plots (IRAF imexamine) of a stack of 100 single simulated exposures

----
{{:projects:micado:singlecoadd.tgz|EPS versions of 6 figures above}}

=== perfect pupil place ADC simulations

  * for all tables below
    * ID 1 is a source 5 arcsec from AO/guide star and ID 2 is a source 10 arcsec from AO/guide star
    * atmospheric sampling is 1 second

**ZD = 0 degree with ADC**
^ ID ^ exposures per stack ^ number of stacks ^ mean number of spurious detections per stack ^ mean X-distance (sigma) ^ mean XWIN-distance (sigma) ^ mean Y-distance (sigma) ^ mean YWIN-distance (sigma) ^ mean X_IMAGE (sigma) ^ mean XWIN_IMAGE (sigma) ^ mean Y_IMAGE (sigma) ^ mean YWIN_IMAGE (sigma) ^ mean FWHM_IMAGE (sigma) ^ mean FLUX_APER (sigma) ^ NPIX ^
^    ^                     ^                  ^                                              ^ [pixel]                 ^ [pixel]                    ^ [pixel]                 ^ [pixel]                    ^ [pixel]              ^ [pixel]                 ^ [pixel]              ^ [pixel]                 ^ [pixel]                 ^ [ADU]                  ^      ^
|  1  | 1 | 100 |      |             |             |             |             | 2049.120 (0.535) | 2049.122 (0.526) | 2049.013 (0.365) | 2048.993 (0.398) | 3.940 (0.254) | 1325284 (110330) | 221 |
|  2  | 1 | 100 |      |             |             |             |             | 3715.675 (0.773) | 3715.675 (0.771) | 2049.028 (0.406) | 2048.995 (0.413) | 3.974 (0.309) | 1025233 (93542) | 208 |
| 2-1 |    |    | 2.7 | 1666.555 (0.900) | 1666.553 (0.919) | 0.015 (0.450) | 0.001 (0.470) |             |             |             |             |             |             |      |
|  1  | 2 | 50 |      |             |             |             |             | 2049.125 (0.469) | 2049.121 (0.404) | 2049.034 (0.250) | 2048.994 (0.271) | 4.027 (0.246) | 1415423 (89167) | 256 |
|  2  | 2 | 50 |      |             |             |             |             | 3715.652 (0.631) | 3715.668 (0.635) | 2049.012 (0.370) | 2048.996 (0.312) | 4.070 (0.244) | 1117526 (76906) | 259 |
| 2-1 |    |    | 1.1 | 1666.526 (0.778) | 1666.547 (0.769) | -0.022 (0.400) | 0.002 (0.345) |             |             |             |             |             |             |      |
|  1  | 4 | 25 |      |             |             |             |             | 2049.094 (0.333) | 2049.122 (0.300) | 2049.047 (0.187) | 2048.994 (0.197) | 3.990 (0.094) | 1463112 (62426) | 283 |
|  2  | 4 | 25 |      |             |             |             |             | 3715.696 (0.509) | 3715.662 (0.472) | 2048.984 (0.304) | 2048.997 (0.251) | 4.174 (0.251) | 1170023 (48370) | 290 |
| 2-1 |    |    | 0.2 | 1666.602 (0.615) | 1666.540 (0.563) | -0.063 (0.282) | 0.002 (0.216) |             |             |             |             |             |             |      |
|  1  | 5 | 20 |      |             |             |             |             | 2049.093 (0.313) | 2049.123 (0.268) | 2049.044 (0.172) | 2048.992 (0.203) | 3.970 (0.098) | 1465182 (68592) | 282 |
|  2  | 5 | 20 |      |             |             |             |             | 3715.650 (0.526) | 3715.664 (0.452) | 2048.991 (0.131) | 2048.998 (0.141) | 4.128 (0.167) | 1172167 (54107) | 294 |
| 2-1 |    |    | 0.2 | 1666.557 (0.599) | 1666.540 (0.507) | -0.053 (0.175) | 0.006 (0.177) |             |             |             |             |             |             |      |
|  1  | 10 | 10 |      |             |             |             |             | 2049.086 (0.185) | 2049.124 (0.176) | 2049.029 (0.101) | 2048.995 (0.079) | 3.941 (0.061) | 1472070 (49882) | 290 |
|  2  | 10 | 10 |      |             |             |             |             | 3715.658 (0.383) | 3715.659 (0.294) | 2048.992 (0.097) | 2048.999 (0.097) | 4.146 (0.054) | 1180005 (38939) | 301 |
| 2-1 |    |    | 0.0 | 1666.572 (0.442) | 1666.536 (0.364) | -0.036 (0.121) | 0.004 (0.096) |             |             |             |             |             |             |      |
|  1  | 20 | 5 |      |             |             |             |             | 2049.085 (0.173) | 2049.124 (0.145) | 2049.028 (0.070) | 2048.995 (0.044) | 3.929 (0.036) | 1471528 (36941) | 290 |
|  2  | 20 | 5 |      |             |             |             |             | 3715.653 (0.129) | 3715.666 (0.164) | 2048.990 (0.047) | 2049.000 (0.040) | 4.090 (0.060) | 1179025 (26043) | 302 |
| 2-1 |    |    | 0.0 | 1666.567 (0.236) | 1666.542 (0.195) | -0.038 (0.100) | 0.005 (0.049) |             |             |             |             |             |             |      |

**ZD = 30 degree with ADC**
^ ID ^ exposures per stack ^ number of stacks ^ mean number of spurious detections per stack ^ mean X-distance (sigma) ^ mean XWIN-distance (sigma) ^ mean Y-distance (sigma) ^ mean YWIN-distance (sigma) ^ mean X_IMAGE (sigma) ^ mean XWIN_IMAGE (sigma) ^ mean Y_IMAGE (sigma) ^ mean YWIN_IMAGE (sigma) ^ mean FWHM_IMAGE (sigma) ^ mean FLUX_APER (sigma) ^ NPIX ^
^    ^                     ^                  ^                                              ^ [pixel]                 ^ [pixel]                    ^ [pixel]                 ^ [pixel]                    ^ [pixel]              ^ [pixel]                 ^ [pixel]              ^ [pixel]                 ^ [pixel]                 ^ [ADU]                  ^      ^
|  1  | 1 | 100 |      |             |             |             |             | 2049.155 (0.519) | 2049.137 (0.526) | 2037.737 (0.379) | 2037.731 (0.398) | 3.925 (0.232) | 1320743 (108038) | 214 |
|  2  | 1 | 100 |      |             |             |             |             | 3715.693 (0.769) | 3715.690 (0.771) | 2037.079 (0.411) | 2037.044 (0.413) | 3.993 (0.317) | 1031447 (95043) | 210 |
| 2-1 |    |    | 2.8 | 1666.537 (0.897) | 1666.553 (0.919) | -0.658 (0.469) | -0.687 (0.470) |             |             |             |             |             |             |      |
|  1  | 2 | 50 |      |             |             |             |             | 2049.145 (0.441) | 2049.137 (0.404) | 2037.762 (0.238) | 2037.732 (0.271) | 4.010 (0.244) | 1425150 (86403) | 256 |
|  2  | 2 | 50 |      |             |             |             |             | 3715.678 (0.628) | 3715.684 (0.634) | 2037.061 (0.365) | 2037.045 (0.311) | 4.065 (0.245) | 1123847 (76615) | 259 |
| 2-1 |    |    | 1.1 | 1666.533 (0.780) | 1666.547 (0.769) | -0.701 (0.401) | -0.687 (0.344) |             |             |             |             |             |             |      |
|  1  | 4 | 25 |      |             |             |             |             | 2049.117 (0.318) | 2049.137 (0.300) | 2037.790 (0.178) | 2037.733 (0.197) | 3.982 (0.092) | 1473584 (54922) | 281 |
|  2  | 4 | 25 |      |             |             |             |             | 3715.704 (0.503) | 3715.677 (0.472) | 2037.027 (0.301) | 2037.046 (0.251) | 4.170 (0.244) | 1177129 (48664) | 291 |
| 2-1 |    |    | 0.3 | 1666.587 (0.603) | 1666.539 (0.563) | -0.762 (0.275) | -0.687 (0.216) |             |             |             |             |             |             |      |
|  1  | 5 | 20 |      |             |             |             |             | 2049.116 (0.313) | 2049.139 (0.268) | 2037.773 (0.175) | 2037.730 (0.203) | 3.973 (0.094) | 1478942 (62801) | 285 |
|  2  | 5 | 20 |      |             |             |             |             | 3715.665 (0.526) | 3715.679 (0.452) | 2037.043 (0.132) | 2037.047 (0.141) | 4.126 (0.163) | 1180548 (53634) | 295 |
| 2-1 |    |    | 0.1 | 1666.548 (0.600) | 1666.540 (0.507) | -0.730 (0.175) | -0.683 (0.177) |             |             |             |             |             |             |      |
|  1  | 10 | 10 |      |             |             |             |             | 2049.101 (0.187) | 2049.139 (0.176) | 2037.767 (0.102) | 2037.733 (0.079) | 3.933 (0.068) | 1490360 (52927) | 289 |
|  2  | 10 | 10 |      |             |             |             |             | 3715.670 (0.383) | 3715.674 (0.294) | 2037.040 (0.097) | 2037.049 (0.097) | 4.145 (0.056) | 1194461 (41288) | 303 |
| 2-1 |    |    | 0.0 | 1666.569 (0.441) | 1666.535 (0.364) | -0.727 (0.124) | -0.685 (0.097) |             |             |             |             |             |             |      |
|  1  | 20 | 5 |      |             |             |             |             | 2049.103 (0.176) | 2049.139 (0.145) | 2037.769 (0.072) | 2037.734 (0.044) | 3.928 (0.059) | 1489307 (36116) | 291 |
|  2  | 20 | 5 |      |             |             |             |             | 3715.670 (0.131) | 3715.681 (0.164) | 2037.043 (0.046) | 2037.049 (0.040) | 4.087 (0.053) | 1193503 (24630) | 305 |
| 2-1 |    |    | 0.0 | 1666.567 (0.243) | 1666.542 (0.195) | -0.726 (0.100) | -0.685 (0.049) |             |             |             |             |             |             |      |

**ZD = 35 degree with ADC**
^ ID ^ exposures per stack ^ number of stacks ^ mean number of spurious detections per stack ^ mean X-distance (sigma) ^ mean XWIN-distance (sigma) ^ mean Y-distance (sigma) ^ mean YWIN-distance (sigma) ^ mean X_IMAGE (sigma) ^ mean XWIN_IMAGE (sigma) ^ mean Y_IMAGE (sigma) ^ mean YWIN_IMAGE (sigma) ^ mean FWHM_IMAGE (sigma) ^ mean FLUX_APER (sigma) ^ NPIX ^
^    ^                     ^                  ^                                              ^ [pixel]                 ^ [pixel]                    ^ [pixel]                 ^ [pixel]                    ^ [pixel]              ^ [pixel]                 ^ [pixel]              ^ [pixel]                 ^ [pixel]                 ^ [ADU]                  ^      ^
|  1  | 1 | 100 |      |             |             |             |             | 2049.141 (0.521) | 2049.133 (0.526) | 2039.717 (0.371) | 2039.705 (0.398) | 3.931 (0.241) | 1323013 (107868) | 216 |
|  2  | 1 | 100 |      |             |             |             |             | 3715.681 (0.776) | 3715.685 (0.771) | 2040.013 (0.410) | 2039.972 (0.413) | 3.988 (0.325) | 1030178 (94907) | 209 |
| 2-1 |    |    | 2.8 | 1666.540 (0.895) | 1666.553 (0.920) | 0.296 (0.464) | 0.267 (0.470) |             |             |             |             |             |             |      |
|  1  | 2 | 50 |      |             |             |             |             | 2049.145 (0.440) | 2049.132 (0.404) | 2039.740 (0.236) | 2039.706 (0.270) | 4.014 (0.244) | 1425881 (86213) | 256 |
|  2  | 2 | 50 |      |             |             |             |             | 3715.672 (0.660) | 3715.679 (0.635) | 2039.976 (0.372) | 2039.973 (0.312) | 4.099 (0.298) | 1127743 (73798) | 260 |
| 2-1 |    |    | 1.1 | 1666.527 (0.800) | 1666.547 (0.770) | 0.236 (0.416) | 0.267 (0.344) |             |             |             |             |             |             |      |
|  1  | 4 | 25 |      |             |             |             |             | 2049.103 (0.332) | 2049.132 (0.300) | 2039.759 (0.184) | 2039.707 (0.197) | 3.986 (0.093) | 1477319 (55055) | 283 |
|  2  | 4 | 25 |      |             |             |             |             | 3715.692 (0.494) | 3715.672 (0.472) | 2039.967 (0.279) | 2039.974 (0.249) | 4.160 (0.230) | 1178620 (41568) | 291 |
| 2-1 |    |    | 0.2 | 1666.589 (0.614) | 1666.540 (0.564) | 0.208 (0.267) | 0.267 (0.216) |             |             |             |             |             |             |      |
|  1  | 5 | 20 |      |             |             |             |             | 2049.111 (0.314) | 2049.134 (0.268) | 2039.747 (0.175) | 2039.704 (0.203) | 3.975 (0.098) | 1480382 (61721) | 285 |
|  2  | 5 | 20 |      |             |             |             |             | 3715.661 (0.526) | 3715.674 (0.452) | 2039.971 (0.132) | 2039.974 (0.141) | 4.128 (0.165) | 1182443 (51613) | 295 |
| 2-1 |    |    | 0.1 | 1666.550 (0.600) | 1666.540 (0.507) | 0.223 (0.177) | 0.270 (0.177) |             |             |             |             |             |             |      |
|  1  | 10 | 10 |      |             |             |             |             | 2049.097 (0.188) | 2049.134 (0.176) | 2039.741 (0.103) | 2039.707 (0.079) | 3.934 (0.067) | 1493535 (53225) | 290 |
|  2  | 10 | 10 |      |             |             |             |             | 3715.664 (0.383) | 3715.670 (0.294) | 2039.970 (0.098) | 2039.976 (0.097) | 4.141 (0.057) | 1197191 (41426) | 303 |
| 2-1 |    |    | 0.0 | 1666.567 (0.442) | 1666.535 (0.363) | 0.229 (0.126) | 0.269 (0.096) |             |             |             |             |             |             |      |
|  1  | 20 | 5 |      |             |             |             |             | 2049.098 (0.174) | 2049.134 (0.145) | 2039.743 (0.070) | 2039.708 (0.044) | 3.927 (0.057) | 1492422 (36189) | 290 |
|  2  | 20 | 5 |      |             |             |             |             | 3715.664 (0.128) | 3715.677 (0.164) | 2039.968 (0.045) | 2039.976 (0.040) | 4.093 (0.068) | 1196253 (24687) | 304 |
| 2-1 |    |    | 0.0 | 1666.566 (0.241) | 1666.542 (0.195) | 0.225 (0.099) | 0.269 (0.049) |             |             |             |             |             |             |      |

**ZD = 60 degree with ADC**
^ ID ^ exposures per stack ^ number of stacks ^ mean number of spurious detections per stack ^ mean X-distance (sigma) ^ mean XWIN-distance (sigma) ^ mean Y-distance (sigma) ^ mean YWIN-distance (sigma) ^ mean X_IMAGE (sigma) ^ mean XWIN_IMAGE (sigma) ^ mean Y_IMAGE (sigma) ^ mean YWIN_IMAGE (sigma) ^ mean FWHM_IMAGE (sigma) ^ mean FLUX_APER (sigma) ^ NPIX ^
^    ^                     ^                  ^                                              ^ [pixel]                 ^ [pixel]                    ^ [pixel]                 ^ [pixel]                    ^ [pixel]              ^ [pixel]                 ^ [pixel]              ^ [pixel]                 ^ [pixel]                 ^ [ADU]                  ^      ^
|  1  | 1 | 100 |      |             |             |             |             | 2049.155 (0.540) | 2049.147 (0.526) | 2021.238 (0.370) | 2021.233 (0.398) | 3.942 (0.288) | 1324296 (109351) | 219 |
|  2  | 1 | 100 |      |             |             |             |             | 3715.695 (0.765) | 3715.699 (0.771) | 2021.299 (0.415) | 2021.271 (0.414) | 3.960 (0.286) | 1022258 (96416) | 207 |
| 2-1 |    |    | 2.7 | 1666.540 (0.899) | 1666.551 (0.920) | 0.062 (0.453) | 0.039 (0.470) |             |             |             |             |             |             |      |
|  1  | 2 | 50 |      |             |             |             |             | 2049.153 (0.438) | 2049.147 (0.404) | 2021.289 (0.246) | 2021.233 (0.271) | 4.024 (0.238) | 1427484 (82806) | 259 |
|  2  | 2 | 50 |      |             |             |             |             | 3715.682 (0.628) | 3715.693 (0.635) | 2021.289 (0.369) | 2021.273 (0.311) | 4.066 (0.242) | 1122864 (78126) | 260 |
| 2-1 |    |    | 1.1 | 1666.528 (0.773) | 1666.545 (0.769) | 0.001 (0.416) | 0.040 (0.344) |             |             |             |             |             |             |      |
|  1  | 4 | 25 |      |             |             |             |             | 2049.122 (0.332) | 2049.148 (0.300) | 2021.285 (0.181) | 2021.234 (0.197) | 3.986 (0.098) | 1474357 (64167) | 281 |
|  2  | 4 | 25 |      |             |             |             |             | 3715.704 (0.498) | 3715.686 (0.472) | 2021.234 (0.288) | 2021.273 (0.251) | 4.129 (0.192) | 1176266 (54692) | 290 |
| 2-1 |    |    | 0.3 | 1666.582 (0.599) | 1666.538 (0.563) | -0.051 (0.275) | 0.039 (0.216) |             |             |             |             |             |             |      |
|  1  | 5 | 20 |      |             |             |             |             | 2049.134 (0.318) | 2049.149 (0.268) | 2021.261 (0.164) | 2021.231 (0.203) | 3.960 (0.096) | 1472308 (65987) | 278 |
|  2  | 5 | 20 |      |             |             |             |             | 3715.670 (0.525) | 3715.688 (0.452) | 2021.263 (0.127) | 2021.275 (0.141) | 4.114 (0.179) | 1181814 (51974) | 295 |
| 2-1 |    |    | 0.2 | 1666.536 (0.594) | 1666.539 (0.507) | 0.002 (0.172) | 0.044 (0.177) |             |             |             |             |             |             |      |
|  1  | 10 | 10 |      |             |             |             |             | 2049.111 (0.187) | 2049.149 (0.176) | 2021.267 (0.106) | 2021.235 (0.079) | 3.930 (0.056) | 1487282 (50494) | 289 |
|  2  | 10 | 10 |      |             |             |             |             | 3715.681 (0.382) | 3715.683 (0.294) | 2021.270 (0.096) | 2021.277 (0.097) | 4.122 (0.063) | 1192211 (39429) | 302 |
| 2-1 |    |    | 0.0 | 1666.570 (0.441) | 1666.534 (0.363) | 0.003 (0.128) | 0.042 (0.097) |             |             |             |             |             |             |      |
|  1  | 20 | 5 |      |             |             |             |             | 2049.113 (0.176) | 2049.149 (0.145) | 2021.264 (0.072) | 2021.235 (0.044) | 3.920 (0.037) | 1486682 (39624) | 291 |
|  2  | 20 | 5 |      |             |             |             |             | 3715.680 (0.129) | 3715.690 (0.164) | 2021.271 (0.049) | 2021.277 (0.040) | 4.073 (0.058) | 1191547 (28934) | 303 |
| 2-1 |    |    | 0.0 | 1666.567 (0.241) | 1666.541 (0.195) | 0.007 (0.103) | 0.042 (0.050) |             |             |             |             |             |             |      |

**ZD = 0 without ADC**
^ ID ^ exposures per stack ^ number of stacks ^ mean number of spurious detections per stack ^ mean X-distance (sigma) ^ mean XWIN-distance (sigma) ^ mean Y-distance (sigma) ^ mean YWIN-distance (sigma) ^ mean X_IMAGE (sigma) ^ mean XWIN_IMAGE (sigma) ^ mean Y_IMAGE (sigma) ^ mean YWIN_IMAGE (sigma) ^ mean FWHM_IMAGE (sigma) ^ mean FLUX_APER (sigma) ^ NPIX ^
^    ^                     ^                  ^                                              ^ [pixel]                 ^ [pixel]                    ^ [pixel]                 ^ [pixel]                    ^ [pixel]              ^ [pixel]                 ^ [pixel]              ^ [pixel]                 ^ [pixel]                 ^ [ADU]                  ^      ^
|  1  | 1 | 100 |      |             |             |             |             | 1459.868 (0.524) | 1459.866 (0.526) | 1459.732 (0.406) | 1459.736 (0.400) | 3.884 (0.141) | 1254071 (119759) | 163 |
|  2  | 1 | 100 |      |             |             |             |             | 2638.264 (0.757) | 2638.265 (0.772) | 2638.255 (0.420) | 2638.249 (0.415) | 3.930 (0.214) | 966125 (110074) | 158 |
| 2-1 |    |    | 5.1 | 1178.396 (0.893) | 1178.399 (0.921) | 1178.524 (0.475) | 1178.513 (0.472) |             |             |             |             |             |             |      |
|  1  | 2 | 50 |      |             |             |             |             | 1459.878 (0.417) | 1459.866 (0.404) | 1459.768 (0.275) | 1459.737 (0.272) | 3.987 (0.177) | 1366792 (128284) | 198 |
|  2  | 2 | 50 |      |             |             |             |             | 2638.258 (0.656) | 2638.258 (0.637) | 2638.276 (0.363) | 2638.252 (0.312) | 4.086 (0.267) | 1095880 (88204) | 224 |
| 2-1 |    |    | 2.5 | 1178.380 (0.805) | 1178.392 (0.772) | 1178.507 (0.435) | 1178.515 (0.346) |             |             |             |             |             |             |      |
|  1  | 4 | 25 |      |             |             |             |             | 1459.863 (0.357) | 1459.867 (0.301) | 1459.769 (0.194) | 1459.738 (0.198) | 3.969 (0.097) | 1439098 (94023) | 231 |
|  2  | 4 | 25 |      |             |             |             |             | 2638.240 (0.454) | 2638.251 (0.472) | 2638.251 (0.306) | 2638.253 (0.252) | 4.188 (0.239) | 1150882 (87792) | 250 |
| 2-1 |    |    | 1.3 | 1178.377 (0.567) | 1178.384 (0.563) | 1178.483 (0.273) | 1178.516 (0.217) |             |             |             |             |             |             |      |
|  1  | 5 | 20 |      |             |             |             |             | 1459.873 (0.362) | 1459.869 (0.268) | 1459.772 (0.193) | 1459.735 (0.204) | 4.023 (0.111) | 1435484 (109567) | 232 |
|  2  | 5 | 20 |      |             |             |             |             | 2638.252 (0.529) | 2638.253 (0.452) | 2638.231 (0.150) | 2638.252 (0.144) | 4.152 (0.241) | 1163133 (65725) | 263 |
| 2-1 |    |    | 1.1 | 1178.379 (0.626) | 1178.385 (0.508) | 1178.459 (0.245) | 1178.517 (0.180) |             |             |             |             |             |             |      |
|  1  | 10 | 10 |      |             |             |             |             | 1459.876 (0.194) | 1459.869 (0.177) | 1459.752 (0.102) | 1459.739 (0.079) | 3.980 (0.073) | 1485786 (62257) | 259 |
|  2  | 10 | 10 |      |             |             |             |             | 2638.218 (0.361) | 2638.248 (0.294) | 2638.290 (0.133) | 2638.255 (0.097) | 4.147 (0.085) | 1191883 (42554) | 278 |
| 2-1 |    |    | 0.4 | 1178.343 (0.385) | 1178.379 (0.364) | 1178.539 (0.200) | 1178.516 (0.096) |             |             |             |             |             |             |      |
|  1  | 20 | 5 |      |             |             |             |             | 1459.840 (0.168) | 1459.869 (0.145) | 1459.779 (0.060) | 1459.740 (0.044) | 3.951 (0.045) | 1446750 (113709) | 241 |
|  2  | 20 | 5 |      |             |             |             |             | 2638.249 (0.130) | 2638.256 (0.165) | 2638.253 (0.051) | 2638.255 (0.040) | 4.083 (0.061) | 1196224 (22430) | 288 |
| 2-1 |    |    | 0.2 | 1178.409 (0.237) | 1178.387 (0.195) | 1178.474 (0.094) | 1178.516 (0.049) |             |             |             |             |             |             |      |

  * summary tables

**ZD = 0,30,35,60 degree with ADC**:
^ ID ^ exposures per stack ^ number of stacks ^ mean number of spurious detections per stack ^ mean X-distance (sigma) ^ mean XWIN-distance (sigma) ^ mean Y-distance (sigma) ^ mean YWIN-distance (sigma) ^ mean X_IMAGE (sigma) ^ mean XWIN_IMAGE (sigma) ^ mean Y_IMAGE (sigma) ^ mean YWIN_IMAGE (sigma) ^ mean FWHM_IMAGE (sigma) ^ mean FLUX_APER (sigma) ^ NPIX ^
^    ^                     ^                  ^                                              ^ [pixel]                 ^ [pixel]                    ^ [pixel]                 ^ [pixel]                    ^ [pixel]              ^ [pixel]                 ^ [pixel]              ^ [pixel]                 ^ [pixel]                 ^ [ADU]                  ^      ^
|  1  | 10 | 10 |      |             |             |             |             | 2049.086 (0.185) | 2049.124 (0.176) | 2049.029 (0.101) | 2048.995 (0.079) | 3.941 (0.061) | 1472070 (49882) | 290 |
|  2  | 10 | 10 |      |             |             |             |             | 3715.658 (0.383) | 3715.659 (0.294) | 2048.992 (0.097) | 2048.999 (0.097) | 4.146 (0.054) | 1180005 (38939) | 301 |
| 2-1 |    |    | 0.0 | 1666.572 (0.442) | 1666.536 (0.364) | -0.036 (0.121) | 0.004 (0.096) |             |             |             |             |             |             |      |
|  1  | 10 | 10 |      |             |             |             |             | 2049.101 (0.187) | 2049.139 (0.176) | 2037.767 (0.102) | 2037.733 (0.079) | 3.933 (0.068) | 1490360 (52927) | 289 |
|  2  | 10 | 10 |      |             |             |             |             | 3715.670 (0.383) | 3715.674 (0.294) | 2037.040 (0.097) | 2037.049 (0.097) | 4.145 (0.056) | 1194461 (41288) | 303 |
| 2-1 |    |    | 0.0 | 1666.569 (0.441) | 1666.535 (0.364) | -0.727 (0.124) | -0.685 (0.097) |             |             |             |             |             |             |      |
|  1  | 10 | 10 |      |             |             |             |             | 2049.097 (0.188) | 2049.134 (0.176) | 2039.741 (0.103) | 2039.707 (0.079) | 3.934 (0.067) | 1493535 (53225) | 290 |
|  2  | 10 | 10 |      |             |             |             |             | 3715.664 (0.383) | 3715.670 (0.294) | 2039.970 (0.098) | 2039.976 (0.097) | 4.141 (0.057) | 1197191 (41426) | 303 |
| 2-1 |    |    | 0.0 | 1666.567 (0.442) | 1666.535 (0.363) | 0.229 (0.126) | 0.269 (0.096) |             |             |             |             |             |             |      |
|  1  | 10 | 10 |      |             |             |             |             | 2049.111 (0.187) | 2049.149 (0.176) | 2021.267 (0.106) | 2021.235 (0.079) | 3.930 (0.056) | 1487282 (50494) | 289 |
|  2  | 10 | 10 |      |             |             |             |             | 3715.681 (0.382) | 3715.683 (0.294) | 2021.270 (0.096) | 2021.277 (0.097) | 4.122 (0.063) | 1192211 (39429) | 302 |
| 2-1 |    |    | 0.0 | 1666.570 (0.441) | 1666.534 (0.363) | 0.003 (0.128) | 0.042 (0.097) |             |             |             |             |             |             |      |

to be compared with:

**ZD = 0 deg without ADC**:
^ ID ^ exposures per stack ^ number of stacks ^ mean number of spurious detections per stack ^ mean X-distance (sigma) ^ mean XWIN-distance (sigma) ^ mean Y-distance (sigma) ^ mean YWIN-distance (sigma) ^ mean X_IMAGE (sigma) ^ mean XWIN_IMAGE (sigma) ^ mean Y_IMAGE (sigma) ^ mean YWIN_IMAGE (sigma) ^ mean FWHM_IMAGE (sigma) ^ mean FLUX_APER (sigma) ^ NPIX ^
^    ^                     ^                  ^                                              ^ [pixel]                 ^ [pixel]                    ^ [pixel]                 ^ [pixel]                    ^ [pixel]              ^ [pixel]                 ^ [pixel]              ^ [pixel]                 ^ [pixel]                 ^ [ADU]                  ^      ^
|  1  | 10 | 10 |      |             |             |             |             | 1459.876 (0.194) | 1459.869 (0.177) | 1459.752 (0.102) | 1459.739 (0.079) | 3.980 (0.073) | 1485786 (62257) | 259 |
|  2  | 10 | 10 |      |             |             |             |             | 2638.218 (0.361) | 2638.248 (0.294) | 2638.290 (0.133) | 2638.255 (0.097) | 4.147 (0.085) | 1191883 (42554) | 278 |
| 2-1 |    |    | 0.4 | 1178.343 (0.385) | 1178.379 (0.364) | 1178.539 (0.200) | 1178.516 (0.096) |             |             |             |             |             |             |      |



== Reference frame construction **editor GVK**

* what makes most sense to determine astrometric error budget for ADC:
   *reference frame = first image: to quantify astrometric rms
   *reference frame = derived from ensemble of images: to quantify sqrt(n) effects
   *reference frame = true positions?: best for simulated images to assess both of the above.


== Relative astrometry **editor GVK**

We then perform relative astrometric analysis.
For each dataset we quantify the relative astrometric quality: rms(position) within a set of dithers using the following approaches:

* residuals: pos(im)-pos(ref_im)
* do usual approach? (e.g., Kuijken/Rich papers) 
* Existing astrometric pipeline in Astro-WISE. It uses sextractor for source extraction and LDAC and SCAMP for astrometric calibration. Configuration is polynomial mapping. Needs fine-tuning to get beyond global astrometry of rms(0.03)/fwhm(0.7)=4e-2
* GEMS method
* astrom solution based on ensemble of (dithered) exposures (use current astrometric pipeline (SCAMP and/or LDAC))
  * output: x,y,ra,dec of stars + parametric fit to systematic distorsions over FoV.
  * Quantify distorsions: parameterize or look-up table?



== Discussion

* Compare MICADO astrometric results to other (MAD, GemS) results in units PSF/pixelscale?

== Conclusions


* outlook on future:
  * source detection on stacked data (light-curve software relevant here?) = master astrometric reference frame
  * can we exploit non-destructive read-outs? (contact Gert Finger)
  * better PSF modelling (starfinder, daophot, inspiration in Anderson + King-2000PASP112,1360 and KiDS PSF anisotropy/homogenization measurements?)
  * possibly astrometric calibration of GemS/MAD data (to get handle on MCAO case)

== References and handy-links

* [[https://www.dropbox.com/home/MICADO_ADC/ADC_documents | Dropbox astrometry documents]] 
* [[http://www.astro-wise.org/Public/Astro-WISEAstrometryReport.pdf | Astro-WISE astrometry report]]
* [[https://docs.google.com/spreadsheet/ccc?key=0Ah9tUSrro9zCdEw5TGFMX1VwTDFhOUhpQy1JVWxRQmc&usp=sharing | Gijs' MICADO spreadsheets]]
* [[http://www.astro-wise.org/viewvc/awe/astro/experimental/micado/ | MICADO code repository at OmegaCEN]]