<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Phys.</journal-id>
<journal-title>Frontiers in Physics</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Phys.</abbrev-journal-title>
<issn pub-type="epub">2296-424X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1091848</article-id>
<article-id pub-id-type="doi">10.3389/fphy.2023.1091848</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Physics</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Polarized Shack-Hartmann wavefront sensor</article-title>
<alt-title alt-title-type="left-running-head">Yang et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fphy.2023.1091848">10.3389/fphy.2023.1091848</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Yang</surname>
<given-names>Yanrong</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2088292/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Huang</surname>
<given-names>Linhai</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Xiao</surname>
<given-names>Yawei</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Gu</surname>
<given-names>Naiting</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="aff" rid="aff5">
<sup>5</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>College of Ophthalmology</institution>, <institution>Chengdu University of Traditional Chinese Medicine</institution>, <addr-line>Chengdu</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Key Laboratory of Sichuan Province Ophthalmopathy Prevention and Cure and Visual Function Protection with TCM</institution>, <addr-line>Chengdu</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Ineye Hospital of Chengdu University of TCM</institution>, <addr-line>Chengdu</addr-line>, <country>China</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>The Key laboratory of Adaptive Optics</institution>, <institution>Chinese Academy of Sciences</institution>, <addr-line>Chengdu</addr-line>, <country>China</country>
</aff>
<aff id="aff5">
<sup>5</sup>
<institution>University of Chinese Academy of Science</institution>, <addr-line>Beijing</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/782153/overview">Yudong Cui</ext-link>, Zhejiang University, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1684842/overview">Rakesh Kumar Singh</ext-link>, Indian Institute of Technology (BHU), India</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2167349/overview">Maham Aftab</ext-link>, University of Arizona, United States</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2187722/overview">Jian Bai</ext-link>, Zhejiang University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Naiting Gu, <email>gnt7328@163.com</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Optics and Photonics, a section of the journal Frontiers in Physics</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>22</day>
<month>02</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1091848</elocation-id>
<history>
<date date-type="received">
<day>07</day>
<month>11</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>08</day>
<month>02</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Yang, Huang, Xiao and Gu.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Yang, Huang, Xiao and Gu</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Shack-Hartmann wavefront sensor (SHWFS) has been widely used in adaptive optics (AO) systems to detect phase distortion characteristics. In laser communication, target detection, vision optics and other application fields, the performance of SHWFS is affected by bright skylight and scattered light, which restricted the working ability of the AO system severely. Therefore, this paper proposes a new polarized SHWFS (p-SHWFS) based on the principle of polarization imaging, which utilize the difference of state of polarization (SoP) between signal and stray light to improve the image contrast. The p-SHWFS can be composed simply by a micro-lens array and a linear polarization camera. The camera uses four directional polarizing filters at 0&#xb0;, 135&#xb0;, 45&#xb0; and 90&#xb0; on every four pixels. Thus, the degree of linear polarization (DoLP) and angle of linear polarization (AoLP) for the incident wavefront can be analyzed, and the signal-to-background ratio (SBR) can be improved in some certain depending on the difference of SoP. In this paper, we introduce the basic principle of the p-SHWFS and validate the feasibility and accuracy improvement by numerical simulation and practical experiments. The experimental results show that the p-SHWFS can improve obviously the measurement accuracy under strong stray light when the difference of SoP exists. That may give us some initial reference to reduce the influence of stray light in laser communication, target detection, vision optics and other application fields.</p>
</abstract>
<kwd-group>
<kwd>adaptive optics</kwd>
<kwd>Shack-Hartmann wavefront sensor (SHWFS)</kwd>
<kwd>polarization detection method</kwd>
<kwd>polarized dimension</kwd>
<kwd>polarization imaging</kwd>
</kwd-group>
<contract-num rid="cn001">12022308 61905252 12073031</contract-num>
<contract-sponsor id="cn001">National Natural Science Foundation of China<named-content content-type="fundref-id">10.13039/501100001809</named-content>
</contract-sponsor>
<contract-sponsor id="cn002">Youth Innovation Promotion Association of the Chinese Academy of Sciences<named-content content-type="fundref-id">10.13039/501100004739</named-content>
</contract-sponsor>
<contract-sponsor id="cn003">National Key Research and Development Program of China<named-content content-type="fundref-id">10.13039/501100012166</named-content>
</contract-sponsor>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Adaptive optics (AO) technology can effectively correct the dynamic wavefront distortion in real time [<xref ref-type="bibr" rid="B1">1</xref>]. It has been widely used in astronomical observation [<xref ref-type="bibr" rid="B2">2</xref>], laser beam purification, retinal high-resolution imaging, laser communication and other fields [<xref ref-type="bibr" rid="B3">3</xref>&#x2013;<xref ref-type="bibr" rid="B5">5</xref>]. However, all AO systems applied in these fields would be affected by a bright light from background, backscatter or straylight. Therefore, in order to avoid the effect of the daytime stray light in astronomical observation, so far, almost of astronomical observation are working at night or the other low background or stray light scene. Similarly, in the field of laser communication, the stray light generated from high-power laser transmission will be received by the high-sensitivity detector. During the measurement of ocular aberrations, the signal light scattered from fundus will be mixed with corneal reflex. The bright stray light or corneal reflex which would affect the sensitivity of AO system. Shack-Hartmann wavefront sensor (SHWFS) [<xref ref-type="bibr" rid="B6">6</xref>, <xref ref-type="bibr" rid="B7">7</xref>] is one of the most important components of an AO system, which is used to detect the wavefront distortion from the target to observation [<xref ref-type="bibr" rid="B8">8</xref>]. Many applications show that the SHWFS can meet a high accuracy and achieve a stable property under the night or the other low background or stray light scene. However, under the bright background light scenes such as the above conditions of daytime stray light and corneal reflex, the signal-to-background ratio (SBR) [<xref ref-type="bibr" rid="B9">9</xref>] or signal-to-noise ratio (SNR) of each sub-aperture of SHWFS will be limited obviously, where SBR usually used to indicate the intensity level of the background light. Under ultra-worse condition, one may not find any signal light from the SHWFS, which limit its measurement ability of wavefront aberration. For solving this problem, some authors have proposed several kinds of methods. Beckers [<xref ref-type="bibr" rid="B10">10</xref>] applied the active wavefront detection technology based on ultra-narrow band filtering principle to daytime AO system. Gonglewski [<xref ref-type="bibr" rid="B11">11</xref>] used the field stop to eliminate the skylight, but this could not effectively solve the impact of bright skylight on wavefront detection in daytime. [<xref ref-type="bibr" rid="B12">12</xref>] proposed the method of field shift SHWFS, which could detect the mixed information and skylight separately through a tilt mirror or beam splitting prism. However, this field shift method requires a time-sharing detection or two-way wavefront detection set-up, which increased the complexity of control. On the anther hand, the method can only eliminate the evenly distributed skylight. [<xref ref-type="bibr" rid="B13">13</xref>] proposed using polarization filtering technology to suppress polarized background light, but the polarization of skylight is relatively weak, and the filtering method also weaken the target signal. Therefore, this polarization filtering technology only have a limited improvement in detection contrast of the target signal. Huang et al. [<xref ref-type="bibr" rid="B14">14</xref>] proposed a new daytime pyramid wavefront sensor (PyWFS) to separate the object signal from background signal to improve the SNR in daylight AO system. The results shown that the new PyWFS approach was practical when SNR is greater than 1.9 but less than 16.9. In fact, these proposed methods have limited improvement on the performance of SHWFS for daytime AO system in terms of hardware. In addition, there are some improvements come from the algorithm of data processing. Li et al. [<xref ref-type="bibr" rid="B15">15</xref>] proposed a SHWFS-Neural Network (SHNN) to calculate centroid in SHWFS from the extreme situations such as strong environmental light and noise pollution. The results shown that a False Rate with 50 hidden layer neurons was 6% under the peak SNR of 3. Xu et. al [<xref ref-type="bibr" rid="B16">16</xref>] propose a Gaussian modeling centroid extraction algorithm which performs real-time daylight AO closed-loop corrections. This method can provide a lower centroid estimation error under bright sky background conditions. Guo et al. [<xref ref-type="bibr" rid="B17">17</xref>] proposed a new wavefront detection method based on linear polarized modulation. The proposed method can reduce the influence of bright skylight, but the rotating polarizer limits the performance of real-time and accuracy. The mentioned methods can improve in some certain the SBR or SNR of the SHWFS, but the performance of the SHWFS cannot meet the requirement under bright stray light.</p>
<p>Based on above methods, we proposed a polarized SHWFS (p-SHWFS) in this paper, which utilize the difference of state of polarization (SoP) between signal and stray light to improve the image contrast. The p-SHWFS can be composed simply by a micro-lens array and a linear polarization camera. The camera uses four directional polarizing filters at 0&#xb0;, 135&#xb0;, 45&#xb0; and 90&#xb0; on every four pixels. Light passing through these four filters can then be linear interpolated to provide a signal pixel value in polarized dimension. Then, the unpolarized and polarized light can be identified and separated. Thus, the degree of linear polarization (DoLP) and angle of linear polarization (AoLP) for each sub-aperture can be analyzed, and the SBR can be improved in some certain depending on the difference of SoP. In this paper, the basic principle of the proposed p-SHWFS is described, and its performance is also validated by numerical simulation and practical experiments. The experimental results show that the p-SHWFS can improve obviously the measurement accuracy under strong stray light when the difference of SoP exists.</p>
</sec>
<sec id="s2">
<title>2 Principle and methods</title>
<sec id="s2-1">
<title>2.1 Principle of p-SHWFS</title>
<p>The p-SHWFS is composed of a lenslet array, a micro linear polarizer array and a pixeled detector, as shown in <xref ref-type="fig" rid="F1">Figure 1</xref>. The wavefront <italic>W</italic> (<italic>x</italic>, <italic>y</italic>) in the circular domain is divided into <italic>N</italic> sub wavefronts by the lenslet array. Each sub wavefront is modulated by the micro linear polarizer array and imaged on the detector. The micro linear polarizer array is composed of linear polarizing filters placed on every four adjacent pixels of the detector at the directions of 90&#xb0;, 45&#xb0;, 135&#xb0; and 0&#xb0;, which is called a polarized dimension calculation unit. Light passing through these four filters can then be linearly interpolated to provide a single intensity pixel value with its associated DoLP and AoLP in polarized dimension.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>P-SHWFS structure and polarization dimension calculation unit [<xref ref-type="bibr" rid="B18">18</xref>].</p>
</caption>
<graphic xlink:href="fphy-11-1091848-g001.tif"/>
</fig>
<p>According to the actual application scenario, <italic>W</italic> (<italic>x</italic>, <italic>y</italic>) is a mixture of object and stray light. In the field of astronomical observation, the object light from the Sun or distant stars will become unpolarized, the stray light from sky background will become partially linear polarized. Similarly, in the application of human eye aberrations measurements, there will also exist differences of SoP between fundus scattered lights and corneal reflections. Fortunately, this stray light or corneal reflections can be identified by the polarization signature it carries. Therefore, the difference of SoP in the mixture can be used to eliminate or suppress stray light. In order to illustrate the principle of p-SHWFS, we assume that the incident intensity of the <italic>k</italic>-th calculation unit is <inline-formula id="inf1">
<mml:math id="m1">
<mml:mrow>
<mml:msup>
<mml:mi>I</mml:mi>
<mml:mi>k</mml:mi>
</mml:msup>
</mml:mrow>
</mml:math>
</inline-formula>, the object light <inline-formula id="inf2">
<mml:math id="m2">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is modulated by a linear polarizer (the calculation unit can only be used to process linear polarization), and the stray light is <inline-formula id="inf3">
<mml:math id="m3">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> or <italic>vice versa</italic>. Thus,<disp-formula id="e1">
<mml:math id="m4">
<mml:mrow>
<mml:msup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msup>
<mml:mo>&#x3d;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">P</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
<mml:mo>&#x2b;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>The mix intensity of every pixel can be captured simultaneously by the p-SHWFS through one exposure. The intensity of the four adjacent pixels in the <italic>k</italic>-th unit of the detector are <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mn>0</mml:mn>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf5">
<mml:math id="m6">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mn>135</mml:mn>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf6">
<mml:math id="m7">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mn>45</mml:mn>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf7">
<mml:math id="m8">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mn>90</mml:mn>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>,<disp-formula id="e2">
<mml:math id="m9">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3b7;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">P</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
<mml:msup>
<mml:mi mathvariant="bold-italic">cos</mml:mi>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msup>
<mml:mo>&#x2061;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3b8;</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3b7;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">P</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
<mml:msup>
<mml:mi mathvariant="bold-italic">cos</mml:mi>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3b8;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3b7;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">P</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
<mml:msup>
<mml:mi mathvariant="bold-italic">cos</mml:mi>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3b8;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3b7;</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">P</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
<mml:msup>
<mml:mi mathvariant="bold-italic">cos</mml:mi>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi mathvariant="bold-italic">&#x3b8;</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mi mathvariant="bold-italic">k</mml:mi>
</mml:msubsup>
</mml:mrow>
<mml:mrow>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>where &#x3b8; is the polarized direction of <inline-formula id="inf8">
<mml:math id="m10">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, &#x3b7; is the quantum efficiency of the detector. <inline-formula id="inf9">
<mml:math id="m11">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> may be slightly changed as the actual positions of the four pixels do not coincide and the exist aberrations. These coincides will reduce the ability of p-SHWFS in wavefront measurement. But such a high resolution in the unit makes the assumption of the same <inline-formula id="inf10">
<mml:math id="m12">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mi>P</mml:mi>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>. Through Eq. <xref ref-type="disp-formula" rid="e2">2</xref>, the intensity difference between polarized components of <inline-formula id="inf11">
<mml:math id="m13">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mn>0</mml:mn>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf12">
<mml:math id="m14">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mn>135</mml:mn>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf13">
<mml:math id="m15">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mn>45</mml:mn>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula>, <inline-formula id="inf14">
<mml:math id="m16">
<mml:mrow>
<mml:msubsup>
<mml:mi>I</mml:mi>
<mml:mn>90</mml:mn>
<mml:mi>k</mml:mi>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> can be used to calculate the four-dimensional Stokes vector. Then the stray light can be removed with the difference of SoP.</p>
</sec>
<sec id="s2-2">
<title>2.2 Data processing method of p-SHWFS</title>
<p>As shown in <xref ref-type="fig" rid="F1">Figure 1</xref>, each pixel value of the polarized image can be calculated from four pixels calculation unit. To simplify, it is assumed that the <italic>j</italic>-th sub aperture on the detector contains a 2 &#xd7; 2 units. The sub aperture intensity on the detector can be represented by matrix <italic>D</italic>
<sub>j</sub>,<disp-formula id="e3">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">D</mml:mi>
<mml:mi mathvariant="bold-italic">j</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
<mml:mtd>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>where each rectangle is the calculation unit indicated with the superscripts 1&#x2013;4. Then the intensities of <italic>I</italic>
<sub>0</sub>, <italic>I</italic>
<sub>135</sub>, <italic>I</italic>
<sub>45</sub> and <italic>I</italic>
<sub>90</sub> at the polarization directions of 0&#xb0;, 135&#xb0;, 45&#xb0; and 90&#xb0; can be obtained as,<disp-formula id="e4">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>,</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
</mml:mtd>
<mml:mtd>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>Substitute the difference between polarized components of 0&#xb0; and 90&#xb0; directions with <inline-formula id="inf15">
<mml:math id="m19">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mo>&#x22a5;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>. Also, <inline-formula id="inf16">
<mml:math id="m20">
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mo>&#x2225;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the difference between 135&#xb0; and 45&#xb0; directions. Which can be expressed as,<disp-formula id="e5">
<mml:math id="m21">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mo>&#x22a5;</mml:mo>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:mo>&#x394;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mo>&#x2225;</mml:mo>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mfenced open="|" close="|" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">3</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
<mml:mtd>
<mml:mrow>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
<mml:mo>&#x2212;</mml:mo>
<mml:msubsup>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
<mml:mn mathvariant="bold">4</mml:mn>
</mml:msubsup>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>Thus, the DoLP is expressed as <italic>I</italic>
<sub>S</sub>,<disp-formula id="e6">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">S</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mo>&#x22a5;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msup>
<mml:mo>&#x2b;</mml:mo>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mo>&#x2225;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
</p>
<p>Since the p-SHWFS employs linear polarizers, it can be used to calculate the DoLP and remove the unpolarized components. The AoLP replaced with <italic>I</italic>
<sub>A</sub> can be expressed as,<disp-formula id="e7">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">A</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mrow>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:mrow>
</mml:mfrac>
<mml:mi mathvariant="bold-italic">arctan</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mo>&#x2225;</mml:mo>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="bold-italic">&#x394;</mml:mi>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mo>&#x22a5;</mml:mo>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>
</p>
<p>In fact, due to the variety of polarization extinction ratio and the quantum efficiency in every pixel, the DoLP errors caused by the polarizer will be exist. However, when calculating the spot centroids, the errors around the spots can be removed by a fixed threshold value. Using <italic>I</italic>
<sub>0</sub>, <italic>I</italic>
<sub>135</sub>, <italic>I</italic>
<sub>45</sub> and <italic>I</italic>
<sub>90</sub>, we can calculate the total intensity <inline-formula id="inf17">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>D</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>, which the result will be the same as measured by the traditional COMS detector,<disp-formula id="e8">
<mml:math id="m25">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">D</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">0</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">135</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">45</mml:mn>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mn mathvariant="bold">90</mml:mn>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn mathvariant="bold">2</mml:mn>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>Thus, the DoLP, AoLP, and <italic>I</italic>
<sub>D</sub> will then be one-quarter of the original resolution, or can be interpolated to provide the full resolution similar to traditional bayer pattern display. Combining Eqs <xref ref-type="disp-formula" rid="e2">2</xref>, <xref ref-type="disp-formula" rid="e8">8</xref>, it can be seen that the traditional method can only extract the total intensity of light, and cannot effectively eliminate the mixed strong stray light. But using Eq. <xref ref-type="disp-formula" rid="e6">6</xref>, the unpolarized light can be completely eliminated in the mixed incident wavefront. Thus, the p-SHWFS use the difference of SoP between signal and stray light to improve the imaging contrast in wavefront measurement. Different from the traditional method in wavefront measurement, the p-SHWFS uses DoLP to calculate the spot centroid in the <italic>j</italic>-th sub aperture as,<disp-formula id="e9">
<mml:math id="m26">
<mml:mrow>
<mml:mfenced open="{" close="" separators="|">
<mml:mrow>
<mml:mtable columnalign="center">
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mi mathvariant="bold-italic">j</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">R</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="bold-italic">q</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">R</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mi mathvariant="bold-italic">q</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">S</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">q</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">R</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="bold-italic">q</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">R</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">S</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">q</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
<mml:mtr>
<mml:mtd>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">y</mml:mi>
<mml:mi mathvariant="bold-italic">j</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">R</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="bold-italic">q</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">R</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">y</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mi mathvariant="bold-italic">q</mml:mi>
</mml:mrow>
</mml:msub>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">S</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">q</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">K</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi mathvariant="bold-italic">q</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">1</mml:mn>
</mml:mrow>
<mml:mi mathvariant="bold-italic">K</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">S</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">q</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mtd>
</mml:mtr>
</mml:mtable>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>where, (<italic>x</italic>
<sub>
<italic>j</italic>
</sub>, <italic>y</italic>
<sub>
<italic>j</italic>
</sub>) is the spot centroid of the <italic>j</italic>-th sub wavefront, <inline-formula id="inf18">
<mml:math id="m27">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>S</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>q</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the intensity of DoLP calculated according to Eq. <xref ref-type="disp-formula" rid="e6">6</xref>. (<italic>x</italic>
<sub>
<italic>mq</italic>
</sub>, <italic>y</italic>
<sub>
<italic>mq</italic>
</sub>) is the pixel coordinates of the detector, <italic>R</italic>&#xd7;<italic>R</italic> is the centroid calculation region. Thus, the spot displacement can be calculated through comparison of (<italic>x</italic>
<sub>
<italic>j</italic>
</sub>, <italic>y</italic>
<sub>
<italic>j</italic>
</sub>) and the calibrated spot position. By sampling the coming continuous wavefront with an array of sub apertures, all of these spot centroid displacements can be measured and the whole wavefront can be approximated [<xref ref-type="bibr" rid="B19">19</xref>]. In addition, as long as meeting the Nyquist sampling theorem, the resolution loss will have no impact on the centroid measurement accuracy.</p>
</sec>
</sec>
<sec id="s3">
<title>3 Numerical simulation of p-SHWFS</title>
<p>Parameters of p-SHWFS used in simulation are shown in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Structure parameters of p-SHWFS used in simulation.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Parameters</th>
<th align="center">Value</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">The&#xa0;space&#xa0;of&#xa0;two&#xa0;adjacent&#xa0;lenses</td>
<td align="center">0.36&#xa0;mm</td>
</tr>
<tr>
<td align="center">Microlens focal length</td>
<td align="center">14&#xa0;mm</td>
</tr>
<tr>
<td align="center">Array of sub apertures</td>
<td align="center">19 &#xd7; 19</td>
</tr>
<tr>
<td align="center">Number of effective sub apertures</td>
<td align="center">253</td>
</tr>
<tr>
<td align="center">Calculation cell size (pixel)</td>
<td align="center">6.9&#xa0;&#x3bc;m</td>
</tr>
<tr>
<td align="center">Detector effective area (pixels)</td>
<td align="center">992 &#xd7; 992</td>
</tr>
<tr>
<td align="center">Wave length (&#x3bb;) with beacon light</td>
<td align="center">0.653&#xa0;&#x3bc;m</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>In order to completely eliminated the effect of unpolarized light on wavefront aberrations measurement, the incident light intensity <italic>I</italic> (<italic>x</italic>, <italic>y</italic>) is generated by numerical method as,<disp-formula id="e10">
<mml:math id="m28">
<mml:mrow>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">P</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">N</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>where, (<italic>x</italic>, <italic>y</italic>) is the sampling point at the detector. The object and stray light may be partially light. Thus, <inline-formula id="inf19">
<mml:math id="m29">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>P</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the polarized parts of the total light and the direction is <italic>&#x3b8;</italic>. <inline-formula id="inf20">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is the unpolarized parts of the total light. To quantify the relative intensity of the stray light, the SBR is defined as,<disp-formula id="e11">
<mml:math id="m31">
<mml:mrow>
<mml:mtext mathvariant="bold">SBR</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:mn mathvariant="bold">10</mml:mn>
<mml:mi mathvariant="bold-italic">lg</mml:mi>
<mml:mfrac>
<mml:mrow>
<mml:munder>
<mml:mo>&#x222c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2208;</mml:mo>
<mml:mi mathvariant="bold-italic">D</mml:mi>
</mml:mrow>
</mml:munder>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">N</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:munder>
<mml:mo>&#x222c;</mml:mo>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2208;</mml:mo>
<mml:mi mathvariant="bold-italic">D</mml:mi>
</mml:mrow>
</mml:munder>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mi mathvariant="bold-italic">P</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi mathvariant="bold-italic">x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi mathvariant="bold-italic">y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
</p>
<p>The performance of p-SHWFS is affected by the difference of SoP between the signal and the stray light. Thus, in order to quantify the polarization characteristics of <inline-formula id="inf21">
<mml:math id="m32">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>, the polarization degree <italic>P</italic> of <inline-formula id="inf22">
<mml:math id="m33">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> is defined as,<disp-formula id="e12">
<mml:math id="m34">
<mml:mrow>
<mml:mi mathvariant="bold-italic">P</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mi mathvariant="bold-italic">a</mml:mi>
<mml:mi mathvariant="bold-italic">x</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi mathvariant="bold-italic">I</mml:mi>
<mml:mrow>
<mml:mi mathvariant="bold-italic">N</mml:mi>
<mml:mi mathvariant="bold-italic">m</mml:mi>
<mml:mi mathvariant="bold-italic">i</mml:mi>
<mml:mi mathvariant="bold-italic">n</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>p-SHWFS can remove the <inline-formula id="inf23">
<mml:math id="m35">
<mml:mrow>
<mml:msub>
<mml:mi>I</mml:mi>
<mml:mi>N</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:mi>x</mml:mi>
<mml:mo>,</mml:mo>
<mml:mi>y</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula> completely when the <italic>P</italic> is zero. By simulation, the total light intensity calculated from <italic>I</italic>
<sub>0</sub>, <italic>I</italic>
<sub>135</sub>, <italic>I</italic>
<sub>45</sub> and <italic>I</italic>
<sub>90</sub> are shown in <xref ref-type="fig" rid="F2">Figure 2</xref>. It can be seen that the signal is completely submerged in the background (such as <xref ref-type="fig" rid="F2">Figure 2C</xref>) with the increment of stray light intensity. Thus, it is difficult to measure wavefront aberrations accurately by the traditional method from the total light.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>The light intensity calculated with the p-SHWFS. <bold>(A)</bold> no stray light and <bold>(B)</bold> mixed light with SBR &#x3d; 6.29&#xa0;dB and <bold>(C)</bold> SBR &#x3d; 15.03&#xa0;dB.</p>
</caption>
<graphic xlink:href="fphy-11-1091848-g002.tif"/>
</fig>
<p>In order to verify the wavefront measurement ability of p-SHWFS, the target wavefront and aberration components (the first 35 Zernike coefficients) are shown in <xref ref-type="fig" rid="F3">Figure 3</xref>.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Simulated target wave <bold>(A)</bold> and its Zernike coefficient components <bold>(B)</bold>. All the color bar of the wavefront in this work are in units: &#x3bc;m.</p>
</caption>
<graphic xlink:href="fphy-11-1091848-g003.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F4">Figures 4A,B</xref> show the difference of wavefront measuring accuracy by traditional method and p-SHWFS under the polarization degree <italic>P</italic> of stray light is zero. It can be seen that the traditional method cannot accurately measure the target wavefront from the bright stray light. While the p-SHWFS can completely eliminate the stray light, as can be derived from Eq. <xref ref-type="disp-formula" rid="e6">6</xref>. Therefore, under the <italic>P</italic> of stray light is zero, the RMSE (root mean square error) of Zernike coefficient restored by p-SHWFS is approximately zero (as in <xref ref-type="fig" rid="F4">Figure 4A</xref>, the value of blue line is about 10<sup>&#x2013;14</sup>&#xa0;&#x3bc;m). The RMSE and PV (peak valley) values of the residual wavefront are also zeros (<xref ref-type="fig" rid="F4">Figure 4B</xref>, the values marked with &#x201c;&#x394;&#x201d; and hexagon are approximately 10<sup>&#x2013;13</sup>&#xa0;&#x3bc;m). However, the ability of p-SHWFS to suppress stray light can be affected by uneven extinction ratio of micro polarizer array, detector quantum noise, readout noise and <italic>I</italic>
<sub>N</sub>(<italic>x</italic>,y). Limited to space, we only analyzed the influence of <italic>I</italic>
<sub>N</sub>(<italic>x</italic>,y) on the measurement accuracy of p-SHWFS.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>The error comparison between traditional method and p-SHWFS. Zernike coefficient error <bold>(A)</bold> and wavefront error <bold>(B)</bold> under the stray light with <italic>p</italic> &#x3d; 0, Zernike coefficient error <bold>(C)</bold> and wavefront error <bold>(D)</bold> under the stray light with <italic>p</italic> &#x3d; 20%.</p>
</caption>
<graphic xlink:href="fphy-11-1091848-g004.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F4">Figures 4C,D</xref> show the wavefront measurement accuracy of the traditional method and p-SHWFS under different SBR when <italic>p</italic> &#x3d; 20%. The results show that the traditional methods cannot effectively eliminate the stray light, which obviously restrict its wavefront measurement ability. However, the p-SHWFS can effectively eliminate the stray light (<italic>p</italic> &#x3d; 20%, SBR &#x2264;12.35&#xa0;dB). RMSE of the residual wavefront is less than 0.0345&#xa0;&#x3bc;m (&#x2264;0.0528 &#x3bb;), PV is less than 0.2327&#xa0;&#x3bc;m (&#x2264;0.3564 &#x3bb;). RMSE of Zernike coefficient &#x2264;0.0251&#xa0;&#x3bc;m (0.0384 &#x3bb;). From the results, we can know that the ability of p-SHWFS to eliminate the stray light will declined just when the SBR &#x2265;13.10&#xa0;dB (<italic>p</italic> &#x3d; 20%).</p>
<p>Thus, the p-SHWFS can take advantage of the difference of SoP between the signal and stray light to effectively eliminate the influence of strong stray light on the wavefront detection accuracy.</p>
</sec>
<sec id="s4">
<title>4 Experimental verification</title>
<p>In order to verify the wavefront detection ability of p-SHWFS, the optical path shown in <xref ref-type="fig" rid="F5">Figure 5</xref> was built. Where LD is the fiber laser with the wavelength &#x3bb;&#x3d; 653&#xa0;nm. The lenslet array of p-SHWFS is 19 &#xd7; 19 (lenslet spacing 0.36&#xa0;mm, the number <italic>N</italic> of effective sub aperture is 253). The resolution of the detector is 2448 &#xd7; 2048 (pixel size 3.45&#xa0;&#x3bc;m, pixel format with polarized mono12, TRI050S-P, LICID Ltd.). The light beam emitted by LD passes through the collimator L and the polarizer Pr. The passing components will become the polarized light <italic>I</italic>
<sub>P</sub>. A strong stray light <italic>I</italic>
<sub>N</sub> from an incandescent lamp is transmitted to the p-SHWFS along the optical path. <italic>I</italic>
<sub>P</sub> and <italic>I</italic>
<sub>N</sub> are similar to the simulation and also can be explained as in Eq. <xref ref-type="disp-formula" rid="e10">10</xref>. After passing through the linear polarizing filters of the p-SHWFS, the DoLP (one-quarter of the original 2448 &#xd7; 2048 pixels) of the incident intensity will only remain the polarized components <italic>I</italic>
<sub>P</sub> and remove the unpolarized <italic>I</italic>
<sub>N</sub>. On the contrary, it is still true for the p-SHWFS to remain the unpolarized components <italic>I</italic>
<sub>N</sub> and remove the polarized components <italic>I</italic>
<sub>P</sub>, as <italic>I</italic>
<sub>N &#x3d;</sub> <italic>I</italic>&#x2014;<italic>I</italic>
<sub>P</sub>.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Schematic diagram of wavefront detection based on p-SHWFS. L: Lens; Pr: Polarizor; LD: Laser diode. I<sub>N</sub>: strong stray light.</p>
</caption>
<graphic xlink:href="fphy-11-1091848-g005.tif"/>
</fig>
<sec id="s4-1">
<title>4.1 Performance measurement of p-SHWFS without stray light</title>
<p>According to the intensities of <italic>I</italic>
<sub>0</sub>, <italic>I</italic>
<sub>135</sub>, <italic>I</italic>
<sub>45</sub> and <italic>I</italic>
<sub>90</sub> at different polarization directions, the total intensity and DoLP can be obtained according to Eqs <xref ref-type="disp-formula" rid="e8">8</xref>, <xref ref-type="disp-formula" rid="e6">6</xref>. The calculated total intensity is shown in <xref ref-type="fig" rid="F6">Figure 6</xref>. The traditional method and p-SHWFS should have the same results under without stray light, as shown in <xref ref-type="fig" rid="F7">Figure 7</xref>.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>The total light intensity obtained by the intensities of <italic>I</italic>
<sub>0</sub>, <italic>I</italic>
<sub>135</sub>, <italic>I</italic>
<sub>45</sub> and <italic>I</italic>
<sub>90</sub> without stray light.</p>
</caption>
<graphic xlink:href="fphy-11-1091848-g006.tif"/>
</fig>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Comparison of wavefront detection results between traditional method <bold>(A)</bold> and p-SHWFS <bold>(B)</bold> under without stray light, <bold>(C)</bold> the wavefront residual error between the two methods.</p>
</caption>
<graphic xlink:href="fphy-11-1091848-g007.tif"/>
</fig>
<p>
<xref ref-type="fig" rid="F7">Figures 7A,B</xref> show the wavefront detected by the traditional method and p-SHWFS, respectively. The residual wavefront error between the two methods as shown in <xref ref-type="fig" rid="F7">Figure 7C</xref>. The results show that the PV and RMSE of residual wavefront are 0.0058&#xa0;&#x3bc;m (0.0089&#x3bb;) and 0.0017&#xa0;&#x3bc;m (0.0026&#x3bb;), respectively.</p>
<p>Taking the wavefront measured by the traditional method as the target (<xref ref-type="fig" rid="F7">Figure 7A</xref>), and the target aberration components are shown in <xref ref-type="fig" rid="F8">Figure 8A</xref>. The error of the first 35 Zernike coefficients between the traditional method and p-SHWFS is shown in <xref ref-type="fig" rid="F8">Figure 8B</xref>. It can be seen that the RMSE is 0.0269&#xa0;&#x3bc;m (0.0412&#xa0;&#x3bb;). The errors are mainly tip-tilt terms (first and second Zernike coefficient).</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>The target aberration components of the first 35 Zernike coefficients restored by traditional method <bold>(A)</bold> and the error between the traditional method and p-SHWFS <bold>(B)</bold>.</p>
</caption>
<graphic xlink:href="fphy-11-1091848-g008.tif"/>
</fig>
<p>Therefore, the wavefront detection results of the p-SHWFS are consistent with the traditional method without stray light. Next, we take <xref ref-type="fig" rid="F7">Figure 7A</xref>, <xref ref-type="fig" rid="F8">Figure 8A</xref> as the target wavefront to study the detection ability of the p-SHWFS under the influence of stray light.</p>
</sec>
<sec id="s4-2">
<title>4.2 Performance measurement of p-SHWFS under strong stray light</title>
<p>The source <italic>I</italic>
<sub>N</sub> (in <xref ref-type="fig" rid="F5">Figure 5</xref>) passes through an attenuator and ground glass to simulate the unpolarized stray light. By changing the intensity of <italic>I</italic>
<sub>N</sub>, different noise levels of <italic>I</italic>
<sub>N1</sub>, <italic>I</italic>
<sub>N2</sub>, <italic>I</italic>
<sub>N3</sub> are obtained with the SBR of 12.51&#xa0;dB, 15.60&#xa0;dB and 17.04&#xa0;dB, respectively. As shown in <xref ref-type="fig" rid="F9">Figures 9A&#x2013;C</xref> are the intensity calculated under the three noise levels. The strong stray light cannot be eliminated by the conventional method based on the total intensity. It would mix or even submerge the signal light. Thus, strong stray light limits the detection accuracy and the ability of traditional SHWFS. For example, <xref ref-type="fig" rid="F9">Figures 9D&#x2013;F</xref> are wavefront measured results based on the total intensity, and the results obviously deviate from the target wavefront. The wavefront measured by p-SHWFS are shown in <xref ref-type="fig" rid="F9">Figures 9G,H</xref>, (I). Obviously, the results are similar to the target wavefront and not affected by the strong stray light.</p>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>
<bold>(A&#x2013;C)</bold> the total intensity calculated by p-SHWFS under the three noise levels, the wavefront measurement results between traditional method <bold>(D&#x2013;F)</bold> and p-SHWFS <bold>(G&#x2013;I)</bold> under different noise levels.</p>
</caption>
<graphic xlink:href="fphy-11-1091848-g009.tif"/>
</fig>
<p>The errors of target wavefront measured by the traditional method and p-SHWFS under different stray light are shown in <xref ref-type="fig" rid="F10">Figure 10</xref>. <xref ref-type="fig" rid="F10">Figure 10A&#x2013;C</xref> are the wavefront residual error of the traditional method, and the PV and RMES of the residual errors are significantly greater than those of the p-SHWFS shown in <xref ref-type="fig" rid="F10">Figure 10D&#x2013;F</xref>. Therefore, under the influence of strong stray light <italic>I</italic>
<sub>N1</sub>, <italic>I</italic>
<sub>N2</sub>, <italic>I</italic>
<sub>N3</sub>, the target wavefront can be accurately measured based on the p-SHWFS. The PV values of the measurement errors are 0.0128&#xa0;&#x3bc;m (0.0196 &#x3bb;), 0.0180&#xa0;&#x3bc;m (0.0276&#xa0;&#x3bb;), 0.0694&#xa0;&#x3bc;m (0.1063&#xa0;&#x3bb;) respectively. RMSE are 0.0031&#xa0;&#x3bc;m (0.0047&#xa0;&#x3bb;), 0.0046&#xa0;&#x3bc;m (0.0070&#xa0;&#x3bb;), 0.0187&#xa0;&#x3bc;m (0.0286&#xa0;&#x3bb;) respectively. It can be seen that the wavefront measurement accuracy of p-SHWFS is the same with the results without stray light (RMSE &#x3d; 0.0017&#xa0;&#x3bc;m). In addition, PV is slightly greater than the result without stray light (0.0058&#xa0;&#x3bc;m). While the RMSE and PV errors of target wavefront measured by the traditional method are obviously large.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>The residual error of wavefront measured by the traditional method <bold>(A&#x2013;C)</bold> and p-SHWFS <bold>(D&#x2013;F)</bold> under different noise levels.</p>
</caption>
<graphic xlink:href="fphy-11-1091848-g010.tif"/>
</fig>
<p>The error of target aberration components restored by the traditional method and p-SHWFS under <italic>I</italic>
<sub>N1</sub>, <italic>I</italic>
<sub>N2</sub>, <italic>I</italic>
<sub>N3</sub> are shown in <xref ref-type="fig" rid="F11">Figure 11</xref>. It can be seen that the RMSE of p-SHWFS shown in <xref ref-type="fig" rid="F11">Figure 11D&#x2013;E</xref> are 0.0269&#xa0;&#x3bc;m, 0.0275&#xa0;&#x3bc;m and 0.0394&#xa0;&#x3bc;m. It is obviously smaller than the result of traditional method <xref ref-type="fig" rid="F11">Figure 11A&#x2013;C</xref>. The restoration error of p-SHWFS is similar to the results without noise (RMSE &#x3d; 0.0269&#xa0;&#x3bc;m). Similarly, the main error sources are tip-tilt terms.</p>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Zernike coefficient error restored by traditional method <bold>(A&#x2013;C)</bold> and polarization detection method <bold>(D&#x2013;F)</bold> under different noise levels.</p>
</caption>
<graphic xlink:href="fphy-11-1091848-g011.tif"/>
</fig>
<p>From the results, we can conclude that the p-SHWFS can effectively suppress the stray light even when the SBR is &#x2265;15.60&#xa0;dB (<italic>I</italic>
<sub>N2</sub>, <italic>I</italic>
<sub>N3</sub>), and the result is better than the simulation (SBR &#x3d; 13.10&#xa0;dB, <italic>p</italic> &#x3d; 20%). This improvement in suppression the stray light can be attributable to the different polarization degree of the stray light in the experimental system where <italic>P</italic> &#x2260; 20%. Therefore, the experimental and simulation results show that the p-SHWFS can take full use of the significant difference of SoP between the stray light and the object light in the polarization dimension. It would eliminate the influence of strong stray light on wavefront measurement and wavefront component restoration. It can help to accurately measure wavefront distortion under extreme conditions such as low SNR.</p>
</sec>
</sec>
<sec id="s5">
<title>5 Conclusions and discussions</title>
<p>The results in our experiment and simulation have shown that the p-SHWFS would be a promising approach to improve the imaging contrast in AO systems to detect phase distortion characteristics. In order to facilitate the analysis, we used a mathematical model as shown in Eq. <xref ref-type="disp-formula" rid="e10">10</xref> to simulate the mixed incident light, where the difference of SoP between the object and the stray light can be modulated by changing the polarization degree of stray light. The simulated results shown that the performance measurement of p-SHWFS can be affected by the polarization degree of the stray light. In the future analysis and experiments, we will further focus on the performance measurement of p-SHWFS under greater polarization degree of partially polarized stray light.</p>
<p>In summary, a new p-SHWFS based on polarization imaging principle is proposed in this paper. The intensity information would be transformed into polarization dimension using micro linear polarizer array at the directions of 0&#xb0;, 135&#xb0;, 45&#xb0; and 90&#xb0;, so as to take advantage of the significant difference of SoP between stray and object light. The incident mixed light of each sub aperture of p-SHWFS is modulated in four directions to eliminate the interference of stray light. It can improve the SBR and wavefront detection ability under extreme conditions. The feasibility and accuracy of this method have been verified by simulation analysis and experimental system. The experiment results show that the RMSE and PV of wavefront measured by p-SHWFS are &#x2264;0.0070&#xa0;&#x3bc;m and &#x2264;0.0694&#xa0;&#x3bc;m under the stray light with SBR are 12.51&#xa0;dB, 15.60&#xa0;dB and 17.04&#xa0;dB. RMSE of Zernike coefficient are less than 0.0394&#xa0;&#x3bc;m. Compared to the traditional method, the ability of wavefront measurement by the p-SHWFS is obviously improved under the strong stray light when the difference of SoP exists. That may give us some initial reference to reduce the influence of background or stray light in laser communication, target detection, vision optics and other application fields.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s7">
<title>Author contributions</title>
<p>Conceptualization, NG and YY; Methodology, YY; Software, YY; Validation, YY, LH, and YX; Writing&#x2014;original draft preparation, YY; Writing&#x2014;review and editing, YY, LH, NG. All authors have read and agreed to the published version of the manuscript.</p>
</sec>
<sec id="s8">
<title>Funding</title>
<p>This research was funded by National Natural Science Foundation of China (NSFC), grant number 12022308, 12293031, 12073031 and 61905252. The National Key R&#x26;D Program of China, grant number 2021YFC2202204 and 2021YFC2202200.</p>
</sec>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<label>1.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guo</surname>
<given-names>YM</given-names>
</name>
<name>
<surname>Zhong</surname>
<given-names>LB</given-names>
</name>
<name>
<surname>Min</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>JY</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>KL</given-names>
</name>
<etal/>
</person-group> <article-title>Adaptive optics based on machine learning: A review</article-title>. <source>Opto-electron Adv</source> (<year>2022</year>) <volume>5</volume>:<fpage>200082</fpage>. <pub-id pub-id-type="doi">10.29026/oea.2022.200082</pub-id>
</citation>
</ref>
<ref id="B2">
<label>2.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kim</surname>
<given-names>D</given-names>
</name>
<name>
<surname>Choi</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Brendel</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Quach</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Esparza</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>H</given-names>
</name>
<etal/>
</person-group> <article-title>Advances in optical engineering for future telescopes</article-title>. <source>Opto-electron Adv</source> (<year>2021</year>) <volume>4</volume>:<fpage>210040</fpage>. <pub-id pub-id-type="doi">10.29026/oea.2021.210040</pub-id>
</citation>
</ref>
<ref id="B3">
<label>3.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Furukawa</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Takaie</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Maeda</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Ohsaki</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Hasegawa</surname>
<given-names>M</given-names>
</name>
</person-group>. <article-title>Development of one-shot aspheric measurement system with a Shack-Hartmann sensor</article-title>. <source>Appl Opt</source> (<year>2016</year>) <volume>55</volume>(<issue>29</issue>):<fpage>8138</fpage>&#x2013;<lpage>44</lpage>. <pub-id pub-id-type="doi">10.1364/ao.55.008138</pub-id>
</citation>
</ref>
<ref id="B4">
<label>4.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname>
<given-names>GM</given-names>
</name>
<name>
<surname>He</surname>
<given-names>YM</given-names>
</name>
<name>
<surname>Liao</surname>
<given-names>Z</given-names>
</name>
</person-group>. <article-title>Radial Hartmann method for measuring large optical system</article-title>. <source>Chin J Lasers</source> (<year>2010</year>) <volume>37</volume>(<issue>3</issue>):<fpage>795</fpage>&#x2013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.3788/cjl20103703.0795</pub-id>
</citation>
</ref>
<ref id="B5">
<label>5.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>He</surname>
<given-names>JC</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>XT</given-names>
</name>
<name>
<surname>Chu</surname>
<given-names>RY</given-names>
</name>
</person-group>. <article-title>A limitation of Hartmann-Shack system in measuring wavefront aberrations for patients received laser refractive surgery</article-title>. <source>PLoS One</source> (<year>2015</year>) <volume>10</volume>(<issue>2</issue>):<fpage>e0117256</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0117256</pub-id>
</citation>
</ref>
<ref id="B6">
<label>6.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Noll</surname>
<given-names>RJ</given-names>
</name>
</person-group>. <article-title>Zernike polynomials and atmospheric turbulence</article-title>. <source>J Opt Soc Am</source> (<year>1976</year>) <volume>66</volume>(<issue>3</issue>):<fpage>207</fpage>. <pub-id pub-id-type="doi">10.1364/josa.66.000207</pub-id>
</citation>
</ref>
<ref id="B7">
<label>7.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dai</surname>
<given-names>FZ</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>YZ</given-names>
</name>
<name>
<surname>Bu</surname>
<given-names>Y</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>XZ</given-names>
</name>
</person-group>. <article-title>Modal wavefront reconstruction based on Zernike polynomials for lateral shearing interferometry</article-title>. <source>Appl Opt</source> (<year>2012</year>) <volume>56</volume>(<issue>1</issue>):<fpage>61</fpage>&#x2013;<lpage>8</lpage>. <pub-id pub-id-type="doi">10.1364/ao.56.000061</pub-id>
</citation>
</ref>
<ref id="B8">
<label>8.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>C</given-names>
</name>
<name>
<surname>Cai</surname>
<given-names>DM</given-names>
</name>
<name>
<surname>Jia</surname>
<given-names>P</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>JX</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>JJ</given-names>
</name>
</person-group>. <article-title>Research on reconstruction of atmospheric turbulence wavefront compressed sensing measurement</article-title>. <source>Opto-Electronic Eng</source> (<year>2018</year>) <volume>45</volume>(<issue>4</issue>):<fpage>170617</fpage>.</citation>
</ref>
<ref id="B9">
<label>9.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Alam</surname>
<given-names>SU</given-names>
</name>
<name>
<surname>Soni</surname>
<given-names>NK</given-names>
</name>
<name>
<surname>Rao</surname>
<given-names>AS</given-names>
</name>
<name>
<surname>He</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>YX</given-names>
</name>
<name>
<surname>Wong</surname>
<given-names>K</given-names>
</name>
</person-group>. <article-title>Two-photon microscopy with enhanced resolution and signal-to-background ratio using hollow Gaussian beam excitation</article-title>. <source>Opt Lett</source> (<year>2022</year>) <volume>47</volume>(<issue>8</issue>):<fpage>2048</fpage>&#x2013;<lpage>51</lpage>. <pub-id pub-id-type="doi">10.1364/ol.454140</pub-id>
</citation>
</ref>
<ref id="B10">
<label>10.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Beckers</surname>
<given-names>JM</given-names>
</name>
<name>
<surname>Cacciani</surname>
<given-names>A</given-names>
</name>
</person-group>. <article-title>Using laser beacons for daytime adaptive optics</article-title>. <source>Exp Astron</source> (<year>2001</year>) <volume>11</volume>(<issue>2</issue>):<fpage>133</fpage>&#x2013;<lpage>43</lpage>. <pub-id pub-id-type="doi">10.1023/a:1011140920850</pub-id>
</citation>
</ref>
<ref id="B11">
<label>11.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gonglewski</surname>
<given-names>JD</given-names>
</name>
<name>
<surname>Highland</surname>
<given-names>RG</given-names>
</name>
<name>
<surname>Dayton</surname>
<given-names>DC</given-names>
</name>
<name>
<surname>Sandven</surname>
<given-names>SS</given-names>
</name>
<name>
<surname>Rogers</surname>
<given-names>SC</given-names>
</name>
<name>
<surname>BroneADONIS</surname>
<given-names>SL</given-names>
</name>
</person-group>. <article-title>Daylight imaging through atmospheric turbulence</article-title>. <source>Proc SPIE</source> (<year>1996</year>) <volume>2827</volume>:<fpage>152</fpage>&#x2013;<lpage>61</lpage>.</citation>
</ref>
<ref id="B12">
<label>12.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>CH</given-names>
</name>
<name>
<surname>Xian</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Rao</surname>
<given-names>CH</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>WH</given-names>
</name>
</person-group>. <article-title>Field-of-view shifted shack-hartmann wavefront sensor for daytime adaptive optics system</article-title>. <source>Opt Lett</source> (<year>2006</year>) <volume>31</volume>(<issue>19</issue>):<fpage>2821</fpage>&#x2013;<lpage>3</lpage>. <pub-id pub-id-type="doi">10.1364/ol.31.002821</pub-id>
</citation>
</ref>
<ref id="B13">
<label>13.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>RJ</given-names>
</name>
<name>
<surname>Xian</surname>
<given-names>H</given-names>
</name>
<name>
<surname>Rao</surname>
<given-names>CH</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>XJ</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>A</given-names>
</name>
</person-group>. <article-title>Study on effect of polarization filter for suppressing sky background light in daytime</article-title>. <source>Acta Opt Sin</source> (<year>2012</year>) <volume>32</volume>(<issue>5</issue>):<fpage>0501003</fpage>. <pub-id pub-id-type="doi">10.3788/aos201232.0501003</pub-id>
</citation>
</ref>
<ref id="B14">
<label>14.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname>
<given-names>LS</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>JL</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>HF</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>HZ</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>KN</given-names>
</name>
</person-group>. <article-title>Visible pyramid wavefront sensing approach for daylight adaptive optics</article-title>. <source>Opt Express</source> (<year>2022</year>) <volume>30</volume>:<fpage>10833</fpage>&#x2013;<lpage>49</lpage>. <pub-id pub-id-type="doi">10.1364/oe.449021</pub-id>
</citation>
</ref>
<ref id="B15">
<label>15.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>ZQ</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>XY</given-names>
</name>
</person-group>. <article-title>Centroid computation for Shack-Hartmann wavefront sensor in extreme situations based on artificial neural networks</article-title>. <source>Opt Express</source> (<year>2018</year>) <volume>26</volume>:<fpage>31675</fpage>&#x2013;<lpage>92</lpage>. <pub-id pub-id-type="doi">10.1364/oe.26.031675</pub-id>
</citation>
</ref>
<ref id="B16">
<label>16.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>L</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>JL</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>KN</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>LQ</given-names>
</name>
</person-group>. <article-title>Application of the Gaussian modeling algorithm to a Shack&#x2013;Hartmann wavefront sensor for daylight adaptive optics</article-title>. <source>Opt Lett</source> (<year>2021</year>) <volume>46</volume>:<fpage>4196</fpage>&#x2013;<lpage>9</lpage>. <pub-id pub-id-type="doi">10.1364/ol.434941</pub-id>
</citation>
</ref>
<ref id="B17">
<label>17.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guo</surname>
<given-names>T</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>B</given-names>
</name>
<name>
<surname>Gu</surname>
<given-names>NT</given-names>
</name>
<name>
<surname>Rao</surname>
<given-names>CH</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>LH</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>D</given-names>
</name>
<etal/>
</person-group> <article-title>Research on polarization Hartmann wavefront detection technology</article-title>. <source>Opto-electron Eng</source> (<year>2021</year>) <volume>48</volume>(<issue>7</issue>):<fpage>210076</fpage>.</citation>
</ref>
<ref id="B18">
<label>18.</label>
<citation citation-type="confproc">
<person-group person-group-type="author">
<name>
<surname>Huynh</surname>
<given-names>TTM</given-names>
</name>
<name>
<surname>Nguyen</surname>
<given-names>TD</given-names>
</name>
<name>
<surname>Vo</surname>
<given-names>M-T</given-names>
</name>
<name>
<surname>Dao</surname>
<given-names>SVT</given-names>
</name>
</person-group>. <article-title>High dynamic range imaging using A 2x2 camera array with polarizing filters</article-title>. In: <conf-name>2019 19th International Symposium on Communications and Information Technologies (ISCIT)</conf-name>; <conf-date>25-27 September 2019</conf-date>; <conf-loc>Ho Chi Minh City, Vietnam</conf-loc> (<year>2019</year>). p. <fpage>183</fpage>&#x2013;<lpage>7</lpage>.</citation>
</ref>
<ref id="B19">
<label>19.</label>
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ares</surname>
<given-names>M</given-names>
</name>
<name>
<surname>Royo</surname>
<given-names>S</given-names>
</name>
</person-group>. <article-title>Comparison of cubic B-spline and Zernike-fitting techniques in complex wavefront reconstruction</article-title>. <source>Appl Opt</source> (<year>2006</year>) <volume>45</volume>(<issue>27</issue>):<fpage>6954</fpage>&#x2013;<lpage>64</lpage>. <pub-id pub-id-type="doi">10.1364/ao.45.006954</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>