Ladar range estimate with range rate compensation

Patent metadata

Field	Value
Publication number	US-10620315-B2
Application number	US-201715490857-A
Country	US
Kind code	B2
Filing date	Apr 18, 2017
Priority date	Apr 18, 2017
Publication date	Apr 14, 2020
Grant date	Apr 14, 2020

Abstract

A system and method for forming a range estimate for a target with a laser detection and ranging system. The system includes a laser transmitter and an array detector. The method includes: transmitting a plurality of laser pulses; for each transmitted laser pulse: detecting, with the array detector, a plurality of ladar return photons from the laser pulse, each detection producing an electrical pulse; identifying, for each of the electrical pulses, a time bin of a plurality of time bins corresponding to the laser pulse, within which the electrical pulse was produced; forming a one dimensional range histogram array having, for each of a subset of the plurality of bins, an element with a value equal to the number of electrical pulses produced in the array detector during a time interval corresponding to the bin; and forming an estimated range rate for the target.

First claim

What is claimed is:

1. A method for forming a range estimate for a target with a laser detection and ranging system comprising a laser transmitter and an array detector, the method comprising:
transmitting a plurality of laser pulses;
for each transmitted laser pulse:
detecting, with the array detector, a plurality of ladar return photons from the laser pulse, each detection producing an electrical pulse;
identifying, for each of the electrical pulses, a time bin of a plurality of time bins corresponding to the laser pulse, within which the electrical pulse was produced; and
forming a one dimensional range histogram array of a plurality of range histogram arrays having, for each of a subset of the plurality of time bins, an element with a value equal to a number of electrical pulses produced in the array detector during a time interval corresponding to the time bin;

forming an estimated range rate for the target;
forming a two dimensional adjusted range versus frame array having a first plurality of columns, each of the first plurality of columns being identified by a column number and including a portion of one of the range histogram arrays, shifted in proportion to a product of the estimated range rate and the column number;
summing rows of the two dimensional adjusted range versus frame array to form a collapsed histogram; and
determining the range from a location of a peak in the collapsed histogram,
wherein the forming of an estimated range rate comprises:
forming a two dimensional range versus frame array having a second plurality of columns, each of the second plurality of columns corresponding to one of the range histogram arrays;
forming a Radon transform of the range versus frame array to form a Radon transform array, the Radon transform array being a two dimensional array having a plurality of rows each corresponding to a respective range rate and a third plurality of columns;
summing the third plurality of columns of element-wise square of the Radon transform array to form a collapsed Radon transform array; and
calculating an initial estimated range rate as a range rate corresponding to a maximum value of the collapsed Radon transform array.

2. The method of claim 1 , wherein the forming of an estimated range rate further comprises:
setting a tentative range rate to be equal to the initial estimated range rate;
iteratively adjusting the tentative range rate, through a sequence of values until a completion criterion is satisfied, and for each of the sequence of values:
generating a two dimensional tentative adjusted range versus frame array having a plurality rows and a fourth plurality of columns, each of the fourth plurality of columns corresponding to one of the range histogram arrays, shifted in proportion to the product of the estimated range rate and the column number;
summing the rows of the tentative adjusted range versus frame array to form a tentative adjusted collapsed histogram; and
calculating a measure of sharpness for a peak in the tentative adjusted collapsed histogram.

3. The method of claim 2 , wherein the completion criterion is satisfied when the measure of sharpness exceeds a threshold value.

4. The method of claim 2 , wherein the calculating of a measure of sharpness comprises:
determining a maximum value of the tentative adjusted collapsed histogram.

5. The method of claim 2 , wherein the calculating of a measure of sharpness comprises:
fitting a Gaussian function to the tentative adjusted collapsed histogram; and
calculating a reciprocal of a width of the Gaussian function.

6. The method of claim 2 , wherein the calculating of a measure of sharpness comprises:
calculating a ratio of:
a numerical integral of a square of the tentative adjusted collapsed histogram and
a numerical integral of the tentative adjusted collapsed histogram.

7. The method of claim 2 , wherein the portion of one of the range histogram arrays is a portion centered on a bin corresponding to an initial estimate of a range to the target.

8. A ladar system for forming a range estimate for a target, the ladar system comprising:
a pulsed laser transmitter configured to transmit a sequence of laser pulses;
a camera comprising an array detector, the camera being configured to:
detect, for each laser pulse of the sequence of laser pulses, a plurality of ladar return photons from the laser pulse, each detection producing an electrical pulse; and
identify, for each of the electrical pulses, a time bin of a plurality of time bins corresponding to the laser pulse, within which the electrical pulse was produced; and

a processing circuit, configured to:
form a one dimensional range histogram array of a plurality of range histogram arrays having, for each of a subset of the plurality of time bins, an element with a value equal to a number of electrical pulses produced in the array detector during a time interval corresponding to the time bin;
form an estimated range rate for the target;
form a two dimensional adjusted range versus frame array having a first plurality of columns, each of the first plurality of columns being identified by a column number and including a portion of one of the range histogram arrays, shifted in proportion to a product of the estimated range rate and the column number;
sum rows of the two dimensional adjusted range versus frame array to form a collapsed histogram; and
determine a range from a location of a peak in the collapsed histogram,

wherein the processing circuit is configured to form an estimated range rate by:
forming a two dimensional range versus frame array having a second plurality of columns, each of the second plurality of columns corresponding to one of the range histogram arrays;
forming a Radon transform of the range versus frame array to form a Radon transform array, the Radon transform array being a two dimensional array having a plurality of rows each corresponding to a respective range rate and a third plurality of columns;
summing the third plurality of columns of element-wise square of the Radon transform array to form a collapsed Radon transform array; and
calculating an initial estimated range rate as a range rate corresponding to a maximum value of the collapsed Radon transform array.

9. The ladar system of claim 8 , wherein the processing circuit is configured to form an estimated range rate by:
setting a tentative range rate to be equal to the initial estimated range rate;
iteratively adjusting the tentative range rate, through a sequence of values until a completion criterion is satisfied, and for each of the sequence of values:
generating a two dimensional tentative adjusted range versus frame array having a plurality rows and a fourth plurality of columns, each of the fourth plurality of columns corresponding to one of the range histogram arrays, shifted in proportion to the product of the estimated range rate and the column number;
summing the rows of the tentative adjusted range versus frame array to form a tentative adjusted collapsed histogram; and
calculating a measure of sharpness for a peak in the tentative adjusted collapsed histogram.

10. The ladar system of claim 9 , wherein the completion criterion is satisfied when the measure of sharpness exceeds a threshold value.

11. The ladar system of claim 9 , wherein the processing circuit is configured to calculate a measure of sharpness by:
determining a maximum value of the tentative adjusted collapsed histogram.

12. The ladar system of claim 9 , wherein the processing circuit is configured to calculate a measure of sharpness by:
fitting a Gaussian function to the tentative adjusted collapsed histo

Assignees

• Raytheon Co

Inventors

• Zellinger Louis
• Kaganovich Vitaliy M
• Marcus Eran

Classifications

• G01S7/4863 Primary

  Detector arrays, e.g. charge-transfer gates · CPC title

• G01S7/4865

  Time delay measurement, e.g. time-of-flight measurement, time of arrival measurement or determining the exact position of a peak (peak detection in noise, signal conditioning G01S7/487) · CPC title

• G01S17/58

  Velocity or trajectory determination systems; Sense-of-movement determination systems · CPC title

• G01S17/50

  Systems of measurement based on relative movement of target · CPC title

• G01S7/497

  Means for monitoring or calibrating · CPC title