Transcript Considering Polarimeters as Modulated Information Systems
Microgrid Modulation Study F
High Resolution Sinusoidal Image High Resolution Image Image FFT
F
Simulated Microgrid Image (Fully Polarized) Fully Polarized Sinusoidal Microgrid Image FFT
Microgrid Demodulation Study F
Reconstructed sinusoidal s 0 image. Orientation ~30 deg, freq~.3 cycles/pixel
F
s
0 FFT using no interpolation Reconstructed s 1 image. Should be zero.
s
1 FFT using no interpolation
Systems Model of Microgrid Stokes Parameters: Ideal Microgrid Analyzer Vectors: H -45 45 V Modulated Analyzer Vector: Analyzer Layout on Focal Plane Analytical Expression for Modulated Intensity:
s 0 Analytic Test Scene s 1 Test scene using Gaussian stokes parameters.
s 2 I Intensity output simulated according to:
Data Distribution in Fourier Plane s 2 -s 1 s 2 +s 1 Band limited data easily recovered with use of simple filters s 0
Application to Real Data
Microgrid Calibrated Output Fourier Transform of Raw Output
In this imagery example, the emissive sphere provides spatially varying angle and degree of polarization. The spatial frequency plane representation of the information clearly shows the base band data (s s
2
).
0
) and the horizontal and vertical side lobes (s
1
and
Edge Artifacts and Polarimetric Aliasing
By properly reconstructing the data in the frequency domain (or by using appropriately constructed spatial kernels) polarimetric aliasing can be completely eliminated, leading to ideal reconstruction of band limited imagery
Gaussian Kernel NLPN
Modulated Polarimeters
Modulated polarimeters
are a second class of device that uses a single detector or detector array to measure an intensity distribution that has been modulated in space, time, and/or wavelength – – –
Spatial Modulation: Temporal Modulation: Spectral Modulation:
Prismatic polarimeter, microgrid polarimeter Rotating retarder, PEM-based devices Channeled spectropolarimeter Modulated polarimeters have several inherent strengths: – Inherently spatiotemporally aligned – Can be mechanically rugged and vibration tolerant – Need for cross-calibration of sensors is eliminated These benefits come at the cost of a reduction in overall system bandwidth
Communications Theory Revisited Consider two signals present in the same interval of time:
s
1 gaus
t t
s
2 ;
s
2 gaus How do we separate the two signals after measurement?
In order to use the available bandwidth of the detector to measure both signals, we have to create multiplexed
channels
. The simplest and most common method is OFDM, but other methods can also be used.
s
2
Frequency Domain Channels After taking a Fourier transform of the measured intensity, we see that the signals have been separated in frequency
S 2 S 1
1 2
S
1
f
f
0
S
1
f
f
0
S
2 Demodulate the
s
1 channel by
homodyning
2 and
LPF.
• Multiply by cos(2
f
0
t
) • Equalize the channels (factor of 2) • Low Pass Filter
Polarimeters as Multiplexed Systems
A
S
* Modulated Analyzer Vector Stokes parameters to be measured Detector response System PSF • • • System PSF is really a Mueller matrix that alters the polarization state Consider continuous systems with ideal sampling (
h
=
(x,y,t,
)
,
d
=
(x,y,t,
)
) Consider LSI systems, even though this is not always true
a
0
a
2
a
3 , , , , , ,
Spatial Modulation H -45 45 V
Spectral Modulation
Channeled Spectropolarimetry – Slide courtesy Julia Craven-Jones
Two high-order, stationary retarders (usually 1:2 or 3:1 thickness ratio) and an analyzer placed in front of a spectrometer.
Modulates the spectrum with the Stokes parameters.
All Stokes parameters encoded in a single spectrometer measurement.
Channels are isolated and Fourier filtered to recover the Stokes vector.
center burst 7 channels
13
Temporal Modulation The UA-JPL MSPI Polarimeter is an examply of a more complicated modulated system. This 2-channel, temporally modulated polarimeter has analyzer vectors given as: 0 1, 1
J
0
2 0
f t
0
Consider the DRM:
W
1
W W
1
W
1
Inverse of the modulator inner product matrix that serves to unmix the Stokes parameter components and equalize their amplitudes Homodyne process that multiplies by the original modulation functions and filters using a rectangular window.
For simplicity we consider temporal modulation only, but these concepts are general.
Homodyning
n
W I
W WS
t
0
t
0
T
T
/2 3 2
m
0
n
This is nothing more than a homodyne plus a LPF with a rectangular window. We can buy ourselves some flexibility in choosing our LPF by considering the two operations separately
W
I
Modulator Inner Product Matrix The nextstep in the DRM is to compute (
W
T
W
) -1 =
Z
-1 . This multiplication also implicitly includes a rectangular-window LPF In the general case, define the this matrix as:
Z ij
i j
0
i j w
is an arbitrary window function that corresponds to the LPF we will use in reconstruction. In the traditional DRM,
w
is a rectangular window that includes
N
sample intervals. ` For modulation schemes with orthogonal modulators,
Z
is diagonal, but for more complicated schemes,
Z
also unmixes the Stokes parameter signals.
Generalized Inversion
Z
1
Filtered Equalizer Filtered Homodyne
As an aside, this process can be used to create a modified DRM that includes the arbitrary window function with more desirable LPF characteristics.
Example Case – Rotating Retarder Consider a rotating retarder polarimeter with retardance rotating at angular frequency 2
f
0 .
A
0 1 2 3 2 2 1 2 sin 8 2 sin sin 4 cos8
f t
0
f t
0
f t
0
Z
1 12 cos 2 cos 16 cos cos 0 0 8cos 1 1 2 16 2 12 16 cos cos 16 1 2 cos 32 1 2 0 0 0 0 cos 32 1 2 0 0 0 0 8 sin 2
Band Limited Input Signal The input signal is fully polarized and band limited:
S
1 1 1 3 sinc 2
Measured Signal by Component
Each quadrant gives the product a n
(
t
) s n (
t
), and the total signal is obtained by summing these four. The Fourier transform is broken up by parameter.
Homodyne Output in Frequency
Unmixed Output in Frequency
Convention DRM Output in Frequency
Reconstructred Signal in Time
We see the clear effect of high frequencies that “leak” through and corrupt the reconstruction. We should note that this error is completely avoidable, and is a result of the implied rectangular window used with the DRM method.
Reconstruction vs Signal Bandwidth
Aliasing and Cross Talk in Time
When the bandwidth exceeds limits set by the modulation frequency and the sampling frequency, we end up with both aliasing (self interference) and cross talk that corrupt the reconstructed signal.
Space-Time Modulated Polarimetry Microgrid Polarimeter modulates in space and creates up to three side bands in spatial frequency space:
Rotating Retarder Polarimeter The rotating retarder modulates the intensity in time, creating two pairs of complex side bands along the temporal frequency axis
Space-Time Modulation A rotating retarder followed by a microcrid now creates side bands along the temporal frequency axis for each of the side bands in the spatial frequency plane
Space-Time Modulation Continued If that retarder is a HWP, we lose the ability to sense
s
3 , but we get the maximum separation of our side bands in the spatio-temporal frequency cube.
Space-Time Simulation
Input Data from DoAmp Polarimeter
Single Pixel Temporal Traces
Simulation
s
0 – Microgrid and RR Microgrid RR
Polarization Data (
s
1 ) – Microgrid vs RR Microgrid RR
Optimized versus Microgrid (
s
0 ) Optimized Microgrid
Optimized versus RR (
s
0 ) Optimized RR
Optimized versus Microgrid (
s
1 ) Optimized Microgrid
Optimized versus RR (
s
1 ) Optimized RR
The Role of the Operator Null Space Let’s consider different window sizes and shapes in order to assess which method is “better” for any given reconstruction task.
Reconstructed Signal at
t
= 0
S
t
0
Z
1
Z
1
f
0
Consider a rectangular window with N
.
S
t
0
Z
1
A
Z
1
A
This is just the inner product of the Fourier transforms
of A and
I
. For periodic modulators, we select
only
the frequencies included in the modulators. All others are in the null space.
Domain of the TraditionalDRM Operator
This is just the inner product of the Fourier transforms
of A and
I
. For periodic modulators, we select
only
the frequencies included in the modulators. All others are in the null space.
Domain of the Band Limited Operator
S
t
0
rect 2
f f
0 *
A
For the band limited inversion operator, the null space includes all frequencies other than the band assigned to a particular Stokes parameter.
Wiener-Helstrum Window In the general case where we have an expected signal PSD and a known noise PSD, we can construct “optimal” inversion operators that use a window that emphasizes the areas of the spectrum with greatest SNR. An example is the W-H window
w n
W pn s n
2
pn
W pn
2 Assuming uncorrelated, band-limited Stokes parameters with Z independent of time.
Constructing the Wiener Filter
Wiener filters for
s
0 and
s
1
Results with the Microgrid
Results with the Microgrid
Results with the Microgrid