Fundamental flow measurement capabilities of optical Doppler and time-of-flight principles

In order to understand the fundamental measurement capabilities of different flow velocity measurement principles based on Mie scattering, a fundamental equation of how to calculate the shot noise limit for a respective signal model is derived. The derivation is based on the well-known rules of uncertainty propagation and yields the Cramér–Rao bound without the necessity to calculate the Fisher information. The derived equation is next applied to compare the shot noise limit for Doppler and time-of-flight principles including laser Doppler anemometry (LDA), planar Doppler velocimetry (PDV), laser-two-focus velocimetry (L2F), particle tracking velocimetry (PTV) and particle image velocimetry (PIV). The comparison is performed for an identical mean laser power, while two cases are studied in detail: measuring on a single seeding particle as well as measuring on multiple seeding particles and averaging. LDA, L2F and PTV/PIV obey a similar shot noise limit. For the case of a measurement on multiple seeding particles, the minimal achievable measurement uncertainty is directly proportional to the absolute value of the measured velocity component and inversely proportional to the spatial resolution. The respective shot noise limit for PDV is almost independent of the measured flow velocity component and the spatial resolution. Since PDV is sensitive with respect to a different flow velocity component depending on the observation direction, a comparison with the other principles is only reasonable to a certain extent. However, all shot noise limits in case of measuring on multiple seeding particles show the expected inverse proportionality to the square root of the total number of detected photons and thus also to the square root of the measurement time. Considering a comparable spatiotemporal resolution, an identical mean light power and typical measurement configurations, the PDV shot noise limit is the largest. As a final result, it is derived that each measurement principle obeys an uncertainty principle between position and the respective component of the wave vector, which is in agreement with Heisenberg’s uncertainty principle. Therefore, a common basis is provided to assess the fundamental measurement capabilities of Doppler and time-of-flight measurement systems on the basis of what is possible within the quantum mechanical constraints.


Motivation
When confronted with a flow measurement task, finding the proper choice from the huge variety of different measurement principles is difficult. Various optical principles exist, which enable contactless and fast flow velocity measurements. Even hybrid optical measurement methods are known and were, for example, proposed by Grosjean et al. (1997), Förster et al. (2000), Wernet (2004) and Willert et al. (2006). However, regarding the capabilities of the different optical measurement methods, the majority of the studies (aside from some exceptions such as Willert et al. 2005;Thacker et al. 2010) deal only with a single measurement approach. Benchmark tests, in particular over more than two different measurement methods, are rare or even missing. Furthermore, a flow measurement system is typically rated or chosen from what is currently available considering technological and economical aspects. In order to understand the fundamental measurement capabilities of the various optical flow measurement principles-even with respect to future technological developments-and to identify the minimal achievable measurement uncertainty, a scientific study about the fundamental limits of measurability is required.

