{"id":11051,"date":"2020-04-03T04:20:32","date_gmt":"2020-04-03T14:20:32","guid":{"rendered":"http:\/\/www.eaobservatory.org\/jcmt\/?p=11051"},"modified":"2020-12-29T11:03:29","modified_gmt":"2020-12-29T21:03:29","slug":"new-de-biasing-method-for-pol2-data","status":"publish","type":"post","link":"https:\/\/www.eaobservatory.org\/jcmt\/2020\/04\/new-de-biasing-method-for-pol2-data\/","title":{"rendered":"New de-biasing method for POL2 data"},"content":{"rendered":"<p>Calculation of the polarised intensity (PI) and percentage polarisation (P) values from measurements of I, Q and U usually need to be corrected to take account of the bias towards larger values caused by the presence of noise (particularly noticeable\u00a0 in areas of low signal-to-noise). This bias occurs because of the squaring of Q and U involved in calculating PI:<\/p>\n<p>PI = \u221a( Q<sup>2<\/sup> + U<sup>2 <\/sup>)<\/p>\n<p>Thus both positive and negative Q and U values give a positive PI value, pushing the mean PI value upwards. A correction for this effect is applied if the parameter setting &#8220;<span style=\"font-family: 'courier new', courier, monospace\">debias=yes<\/span>&#8221; is specified when running the <span style=\"font-family: 'courier new', courier, monospace\">SMURF<\/span> <span style=\"font-family: 'courier new', courier, monospace\">pol2map<\/span> command (note, the default is &#8220;<span style=\"font-family: 'courier new', courier, monospace\">debias=no<\/span>&#8220;). To date, this correction has used the so-called &#8220;asymptotic estimator&#8221; (AS):<\/p>\n<p>PI<sub>AS<\/sub> = \u221a( Q<sup>2<\/sup> + U<sup>2 <\/sup> &#8211; \u03c3<sup>2 <\/sup>)<\/p>\n<p>where \u03c3<sup>2 <\/sup>is the weighted mean of the variances on Q and U:<\/p>\n<p>\u03c3<sup>2<\/sup> = ( Q<sup>2<\/sup>*\u03c3<sub>Q<\/sub><sup>2<\/sup> + U<sup>2<\/sup>*\u03c3<sub>U<\/sub><sup>2<\/sup> ) \/ ( Q<sup>2<\/sup> + U<sup>2<\/sup> )<\/p>\n<p>However, this estimator has the problem that it becomes undefined when<\/p>\n<p>Q<sup>2<\/sup> + U<sup>2 <\/sup>&lt; \u03c3<sup>2\u00a0<\/sup><\/p>\n<p>The usual practice is to use an exact value of zero for the PI at such points. However, this can upset the PI noise statistics in very low SNR regions.<\/p>\n<p>However, there are other forms of de-biasing that avoid this problem, such as the &#8220;modified asymptotic estimator&#8221; (MAS) described by\u00a0<a href=\"https:\/\/academic.oup.com\/mnras\/article\/439\/4\/4048\/1175299\">Plaszczynski et al\u00a0 (MNRAS 2014).<\/a>\u00a0 \u00a0The de-biased MAS estimate of PI is defined by<\/p>\n<p>PI<sub>MAS<\/sub> = PI &#8211; 0.5*\u03c3<sup>2<\/sup>*( 1 &#8211; e<sup>-(PI\/\u03c3)<sup>2<\/sup><\/sup>)\/PI<\/p>\n<p>where PI is \u221a( Q<sup>2<\/sup> + U<sup>2 <\/sup>). The PI<sub>MAS <\/sub>value is defined for all non-zero PI. The following figure shows PI<sub>AS <\/sub>\u00a0(blue) and PI<sub>MAS<\/sub> (green) plotted against uncorrected PI, assuming \u03c3<sup>2\u00a0<\/sup>= 1:<\/p>\n<p><a href=\"http:\/\/www.eaobservatory.org\/jcmt\/wp-content\/uploads\/sites\/2\/2020\/04\/mas.png\"><img loading=\"lazy\" class=\"wp-image-11054 aligncenter\" src=\"http:\/\/www.eaobservatory.org\/jcmt\/wp-content\/uploads\/sites\/2\/2020\/04\/mas-300x224.png\" alt=\"\" width=\"447\" height=\"334\" srcset=\"https:\/\/www.eaobservatory.org\/jcmt\/wp-content\/uploads\/sites\/2\/2020\/04\/mas-300x224.png 300w, https:\/\/www.eaobservatory.org\/jcmt\/wp-content\/uploads\/sites\/2\/2020\/04\/mas-201x150.png 201w, https:\/\/www.eaobservatory.org\/jcmt\/wp-content\/uploads\/sites\/2\/2020\/04\/mas-150x112.png 150w, https:\/\/www.eaobservatory.org\/jcmt\/wp-content\/uploads\/sites\/2\/2020\/04\/mas.png 535w\" sizes=\"(max-width: 447px) 100vw, 447px\" \/><\/a><\/p>\n<p>The MAS estimator may now be used in <span style=\"font-family: 'courier new', courier, monospace\">pol2map<\/span> by adding &#8220;<span style=\"font-family: 'courier new', courier, monospace\">debiastype=mas<\/span>&#8221; to the <span style=\"font-family: 'courier new', courier, monospace\">pol2map<\/span> command line. Note, the default de-biasing method is still the asymptotic estimator. Using the MAS estimator may facilitate investigation of PI noise statistics in low SNR areas.<\/p>\n<p>&nbsp;<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Calculation of the polarised intensity (PI) and percentage polarisation (P) values from measurements of I, Q and U usually need to be corrected to take account of the bias towards larger values caused by the presence of noise (particularly noticeable\u00a0 in areas of low signal-to-noise). This bias occurs because of\u2026 <a class=\"continue-reading-link\" href=\"https:\/\/www.eaobservatory.org\/jcmt\/2020\/04\/new-de-biasing-method-for-pol2-data\/\">Continue reading<\/a><\/p>\n","protected":false},"author":13,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":[],"categories":[36],"tags":[37],"_links":{"self":[{"href":"https:\/\/www.eaobservatory.org\/jcmt\/wp-json\/wp\/v2\/posts\/11051"}],"collection":[{"href":"https:\/\/www.eaobservatory.org\/jcmt\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.eaobservatory.org\/jcmt\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.eaobservatory.org\/jcmt\/wp-json\/wp\/v2\/users\/13"}],"replies":[{"embeddable":true,"href":"https:\/\/www.eaobservatory.org\/jcmt\/wp-json\/wp\/v2\/comments?post=11051"}],"version-history":[{"count":4,"href":"https:\/\/www.eaobservatory.org\/jcmt\/wp-json\/wp\/v2\/posts\/11051\/revisions"}],"predecessor-version":[{"id":11056,"href":"https:\/\/www.eaobservatory.org\/jcmt\/wp-json\/wp\/v2\/posts\/11051\/revisions\/11056"}],"wp:attachment":[{"href":"https:\/\/www.eaobservatory.org\/jcmt\/wp-json\/wp\/v2\/media?parent=11051"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.eaobservatory.org\/jcmt\/wp-json\/wp\/v2\/categories?post=11051"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.eaobservatory.org\/jcmt\/wp-json\/wp\/v2\/tags?post=11051"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}