State of the art
Before reviewing the scientific studies regarding a physical lower limit of the measurement uncertainty, the measurement principles are introduced. The present article is focused on optical flow measurements based on laser light and Mie scattering, because these measurements provide the highest signal-to-noise ratio that for instance allows high-speed measurements even with MHz rate and imaging capability (Thurow et al. 2013). The so-called seeding particles are thus added to the fluid (if not naturally present), which are sufficiently small to follow the flow with negligible slip but at the same time sufficiently large to provide Mie scattering (Albrecht et al. 2003). The scattered light intensity for the Mie scattering on seeding particles is about 10 to 15 orders of magnitude larger than what can be achieved for Rayleigh scattering on the fluid molecules (Bohren and Huffman 2004). As a result, the particle velocity is actually measured and, by assuming a negligible slip, is then interpreted as the flow velocity. According to the review article (Fischer 2017a) from 2017, the flow velocity measurement principles using Mie scattering can by categorized into Doppler and time-of-flight principles. For each category, the developed measurement methods are subsequently summarized.
Doppler principles evaluate the frequency change of the scattered light from a particle that is caused by a nonzero flow velocity relative to the measurement system. The Doppler frequency is the light frequency difference between illumination and observation. One technique to measure the Doppler frequency directly is to mix the light from illumination and observation, since the intensity of the mixed light signal oscillates with the Doppler frequency. The latter is thus extracted from the detected intensity of the mixed signal using a frequency analysis. An early but still state-of-the-art measurement method that is based on the direct Doppler frequency measurement is Laser Doppler Anemometry (LDA) (Yeh and Cummins 1964). Originally introduced as a point measurement device, enhanced LDA techniques enabling profile (Czarske 2001) and field measurements (Coupland 2000;Voigt et al. 2008;Meier and Rösgen 2012) are also known today. Another Doppler technique is to separately measure the frequency of the illuminated light and the observed scattered light, and evaluating the frequency difference in a second step. The required light frequency measures are enabled using an optical filter whose spectral transmission has a steep edge that converts changes in light frequency to changes in transmitted light intensity. The latter is measurable with a photodetector or a camera. Respective Doppler measurement methods have been developed with atomic or molecular filters (Meyers 1995;Charrett et al. 2004;Fischer et al. 2007;Müller et al. 2007;Cadel and Lowe 2015) and with interferometric filters such as a Fabry-Pérot interferometer (Jackson and Paul 1970;Büttner et al. 2013), a Michelson interferometer (Smeets and George 1978;Seiler and Oertel 1983;Landolt and Rösgen 2009) and a Mach-Zehnder interferometer (Lu et al. 2009). These methods are mostly denominated as Planar Doppler Velocimetry (PDV) or Doppler Global Velocimetry (DGV), since planar measurements are easily obtained with a light sheet illumination and cameras. However, scanning or multiplexing techniques also enable volumetric PDV measurements (Fischer et al. , 2015. Time-of-flight principles directly apply the definition of velocity as the derivative of the position with respect to the time. Therefore, the change in particle position for a priori determined time steps is evaluated, or vice versa, the elapsed time when the particle moves over a priori-determined spatial distances. In both cases, position as well as time differences need to be measured, but the uncertainty of the a priori-determined quantity is usually negligibly small. Therefore, one distinguishes between time-of-flight methods with position measurements and time measurements. Note that the differential quotient of the position over time is approximated for most of the applications by the difference quotient. The measurement method that is based on two time measurements is Laser-2-Focus Anemometry (L2F) (Tanner 1973), where the particle subsequently crosses two parallel laser beams and the two particle arrival times are obtained from the detected scattered light pulses. The method is established for point measurements, but can in principle be enhanced toward profile and field measurements with an increased number of laser beam pairs and a camera. A more popular field measurement method that directly measures the time differences of predefined spatial distances is Spatial Filter Velocimetry (SFV) (Aizu and Asakura 1987). Instead of a spatially structured illumination, SFV uses a receiving aperture with a spatially structured sensitivity by applying an optical fringe-type absorption filter. Volumetric SFV measurements are also possible (Hosokawa et al. 2013). However, the time-of-flight measurement methods that are based on position measurements are currently the most popular methods for optical flow velocity field measurements. Particles in a laser light sheet are imaged with a camera for two short laser pulses with a predefined pulse separation time. The position change of each single particle or a pattern of multiple particle images is then obtained by the cross-correlation of sub-images. The method of tracking single particles is Particle Tracking Velocimetry (PTV) (Adrian 1991;Maas et al. 1993). Tracking a pattern of multiple particle images in the interrogation window is Particle Imaging Velocimetry (PIV) (Adrian 2005). Several volumetric PTV/PIV measurement methods are known aside from scanning: for instance, by using tomography (Scarano 2013), holography (Hinsch 2002), light field technique (Cenedese et al. 2012;Fahringer et al. 2015) and synthetic apertures (Cierpka et al. 2010;Belden et al. 2010).
As a result, different Doppler and time-of-flight principles for flow velocity field measurements exist, which are in the following represented by LDA as Doppler principle using light mixing, PDV as Doppler principle using light filtering, L2F as time-of-flight principle based on time measurements and PTV/PIV as time-of-flight principle based on position measurements.
Regarding a physical lower limit of the measurement uncertainty, the error propagation from the naturally fluctuating number of scattered photons was studied for each measurement method. The effect of fluctuating number of photons is known as photon shot noise, which is a fundamental, ultimate and inevitable source of uncertainty. Considering non-squeezed monochromatic laser light, the number of scattered photons that occur in a certain time interval follows a Poisson distribution. Therefore, an uncertainty propagation calculation considering Poissonian noise was accomplished for LDA and L2F by Oliver (1980) in 1980, and Lading and Edwards (1993 studied the extraction of the velocity information from LDA and L2F signals with noise in 1993. McKenzie (1996) investigated the PDV measurement uncertainty in 1996 and included an uncertainty propagation calculation for Poissonian noise. Using an information theoretic approach, the Cramér-Rao bound (CRB) was further calculated for all variants of PDV with and without laser frequency modulation by  in 2017. The Cramér-Rao bound is obtained from the inverse of the Fisher information and equals the minimal achievable measurement uncertainty squared for any unbiased estimator. The same information theoretic approach was used to derive the lower uncertainty limit for PTV by Wernet and Pline (1993) in 1993. Note that the PTV solution is also applicable for PIV when the particle images do not overlap. The measurement uncertainty then reduces inversely proportional to the square root of the number of particle images in the interrogation window. Furthermore, a first attempt to compare the photon shot noise limits of LDA, PDV and L2F using a uniform information theoretic approach took place in 2010 (Fischer et al. 2010) and was briefly complemented with PTV in 2017 (Fischer 2017a). As a result, the minimal achievable measurement uncertainty for LDA, L2F and PTV is similar while it is comparably high for PDV because the scattered light and the information obtained is attenuated by the light filtering. However, the relation between the applied Fisher information calculation (approach by Rao 1945 andCramér 1946, see also Casella and Berger 1990) and a classical uncertainty propagation calculation using linear approximation (approach by C. F. Gauß that is included in the international guide to the expression of uncertainty in measurement Joint Committee for Guides in Metrology 2008) was not addressed. In particular, it is worth studying under which circumstances an uncertainty propagation calculation equals the same flow velocity uncertainty limit as the more cumbersome calculation of the Fisher information and the Cramér-Rao bound.
Furthermore, the existing comparative studies do not focus on the dependency of the shot noise limit on the spatiotemporal resolution. This dependency is important to assess the fundamental measurement capabilities.
Since photons follow the natural laws of quantum mechanics, Heisenberg's uncertainty principle (Heisenberg 1927) is another approach to calculate the fundamental limit of measurement uncertainty due to photon shot noise. This physical approach by W. Heisenberg leads to the same results as with Poissonian photon distribution and information theory (Teich and Saleh 1989), which was explicitly proven for LDA and L2F (Fischer 2016). The surprise of this comparative study was that aside from proportionality factors, an equal limit for the flow velocity measurement uncertainty is obtained for L2F and LDA. This seems remarkable because LDA as Doppler principle evaluates the momentum information of the scattered photons while L2F as time-of-flight principle evaluates the position information of the scattered light photons, so from the physical perspective two completely different, complementary measurement principles. A broader respective study of Heisenberg's uncertainty principle that additionally covers PDV and PTV/PIV is an open task. However, the results from Heisenberg's uncertainty principle will be identical with the consideration of a Poissonian photon distribution. The interesting statement of Heisenberg's uncertainty principle is that the product of the position uncertainty and the momentum uncertainty of a quantum mechanical object (such as the photon) cannot be reduced below a certain quantum limit. A reduced position measurement uncertainty is thus only enabled by an increased momentum uncertainty and vice versa. This quantum mechanical perspective initiated a uniform physical understanding and review of the various optical position and shape measurement methods (Pavliček and Häusler 2014;Pavliček and Pech 2016;Fischer 2019), but the aspect of a fundamental uncertainty relation that enables the minimization of the flow velocity measurement uncertainty at the cost of a maximized uncertainty of a complementary quantity has not yet been fully considered in a comparative study of LDA, PDV, L2F and PTV/PIV.

Aim and structure of the article
The first aim of the article is to replace the calculation of the Cramér-Rao bound by the more common and simple rules of uncertainty propagation and to derive the photon shot noise limits for LDA, PDV, L2F and PTV/PIV. The second aim is to compare the derived fundamental measurement limits and to understand the respective dependency on the spatiotemporal resolution. The third aim is to understand the uncertainty relations for the optical flow velocity measurements based on laser Mie scattering, i.e., what uncertainty has to be maximized to minimize the flow velocity measurement uncertainty?
The considered measurement methods, which represent Doppler principles with light mixing and light filtering as well as time-of-flight principles with position and time measurements, respectively, are described in Sect. 2. In the description of the principles, the comparability of the different measurement methods is addressed, i.e., fundamentally different measurement features are emphasized that indicate the uniqueness of each measurement principle. As a result, signal models for LDA, PDV, L2F and PTV/PIV are derived that are applicable to an uncertainty propagation calculation. The uncertainty propagation calculation that yields the minimal achievable measurement uncertainty for photon shot noise is presented in Sect. 3. The simplified consideration of a single measurand as the only unknown quantity allows to apply the uncertainty propagation calculation in order to obtain the square root of the Cramér-Rao bound as the minimal achievable measurement uncertainty. This finding is applied to calculate the shot noise limits of LDA, PDV, L2F and PTV/PIV in Sect. 4. In addition, uncertainty relations are identified in the results and compared between the different measurement methods. The conclusions are drawn in Sect. 5 together with a discussion concerning remaining research questions.
Page 5 of 19 37 2 Measurement principles A physical signal model for each measurement principle is required to calculate the respective shot noise limit. The signal modeling is based on several assumptions. Particle slip is not part of the following investigation, so that the actually measured particle velocity can be assumed to equal the flow velocity. Furthermore, the flow velocity components are considered as constant during the measurement time, because flow acceleration effects on the measurement result are not studied here. Flow turbulence is also excluded here, i.e., all particles that occur in the measurement volume have an identical velocity. For the sake of simplicity, monodisperse particles and a flow velocity vector with no out-of-plane component are further considered. Out-of-plane means out of the measurement plane, which is defined by the illumination. In order to focus on the fundamental physics of each measurement principle, technical aspects regarding the photo detection unit (photo detector or camera) and the signal digitization are not included in the signal modeling. In particular, the discretization of the detected light signal (or image) over time (due to sampling) and over space (due to camera pixels) is not part of the present study. Finally, the signal processing algorithm and its efficiency to precisely estimate the flow velocity are not discussed, because the article is focused on the comparison of different measurement principles using the Cramér-Rao bound. The Cramér-Rao bound for any unbiased estimator is independent from the signal processing algorithm.

LDA
A state-of-the-art differential LDA setup is considered as depicted in Fig. 1. The original laser beam with the wavelength is split into two beams, which are here modeled as symmetric Gaussian beams with the 1∕e 2 -radius 2w beam . Both beams are brought to intersection with the crossing angle 2 . The intersection region is the measurement volume, in which an interference fringe system occurs that is observed with a photodetector. The fringe distance is d fringe = 2 sin . Due to the interference fringe system, the intensity distribution is a cosine-type modulation with a Gaussian envelope. Considering one seeding particle that crosses the fringe system at the arrival time t 0 , the detected scattered light signal over time t in unit number of photons (accumulated over one sampling period t ) reads with The amplitude A LDA is the total number of detected photons, which is here defined by a mean photon rate Ṅ photon,LDA multiplied by the transit time where 4w fringes = 4w beam ∕ cos is the (1∕e 2 )-extension of the measurement volume in x-direction and v x denotes the respective flow velocity component. The temporal width of the Gaussian envelope is thus described by the lateral width of the beam intersection region divided by the flow velocity. The beat frequency of the signal due to the Doppler effect is the Doppler frequency f D . This LDA signal model is shown in Fig. 2.
The signal from Eq. (1) is evaluated for instance using harmonic analysis to estimate the Doppler frequency. Together with the previously determined fringe distance d fringe during calibration, the flow velocity component perpendicular to the fringes is finally obtained by calculating Sensitivity: The basic LDA measurement configuration yields one flow velocity component at a single point. The (1) sensitivity direction is perpendicular to fringes, which is in the same plane as the two laser beams and perpendicular to the optical axis of the illumination. The observation direction has no effect on the sensitivity direction. Spatial resolution: The spatial resolution in all three dimensions is typically determined by the spatial extension of the fringe system (illumination) but might be additionally limited by the optical imaging on the photodetection unit (observation).
Temporal resolution: The temporal resolution for a single particle amounts to T LDA . It depends on the flow velocity and the width of the laser beams as well as the beam crossing angle. Note that varying the laser beam width or the beam crossing angle will affect the temporal as well as the spatial resolution. For a given measurement time T ≫ T LDA that allows to observe multiple particles crossing the fringe system, an averaging over floor(T∕T LDA ) consecutive single particle measurements is considered here as the theoretic maximum.

PDV
The flow is continuously illuminated with a laser light sheet which represents the measurement plane. Illumination intensity variations in the observed region of interest in the measurement plane are neglected here. The laser has a narrow bandwidth and has a stabilized center frequency f c = c∕ according to the laser wavelength and the light speed c. To detect the flow velocity-dependent frequency change of the light scattered on seeding particles, an imaging system including a beam splitter (here 50 %/50 %) and two digital cameras as well as a narrow-band spectral filter is applied, see Fig. 3. Half of the scattered photons from the illuminated seeding particles in the measurement plane are imaged with the reference camera. The other half of the scattered photons pass a spectral filter with the spectral transmission (f c + f D ) , and the transmitted photons are imaged with the signal camera. The symbol f D denotes the Doppler frequency shift. Note that the filter bandwidth is significantly larger than the bandwidth of the continuouswave laser, i.e., the illuminating laser can be considered as an ideal laser source with single frequency. As a result, the signals from the signal camera and the reference camera over time t in unit number of photons (accumulated over one sampling period t ) in the corresponding interrogation window (single pixel or superpixel) read for the presence of N particle,PDV seeding particles with A PDV as the total number of scattered photons per seeding particle during the particle transit time T PDV through the interrogation window. Note that the number of observed seeding particles increases with an increasing size of the interrogation window and an increasing light sheet thickness. The mean photon rate Ṅ photon,PDV per seeding particle is defined by the relation with This PDV signal model is shown in Fig. 4. The symbol w transit denotes the width of the illuminated interrogation region with respect to the PDV sensitivity direction ⃗ , and v denotes the measured velocity component along the PDV sensitivity direction.
The Doppler frequency f D is determined by evaluating the ratio , by subtracting the laser working point (f c ) and by applying the inverse transmission function that was previously calibrated. The flow is then calculated according to the relation  The unit vectors ⃗ i and ⃗ o describe the direction of illumination and observation, respectively.
Note that only a classical PDV without laser frequency modulation is considered here, since it was already shown in ) that all variants of PDV with laser frequency modulation and without laser frequency modulation have similar measurement capabilities with respect to the shot noise limit.
Sensitivity: The PDV principle provides sensitivity with respect to the out-of-plane velocity component.
and thus depends on the illumination direction and the observation direction. For a typical perpendicular arrangement of the light sheet plane and the observation, a mixture of in-plane and out-of-plane velocity component is measured. The in-plane velocity component along the light sheet propagation direction is measurable independent of the observation direction when two consecutive Doppler frequency measurements are performed with opposite illumination directions and the same observation direction (Roehle and Willert 2001). Due to the use of cameras, imaging capability is provided.
Spatial resolution: The out-of-plane spatial resolution is usually determined by the illumination, i.e., the light sheet thickness. This holds in particular for an observation direction perpendicular to the light sheet plane, which is a typical measurement configuration. The in-plane spatial resolutions are then determined by the chosen region of interest in the captured frame, which is either a single pixel or an area of binned pixels. A perpendicular arrangement between the light sheet plane and the observation direction allows the highest spatial resolution.
Temporal resolution: Note that multiple particles are allowed to be present in the interrogation window of the illuminated region and are measured simultaneously. The PDV measurement is possible even for overlapping particle images so that very high seeding particle concentrations can be applied. The temporal resolution is thus typically determined and adjustable by the camera exposure time T. Considering T PDV as the time for a single measurement and T > T PDV , an averaging over (T∕T PDV ) measurements takes place implicitly.

L2F
The flow is illuminated by two parallel laser beams, see Fig. 5. A photodetector is applied to detect the scattered light from a seeding particle that crosses both laser beams in the region of interest. Assuming identical Gaussian beam intensity profiles, the sensor signal over the time t contains two identical Gaussian pulses at time t 1 and t 2 and reads in unit number of photons (accumulated over the sampling period t) with The temporal width of the pulses as well as the energy content of the pulses scale with the inverse of the soughtafter velocity component v x . Therefore, the total number of detected photons A L2F is written as a mean photon rate Ṅ photon,L2F multiplied by the transit time where d beam is the distance of the two laser beams. The temporal pulse half-width is (2w beam ∕|v x |) , where 2w beam denotes the 1∕e 2 -radius of each laser beam. This L2F signal model is shown in Fig. 6. Since the overlap of the pulses is negligible, the sensor signal to be solved for the two parameters t 1 , t 2 can also be written in the form of a combination of two signals i = 1, 2 with only one parameter t i each: In practice, the signal s t (t 1 , t 2 ) is evaluated (for instance by using an auto-correlation) to directly determine the temporal distance (t 2 − t 1 ) of the two scattered light pulses. Together with the previously calibrated beam distance d beam , the flow velocity component v x perpendicular to the direction of the laser beams is then obtained according to the relation Sensitivity: The described L2F measurement principle in its original configuration allows to measure a flow velocity component at a single point. The sensitivity direction is perpendicular to the propagation direction of the illuminating laser beams and lies in the same plane as the two laser beams. The observation direction has no effect on the sensitivity direction.
Spatial resolution: The spatial resolution along the sensitivity direction (x-direction) is mainly determined by the beam distance, while the spatial resolution along the beam propagation (y-direction) is determined by the beam divergence or the field of view for the photodetector with respective optics together with the observation angle. The spatial resolution in z-direction, i.e., perpendicular to the measurement plane which is defined by the two laser beams, follows from the beam diameter.
Temporal resolution: The temporal resolution for a single particle is the transit time T L2F . Thus, it depends on the flow velocity and the beam distance. A change of the beam distance will affect the temporal and the spatial resolution. In case of a given measurement time T that is larger than T L2F , an averaging over floor(T∕T L2F ) consecutive seeding particles is considered here as the theoretic maximal number of measurements, which is analogous to LDA.

PTV/PIV
The flow is illuminated with two sequential light sheet pulses . The illuminated plane is the measurement plane and the illuminated seeding particles are imaged using an imaging system with a camera, see Fig. 7. Note that spatial variations of the illumination intensity can occur (Nobach 2011), in particular when a Gaussian intensity profile is used for illumination. For the sake of simplicity, a constant illumination intensity profile is subsequently assumed in the observation region. Furthermore, the image quantization due to the limited pixel size of the camera is neglected in the present study, because it is considered a technological limitation depending on the available cameras. If the pixel size is near or larger than the size of the particle image and cannot be neglected, the resulting Cramér-Rao bound is higher than in the present study, see (Wernet and Pline 1993) for Poisson noise and (Westerweel 1997(Westerweel , 2000 for Gaussian noise. Assuming a negligibly small pixel size thus leads to the physically lowest possible Cramér-Rao bound and a signal model with a continuous space variable. As a result, each image of a particle in unit number of detected scattered photons (accumulated over the area x ⋅ y of one pixel) as a function of the inplane coordinates x and y can be approximately described by the two-dimensional Gaussian function The position (x i , y i ) is the particle position at the time t i of the i-th light pulse with i = 1, 2 . The 1∕e 2 -radius of the particle image amounts to 2w particle . The particle image size in the units of the original space coordinates is either the true particle size or it is increased due to diffraction. The total number of detected scattered photons per seeding particle is A PTV ∕2 for each of both images and does not depend on the flow velocity. Although A PTV can be adjusted directly by varying the pulse energy of the laser, the definition with a mean photon rate Ṅ photon,PTV over the pulse interval time Fig. 7 PTV/PIV measurement principle illustrated with one (PTV) and for instance three (PIV) seeding particles in the interrogation region. The light sheet is shown in side-view, which is the camera view when the observation direction is chosen perpendicular to the light sheet plane is in principle also possible analogous to a continuous illumination, cf. LDA and L2F.
Note that the task of solving Eq. (14) for the two unknown particle coordinates x i , y i can be transformed in principle to the task of separately solving one equation for each coordinate. By integrating over x or y and multiplying with the respective pixel size x or y , the separate signals in unit number of photons (accumulated over one sample size) that do solely depend on x i or y i , respectively, read The resulting PTV/PIV signal model for the x-coordinate is shown in Fig. 8.
In practice, the two sequential particle images according to Eq. (14) are cross-correlated to directly determine the particle position change (x 2 − x 1 , y 2 − y 1 ) during the interval time T PTV of the two pulses. Assuming a constant flow velocity, the in-plane velocity components then follow from the calculation Having more than one particle image in the interrogation region, which is the case for PIV in contrast to PTV, the measurement procedure remains the same. However, the resulting velocity measurement is an average over all N particle,PIV particles in the interrogation region according to the seeding particle concentration.
Sensitivity: The described standard PTV/PIV measurement arrangement has imaging capability and allows to measure simultaneously the axial velocity component and one lateral velocity component with respect to the illumination. Both sensitivity directions result from the chosen illumination. Spatial resolution: A perpendicular arrangement between the light sheet plane and the camera direction is a preferred measurement configuration, since it allows to achieve a high spatial resolution. The spatial resolution along the out-ofplane direction is then mostly determined by the light sheet thickness (illumination path) and typically not by the depth of focus of the imaging system (observation path). On the contrary, the spatial resolution along the in-plane directions is usually determined by the chosen interrogation region (observation path) and not by the illumination (PIV). However, the minimal achievable spatial resolution is determined by the particle motion between the two light pulses (PTV).
Temporal resolution: The temporal resolution is given by the pulse interval time T PTV and is adjustable by varying the temporal distance between both light pulses. Note, however, that the pulse interval time is also connected with the minimal achievable spatial resolution that is the particle motion. If a measurement time T is possible that is larger than a pulse interval time constraint, which follows from a spatial resolution requirement, a theoretic maximum of floor(T∕T PTV ) measurements can be performed. Furthermore, multiple particles in the interrogation region can be evaluated simultaneously for PIV. Note, however, that the respective averaging over multiple measurements according to the derived signal model reflects only the case of non-overlapping particle images, which is the optimal case.

Comparison
A key difference between the different measurement principles is that PDV measures a mixture of in-plane and out-ofplane velocity component (with respect to the measurement plane defined by the illumination), LDA and L2F measure the in-plane velocity component transverse to optical axis of the illumination and PTV/PIV allows to measure both inplane velocity components. Note that this statement holds when only a single illumination direction is allowed. On the contrary, PDV is the only principle where the sensitivity depends on the observation direction. It seems important to remind that a reference LDA, in contrast to the differential LDA that is considered here, actually has the same capability. The reason is that PDV and LDA are based on the Doppler effect, which is affected by the illumination and the observation direction. With the differential LDA setup, which is nowadays more common and therefore considered here, the influence of the observation direction is eliminated by evaluating the difference of two Doppler effects resulting from two different illumination directions but one common observation direction. Due to the different sensitivity direction, the PDV measurement capabilities are not one-to-one The relevant feature for the velocity measurement with L2F and LDA is a characteristic size of the illumination profile. For L2F, it is the radius of the laser beams, and for LDA, it is the fringe distance. The size of the respective feature is diffraction-limited according to the illumination. The particle size is typically smaller than this feature, but the particle size can limit the feature size if it is larger than the illumination diffraction limit. The opposite is true for PTV/PIV, where the relevant information feature for the velocity measurement is the particle image size, which is a diffraction-limited observation feature. If the particle size is smaller than the diffraction limit, the feature size is limited according to the diffraction limit. For PDV, the relevant information feature is the light wavelength (or the respective mean light frequency), which is a characteristic size in axial direction, i.e., in the light propagation direction. Therefore, both the illumination direction (particle motion with respect to the illumination direction) and the observation direction (particle motion with respect to the observation direction) have an influence on the resulting Doppler effect.
While LDA and L2F optimally operate for a single particle in the interrogation region, PIV and PDV allow measurements with multiple particles in the interrogation region without degrading the sensitivity. In case of PIV, however, the particles should not overlap to not degrade the sensitivity. Therefore, the highest number of simultaneous particle velocity measurements (which means the highest seeding particle concentration) is possible with PDV.
Regarding the spatiotemporal resolution (in the sensitivity direction), the spatial resolution for LDA and L2F is defined by the chosen spatial distribution of the illumination while the temporal resolution then follows according to the seeding particle motion (transit time). For PTV, the temporal illumination characteristic sets the temporal resolution, and the spatial resolution then follows from the respective particle movement. For PDV and PIV, the spatial and the temporal resolutions are defined by the illumination and observation in space and time.

Method of uncertainty propagation that yields the Cramér-Rao bound
A measurement signal s m in unit number of photons is considered that depends on a single measurand and contains M uncorrelated samples: The symbol w m denotes the signal noise. In order to investigate the shot noise limit, each measurement signal sample is assumed to follow a Poissonian distribution, i.e., the variance of the signal sample s m equals the expectation value: Var s m = E s m . The expectation value is subsequently denoted by the symbol s m , i.e., E s m = s m . Note that a Poissonian photon distribution minimizes Heisenberg's uncertainty relation (Saleh and Teich 2007), so that the considered case is a quantum mechanical lower limit. According to the estimation theory of Cramér and Rao, the minimal achievable variance for any unbiased estimator ̂ of the unknown measurand reads for the given signal in Eq. (20) with uncorrelated noise (Casella and Berger 1990). The right side of the equation is the Cramér-Rao bound, which is identical to the minimal achievable measurement uncertainty squared . Note that the Cramér-Rao bound is the information theoretic limit of uncertainty since it follows from the inverse of the Fisher information.
In order to derive Eq. (21) using the uncertainty propagation calculation, the weighted average over the estimations ̂m from each signal sample is considered as the measurement result: with 1 Var(̂m) as the optimal weighting factor according to the Gauß-Markow theorem. Since the noise w m of the measurement signal is uncorrelated, the estimators ̂m are uncorrelated as well. Therefore, the uncertainty propagation calculation at first yields Estimating Var ̂m with a second uncertainty propagation calculation based on Eq. (20), the relation then reads Var(s m) .
Comparing Eqs. (21) and (24), the uncertainty propagation calculation applied to the optimal estimator according to the Gauß-Markow theorem directly leads to the Cramér-Rao bound.
Since the variance of each signal sample is Var s m = E s m = s m according to the Poissonian noise and Eqs. (20), (24) finally becomes as the shot noise limit of the measurand .
For a sufficiently small sampling period, the signal samples s m can be approximated by a continuous signal (for instance over time or space). Considering the continuous signal s t ≈ s m over the variable t = m ⋅ t with the sampling period t → dt , the shot noise limit solution in Eq. (25) is writable in the integral form with the measurement interval T = M ⋅ t . Note that s t denotes the number of detected photons accumulated over one sampling period t , so that integrating s t over t and dividing by t yields the same total number of detected photons as for summing up s m over m. The result is interpretable as the quotient of the noise power spectral density 1 1 t for one photon and an accumulated (optimal weighted) sensitivity. Enhancing denominator and numerator of the quotient with as weighting normalization, the quotient can also be considered as an effective noise power spectral density divided by the average sensitivity with an optimal weighting. Note that the shot noise limit in Eqs. (25) or (26), respectively, is valid when the measurand is the only unknown quantity. Additional unknown quantities can increase the minimal achievable measurement uncertainty, and the calculation of the Cramér-Rao lower bound needs to be enhanced using covariances and matrix calculus (Casella and Berger 1990;Kay 1993). The increased calculation effort is not necessary here, since all enhanced solutions remain bounded by the shot noise limit derived in Eqs. (25) and (26). Both formulas thus represent the sought-after ultimate lower bound of the minimal achievable measurement uncertainty. (25) s t dt 4 Photon shot noise limit results The shot noise limit for each measurement principle is calculated according to Eq. (26) and the respective signal model derived in Sect. 2. Here, the integral is evaluated from minus infinity to plus infinity, which implicitly implies a sufficiently large measurement time, in particular for LDA and L2F. The following discussions are thus always focused on the case of a measurement time that is equal or larger than the particle transit time. The calculations for LDA, L2F and PTV/PIV make use of the integral solution ∫ ∞ −∞ a ⋅ exp(−z 2 ∕(2b 2 )) ⋅ z 2 b 2 dz = √ 2 ab and all calculations are straightforward.
For LDA, the signal model from Eq. (1) is applied and the resulting variance of the Doppler frequency is converted to the corresponding variance of the measured flow velocity component by using Eq. (4). Note that a fringe spacing is assumed that is significantly smaller than the width of the Gaussian envelope. The assumption is plausible, and it is further valid for heterodyne LDA setups (Czarske 2006), where the cosine signal has an offset frequency. An enhanced study for LDA without this assumption is presented in Sobolev and Feshenko (2006), but note that Gaussian noise is considered there. The resulting Cramér-Rao bound then only differs by a different pre-factor and a marginal deviation below seven percent.
For PDV, the signal model from Eq. (5) is applied, where both signals are treated as one concatenated signal with double length. The resulting variance of the Doppler frequency is then used together with Eq. (8) to determine the corresponding variance of the measurement flow velocity component. For L2F, the shot noise limit for both pulse occurrence times is determined at first using the signal model from Eq. (12). The resulting shot noise limit of the measurement flow velocity component is then obtained by an uncertainty propagation calculation of the uncorrelated estimations of the two occurrence times on the basis of Eq. (13). For PTV/ PIV, the signal model from Eq. (17) is used, and the calculated variance of the particle position is applied to Eq. (19) to yield the shot noise limit for the measured flow velocity component. As a result, the calculated shot noise limits read in summary Var �v where the relations from Eqs. (3), (11) and (16) were applied to demonstrate the similar structure of the shot noise limits. It seems remarkable that all identified shot noise limits in Eq. (27) have an identical structure although the measurement principles are fundamentally different. The common inverse proportionality to the total number of detected photons A . is no surprise, because this is a well-known characteristic feature for Poisson noise. Note that the expression N ⋅ ⋅ A PDV for PDV is the total number of detected photons for the signal camera. The numerical pre-factors are one for L2F, PDV and PTV/PIV. The pre-factor 4 2 = 0.41 for LDA is also near one, so that no significant difference between the different measurement principles exists with respect to the pre-factor. The chosen definition of the particle transit time with respect to the width of the Gaussian envelope (here: the time T LDA that corresponds to the full 1∕e 2 -width of the Gaussian envelope) and the definition of a characteristic length (here: the fringe distance d fringe ) are finally responsible for the resulting value of the pre-factor.
All shot noise limits are directly proportional to the square of a characteristic length scale, but the characteristic length scale is different for each measurement principle. It is the fringe distance d fringe for LDA, the laser wavelength scaled with the inverse of the normalized sensitivity vector length �⃗ o − ⃗ i�∕ √ 2 (normalization is with respect to the typical measurement configuration ⃗ o ⟂ ⃗ i ) for PDV, the beam radius 2w beam for L2F and the particle image radius 2w particle for PTV/PIV. Furthermore, the shot noise limits are inversely proportional to the square of a characteristic time scale. For LDA, L2F and PTV/PIV, this time scale is the transit time of the particle to pass the measurement volume, whereas for PDV, it is the derivative � = d (f ) df of the spectral filter transmission divided by the spectral transmission both evaluated at the mean frequency of the scattered laser light.

Comparison for measuring a single particle
To determine and compare the shot noise limits for the flow velocity measurement with a single seeding particle, the following common conditions are applied to Eq. (27) for the sake of comparability: • A single particle is measured, i.e., N particle,PDV = 1. • The observation direction is perpendicular to the illumination direction, so that the scattering angle is the same for each measurement principle. Note that for PDV, this condition ⃗ o ⟂ ⃗ i means �⃗ o − ⃗ i� = √ 2. • All light sources provide the same power on average in the measurement region and the working distance for illumination and observation, respectively, and the numerical aperture is equal. As a result, the mean detected photon rate during the measurement time is the same for each measurement principle. Therefore, the total number of detected photons A . are expressed by using Eqs.
(2), (6), (10) and (15). • A comparable spatial resolution (in all three dimensions) and a comparable temporal resolution are considered. However, the temporal resolution is inversely proportional to the flow velocity for a given spatial resolution. To describe this behavior, the particle transit times T . are expressed by Eqs.
(3), (7), (11) and (16). Note that this implies an infinitely large measurement time if the measured velocity component becomes zero. Since this particular case is not meaningful for a practical measurement, the subsequent discussions are focused on a nonzero flow velocity component to be measured.
The shot noise limits for LDA, PDV, L2F and PTV/PIV then read Note that for LDA, the ratio 4 w fringes ∕d fringe is the number of fringes. The LDA shot noise limit is thus inversely proportional to the squared number of fringes, which was also found in earlier studies regarding Gaussian noise (Sobolev and Feshenko 2006). As a result, each shot noise limit of the measurement uncertainty increases with the flow velocity to be measured. For LDA, L2F and PTV/PIV, the variance lower bound is directly proportional to |v x | 3 , i.e., to the measurand to the power of three, whereas for PDV, it is only directly proportional to |v| if the slope ′ of the spectral filter transmission and the transmission value can be assumed as constant. The condition � = const is indeed a plausible approximation for a large range of the Doppler frequency (as long as the scattered light mean frequency remains at the edge of the spectral filter), while = const holds for PDV setups where the laser frequency is adjusted accordingly.
The direct proportionality to the square of the respective characteristic length scale d fringe , , 2w beam and 2w particle is unchanged, while the inverse proportionality to the square of the characteristic time scale only remains for PDV, when the dependency of ′ and with respect to the flow velocity is neglected. The respective spatial resolution now occurs in the denominator to the power of three for LDA, L2F, PTV/PIV and to the power of one for PDV. In the following comparison, an equal spatial resolution w = 4 w fringes = w transit = d beam = x 2 − x 1 is considered. Note that the different measurement principles are compared for a similar spatial resolution regarding each dimension. Having a radially symmetrical beam profile for LDA, for instance, this means to consider L2F beams that are not circular but elliptical. The beam dimension perpendicular to the sensitivity direction is then larger than the beam dimension in the sensitivity direction. Although such an L2F setup is not a typical setup today, it can be realized using cylinder lenses and is a fair way to compare the potential of the L2F measurement principle with the LDA measurement principle for one common spatial resolution (and a comparable data rate). Without this consideration of a fully comparable spatial resolution, L2F is known to enable lower measurement uncertainties than LDA due to the higher illumination intensity for the case of a comparable illumination light power (Beversdorff et al. 1997).
All shot noise limits are inversely proportional to the mean photon rate, and an equal mean photon rate Ṅ photon is considered in the following comparison. Note that the mean photon rate equality is below the PTV capabilities, because the laser pulse energy can be accumulated over a larger period than the actual measurement time. Thus, the pulsed laser at PTV/PIV in principle allows to achieve a higher mean photon rate over the measurement time than the continuous-wave laser at LDA and L2F. This also holds for PDV when a pulsed laser is used for illumination. Although below theoretical capabilities, the subsequent comparison with an equal photon rate in combination with an equal temporal resolution (measurement time) means a comparison with an equal number of detected photons on average, which is of high interest. The question is: Which measurement principle uses the amount of energy in the most effective way to minimize the achievable measurement uncertainty?
To answer this question, Eq. (28) is evaluated for the parameter values listed in Table 1. A spatial resolution of w = 100 μm is considered as an example, so that the temporal resolution depends on the flow velocity and is w∕|v x | or w/|v|, respectively. The mean photon rate is set to Ṅ photon = 1×10 9 s −1 , which corresponds for instance to a mean power of 373 pW when the laser wavelength amounts to 532 nm. The mean power value is plausible since it is the detected light power that was scattered on a single particle with micrometer or even sub-micrometer dimension. Thus, 100,000 |v x |∕(m∕s) photons are detected in total for the given spatial and temporal resolution in dependence on the flow velocity component v x . The assumed typical characteristic length scales are the fringe distance d fringe = 10 μm (LDA), the scaled laser wavelength √ 2 ∕�⃗ o − ⃗ i� = = 532 nm (PDV), the beam radius 2w beam = 10 μm (L2F), the particle image radius 2w particle = 5 μm (PTV/PIV). Regarding the spectral filter at PDV, the transmission = 0.5 and the slope � = 2.5×10 −9 Hz −1 of the transmission curve are applied, Table 1 Parameter values for the numerical evaluation of Eq. (28) in Fig. 9 and the numerical evaluation of Eq. (31) in Fig. 10 Principle ( which are realistic values for a molecular absorption filter filled with iodine gas . The finally resulting shot noise limits for LDA, PDV, L2F and PTV/PIV are shown in Fig. 9 as a function of the respective measurand. L2F in principle attains the same shot noise limit as PTV/ PIV, but due to the chosen characteristic length scales, the shot noise limit of L2F is here a factor of 2 larger. L2F and LDA have similar shot noise limits, which is considered a remarkable result since the Doppler principle and the time-of-flight measurement principles are complementary principles from a physical perspective. Contrary to the similar shot noise limits of L2F, LDA and PTV/PIV, however, the shot noise limit of PDV is larger for flow velocities below 1.5 km/s (L2F), 2.4 km/s (LDA) and 3 km/s (PTV/PIV), respectively. The sensitivity, which is mainly determined by the slope of the transmission curve of the spectral filter to determine the Doppler frequency, is small in comparison with the LDA light mixing approach. However, a direct comparison of PDV with LDA, L2F and PTV/PIV is not possible because PDV measures (at least in parts) the out-of-plane velocity component, whereas the other measurement principles are only sensitive with respect to the in-plane velocity. Reformulating Eq. (28) by denoting the temporal resolution with the symbol T and the characteristic length scale for LDA, L2F, PTV/PIV with the symbol L x , the uncertainty relation connecting the measurement uncertainty with the temporal resolution reads for LDA, L2F and PTV/PIV and for PDV As a result, an improved temporal resolution means an increase in the velocity measurement uncertainty for a constant mean photon rate and a constant characteristic length. To minimize the uncertainty products on the left side of the formulas, the mean photon rate must be maximized and the characteristic length scale must be minimized. In addition, for PDV, the slope of the filter transmission curve needs to be maximized. Considering the diffraction limit of monochrome light, the characteristic length scale for all measurement principles can be minimized by minimizing the laser wavelength. Furthermore, the characteristic length scale can be minimized by maximizing the numerical aperture of the illumination (LDA, L2F) or observation (PTV/PIV) or between illumination and observation (PDV), which represents a significant difference between the different measurement principles. PTV/PIV can thus cope better than LDA, L2F with large working distances for the illumination, while LDA and L2F seem in principle more suitable to operate with large working distances for the observation.

Comparison for measuring multiple particles with averaging
Two kinds of averaging are possible for the case of measuring on multiple particles. Temporal averaging means that an averaging over consecutive measurements on single particles takes place during a given measurement time T. As a maximum, the averaging is possible over T∕T LDA , T∕T PDV , T∕T L2F , T∕T PTV measurement samples, respectively. In addition, PIV and PDV enable the simultaneous measurement on multiple particles in the interrogation region without disturbance. The latter is true for PIV as long as the particle images do not overlap, whereas no such constraint exists for PDV. Thus, significantly higher seeding particle concentrations can be applied for PDV than for PIV. As a result, an average over N particle,PIV < N particle,PDV particles is possible for PIV and PDV. By applying both kinds of averaging on the shot noise limits from Eq. (28) and using the definitions (3), (7), (11), (16), the shot noise limits for LDA, PDV, L2F and PTV/PIV read Calculated measurement uncertainty that is minimal achievable with LDA, PDV, L2F and PTV/PIV due to the shot noise limit as a function of the measurand v x and v, respectively. The calculation is based on the square root of Eq. (28) and is valid for the measurement with one single particle, with a comparable mean photon rate Ṅ photon = 1×10 9 s −1 as well as with a comparable spatial resolution of w = 100 μ m and temporal resolution w∕|v x | or w/|v|, respectively. A complete list of the applied parameter values is given in Table 1 Remember that the previous calculations are based on the assumption that the measurement time is equal or larger than the particle transit time through the measurement volume. Thus, the derived shot noise limits in Eq. (31) are only applicable for sufficiently large flow velocities. The minimal flow velocity follows from the respective spatial resolution 4 w fringes , w transit , d beam or x 2 − x 1 divided by the measurement time T. All shot noise limits are inversely proportional to the temporal resolution or the measurement time T. Except for PDV, which seems to be velocity-independent, all other shot noise limits now show a proportionality to the flow velocity squared. Therefore, the shot noise limit for LDA, L2F and PTV/PIV can be expressed in the form of a relative measurement uncertainty, whereas it is an absolute measurement uncertainty for PDV. However, note that the PDV result does depend on the velocity, if the linear range of the spectral filter's transmission curve is left and the laser frequency is not adjusted accordingly. Furthermore, all limits are directly proportional to the squared ratio of the characteristic length scale to the spatial resolution.
Inserting the same numerical values from Table 1 already studied in Sect. 4.1 and additionally assuming N particle,PIV = 5 , N particle,PDV = 50 as an example, the shot noise limits for the measurement time T = 1 ms are shown in Fig. 10 as a function of the measurand v x and v, respectively. According to the considered spatial resolution and measurement time, the calculation only applies for velocities larger than 0.1 m/s. Furthermore, a best case scenario is considered for PDV. The laser frequency is always adjusted to measure with the maximal available sensitivity at the edge of the spectral filter's transmission. Therefore, the calculated shot noise limit for PDV is velocity independent.
In comparison with L2F and LDA, the shot noise limit of PTV/PIV is the lowest due to the averaging over simultaneously measured particles. Being independent on the flow Var �v velocity and coping well with an even higher seeding concentration, the shot noise limit of PDV becomes lower than the shot noise limits of L2F, LDA and PTV/PIV at flow velocities > 210 m/s, > 330 m/s and > 950 m/s, respectively. Reformulating Eq. (31) by denoting each characteristic length scale with the symbol L x (LDA, L2F, PTV/PIV) or √ 2 ∕�⃗ o − ⃗ i� (PDV), respectively, the uncertainty relation connecting the measurement uncertainty with the spatial resolution w x and the temporal resolution T reads for LDA, L2F and PTV/PIV and for PDV Eq. (33) is identical with Eq. (30) for N particle = 1 , i.e., the uncertainty relation is the same in both cases, single and multiple particle measurement. Contrary to this, the uncertainty relation for LDA, L2F and PTV/PIV now includes the spatial resolution on the left side of Eq. (32) and also indicates an increasing uncertainty product with an increase in the flow velocity to be measured. Calculated measurement uncertainty that is minimal achievable with LDA, PDV, L2F and PTV/PIV due to the shot noise limit as a function of the measurand v x and v, respectively. The calculation is based on the square root of Eq. (31) and is valid for the measurement with multiple particles and averaging. Note that a common mean photon rate ( Ṅ photon = 1×10 9 s −1 ) and a comparable spatial resolution ( w = 100 μm ) and temporal resolution ( T = 1 ms ) are considered here. Note further that all parameter values used for the calculation are listed in Table 1 In order to show that the shot noise limit is in agreement with Heisenberg's uncertainty principle, limitations due to diffraction are now included in the consideration for LDA, L2F and PTV/PIV. The characteristic length scale L x for LDA results from the fringe distance and for L2F from the beam radius. Both are illumination features that are limited in size according to the Gaussian beam characteristics (illumination diffraction limit): L x ∼ ∕ tan(Θ x,illumination ) with the opening angle Θ x,illumination of the illumination. For PTV/PIV, L is the particle image size, and because of being an imaging technique the minimal resolvable size follows from Abbe's law (observation diffraction limit): L x ∼ ∕ sin(Θ x,observation ) with the opening angle Θ x,observation of the observation. Dividing Eq. (32) by T ⋅ |v x | 2 and applying the diffraction limit for L x in the direction of sensitivity with the approximation tan(Θ x,illumination ) ≈ sin(Θ x,illumination ) for sufficiently small respective opening angles Θ x yields for LDA, L2F and PTV/PIV In other words (after calculating the square root): The position uncertainty (x) of the measured particles times the uncertainty (k x ) of the wave vector component in the direction of sensitivity is inversely proportional to the square root of the total number of detected photons. The uncertainty relation holds for PTV/PIV also for (y) and (k y ) , respectively. The derived result is thus in agreement with Heisenberg's uncertainty principle and other shot noise studies with respect to position measurements (Pavliček and Pech 2016).
Regarding PDV with a molecular absorption filter as spectral filter, the Doppler broadening of the spectral absorption due to the moving absorbing molecules has a Gaussian shape with a full width at half-maximum that is directly proportional to the molecular resonance frequency (Svelto 2010). Furthermore, the resonance frequency is close to the scattered light frequency f. Therefore, it is assumed that the maximal slope of the spectral transmission is approximately � ∼ 1∕f and the working point is at = 0.5 . Inserting both assumptions in Eq. (33), dividing by T and enhancing the expression on the left side by the square of the particle transit time T PDV = w transit ∕|v| (definition see Eq. (7)) results in The derived uncertainty principle now indicates a relationship between the position uncertainty (x transit ) and the (34) Var Var uncertainty (k) of the wave vector component in the direction of sensitivity. The latter follows since the inverse of the particle transit time 1∕T PDV represents a resolution limit (f ) of the light wave frequency f according to the time-frequency uncertainty principle. Note that this is not a transverse but an axial light property. Furthermore, the relation (f ) ∼ c ⋅ (k) holds due to the wave number definition k = 2 and = c∕f with the light speed c. Since the observed light frequency depends on the difference between illumination and observation direction, a scaling factor of |⃗ o − ⃗ i| = 2 sin Θ 2 additionally occurs with Θ = ∠(⃗ o, ⃗ i) as the angle between the observation vector and the illumination vector (McKenzie 1996). As a result, the identified physical uncertainty principle for PDV in Eq. (35) has the same structure as for LDA, L2F and PTV/PIV in Eq. (34). Both identified uncertainty principles obey the same physical law, which is applied to different sensitivity directions.

Conclusion and outlook
The lower limit of the measurement uncertainty for a signal with a single unknown measurand was derived on the basis of an efficient estimator according to the Gauß-Markow theorem and the well-known rules of uncertainty propagation calculation. Inserting Poissonian noise, the shot noise limit was then provided for discrete signal samples and continuous signals. Note that the restriction to one unknown parameter (the measurand) means that an ultimate lower bound for the minimal achievable measurement uncertainty is obtained.
The general solution to determine the shot noise limit was further applied for LDA, PDV, L2F and PTV/PIV. As a result, the shot noise limits for LDA, L2F and PTV/PIV have an identical structure. PDV is different because it is the only measurement principle that is sensitive with respect to the out-of-plane velocity component. Due to the different sensitivity direction and the different measurand, respectively, a one-to-one comparison with the other measurement principles is difficult. However, a numerical example revealed that the minimal achievable measurement uncertainty due to photon shot noise is comparably large although the largest number of particles can be measured simultaneously. Note that the numerical comparison was performed for an identical spatiotemporal resolution as well as an identical mean photon rate.
The shot noise limits were studied for the case of measuring a single seeding particle and the case of measuring multiple seeding particles and averaging. Averaging multiple particle measurements over a given measurement time leads to a minimal achievable velocity measurement uncertainty that is directly proportional to the absolute value of the velocity (LDA, L2F, PTV/PIV) and that is independent on the velocity value (PDV), respectively. Studying the dependency on the spatiotemporal resolution, the measurement uncertainty is always inversely proportional to the square root of the measurement time that here is the temporal resolution. The measurement uncertainty of the flow velocity is further inversely proportional to the spatial resolution for LDA, L2F, PTV/PIV and it does not depend on the spatial resolution for PDV, when an identical mean photon rate is considered.
For all measurement principles, the shot noise limit could be shown to obey an uncertainty principle that is in agreement with Heisenberg's uncertainty principle and Poisson statistics: The product of the position uncertainty (that is the velocity uncertainty times the particle observation time in the interrogation region) and the uncertainty of the respective wave vector component along the sensitivity direction is inversely proportional to the square root of the total number of detected photons. As a result, all studied Doppler and time-of-flight measurement principles obey the same fundamental limit from the perspective of quantum mechanics.
The theoretically identified limits provide a common solid basis to assess different laser-based flow velocity measurement systems with respect to what is feasible from the perspective of photon shot noise (fundamental physical limitation). Experimental studies must follow in the next step, in order to demonstrate shot noise-limited measurements for each measurement principle and to validate the theoretical findings. Furthermore, the dependency of the scattering light power and the signal's sensitivity on the particle size as well as the angular scattering characteristic should be included in the estimation of the minimal achievable measurement uncertainty in future. Finally, shot noise is typically not the only limiting factor for laser-based flow velocity measurements. Another fundamental physical aspect would be to compare the measurement capabilities for a non-constant flow velocity, which means to consider the effect of an accelerated flow and flow turbulence on the flow velocity measurement. Furthermore, the signal processing efficiency should be included in future comparison studies, because the velocity estimation algorithm might not be able to attain the Cramér-Rao bound (signal processing limitation). Future comparison studies should also include the effect of signal quantization over time or space (for PTV see Wernet and Pline 1993 for instance), respectively, as well as other unavoidable error sources such as background light and refractive index changes in flame flows and two-or multi-phase flows. This means to compare the measurement principles with respect to current technological as well as applicationdriven physical limitations.