WO2021251161A1 - 測距装置 - Google Patents

測距装置 Download PDF

Info

Publication number
WO2021251161A1
WO2021251161A1 PCT/JP2021/020214 JP2021020214W WO2021251161A1 WO 2021251161 A1 WO2021251161 A1 WO 2021251161A1 JP 2021020214 W JP2021020214 W JP 2021020214W WO 2021251161 A1 WO2021251161 A1 WO 2021251161A1
Authority
WO
WIPO (PCT)
Prior art keywords
time
light
distance
distance measuring
light receiving
Prior art date
Application number
PCT/JP2021/020214
Other languages
English (en)
French (fr)
Inventor
謙太 東
Original Assignee
株式会社デンソー
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from JP2021070064A external-priority patent/JP7375789B2/ja
Application filed by 株式会社デンソー filed Critical 株式会社デンソー
Priority to CN202180042031.1A priority Critical patent/CN115698756A/zh
Publication of WO2021251161A1 publication Critical patent/WO2021251161A1/ja
Priority to US18/064,010 priority patent/US20230108583A1/en

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88Lidar systems specially adapted for specific applications
    • G01S17/93Lidar systems specially adapted for specific applications for anti-collision purposes
    • G01S17/931Lidar systems specially adapted for specific applications for anti-collision purposes of land vehicles
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C3/00Measuring distances in line of sight; Optical rangefinders
    • G01C3/02Details
    • G01C3/06Use of electric means to obtain final indication
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/06Systems determining position data of a target
    • G01S17/08Systems determining position data of a target for measuring distance only
    • G01S17/10Systems determining position data of a target for measuring distance only using transmission of interrupted, pulse-modulated waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/06Systems determining position data of a target
    • G01S17/42Simultaneous measurement of distance and other co-ordinates
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/481Constructional features, e.g. arrangements of optical elements
    • G01S7/4811Constructional features, e.g. arrangements of optical elements common to transmitter and receiver
    • G01S7/4813Housing arrangements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/481Constructional features, e.g. arrangements of optical elements
    • G01S7/4817Constructional features, e.g. arrangements of optical elements relating to scanning
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/483Details of pulse systems
    • G01S7/486Receivers
    • G01S7/4861Circuits for detection, sampling, integration or read-out
    • G01S7/4863Detector arrays, e.g. charge-transfer gates
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/483Details of pulse systems
    • G01S7/486Receivers
    • G01S7/4865Time delay measurement, e.g. time-of-flight measurement, time of arrival measurement or determining the exact position of a peak
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/483Details of pulse systems
    • G01S7/486Receivers
    • G01S7/487Extracting wanted echo signals, e.g. pulse detection
    • G01S7/4873Extracting wanted echo signals, e.g. pulse detection by deriving and controlling a threshold value

Definitions

  • This disclosure relates to a distance measuring device.
  • the target Distance measuring devices that detect the presence or absence of an object and measure the distance to an object are known.
  • Light the so-called clutter, can be an obstacle in measuring the distance to an object.
  • the reflected light from the object and the clutter are combined and detected as one light receiving pulse, which may make it impossible to measure the flight time correctly and reduce the accuracy of the measured distance. Therefore, in the distance measuring device of Patent Document 1 below, a value larger than the amplitude value (intensity) of the clutter is set as a threshold value (High threshold value) as a threshold value for the light receiving intensity, and a High threshold value is set for the light receiving intensity of the combined light receiving pulse.
  • the cause of clutter is fog
  • the thicker the fog the larger the amplitude value of the clutter
  • the cause of the clutter is a window or a housing for emitting pulsed light
  • the larger the reflectance of the window or the housing surface the larger the amplitude of the clutter, and the amplitude of the clutter and the desired wave.
  • the difference becomes smaller.
  • the difference in amplitude between the clutter and the desired wave becomes small in this way, the time change of the light receiving intensity of the combined wave becomes a shape like a mountain.
  • Such a phenomenon is not limited to the case where a clutter occurs, and for example, when there are a plurality of objects having similar directions and different distances with respect to the distance measuring device, the reflected light from the plurality of objects is present. Can occur together. In this way, when the light receiving intensities of a plurality of reflected lights are combined into a single mountain-like shape, the distance measuring device of Patent Document 1 cannot separate the desired wave from the received wave, and the distance measuring accuracy is high. Cannot be improved. In Patent Document 1, for example, in a bad environment such as heavy fog or heavy rain, the distance measurement is not performed to suppress the deterioration of the distance measurement accuracy. However, even in a situation where time changes in the light receiving intensity of a plurality of reflected lights are combined to form a single mountain-like shape, such as in a bad environment, a technique capable of accurately measuring the distance is desired.
  • a ranging device utilizes the light emitting unit that emits pulsed light, the light receiving unit that receives light including the reflected light of the pulsed light, and the flight time of the light received by the light receiving unit to obtain the pulsed light. It is provided with a calculation unit for calculating the object distance, which is the distance to the object that reflects the reflected light and outputs the reflected light.
  • the calculation unit includes a histogram generating unit that generates a histogram representing the light receiving intensity in each of the light receiving units in the plurality of flight times, and the change in the light receiving intensity with respect to the flight time in the histogram with the peak flight time as a boundary.
  • Two flight times, which are the same as the first threshold light receiving intensity, which is larger than the light receiving intensity of the noise light including the background light, in the compound mountain part estimating part for estimating the existence of A flight time specifying part that specifies a rise time that is a time and a fall time that is a longer flight time, and the rise that is specified for the composite mountain portion when it is estimated that the composite mountain portion exists.
  • the base flight time determination unit that determines the base flight time that is the basis for calculating the object distance based on either the time or the fall time, and the base flight time are used to determine the base flight time.
  • a distance calculation unit for calculating the distance of one of a plurality of objects as the object distance is provided.
  • the base flight time which is the basis for calculating the object distance, is determined based on either the rise time or the fall time specified for the compound mountain portion. Since the object distance is calculated using the determined base flight time, the object distance can be calculated even in a situation where the time changes of the light receiving intensity of multiple reflected lights are combined to form a single mountain-like shape. It can be measured with high accuracy.
  • FIG. 1 is a schematic configuration diagram of a distance measuring device as an embodiment of the present disclosure.
  • FIG. 2 is an explanatory diagram schematically showing the configuration of the light receiving array.
  • FIG. 3 is a circuit diagram schematically showing the configuration of the SPAD circuit.
  • FIG. 4 is a block diagram showing a functional configuration of the distance measuring device according to the first embodiment.
  • FIG. 5 is an explanatory diagram showing an example of a histogram having a composite mountain portion.
  • FIG. 6 is a flowchart showing the procedure of the distance measuring process according to the first embodiment.
  • FIG. 1 is a schematic configuration diagram of a distance measuring device as an embodiment of the present disclosure.
  • FIG. 2 is an explanatory diagram schematically showing the configuration of the light receiving array.
  • FIG. 3 is a circuit diagram schematically showing the configuration of the SPAD circuit.
  • FIG. 4 is a block diagram showing a functional configuration of the distance measuring device according to the first embodiment.
  • FIG. 5 is an explanatory diagram showing an example of
  • FIG. 7 is an explanatory diagram showing an example of a histogram having no compound mountain portion.
  • FIG. 8 is a block diagram showing a functional configuration of the distance measuring device according to the second embodiment.
  • FIG. 9 is a flowchart showing the procedure of the distance measuring process in the second embodiment.
  • FIG. 10 is a plan view schematically showing the positional relationship between the plurality of objects and the distance measuring device in the second embodiment.
  • FIG. 11 is an explanatory diagram showing an example of a histogram having a composite mountain portion in the second embodiment.
  • FIG. 12 is an explanatory diagram showing an example of a histogram having a composite mountain portion in the second embodiment.
  • FIG. 13 is an explanatory diagram showing an example of a histogram having a composite mountain portion in the second embodiment.
  • FIG. 14 is an explanatory diagram showing an example of a histogram having a composite mountain portion in the third embodiment.
  • FIG. 15 is a flowchart showing the procedure of the clutter position setting process in the fourth embodiment.
  • FIG. 16 is an explanatory diagram showing an example of a distance histogram in the fourth embodiment.
  • FIG. 17 is a flowchart showing the procedure of the clutter position setting process in the fifth embodiment.
  • FIG. 18 is an explanatory diagram showing an example of a histogram in the fifth embodiment.
  • FIG. 19 is a flowchart showing the procedure of the distance measuring process according to the sixth embodiment.
  • FIG. 20 is an explanatory diagram showing an example of the histogram in the sixth embodiment.
  • FIG. 21 is a flowchart showing a procedure of distance measurement processing in the seventh embodiment.
  • the distance measuring device 10 shown in FIG. 1 includes an optical system 30 that emits pulsed light for distance measurement and receives reflected light from an external object, and a calculation unit 20 that processes a signal obtained from the optical system 30. Be prepared. External objects are also called "reflectors".
  • the optical system 30 includes a light emitting unit 40 that emits laser light as pulsed light, a scanning unit 50 that scans the laser light within a predetermined viewing range 80, and incident light including reflected light and ambient light from an external object.
  • a light receiving unit 60 for receiving light is provided.
  • the distance measuring device 10 is housed in a housing 90 having a window 92 on the front surface. The window 92 transmits most of the pulsed light emitted from the light emitting unit 40 and reflects a part of the pulsed light.
  • the distance measuring device 10 is, for example, an in-vehicle LiDAR (Laser Imaging Detection and Ringing) mounted on a vehicle such as an automobile.
  • LiDAR Laser Imaging Detection and Ringing
  • the horizontal direction of the visual field range 80 coincides with the horizontal direction X
  • the vertical direction coincides with the vertical direction Y.
  • the light emitting unit 40 was emitted from a semiconductor laser element (hereinafter, also simply referred to as a laser element) 41 that emits laser light including pulsed light, a circuit board 43 incorporating a drive circuit of the laser element 41, and a laser element 41.
  • a collimating lens 45 for converting a laser beam into parallel light is provided.
  • the laser element 41 is a laser diode capable of oscillating a so-called short pulse laser.
  • the laser element 41 constitutes a rectangular laser emission region by arranging a plurality of laser diodes along the vertical direction.
  • the intensity of the laser beam output by the laser element 41 is configured to be adjustable according to the voltage supplied to the laser element 41.
  • the scanning unit 50 is composed of a so-called one-dimensional scanner.
  • the scanning unit 50 includes a mirror 54, a rotary solenoid 58, and a rotating unit 56.
  • the mirror 54 reflects a laser beam that is collimated by the collimated lens 45.
  • the rotary solenoid 58 receives a control signal from the calculation unit 20 and repeats forward rotation and reverse rotation within a predetermined angle range.
  • the rotating portion 56 is driven by the rotary solenoid 58, repeats forward rotation and reverse rotation on a rotation axis whose axial direction is the vertical direction, and scans the mirror 54 in one direction along the horizontal direction.
  • the laser beam emitted from the laser element 41 via the collimating lens 45 is reflected by the mirror 54 and scanned along the horizontal direction by the rotation of the mirror 54.
  • the field of view range 80 shown in FIG. 1 corresponds to the entire scanning range of this laser beam. Since the light receiving intensity is obtained at each pixel position in the visual field range 80, the distribution of the light receiving intensity in the visual field range 80 constitutes a kind of image.
  • the scanning unit 50 may be omitted, and the light emitting unit 40 may emit pulsed light over the entire visual field range 80, and the light receiving unit 60 may receive the reflected light over the entire visual field range 80.
  • pulsed light is irradiated to each position in the scanning range, in other words, each pixel position in the field of view 80. Then, the irradiation of the pulsed light and the distance measuring process described later based on the reflected light from each pixel position are executed at predetermined time intervals for each pixel position.
  • the laser beam output from the light emitting unit 40 is diffusely reflected on the surface, and a part of the laser light returns to the mirror 54 of the scanning unit 50 as reflected light.
  • This reflected light is reflected by the mirror 54, is incident on the light receiving lens 61 of the light receiving unit 60 as incident light together with the ambient light, is condensed by the light receiving lens 61, and is incident on the light receiving array 65.
  • the laser beam output from the ranging device 10 is diffusely reflected not only by an external object but also by an object inside the ranging device 10, for example, a window 92, and a part of the reflected light is incident on the light receiving array 65. do.
  • the light receiving array 65 is composed of a plurality of pixels 66 arranged two-dimensionally.
  • One pixel 66 is composed of a plurality of SPAD (Single Photon Avalanche Diode) circuits 68 arranged so as to have H in the horizontal direction and V in the vertical direction.
  • H and V are integers of 1 or more, respectively.
  • each of the five SPAD circuits 68 is configured in the horizontal direction and the vertical direction.
  • the light receiving result of one pixel 66 is the light receiving intensity at one pixel position within the visual field range 80.
  • the SPAD circuit 68 connects an avalanche diode Da and a quench resistor Rq in series between the power supply Vcc and the ground line, and the voltage at the connection point is one of the logic calculation elements. It is input to the inverting element INV and converted into a digital signal with the voltage level inverted. The output signal Sout of the inverting element INV is output as it is to the outside.
  • the quenching resistor Rq is configured as a FET, and if the selection signal SC is active, the on-resistance thereof acts as the quenching resistor Rq.
  • the quench resistor Rq When the selection signal SC becomes inactive, the quench resistor Rq is in a high impedance state, so that even if light is incident on the avalanche diode Da, the quench current does not flow, and as a result, the SPAD circuit 68 does not operate.
  • the selection signal SC is collectively output to the 5 ⁇ 5 SPAD circuits 68 in the pixel 66, and is used to specify whether to read the signal from each pixel 66 or not.
  • the avalanche diode Da is operated in the Geiger mode, but the avalanche diode Da may be used in the linear mode and the output thereof may be treated as an analog signal. Further, a PIN photodiode may be used instead of the avalanche diode Da.
  • the avalanche diode Da If no light is incident on the SPAD circuit 68, the avalanche diode Da is kept in a non-conducting state. Therefore, the input side of the inverting element INV is maintained in a state of being pulled up via the quench resistor Rq, that is, at a high level H. Therefore, the output of the inversion element INV is maintained at the low level L.
  • the avalanche diode Da is energized by the incident light (photons). As a result, a large current flows through the quench resistor Rq, the input side of the inverting element INV temporarily becomes the low level L, and the output of the inverting element INV is inverted to the high level H.
  • the inverting element INV outputs a high-level pulse signal for a very short time.
  • the output signal of the inverting element INV that is, the output signal Sout from each SPAD circuit 68 is the avalanche diode Da. It becomes a digital signal that reflects the state.
  • the output signal Sout corresponds to a pulse signal generated by receiving incident light including reflected light and ambient light in which the irradiation light is reflected by an external object existing in the scanning range, the window 92, and the like and returned.
  • the calculation unit 20 uses the flight time of the reflected light received by the light receiving unit 60 to determine the distance to the object that reflects the pulsed light and outputs the reflected light (hereinafter referred to as “object distance”). Calculate.
  • object distance the distance to the object that reflects the pulsed light and outputs the reflected light (hereinafter referred to as “object distance”). Calculate.
  • the outline of the calculation method of such a distance is as follows. As shown in FIG. 4, the pulsed light P1 emitted from the light emitting unit 40 is reflected by the reflecting object OBJ which is an external object. In other words, the reflecting object OBJ outputs the reflected light P2 of the pulsed light P1. Further, the pulsed light P1 is also reflected on the inner surface of the window 92, and the reflected light P3 is output. As a result, the reflected light P2 and P3 reach the light receiving unit 60.
  • the time from the emission of the pulsed light P1 to the light reception of the reflected lights P2 and P3 is specified as the flight time Tf of the light.
  • the calculation unit 200 calculates the distance from the distance measuring device 10 (light emitting unit 40 and the light receiving unit 60) to the reflecting object OJB using this flight time Tf.
  • the reflected light obtained by reflecting the pulsed light P1 on the inner surface of the window 92 and the light reflected by the reflected light several times inside the housing of the distance measuring device 10 and received by the light receiving unit 60 are , Also known as "clutter".
  • the calculation unit 20 includes an addition unit 220, a histogram generation unit 230, a composite mountain portion estimation unit 240, a flight time identification unit 250, a base flight time determination unit 260, and a distance calculation unit 270. And a control unit 280 and a memory 290.
  • the addition unit 220 adds the outputs of each SPAD circuit 68 included in the pixels 66 constituting the light receiving array 65.
  • the SPAD circuit 68 included in the pixel 66 operates.
  • the SPAD circuit 68 can detect only one photon incident.
  • the detection of the limited light output from the reflector OBJ must be probabilistic. Therefore, the addition unit 220 adds the output signal Sout from the SPAD circuit 68, which cannot detect the light incident only probabilistically, for all the SPAD circuits 68 included in each pixel 66, thereby adding each pixel 66. It is configured to more reliably detect the reflected light from the reflecting object OBJ in the above.
  • the histogram generation unit 230 generates a histogram of the light receiving intensity by acquiring the addition result of the addition unit 220 in time series, and stores it in the memory 290.
  • the histogram generated by the histogram generation unit 230 can be said to be a graph showing the light receiving intensity at each of the plurality of flight times.
  • the light receiving intensity is the total number of light received SPAD circuits 68 in one pixel 66.
  • the composite mountain portion estimation unit 240 estimates whether or not the composite mountain portion exists in the histogram generated by the histogram generation unit 230.
  • the compound mountain portion is a mountain portion in which the change in the light receiving intensity with respect to the flight time changes from an increase to a decrease with respect to the peak flight time in the histogram, and the difference between the distances from the distance measuring device 10 is within a predetermined range. It means a mountain part obtained by synthesizing the distribution of light intensity from a plurality of objects.
  • the horizontal axis represents the flight time Tf
  • the vertical axis represents the light receiving intensity I.
  • the histogram hr1 shown by the thick solid line is a combination of the histogram hr11 of the light receiving intensity by the clutter shown by the thin alternate long and short dash line and the histogram hr12 of the light receiving intensity of the reflected light by the object outside the ranging device 10 with the thin solid line. It is a histogram obtained by.
  • Each histogram shown in FIG. 5 represents a histogram when an object outside the distance measuring device 10 is located in the vicinity of the distance measuring device 10 (window 92).
  • the composite mountain portion mp0 appears in the histogram hr1.
  • the histogram hr11 and the histogram hr12 overlap each other in time. Then, the light receiving intensity (frequency) of the peak of the histogram hr11 and the light receiving intensity (frequency) of the peak of the histogram hr12 are substantially equal to each other. Therefore, as shown in FIG. 5, the composite mountain portion mp0 appears as a single mountain-like shape.
  • the composite mountain portion estimation unit 240 is used when the time between the rise time and the fall time in the mountain portion, that is, the time when the light receiving intensity becomes equal to or higher than the first threshold value described later in the histogram, is continuous for a predetermined time or longer. It is estimated that there is a compound mountain in the histogram.
  • each value of the light receiving intensity I shown in FIG. 5 that is, the noise intensity I1, the first threshold light receiving intensity I3, and the peak intensity I4 will be described later.
  • the details of each value of the flight time Tf shown in FIG. 5, that is, the rise time Tu1, the fall time Td1, the peak time Tp1, and the intermediate time Tc1 will be described later.
  • the flight time specifying unit 250 specifies two flight times (rise time and fall time, which will be described later) whose light receiving intensity coincides with the first threshold light receiving intensity in the histogram generated by the histogram generation unit 230.
  • the first threshold light receiving intensity is a threshold value of the light receiving intensity used when calculating the target distance, and is used to determine the rising time and the falling time in the mountain portion of the histogram.
  • the first threshold light receiving intensity I3 is set.
  • the first threshold light receiving intensity is set as a value represented by the following equation (1).
  • First threshold light reception intensity (peak intensity-noise intensity) x 0.4 + noise intensity ... (1)
  • the peak intensity means the maximum light receiving intensity in the mountain portion
  • the peak intensity I4 corresponds to the composite mountain portion mp0 in FIG.
  • the noise intensity means the light receiving intensity of light other than the reflected light of the pulsed light (hereinafter, referred to as "noise light").
  • Examples of such noise light include sunlight, reflected light obtained by reflecting sunlight on an external object, and light from a street lamp.
  • the noise intensity can be obtained as an average value of the light receiving intensity measured at the timing when the light emitting unit 40 does not emit the pulsed light in a predetermined time.
  • the (peak intensity-noise intensity) in the formula (1) is also called the reflected light intensity.
  • the equation (1) means that the light receiving intensity obtained by adding the value of 40% of the reflected light intensity to the noise intensity is set as the first threshold light receiving intensity.
  • an arbitrary ratio smaller than 40% or larger than 40% may be used.
  • the flight time specifying unit 250 specifies the time Tu1 earlier in time as the rise time Tu1 among the two time Tu1 and Td1 in which the histogram hr1 and the first threshold light receiving intensity I3 match. Further, the later time Td1 is specified as the fall time Td1.
  • the basic flight time determination unit 260 determines the basic flight time (hereinafter referred to as "basic flight time") for calculating the target distance.
  • the base flight time determination unit 260 determines the fall time as the base flight time, as will be described later. The reason for determining the fall time as the base flight time will be described later.
  • the distance calculation unit 270 calculates the object distance using the base flight time determined by the base flight time determination unit 260. The details of the method of calculating the object distance will be described later.
  • the calculated object distance is used, for example, in an estimation process of whether or not a collision with an object occurs in a vehicle equipped with a distance measuring device 10, and in steering and brake control for collision avoidance. obtain.
  • the distance measuring process shown in FIG. 6 is periodically executed at predetermined time intervals for each pixel position.
  • the light emitting unit 40 emits pulsed light, the reflected light is received by the light receiving unit 60, and the distance measuring process is started at a timing after the histogram is generated by the histogram generation unit 230.
  • the composite mountain portion estimation unit 240 specifies the composite mountain portion in the generated histogram, and the flight time specifying unit 250 specifies the rise time and the fall time (step S105). Specifically, the composite mountain portion estimation unit 240 identifies the composite mountain portion by specifying a portion in the histogram in which the flight time exceeding the first threshold light receiving intensity is continuous for a predetermined time or longer. Further, the flight time specifying unit 250 specifies two flight times that match the first threshold value in the compound mountain portion as the rise time and the fall time.
  • the composite mountain portion mp0 is specified, and of the two time Tu1 and Td1 in which the composite mountain portion mp0 and the first threshold light receiving intensity I3 match, the earlier time Tu1 is the rise time Tu1.
  • the later time Td1 is specified as the fall time Td1.
  • the base flight time determination unit 260 determines whether or not the light receiving position of the reflected light coincides with the position of the clutter (step S110).
  • the clutter that is, the reflected light of the pulsed light by the window 92 or the inner wall surface of the housing always enters the light receiving unit 60 from the same position (direction). Further, the distance to the window 92 and the inner wall surface of the housing is always constant. Therefore, it is possible to specify in advance by experiment or simulation which direction the pixel position and how long the flight time the reflected light of the pulsed light from the window 92 or the inner wall surface of the housing is specified.
  • the pixel position and flight time of the reflected light of the pulsed light from the window 92 and the inner wall surface of the housing are stored in the memory 290 in advance, and the pixel position and the composite mountain portion when the composite mountain portion is detected are stored in the memory 290.
  • the flight time of the light matches the pixel position and flight time stored in the memory 290, it is determined whether or not the light receiving position of the reflected light matches the position of the clutter.
  • the distance calculation unit 270 calculates the distance by a normal method and specifies it as the object distance (step S115). The distance measuring method in step S115 will be described with reference to FIG. 7.
  • the vertical axis and the horizontal axis in FIG. 7 are the same as those in FIG. 5, detailed description thereof will be omitted.
  • the mountain portion mp1 shown in FIG. 7 for example, an object exists outside the distance measuring device 10 in a direction different from the direction in which the clutter is generated, and the reflected light of the pulsed light is received by the light receiving unit 60 from the object. Can occur in such cases.
  • the intermediate time between the rise time Tu2 and the fall time Td2 is the peak time Tp2 of the mountain part mp1 by the flight time specifying part 250. Specified as.
  • the basic flight time determination unit 260 determines this peak time Tp2 as the basic flight time
  • the distance calculation unit 270 obtains the distance from the distance measuring device 10 to these a plurality of objects by using the peak time Tp2. ..
  • the base flight time determination unit 260 uses the fall time specified in step S105 as the basis. Determined as the flight time (step S120).
  • the distance calculation unit 270 calculates the distance based on the base flight time determined in step S120 (step S125). Specifically, for example, when the fall time Td1 of the composite mountain portion mp0 shown in FIG. 5 is determined as the basic flight time, the distance calculation unit 270 traces back from the fall time Td1 by a predetermined time ⁇ t. The time tp1 is specified as the peak time tp1. Then, the distance calculation unit 270 calculates the target distance using this peak time tp1. More specifically, the target distance is calculated assuming that the peak time tp1 is the flight time Tf of the pulsed light and the reflected light thereof.
  • the time (peak time Tp1) that goes back by the time ⁇ t from the time Td1 almost coincides with the peak time of the histogram hr12.
  • the intermediate time Tc1 between the rise time Tu1 and the fall time Td1 is largely deviated from the peak time of the histogram hr1. Therefore, it can be seen that by calculating the distance with the peak time Tp1 as the base flight time as described above, the distance to the object corresponding to the histogram hr1 can be calculated more accurately than the distance calculated by the usual method. .. As described above, the reason why the intermediate time Tc1 deviates significantly from the peak time of the histogram hr1 will be described.
  • the composite mountain portion mp0 is generated by synthesizing the histogram hr12 by the reflected light output from the reflector OBJ located farther away and the histogram h12 by the clutter output from the window 92 closer to it.
  • the rising position Tu1 substantially coincides with the rising position of the histogram hr11 whose peak position is earlier in time.
  • the falling position of the histogram hr11 substantially coincides with the falling position of the histogram h12 due to the reflected light from the object existing farther (longer flight time).
  • the falling position of the histogram hr11 that is ahead in time is largely deviated from the falling position of the histogram hr1 due to the presence of the histogram hr12 that is later in time.
  • the rising position of the histogram hr12, which is later in time is largely deviated from the rising position of the histogram hr1 due to the presence of the histogram hr11, which is earlier in time. Therefore, even if the intermediate time Tc1 which is the intermediate time between the rising time Tu1 and the falling time Td1 is specified in order to calculate the target distance according to a usual method, the intermediate time Tc1 is the peak time and the histogram of the histogram hr11.
  • the peak time of the object is set by setting the time retroactive by the time ⁇ t from the fall time as the reference, based on the fall time of the histogram hr1, which almost coincides with the fall time of the object. It is possible to obtain a distance that is almost the same as the distance calculated when is used as the base flight time.
  • the distance calculation unit 270 determines whether or not the distance calculated by step S125 is smaller than the predetermined threshold distance (step S130).
  • the threshold distance in step S130 is set as a value obtained by adding a predetermined value to the distance from the light emitting unit 40 and the light receiving unit 60 to the window 92 as a distance recognized as not being a clutter if it is larger than such a distance.
  • This step S130 is executed to confirm that the received reflected light is certainly not due to the clutter because the distance calculated based on the base flight time is larger than the distance calculated in the case of the clutter.
  • step S130: NO When it is determined that the distance calculated by step S125 is not smaller than the threshold distance (step S130: NO), the above-mentioned step S115 is executed, and the object distance is calculated by a normal method.
  • step S130: YES when it is determined that the distance calculated by step S120 is smaller than the threshold distance (step S130: YES), the distance calculation unit 270 specifies the distance calculated by step S125 as the object distance. (Step S135). After the completion of step S115 or step S135 described above, the distance measuring process at the corresponding pixel position is completed, and the distance measuring process for the next pixel position is started.
  • the distance measuring device 10 of the first embodiment described above when it is estimated that there is a composite mountain portion including a clutter which is reflected light of pulsed light inside the housing such as a window 92 and an inner wall surface of the housing. Since the base flight time is determined based on the fall time, the object distance can be measured accurately when a clutter occurs. Generally, the reflected light of the pulsed light inside the window 92 or the housing has a shorter flight time than the object existing outside the distance measuring device. In other words, the reflected light of the pulsed light by the object existing outside has a longer flight time than the reflected light of the pulsed light inside the housing.
  • the distance measuring device 10 of the present embodiment it is possible to improve the distance measuring accuracy when it is estimated that the compound mountain portion exists.
  • the intermediate time between the rise time and the fall time is determined as the base flight time. Therefore, even in such a case, either the rise time or the fall time is determined.
  • the object distance can be measured more accurately than the configuration in which the base flight time is determined based on the sword.
  • Second embodiment The configuration of the distance measuring device 10a of the second embodiment shown in FIG. 8 is different from the configuration of the distance measuring device 10 of the first embodiment in that the calculation unit 20 additionally includes the intermediate time specifying unit 235. Since the other configurations of the distance measuring device 10a are the same as those of the distance measuring device 10 of the first embodiment, the same components are designated by the same reference numerals, and detailed description thereof will be omitted.
  • the intermediate time specifying unit 235 coincides with the intermediate time of two flight times (hereinafter referred to as "first intermediate time”) that coincides with the first threshold light receiving intensity and the second threshold light receiving intensity in the compound mountain portion.
  • the intermediate time of one flight time (hereinafter referred to as "second intermediate time") is specified.
  • the first threshold light receiving intensity is the same as the first threshold light receiving intensity of the first embodiment.
  • the second threshold light receiving intensity is a light receiving intensity smaller than that of the first threshold light receiving intensity, and the noise intensity is also larger.
  • the second threshold light receiving intensity is set as a value represented by the following equation (2).
  • Second threshold light reception intensity (peak intensity-noise intensity) x 0.1 + noise intensity ... (2)
  • the distance measurement process of the second embodiment shown in FIG. 9 is the first embodiment shown in FIG. 6 at the point where step S110 is omitted and the point where steps S106, S107, S108, and S123 are additionally executed. It is different from the distance measurement process of. Since the other procedures in the distance measuring process of the second embodiment are the same as those of the distance measuring process of the first embodiment, the same procedures are designated by the same reference numerals, and detailed description thereof will be omitted.
  • the intermediate time specifying unit 235 calculates the first intermediate time and the second intermediate time (step S106).
  • each object OBJ1 and OBJ2 have regions Ar1 and Ar2 exposed from the distance measuring device 10a.
  • the reflected light output from each of these two regions Ar1 and Ar2 enters the light receiving unit 60 of the distance measuring device 10a.
  • the object OBJ1 is located on the front side when viewed from the distance measuring device 10a
  • the object OBJ2 is located on the back side when viewed from the distance measuring device 10a.
  • the exposed area of the region Ar1 seen from the distance measuring device 10a is larger than the exposed area of the region Ar2.
  • the histogram hr3 generated in the situation of FIG. 10 has a complex mountain portion mp3 formed by synthesizing a large mountain portion and a small mountain portion having a peak later in time.
  • the composite mountain portion mp3 is a histogram obtained by synthesizing the histogram hr31 by the reflected light from the region Ar1 shown by the thin alternate long and short dash line and the histogram hr32 by the reflected light from the region Ar2 shown by the thin solid line.
  • the region Ar1 is located closer to the surface and has a larger exposed area.
  • the region Ar2 is located deeper and has a smaller exposed area. Therefore, the mountain portion of the histogram hr31 is located in front of the mountain portion of the histogram hr32 in time, and the peak light receiving intensity is large. Therefore, the composite mountain portion mp3 has a distorted mountain shape having a relatively large peak at the front position in time and a relatively small peak at the rear position in time.
  • the intermediate time specifying unit 235 has two times (rise time Tu3 and fall time) that coincide with the first threshold light receiving intensity I11 in the composite mountain portion mp3. Td3) is specified, and the intermediate time Tc31 between these two times is specified. Further, the intermediate time specifying unit 235 identifies two times (rising time Tu4 and falling time Td4) that match the second threshold light receiving intensity I12 in the composite mountain portion mp3, and the intermediate time Tc32 of these three times. To identify.
  • the intermediate time of the two times matching the first threshold light receiving intensity is referred to as the first intermediate time
  • the intermediate time of the two times matching the second threshold light receiving intensity is referred to as the second intermediate time.
  • the intermediate time Tc31 is also referred to as the first intermediate time Tc31
  • the intermediate time Tc32 is also referred to as the second intermediate time Tc32.
  • the composite mountain portion estimation unit 240 determines whether or not the width (hereinafter referred to as “root width”) between the two times corresponding to the second threshold light receiving intensity is larger than the threshold width. Determination (step S107).
  • the width between the two times Tu4 and Td4 that coincide with the second threshold light receiving intensity I12 corresponds to the root width dt31.
  • the root width is large in the composite mountain portion obtained by synthesizing the mountain portions included in the histogram of the reflected light from a plurality of objects, such as the composite mountain portion mp3.
  • the root width of the mountain part due to the reflected light from a single object is small.
  • the threshold width for identifying the complex mountain portion is specified and set in advance by an experiment or the like.
  • the above-mentioned step S115 is executed. In this case, since the reflected light from a single object is incident, the distance is calculated by a usual method and specified as the object distance.
  • step S107 When it is determined that the root width is larger than the threshold width (step S107: YES), in this case, there is a high possibility that the mountain part is a compound mountain part. Then, in this case, the basic flight time determination unit 260 determines whether or not the first intermediate time is equal to or less than the second intermediate time (step S108).
  • the first intermediate time Tc31 is smaller (faster) than the second intermediate time Tc32. Further, in the example of FIG. 13, the first intermediate time and the second intermediate time coincide with each other.
  • the histogram hr5 shown in FIG. 13 is a histogram obtained by synthesizing the histogram hr51 and the histogram hr52.
  • the histogram hr51 is due to the reflected light from the object OBJ1 shown in FIG. 10, and the histogram hr52 is due to the reflected light from the object OBJ2 shown in FIG.
  • unlike the example of FIG. 13 unlike the example of FIG.
  • the exposed area of the region Ar1 seen from the distance measuring device 10a is equal to the exposed area of the region Ar2.
  • the peak light receiving intensity of the histogram hr51 due to the reflected light from the object OBJ1 becomes equal to the peak light receiving intensity of the histogram hr52 due to the reflected light from the object OBJ2. Therefore, the composite mountain portion mp5 has a large mountain shape with one peak. Then, in this example, both the first intermediate time and the second intermediate time are the intermediate time Tc51.
  • the first intermediate time Tc31 is determined to be smaller (earlier) than the second intermediate time Tc32.
  • the base flight time determination unit 260 determines the rise time specified in step S105. Is determined as the base flight time (step S123). For example, in the example of FIG. 11, the rise time Tu3 is determined as the base flight time. Further, in the example of FIG. 13, the rise time Tu7 is determined as the base flight time.
  • the histogram hr4 shown in FIG. 12 is a histogram obtained by synthesizing the histogram hr41 and the histogram hr42.
  • the histogram hr41 is due to the reflected light from the object OBJ1 shown in FIG. 10, and the histogram hr42 is due to the reflected light from the object OBJ2 shown in FIG.
  • the exposed area of the region Ar1 seen from the distance measuring device 10a at the pixel position 81 is smaller than the exposed area of the region Ar2.
  • the composite mountain portion mp4 has a distorted mountain shape having a relatively small peak at the front position in time and a relatively large peak at the rear position in time.
  • the composite mountain portion mp4 in the composite mountain portion mp4, two times (rising time Tu5 and falling time Td5) that coincide with the first threshold light receiving intensity I21 are specified, and an intermediate time between these two times is specified.
  • the second intermediate time Tc41 is specified.
  • two times (rising time Tu6 and falling time Td6) corresponding to the second threshold light receiving intensity I22 are specified, and the second intermediate time Tc42 which is an intermediate time between these two times is specified. Is identified.
  • the first intermediate time Tc41 is larger (slower) than the second intermediate time Tc42.
  • step S108 when it is determined that the first intermediate time is not equal to or less than the second intermediate time (step S108: NO), the above-mentioned step S120 is executed, and the fall time specified in step S105 is the base flight. Determined as time. In the example of FIG. 12, the fall time Td5 is determined as the base flight time.
  • step S120 or step S123 the above-mentioned steps S125 and S135 are executed.
  • the exposed area of OBJ2 is large. Then, in this case, since the object distance is calculated using the flight time of the reflected light of the pulsed light on the surface of the object OBJ2 as the base flight time, this case also has a greater effect on the vehicle equipped with the distance measuring device 10a. The object distance is calculated based on the large object of.
  • the distance measuring device 10a of the second embodiment described above has the same effect as the distance measuring device 10 of the first embodiment.
  • the base flight time is determined based on the rise time, so that a plurality of objects having the same direction and different distances with respect to the distance measuring device 10a are used as a reference. Is present, and when the irradiation area of the pulsed light in that direction is larger than that of the object on the front side, the distance to the object on the front side can be obtained as the object distance.
  • the base flight time is determined based on the fall time, so that a plurality of objects having different distances in the same direction with respect to the distance measuring device 10a are used as a reference. Is present, and when the irradiation area of the pulsed light in that direction is larger than that of the object on the rear side, the distance to the object on the rear side can be obtained as the object distance.
  • the distance measuring device 10a of the second embodiment when there are a plurality of objects having the same direction and different distances with respect to the distance measuring device 10a, either the front or the rear object. The distance to the object, that is, the distance to the position where the object actually exists can be obtained as the object distance.
  • the distance to the object having a larger pulsed light irradiation area is obtained as the object distance
  • the distance to the object having a larger exposed area when viewed from the distance measuring device 10a can be obtained as the object distance. .. Therefore, for example, in a configuration in which the safety function of the vehicle is exerted by utilizing the object distance, such a function can be exerted so as to ensure safety.
  • the base flight time is determined based on the rise time, so that a plurality of objects having the same direction and different distances with respect to the distance measuring device 10a are present.
  • the distance to any of the front and rear objects can be obtained as the object distance. Therefore, it is possible to suppress that the distance to the position where the object does not exist is obtained as the object distance.
  • the basic flight time determination unit 260 is specified based on the first threshold light receiving intensity for the composite mountain portion when the composite mountain portion is present.
  • the base flight time can be determined based on either the rise time or the fall time.
  • the same components are designated by the same reference numerals, and detailed description thereof will be given. Omit.
  • the threshold light receiving intensity used for specifying the rise time Tu3 and the time Td3 and the intermediate time (first intermediate time) Tc31 is common to the first threshold light receiving intensity I11.
  • the threshold light receiving intensity used to specify the rise time Tu3 and the time Td3 and the threshold light receiving intensity used to specify the intermediate time (first intermediate time) are defined as the threshold light receiving intensity. Different from each other.
  • a specific description will be given with reference to FIG.
  • the rise time Tu3 and the fall time Td3 are specified by using the first threshold light receiving intensity I11, as in the second embodiment.
  • the threshold light receiving intensity I14 for specifying the intermediate time Tc33 is set separately from the above-mentioned first threshold light receiving intensity I11. Then, two times corresponding to the threshold light receiving intensity I14 are specified in the composite mountain portion mp3, and the intermediate time between these two times is specified as the intermediate time Tc33.
  • This intermediate time Tc33 is used as the first intermediate time instead of the intermediate time Tc31 of the second embodiment. Therefore, for example, in step S108 of the distance measuring process shown in FIG.
  • the threshold light receiving intensity I14 is set in advance as a value larger than the first threshold light receiving intensity I11 and equal to or less than the peak light receiving intensity I13 of the composite mountain portion mp3.
  • the threshold light receiving intensity I14 may be set as a value smaller than the first threshold light receiving intensity I11 and larger than the second threshold light receiving intensity I12.
  • the distance measuring device 10a of the third embodiment described above has the same effect as the distance measuring device 10a of the second embodiment.
  • the threshold light receiving intensity I14 in the third embodiment corresponds to the third threshold light receiving intensity in the present disclosure.
  • D. Fourth Embodiment D1. Device configuration: Since the configuration of the distance measuring device 10 of the fourth embodiment is the same as that of the distance measuring device 10 of the first embodiment, the same reference numerals are given to the same configurations, and detailed description thereof will be omitted. As described in the first embodiment, regarding the clutter generated due to the reflected light of the pulsed light by the window 92 and the inner wall surface of the housing, it is determined in which direction and at what flight time the clutter is specified. It can be specified in advance by experiments and simulations. On the other hand, when foreign matter such as sand or mud adheres to the outer surface of the window 92, the reflected light of the pulsed light at such a position becomes very large and may be newly generated as a clutter.
  • the pixel position where the clutter is generated differs depending on the position where the foreign matter adheres, it cannot be specified in advance.
  • the position of the clutter (the pixel position in which direction) caused by the adhesion of foreign matter or the like is specified and set as described above. do.
  • Clutter position setting process The clutter position setting process shown in FIG. 15 is executed when the power of the distance measuring device 10 is turned on.
  • the histogram generation unit 230 determines whether or not the identification of the object distance is completed by the distance measurement process shown in FIG. 6 (step S205). When the identification of the object distance is completed, the distance from the distance measuring device 10 to the object for each pixel is specified.
  • the histogram generation unit 230 generates a histogram of the distance for each pixel including the distance specified by the current distance measurement process (step S210).
  • the memory 290 stores the specified object distance for a predetermined number of times for each pixel. The predetermined number of times depends on the size of the memory 290, but for example, a maximum of 100 times may be stored. Then, when the object distance is newly specified in the state where the object distance for 100 times is stored, the newly specified object distance is stored by overwriting the oldest stored object distance. You may.
  • the distance histogram for each pixel generated in step S210 corresponds to the "statistical value of the object distance calculated for each pixel" in the present disclosure.
  • two histograms hd1 and hd2 are shown as an example of the histogram generated in step S210.
  • the histogram hd1 shown by a thick solid line is a histogram for pixels different from the positions where clutter is known to occur in advance, and shows histrum when foreign matter such as mud adheres to the window 92.
  • the thin broken line histogram hd2 is a histogram for pixels different from the positions where clutter is known to occur in advance, and shows a histogram when foreign matter such as mud does not adhere to the window 92.
  • a frequency peak occurs at a distance corresponding to the position of the window 92.
  • the distance based on the reflected light blocked by the foreign matter is not specified, and the frequency at such distance remains low. It becomes. Therefore, as shown in the histogram hd1, a peak occurs at a relatively small distance d1.
  • the frequency becomes greater than 0 over a wide range of distances, and for example, a gentle hill-shaped histogram is obtained.
  • the histogram generation unit 230 determines whether or not there are pixels having a frequency exceeding the threshold frequency within the range within the threshold distance (step S215).
  • the frequency exceeding the threshold frequency THn is recorded in the range within the threshold distance THd, so that the pixel corresponding to the histogram hd1 sets the threshold frequency within the range within the threshold distance.
  • the histogram hd2 of FIG. 16 is obtained, the frequency exceeding the threshold frequency THn is not recorded in the range within the threshold distance THd, so that the pixel corresponding to the histogram hd2 has the threshold within the range within the threshold distance. It does not correspond to a pixel whose frequency exceeds the frequency.
  • step S215 when it is determined that there is a pixel having a frequency exceeding the threshold frequency in the range within the threshold distance (step S215: YES), the histogram generation unit 230 is the pixel, that is, the range within the threshold distance.
  • the position of the pixel whose frequency exceeds the threshold frequency is additionally set as the clutter position (step S220).
  • the histogram generation unit 230 reads out the information on the clutter position, adds the pixel position specified in step S215, and stores it in the memory 290. .. After the completion of step S220, the process returns to step S205.
  • step S215 when it is determined that there are no pixels having a frequency exceeding the threshold frequency within the range within the threshold distance (step S215: NO), the histogram generation unit 230 does not position all the pixels in the clutter position. Additional settings (step S225). After the completion of step S225, the process returns to step S205.
  • a wiper that performs an operation of wiping the window 92 is provided, and when it is determined in step S215 that there are pixels having a frequency exceeding the threshold frequency within the range within the threshold distance, the window 92 is wiped by the wiper. May be.
  • step In S215 it is determined that there are no pixels having a frequency exceeding the threshold frequency within the range within the threshold distance, and by executing step S225, the clutter position derived from the foreign matter is updated and set as not the clutter position.
  • the distance calculation unit 270 or another functional unit may execute the above steps S205 to S220.
  • a new functional unit other than the functional unit shown in FIG. 4 may execute the above-mentioned steps S205 to S220.
  • the functional unit that executes the above steps S205 to S220 corresponds to the "clutter position setting unit" of the present disclosure.
  • the distance measuring device 10 of the fourth embodiment described above has the same effect as the distance measuring device 10 of the first embodiment.
  • the distance measuring device 10 of the fourth embodiment generates a histogram of the distance of each pixel, and when it is determined that there are pixels having a frequency exceeding the threshold frequency THn within the range within the threshold distance THd, such pixels. Since the position of is set as the clutter position, it is possible to specify the position of the clutter generated by the foreign matter such as mud adhering to the window 92. Therefore, in step S110 of the distance measuring process shown in FIG. 6, it is possible to accurately determine whether or not the position of the reflected light and the clutter position match.
  • E. Fifth Embodiment Since the configuration of the distance measuring device 10 of the fifth embodiment is the same as that of the distance measuring device 10 of the fourth embodiment, the same reference numerals are given to the same configurations, and detailed description thereof will be omitted.
  • the distance measuring device 10 of the fifth embodiment is different from the distance measuring device 10 of the fourth embodiment in the procedure of the clutter position setting process.
  • step S205a is executed instead of step S205
  • step S210a is executed instead of step S210
  • step S215a are executed instead of step S215.
  • the histogram generation unit 230 determines whether or not the histogram has been generated (step S205a). If it is determined that the histogram has not been generated (step S205a: NO), the process returns to step S205a. That is, the histogram generation unit 230 waits until the histogram has been generated.
  • the histogram generation unit 230 When it is determined that the histogram has been generated (step S205a: YES), the histogram generation unit 230 generates a histogram of the flight time and the average intensity for each pixel including the histogram calculated this time (step S210a). ).
  • the histogram for each pixel generated by the histogram generation unit 230 is stored in the memory 290. Then, in step S210a, the histogram generation unit 230 obtains the average value of the intensities for each flight time Tg for the plurality of histograms stored in the memory 290, and thereby generates the histogram of the average intensities.
  • the histogram of the average intensity generated for each pixel corresponds to the "statistical value of the light receiving intensity of each pixel" in the present disclosure.
  • two histograms hr61 and hr62 are shown as an example of the histogram generated in step S210a.
  • the histogram hr61 shown by a thick solid line is a histogram for pixels different from the positions where clutter is known to occur in advance, and shows histrum when foreign matter such as mud adheres to the window 92.
  • the thin broken line histogram hr62 is a histogram for pixels different from the positions where clutter is known to occur in advance, and shows a histogram when foreign matter such as mud does not adhere to the window 92.
  • the object when no foreign matter adheres to the window 92, if the vehicle is moving, the object may exist at various distances in the direction corresponding to the pixel. Therefore, as shown in the histogram hr62, the light receiving intensity I becomes larger than 0 in a wide range of flight time Tf, and for example, a gentle hill-shaped histogram is obtained.
  • the histogram generation unit 230 determines whether or not there is a pixel having a peak whose average intensity exceeds the threshold intensity within the range within the threshold flight time (step S215a).
  • the histogram hr61 of FIG. 18 is obtained, the flight time T6 exceeding the threshold intensity It1 is recorded in the range within the threshold flight time THt, so that the pixel corresponding to the histogram hr61 is within the threshold flight time THt. It corresponds to a pixel having a frequency exceeding the threshold intensity It1 in the range.
  • the frequency of exceeding the threshold intensity It1 is not recorded in the range within the threshold flight time THt, so that the pixel corresponding to the histogram hr62 is within the threshold flight time THt. It does not correspond to a pixel having an intensity exceeding the threshold intensity It1 in the range.
  • the histogram generation unit 230 determines the pixel, that is, within the threshold flight time THt.
  • the position of the pixel having the peak whose average intensity exceeds the threshold intensity It1 is additionally set as the clutter position in the range of (step S220).
  • step S215a determines whether there is no pixel having a peak whose average intensity exceeds the threshold intensity It1 within the range within the threshold flight time THt. If it is determined in step S215a that there is no pixel having a peak whose average intensity exceeds the threshold intensity It1 within the range within the threshold flight time THt (step S215a: NO), the process returns to step S205a. Therefore, in this case, the clutter position stored in the memory 290 is not updated.
  • the distance measuring device 10 of the fifth embodiment described above has the same effect as the distance measuring device 10 of the fourth embodiment.
  • the ranging device 10 of the fifth embodiment generates a histogram of the flight time and the average intensity for each pixel, and it is determined that there is a pixel having an intensity exceeding the threshold intensity It1 within the range within the threshold flight time THt.
  • the position of the pixel is set as the clutter position, the position of the clutter generated by the foreign matter such as mud adhering to the window 92 can be specified. Therefore, in step S110 of the distance measuring process shown in FIG. 6, it is possible to accurately determine whether or not the position of the reflected light and the clutter position match.
  • the pulse width of the compound mountain portion is preset in the distance calculation unit 270. It is determined whether or not it is smaller than the threshold width (step S132).
  • the pulse width of step S132 means the first threshold light receiving intensity, that is, the width of the flight time used to determine the rise time and the fall time in the mountain portion.
  • the threshold width of step S132 is either a mountain portion formed by the clutter and the reflected light from the target object or a mountain portion formed by multiple reflection described later. It means the range of flight time set to discriminate between.
  • the width of the mountain portion formed by multiple reflections is larger than the width of the mountain portion formed by the clutter and the reflected light from the target object. Then, the discriminable threshold width is determined and set by, for example, an experiment or a simulation.
  • the above-mentioned “multiple reflection” means pulsed light and applied between a distance measuring device 10 and an object having high reflectance (hereinafter referred to as "high reflectance object”) such as a reflector provided at the rear of the vehicle. It means that the reflected light reciprocates multiple times and reflects each other. When a high-reflectance object is irradiated with pulsed light, the intensity of the reflected light is high.
  • such reflected light may be reflected inside the window 92 or the housing 90, and the reflected light may irradiate the high-reflectance object again to generate reflected light toward the ranging device 10 again.
  • a plurality of histograms are formed so that the flight time Tf becomes longer and the light receiving intensity becomes smaller as the reflection occurs back and forth between the high reflectance object and the distance measuring device 10. An integrated histogram will be generated.
  • the histogram hr71 obtained by the first reflected light, the histogram hr72 obtained by the second reflected light, and the histogram hr73 obtained by the third reflected light are combined to form a composite mountain.
  • a wide histogram hr7 showing the part mp7 is obtained.
  • the timely small flight time Tu7 is specified as the rise time
  • the timely large flight time Td7 is specified. Is specified as the fall time.
  • step S132 shown in FIG. 19 when it is determined that the pulse width of the compound mountain portion is smaller than the preset threshold width (step S132: YES), the above-mentioned step S135 is executed.
  • the pulse width of the compound mountain portion is smaller than the preset threshold width, it is highly possible that the compound mountain portion is composed of the clutter and the reflected light from the target object. Therefore, in this case, as in the first embodiment, the distance calculated by step S125 is specified as the object distance.
  • step S132 when it is determined that the pulse width of the compound mountain portion is not smaller than the preset threshold width (step S132: NO), the above-mentioned step S115 is executed, and the object distance is obtained by a normal method. Is calculated. That is, the intermediate time between the rise time Tu7 and the fall time Td7 is specified as the peak time, and the object distance is specified using the peak time.
  • the pulse width of the compound mountain portion is not smaller than the preset threshold width, that is, when the pulse width of the compound mountain portion is large and there is a high possibility that multiple reflections are occurring, the first step is made. It is desirable to find the distance based on the reflected light.
  • step S135 since the distance calculated with the falling position as the base flight time is specified as the object distance, the identification accuracy may be lowered. This is because the fall time, which is increased due to the long flight time due to repeated multiple reflections, becomes the basic flight time. Therefore, in the present embodiment, in such a case, step S135 is not executed, and the distance is calculated by a normal method and treated as the object distance.
  • the distance measuring device 10 of the sixth embodiment described above has the same effect as the distance measuring device 10 of the first embodiment.
  • the object distance can be specified more accurately than the configuration in which the object distance is calculated and specified using the fall time as the base flight time.
  • the configuration of the distance measuring device 10a of the sixth embodiment is the same as that of the distance measuring device 10a of the second embodiment except for the detailed procedure of the distance measuring process, the same components are designated by the same reference numerals. A detailed description will be omitted.
  • the distance measuring process of the seventh embodiment shown in FIG. 21 is different from the distance measuring process of the second embodiment shown in FIG. 9 in that step S117 is executed instead of steps S108, S120, and S123. Since the other procedures of the distance measuring process of the seventh embodiment are the same as those of the second embodiment, the same procedures are designated by the same reference numerals and detailed description thereof will be omitted.
  • step S107 when it is determined in step S107 that the root width is larger than the threshold width (step S107: YES), the base flight time determination unit 260 bases the rise time and the fall time, respectively.
  • the flight time is determined (step S117).
  • step S117 the rise time and the fall time are each determined as the base flight time. Therefore, in step S125, the distance is calculated based on each base flight time, and in step S135, the calculated two distances are calculated. Both will be specified as the object distance. As shown in FIG.
  • the object distance for the object located on the front side can be accurately obtained. .. Further, since the distance is obtained based on the fall time, that is, the time derived from the object located on the back side, the object distance for the object located on the back side can be obtained accurately.
  • the distance measuring device 10a of the seventh embodiment described above has the same effect as the distance measuring device 10a of the second embodiment.
  • the root width is larger than the threshold width, that is, there is a high possibility that the mountain portions included in the histogram of the reflected light from a plurality of objects are combined to obtain a composite mountain portion.
  • the rise time and the fall time are each determined by the base flight time, the distances to the two objects existing in front and behind can be accurately specified.
  • step S110: NO When it is presumed that the compound mountain portion does not exist, that is, when it is determined that the light receiving position of the reflected light does not match the position of the clutter in the first embodiment (step S110: NO), the first embodiment is also performed.
  • step S115 when it is determined that the root width is not larger than the threshold width (step S107: NO), step S115 is executed based on the intermediate time between the rise time and the fall time of the mountain portion.
  • the object distance was calculated as the flight time, but the present disclosure is not limited to this. In such a case, the object distance may be calculated by using either the rising time or the falling time of the mountain part as the base flight time.
  • the base flight time is set according to the comparison result between the first intermediate time and the second intermediate time.
  • the rise time or fall time was determined, but the present disclosure is not limited to this.
  • the rise time may be determined as the base flight time without performing the comparison between the first intermediate time and the second intermediate time. According to such a configuration, the distance to the object located on the front side of the distance measuring device 10a can be easily obtained as the object distance, and the object located closer to the object having a greater influence on the vehicle can be obtained. You can find the distance.
  • the base flight time determination unit 260 determines the rise time as the base flight time when the first intermediate time and the second intermediate time are equal to each other. Instead of time, the fall time may be determined as the base flight time.
  • step S130 may be omitted.
  • step S135 the distance may be calculated again based on the base flight time and specified as the object distance.
  • the distance measuring devices 10 and 10a may be mounted not only on a vehicle but also on an arbitrary moving body such as an airplane or a ship. Alternatively, it may be fixedly installed and used for purposes such as security. Further, the housings of the distance measuring devices 10 and 10a may have no window member and may be simply provided with an opening.
  • the arithmetic unit 20 and methods thereof described in the present disclosure are provided by configuring a processor and memory programmed to perform one or more functions embodied by a computer program. It may be realized by a dedicated computer. Alternatively, the arithmetic unit 20 described in the present disclosure and these methods may be realized by a dedicated computer provided by configuring a processor with one or more dedicated hardware logic circuits. Alternatively, the arithmetic unit 20 and these methods described in the present disclosure include a processor and memory programmed to perform one or more functions and a processor configured by one or more hardware logic circuits. It may be realized by one or more dedicated computers configured by a combination. Further, the computer program may be stored in a computer-readable non-transitional tangible recording medium as an instruction executed by the computer.
  • the distance of the base object is specified according to whether or not the vehicle is running. You may weight the frequency you want to do. Specifically, when the vehicle is running, a larger frequency may be accumulated as compared with the case where the vehicle is not running (that is, when the vehicle is stopped). When the vehicle is stopped, the distance to the surrounding objects changes little over time, and when the object is located at a short distance, the histogram is high in the short distance range. Is likely to be created.
  • a larger frequency may be accumulated when the vehicle is running than when the vehicle is not running (that is, when the vehicle is stopped).
  • the clutter position setting process is performed only when the object distance is specified during the vehicle, and the vehicle is stopped. The clutter position setting process may not be performed.
  • the "threshold width" used in step S132 of the sixth embodiment is one fixed value specified in advance by an experiment or the like, but the present disclosure is not limited to this.
  • the threshold width may be changed according to the fall time of the compound mountain portion. Specifically, the larger (longer) the falling time of the compound mountain portion, the larger the threshold width may be set.
  • the fall time of the compound mountain part is long (long)
  • it is assumed that the number of multiple reflections is large. Therefore, in such a case, by setting a larger threshold width, multiple reflections occur when a plurality of objects having the same direction but different distances exist in a relatively long range in the depth direction. It is possible to accurately discriminate from the case where there is.
  • step S115 when it is determined in step S132 that the pulse width of the composite mountain portion is not smaller than the preset threshold width (step S132: NO), step S115 is executed and is normally performed.
  • the distance was calculated with the intermediate time as the base flight time, but the present disclosure is not limited to this. In this case, unlike the usual method, the distance may be calculated using the rise time as the base time.
  • the distance calculated based on the first reflected light is obtained as a value closer to the object distance. This is because the round-trip time of the reflected light is added as the number of times increases from the second time to the third time, and the flight time deviates from the accurate flight time (round-trip time to the reflective object). Therefore, by adopting the above-mentioned configuration, the object distance can be accurately specified in the situation where multiple reflections occur.
  • This disclosure can also be realized in various forms. For example, it is realized in the form of a distance measuring system, a mobile body equipped with a distance measuring device, a distance measuring method, a computer program for realizing these devices and methods, a non-temporary recording medium on which such a computer program is recorded, and the like. Can be done.
  • the present disclosure is not limited to the above-described embodiment, and can be realized by various configurations within a range not deviating from the purpose.
  • the technical features in each embodiment corresponding to the technical features in the embodiments described in the column of the outline of the invention may be used to solve some or all of the above-mentioned problems, or one of the above-mentioned effects. It is possible to replace or combine as appropriate to achieve the part or all. Further, if the technical feature is not described as essential in the present specification, it can be appropriately deleted.

Abstract

測距装置(10、10a)は、発光部(40)、受光部(60)、光の飛行時間を利用して対象物距離を演算する演算部(20)を備える。演算部は、ヒストグラム生成部(230)と、ヒストグラムにおいてピーク飛行時間を境に増加から減少に変わる山部であって測距装置からの距離の互いの差が所定範囲内である複数の物体からの光の強度の分布が合成されて得られる山部である複合山部が存在するかを推定する複合山部推定部(240)と、山部において、立ち上がり時間と立ち下がり時間とを特定する飛行時間特定部(250)と、複合山部が存在すると推定された場合に、立ち上がり時間と立ち下がり時間とのうちのいずれかに基づき基飛行時間を決定する基飛行時間決定部(260)と、基飛行時間を利用して対象物距離として算出する距離算出部(270)と、を備える。

Description

測距装置 関連出願の相互参照
 本出願は、2020年6月12日に出願された日本出願番号2020-102143号、および2021年4月19日に出願された日本出願番号2021-070064号に基づくもので、ここにその記載内容を援用する。
 本開示は、測距装置に関する。
 レーザ光などのパルス光を発光部から射出して、対象物からの反射光を受光部で検出し、照射から受光までの光の飛行時間(TOF:Time Of Flight)を計測することで、対象物の存否の検出や対象物までの距離を測定する測距装置が知られている。このような測距装置では、パルス光を射出するための窓や筐体内部によるパルス光の反射光や、測距装置から対象物までに至る空間に存在する雨や霧などによるパルス光の反射光は、いわゆるクラッタと呼ばれ、対象物までの距離の測定において障害となり得る。具体的には、対象物による反射光とクラッタとが合成されて1つの受光パルスとして検出されてしまい、これにより飛行時間が正しく計測できず、測定される距離の精度が低下するおそれがある。そこで、下記特許文献1の測距装置では、受光強度に対する閾値として、クラッタの振幅値(強度)よりも大きな値を閾値(High閾値)として設定し、合成受光パルスの受光強度に対してHigh閾値を適用することにより、合成受光パルスから対象物による反射波(所望波)を分離し、これにより、対象物までの距離の測定精度を向上させている。
特開2015-194356号公報
 例えば、クラッタの原因が霧である場合には、霧が濃いほどクラッタの振幅値は大きくなり、クラッタと所望波との振幅の差が小さくなる。また、例えば、クラッタの原因がパルス光を射出するための窓や筐体である場合、かかる窓や筐体表面の反射率が大きいほどクラッタの振幅は大きくなり、クラッタと所望波との振幅の差が小さくなる。このようにクラッタと所望波との振幅の差が小さくなる場合には、合成波の受光強度の時間変化は、あたかも1つの山のような形状となる。このような現象は、クラッタが生じる場合に限らず、例えば、測距装置を基準として互いに同様な方角で距離が異なる複数の対象物が存在する場合にも、これら複数の対象物からの反射光が合わさって起こり得る。このように、複数の反射光の受光強度が合成されて1つの山のような形状となる場合においては、特許文献1の測距装置では、受光波から所望波を分離できず、測距精度を向上させることができない。特許文献1では、例えば、濃霧や激しい雨が降っているような悪環境下では、測距を行わないようにすることで、測距精度の低下を抑制している。しかし、悪環境下のように、複数の反射光の受光強度の時間変化が合成されて1つの山のような形状となる状況においても、精度良く測距可能な技術が望まれる。
 本開示の一形態として、測距装置が提供される。この測距装置は、パルス光を射出する発光部と、前記パルス光の反射光を含む光を受光する受光部と、前記受光部で受光された光の飛行時間を利用して、前記パルス光を反射して前記反射光を出力する対象物までの距離である対象物距離を演算する演算部と、を備える。前記演算部は、複数の前記飛行時間における各々の前記受光部における受光強度を表すヒストグラムを生成するヒストグラム生成部と、前記ヒストグラムにおいて、ピーク飛行時間を境に前記飛行時間に対する前記受光強度の変化が増加から減少に変わる山部であって前記測距装置からの距離の互いの差が所定範囲内である複数の物体からの光の強度の分布が合成されて得られる山部である複合山部が存在するかを推定する複合山部推定部と、前記山部において、背景光を含むノイズ光の受光強度よりも大きな第1閾値受光強度と一致する2つの飛行時間であって、より短い飛行時間である立ち上がり時間と、より長い飛行時間である立ち下がり時間と、を特定する飛行時間特定部と、前記複合山部が存在すると推定された場合に、前記複合山部について特定された前記立ち上がり時間と前記立ち下がり時間とのうちのいずれかに基づき、前記対象物距離を算出するための基となる基飛行時間を決定する基飛行時間決定部と、前記基飛行時間を利用して、前記複数の物体のうちの一の物体の距離を、前記対象物距離として算出する距離算出部と、を備える。
 この形態の測距装置によれば、複合山部について特定された立ち上がり時間と立ち下がり時間とのうちのいずれかに基づき、対象物距離を算出するための基となる基飛行時間を決定し、決定された基飛行時間を利用して対象物距離が演算されるので、複数の反射光の受光強度の時間変化が合成されて1つの山のような形状となる状況においても、対象物距離を精度良く測定できる。
 本開示についての上記目的およびその他の目的、特徴や利点は、添付の図面を参照しながら下記の詳細な記述により、より明確になる。その図面は、
図1は、本開示の一実施形態としての測距装置の概略構成図、 図2は、受光アレイの構成を模式的に示す説明図、 図3は、SPAD回路の構成を模式的に示す回路図、 図4は、第1実施形態の測距装置の機能構成を示すブロック図、 図5は、複合山部を有するヒストグラムの一例を示す説明図、 図6は、第1実施形態における測距処理の手順を示すフローチャート、 図7は、複合山部を有しないヒストグラムの一例を示す説明図、 図8は、第2実施形態の測距装置の機能構成を示すブロック図、 図9は、第2実施形態における測距処理の手順を示すフローチャート、 図10は、第2実施形態における複数の対象物と測距装置との位置関係を模式的に示す平面図、 図11は、第2実施形態における複合山部を有するヒストグラムの一例を示す説明図、 図12は、第2実施形態における複合山部を有するヒストグラムの一例を示す説明図、 図13は、第2実施形態における複合山部を有するヒストグラムの一例を示す説明図、 図14は、第3実施形態における複合山部を有するヒストグラムの一例を示す説明図、 図15は、第4実施形態におけるクラッタ位置設定処理の手順を示すフローチャート、 図16は、第4実施形態における距離のヒストグラムの一例を示す説明図、 図17は、第5実施形態におけるクラッタ位置設定処理の手順を示すフローチャート、 図18は、第5実施形態におけるヒストグラムの例を示す説明図、 図19は、第6実施形態における測距処理の手順を示すフローチャート、 図20は、第6実施形態におけるヒストグラムの一例を示す説明図、 図21は、第7実施形態における測距処理の手順を示すフローチャート。
A.第1実施形態:
A1.装置構成:
 図1に示す測距装置10は、測距のためのパルス光を射出して外部物体からの反射光を受ける光学系30と、光学系30から得られた信号を処理する演算部20とを備える。外部物体を、「反射物」とも呼ぶ。光学系30は、パルス光としてのレーザ光を射出する発光部40と、レーザ光を予め定められた視野範囲80内で走査させる走査部50と、外部物体からの反射光や外乱光を含む入射光を受光するための受光部60とを備える。測距装置10は、前面に窓92を有する筐体90に収容されている。窓92は、発光部40から射出されるパルス光の多くを透過し、一部を反射する。
 測距装置10は、例えば、自動車などの車両に搭載される車載用のLiDAR(Laser Imaging Detection and Ranging)である。車両が水平な路面を走行している場合に、視野範囲80の横方向は水平方向Xと一致し、縦方向は鉛直方向Yと一致する。
 発光部40は、パルス光を含むレーザ光を射出する半導体レーザ素子(以下、単にレーザ素子とも呼ぶ)41と、レーザ素子41の駆動回路を組み込んだ回路基板43と、レーザ素子41から射出されたレーザ光を平行光にするコリメートレンズ45とを備える。レーザ素子41は、いわゆる短パルスレーザを発振可能なレーザダイオードである。本実施形態において、レーザ素子41は、複数のレーザダイオードを鉛直方向に沿って配列させることにより矩形状のレーザ発光領域を構成する。レーザ素子41が出力するレーザ光の強度は、レーザ素子41に供給される電圧に応じて調整可能に構成されている。
 走査部50は、いわゆる一次元スキャナによって構成される。走査部50は、ミラー54と、ロータリソレノイド58と、回転部56とによって構成される。ミラー54は、コリメートレンズ45により平行光とされたレーザ光を反射する。ロータリソレノイド58は、演算部20からの制御信号を受けて、予め定められた角度範囲内で正転および逆転を繰り返す。回転部56は、ロータリソレノイド58によって駆動し、鉛直方向を軸方向とする回転軸で正転および逆転を繰り返し、ミラー54を水平方向に沿った一方向に走査させる。コリメートレンズ45を介してレーザ素子41から射出されたレーザ光は、ミラー54によって反射され、ミラー54の回転により水平方向に沿って走査される。図1に示す視野範囲80は、このレーザ光の全走査範囲に相当する。視野範囲80内の各画素位置で受光強度が得られるので、視野範囲80内の受光強度の分布は一種の画像を構成する。なお、走査部50を省略して、発光部40から視野範囲80内の全体にわたってパルス光を射出するとともに、受光部60で視野範囲80内の全体にわたる反射光を受光するようにしてもよい。本実施形態では、走査範囲内における各位置、換言すると視野範囲80内の各画素位置に対してパルス光が照射される。そして、かかるパルス光の照射と各画素位置からの反射光に基づく後述の測距処理とが、各画素位置について所定の時間間隔で実行される。
 発光部40から出力されるレーザ光は、人や車などの外部物体(反射物)があると、その表面で乱反射し、その一部は反射光として走査部50のミラー54に戻ってくる。この反射光は、ミラー54で反射されて、外乱光とともに入射光として受光部60の受光レンズ61に入射し、受光レンズ61で集光されて受光アレイ65に入射する。なお、測距装置10から出力されるレーザ光は、外部物体に限らず、測距装置10内部の物体、例えば、窓92においても乱反射し、その反射光の一部は、受光アレイ65に入射する。
 図2に示すように、受光アレイ65は、二次元配列された複数の画素66で構成される。1つの画素66は、水平方向にH個、鉛直方向にV個となるように配列された複数のSPAD(Single Photon Avalanche Diode)回路68で構成されている。H及びVはそれぞれ1以上の整数である。本実施形態ではH=V=5であり、水平方向および鉛直方向においてそれぞれ5個のSPAD回路68で構成されている。但し、任意の数のSPAD回路68で画素66を構成することが可能であり、一つのSPAD回路68で画素66を構成してもよい。1つの画素66の受光結果は、視野範囲80内の1つの画素位置における受光強度となる。
 図3に示すように、SPAD回路68は、電源Vccと接地ラインとの間に直列にアバランシェダイオードDaとクエンチ抵抗器Rqとを接続し、その接続点の電圧を論理演算素子の一つである反転素子INVに入力し、電圧レベルが反転したデジタル信号に変換している。反転素子INVの出力信号Sout は、外部にそのまま出力される。本実施形態においてクエンチ抵抗器RqはFETとして構成されており、選択信号SCがアクティプとなっていれば、そのオン抵抗がクエンチ抵抗器Rqとして働く。選択信号SCがノンアクティブとなれば、クエンチ抵抗器Rqはハイインピーダンス状態となるので、光がアバランシェダイオードDaに入射しても、クエンチ電流は流れず、結果的にSPAD回路68は、動作しない。選択信号SCは、画素66内の5×5個のSPAD回路68に対しては、一括して出力され、各画素66からの信号を読み出すか読み出さないかを指定するのに用いられる。本実施形態では、アバランシェダイオードDaをガイガーモードで動作させているが、アバランシェダイオードDaをリニアモードで用い、その出力をアナログ信号のまま扱ってもよい。また、アバランシェダイオードDaに代えて、PINフォトダイオードを用いてもよい。
 SPAD回路68に光が入射していなければ、アバランシェダイオードDaは、非導通状態に保たれる。このため、反転素子INVの入力側は、クエンチ抵抗器Rqを介してプルアップされた状態、つまりハイレベルHに保たれている。従って、反転素子INVの出力はロウレベルLに保たれる。各SPAD回路68に外部から光が入射すると、アバランシェダイオードDaは、入射した光(フォトン)により通電状態となる。この結果、クエンチ抵抗器Rqを介して大きな電流が流れ、反転素子INVの入力側は一旦ロウレベルLとなり、反転素子INVの出力はハイレベルHに反転する。クエンチ抵抗器Rqを介して大きな電流が流れた結果、アバランシェダイオードDaに印加される電圧は低下するから、アバランシェダイオードDaへの電力供給は止り、アバランシェダイオードDaは、非導通状態に戻る。この結果、反転素子INVの出力信号も反転してロウレベルLに戻る。結果的に、反転素子INVは、各SPAD回路68に光(フォトン)が入射すると、ごく短時間、ハイレベルとなるパルス信号を出力することになる。そこで、各SPAD回路68が光を受光するタイミングに合わせて、選択信号SCをハイレベルHにすれば、反転素子INVの出力信号、つまり各SPAD回路68からの出力信号Sout は、アバランシェダイオードDaの状態を反映したデジタル信号となる。そして、この出力信号Soutは、照射光が走査範囲に存在する外部物体や窓92等に反射して戻ってくる反射光や外乱光を含む入射光の受光により生じるパルス信号に相当する。
 演算部20は、受光部60で受光された反射光の飛行時間を利用して、パルス光を反射して反射光を出力する対象物までの距離(以下、「対象物距離」と呼ぶ)を演算する。かかる距離の演算方法の概要は以下の通りである。図4に示すように、発光部40から射出されたパルス光P1は、外部物体である反射物OBJにおいて反射される。換言すると、反射物OBJは、パルス光P1の反射光P2を出力する。また、窓92の内側表面においても、パルス光P1は反射され、反射光P3が出力される。その結果、受光部60には、反射光P2、P3が届く。このとき、パルス光P1の射出から反射光P2、P3の受光までの時間が光の飛行時間Tfとして特定される。演算部200は、この飛行時間Tfを利用して測距装置10(発光部40および受光部60)から反射物OJBまでの距離を演算する。なお、窓92の内側表面においてパルス光P1が反射されて得られた反射光や、その反射光が測距装置10の筐体内部で何度か反射されて受光部60に受光される光は、「クラッタ」とも呼ばれる。
 図4に示すように、演算部20は、加算部220と、ヒストグラム生成部230と、複合山部推定部240と、飛行時間特定部250と、基飛行時間決定部260と、距離算出部270と、制御部280と、メモリ290とを備える。
 加算部220は、受光アレイ65を構成する画素66に含まれる各SPAD回路68の出力を加算する。入射する光パルスが一つの画素66に入射すると、画素66に含まれるSPAD回路68が動作する。SPAD回路68は、一つのフォトンが入射しただけでこれを検出することが可能である。しかし、SPAD回路68において、反射物OBJから出力される限られた光の検出は確率的なものにならざるを得ない。そこで、加算部220は、確率的にしか入射しない光を検出し得ないSPAD回路68からの出力信号Soutを、各画素66に含まれるすべてのSPAD回路68分だけ加算することにより、各画素66における反射物OBJからの反射光をより確実に検出するように構成されている。
 ヒストグラム生成部230は、加算部220の加算結果を時系列的に取得することにより、受光強度のヒストグラムを生成し、メモリ290に記憶させる。ヒストグラム生成部230が生成するヒストグラムは、複数の飛行時間の各々における受光強度を表すグラフといえる。受光強度は、1画素66内における受光したSPAD回路68の合計数である。
 複合山部推定部240は、ヒストグラム生成部230により生成されたヒストグラムにおいて、複合山部が存在するかを推定する。複合山部とは、ヒストグラムにおいて、ピーク飛行時間を境に飛行時間に対する受光強度の変化が増加から減少に変わる山部であって、測距装置10からの距離の互いの差が所定範囲内である複数の物体からの光の強度の分布が合成されて得られる山部を意味する。
 図5において、横軸は飛行時間Tfを示し、縦軸は受光強度Iを示す。図5において太い実線で示すヒストグラムhr1は、細い一点鎖線で示すクラッタによる受光強度のヒストグラムhr11と、細い実線で測距装置10の外部の対象物による反射光の受光強度のヒストグラムhr12とが合成されて得られるヒストグラムである。図5に示す各ヒストグラムは、測距装置10の外部の対象物が測距装置10(窓92)の近傍に位置している場合のヒストグラムを表している。図5の例では、ヒストグラムhr1には、複合山部mp0が現れている。ここで、ヒストグラムhr11とヒストグラムhr12とは、互いに時間的に重複している。そして、ヒストグラムhr11のピークの受光強度(度数)と、ヒストグラムhr12のピークの受光強度(度数)とは、互いにほぼ等しい。このため、図5に示すように、複合山部mp0は、1つの山のような形状として現れている。複合山部推定部240は、山部における立ち上がり時間から立ち下がり時間までの間の時間、換言すると、ヒストグラムにおいて、受光強度が後述する第1閾値以上となる時間が所定時間以上連続する場合に、ヒストグラムにおいて複合山部が存在すると推定する。なお、図5に示す受光強度Iの各値、すなわち、ノイズ強度I1、第1閾値受光強度I3およびピーク強度I4の詳細については後述する。また、図5に示す飛行時間Tfの各値、すなわち、立ち上がり時間Tu1、立ち下がり時間Td1、ピーク時間Tp1、および中間時間Tc1の詳細については後述する。
 飛行時間特定部250は、ヒストグラム生成部230により生成されたヒストグラムにおいて、受光強度が、第1閾値受光強度と一致する2つの飛行時間(後述する立ち上がり時間および立ち下がり時間)を特定する。第1閾値受光強度とは、対象距離を算出する際に用いられる受光強度の閾値であり、ヒストグラムの山部における立ち上がり時間と立ち下がり時間とを決定するために用いられる。図5の例では、第1閾値受光強度I3が設定されている。本実施形態において、第1閾値受光強度は、下記式(1)に示す値として設定される。
 第1閾値受光強度=(ピーク強度-ノイズ強度)×0.4+ノイズ強度・・・(1)
 ここで、ピーク強度とは、山部における最大受光強度を意味し、図5の複合山部mp0では、ピーク強度I4が該当する。また、ノイズ強度とは、パルス光の反射光以外の光(以下、「ノイズ光」と呼ぶ)の受光強度を意味する。このようなノイズ光としては、太陽光や、太陽光が外部の対象物において反射して得られる反射光や、街灯の光などが該当する。ノイズ強度は、発光部40がパルス光を射出していないタイミングで測定される受光強度の所定時間における平均値として求めることができる。式(1)における(ピーク強度-ノイズ強度)は、反射光強度とも呼ばれる。したがって、式(1)は、反射光強度の4割の値を、ノイズ強度に加えて得られた受光強度が、第1閾値受光強度として設定されることを意味する。なお、4割に代えて、4割よりも小さい又は4割よりも大きな任意の割合を用いてもよい。図5の例では、飛行時間特定部250は、ヒストグラムhr1と、第1閾値受光強度I3とが一致する2つの時間Tu1、Td1のうち、時間的により前の時間Tu1を、立ち上がり時間Tu1として特定し、また、時間的により後の時間Td1を、立ち下がり時間Td1として特定する。
 基飛行時間決定部260は、対象距離を算出するための基となる飛行時間(以下、「基飛行時間」と呼ぶ)を決定する。本実施形態では、基飛行時間決定部260は、後述するように、立ち下がり時間を基飛行時間として決定する。立ち下がり時間を基飛行時間として決定する理由については、後述する。
 距離算出部270は、基飛行時間決定部260により決定された基飛行時間を利用して、対象物距離を算出する。かかる対象物距離の算出方法の詳細については、後述する。算出された対象物距離は、例えば、測距装置10が搭載された車両において、対象物との衝突が発生するか否かの推定処理や、衝突回避のための操舵やブレーキの制御において用いられ得る。
A2.測距処理:
 図6に示す測距処理は、各画素位置に対して所定の時間間隔で定期的に実行される。発光部40がパルス光を射出して受光部60においてその反射光が受光され、ヒストグラム生成部230によりヒストグラムが生成された後のタイミングで、測距処理は開始される。
 複合山部推定部240は、生成されたヒストグラムにおいて複合山部を特定し、飛行時間特定部250は、立ち上がり時間および立ち下がり時間を特定する(ステップS105)。具体的には、複合山部推定部240は、ヒストグラムにおいて、第1閾値受光強度を超える飛行時間が所定時間以上連続する部位を特定することにより複合山部を特定する。また、飛行時間特定部250は、複合山部において第1閾値と一致する2つの飛行時間を、立ち上がり時間および立ち下がり時間として特定する。
 例えば、図5の例では、複合山部mp0が特定され、複合山部mp0と第1閾値受光強度I3とが一致する2つの時間Tu1、Td1のうち、時間的により早い時間Tu1が立ち上がり時間Tu1として特定され、時間的により遅い時間Td1が立ち下がり時間Td1として特定される。
 図6に示すように、基飛行時間決定部260は、反射光の受光位置はクラッタの位置と一致するか否かを判定する(ステップS110)。クラッタ、すなわち、窓92や筐体内壁面によるパルス光の反射光は、常に同じ位置(方角)から受光部60に入光する。また、窓92や筐体内壁面までに距離は常に一定である。このため、窓92や筐体内壁面によるパルス光の反射光はどの方角の画素位置において、どの程度の飛行時間で特定されるかは、予め実験やシミュレーションにより特定できる。そこで、本実施形態では、窓92や筐体内壁面によるパルス光の反射光の画素位置および飛行時間が予めメモリ290に記憶されており、複合山部が検出された際の画素位置および複合山部の飛行時間が、メモリ290に記憶されている画素位置および飛行時間に一致するか否かを判定することにより、反射光の受光位置はクラッタの位置と一致するか否かを判定するようにしている。
 反射光の受光位置はクラッタの位置と一致しないと判定された場合(ステップS110:NO)、距離算出部270は、通常の方法により距離を算出し、対象物距離として特定する(ステップS115)。かかるステップS115における測距方法を、図7を用いて説明する。
 図7における縦軸および横軸は、図5と同じであるので、その詳細な説明を省略する。図7に示す山部mp1は、例えば、クラッタが生じる方角とは異なる方角に測距装置10の外部において対象物が存在し、かかる対象物からパルス光の反射光を受光部60で受光したような場合に生じ得る。この場合、山部mp1の画素位置(方角)がクラッタの位置と異なる(含まない)ので、飛行時間特定部250により、立ち上がり時間Tu2と立下り時間Td2の中間時間が山部mp1のピーク時間Tp2として特定される。そして、基飛行時間決定部260は、このピーク時間Tp2を、基飛行時間として決定し、距離算出部270は、ピーク時間Tp2を利用して測距装置10からこれら複数の物体までの距離を求める。
 図6に示すように、反射光の受光位置はクラッタの位置と一致すると判定された場合(ステップS110:YES)、基飛行時間決定部260は、ステップS105で特定された立ち下り時間を、基飛行時間として決定する(ステップS120)。距離算出部270は、ステップS120で決定された基飛行時間に基づき、距離を算出する(ステップS125)。具体的には、例えば、図5に示す複合山部mp0の立ち下がり時間Td1が基飛行時間として決定されると、距離算出部270は、立ち下がり時間Td1から予め定められた時間Δtだけ遡った時間tp1をピーク時間tp1として特定する。そして、距離算出部270は、このピーク時間tp1を利用して対象距離を算出する。より具体的には、ピーク時間tp1が、パルス光およびその反射光の飛行時間Tfであるものとして、対象距離を算出する。
 図5に示すように、時間Td1から時間Δtだけ遡った時間(ピーク時間Tp1)は、ヒストグラムhr12のピーク時間とほぼ一致している。他方、立ち上がり時間Tu1と立ち下がり時間Td1の中間時間Tc1は、ヒストグラムhr1のピーク時間から大きくずれている。したがって、上述のようにピーク時間Tp1を基飛行時間として距離を算出することにより、通常の方法で距離を算出するよりも、ヒストグラムhr1に対応する対象物までの距離をより正確に算出できることが分かる。上述のように、中間時間Tc1がヒストグラムhr1のピーク時間から大きくずれる理由について説明する。複合山部mp0は、より遠くにある反射物OBJから出力される反射光によるヒストグラムhr12と、より近くにある窓92から出力されるクラッタによるヒストグラムh12とが合成されて生じたものであるため、その立ち上がり位置Tu1は、ピーク位置がより時間的に前となるヒストグラムhr11の立ち上がり位置とほぼ一致することとなる。また、ヒストグラムhr11の立ち下がり位置は、より遠くに存在する(より飛行時間が長い)対象物からの反射光によるヒストグラムh12の立ち下がり位置とほぼ一致することとなる。他方、時間的に前となるヒストグラムhr11の立ち下がり位置は、時間的により後となるヒストグラムhr12の存在により、ヒストグラムhr1の立ち下がり位置から大きくずれている。同様に、時間的に後となるヒストグラムhr12の立ち上がり位置は、時間的に前となるヒストグラムhr11の存在により、ヒストグラムhr1の立ち上がり位置から大きくずれている。したがって、通常の方法にしたがって対象距離を算出するために、立ち上がり時間Tu1と立ち下がり時間Td1の中間時間である中間時間Tc1を特定しても、かかる中間時間Tc1は、ヒストグラムhr11のピーク時間とヒストグラムhr12のピーク時間とのいずれからも大きく異なるずれた時間となる。これに対して、対象物の立ち下がり時間とほぼ一致するヒストグラムhr1の立ち下がり時間を基準として、かかる立ち下がり時間から時間Δtだけ遡った時間を基飛行時間とすることにより、対象物のピーク時間を基飛行時間とする場合に算出される距離とほぼ同じ距離を求めることができる。
 図6に示すように、距離算出部270は、ステップS125により算出された距離が予め定められた閾値距離よりも小さいか否かを判定する(ステップS130)。ステップS130における閾値距離は、かかる距離よりも大きければ、クラッタではないと認められる距離として、発光部40および受光部60から窓92までの距離に所定値を加えた値として設定されている。このステップS130は、基飛行時間に基づき算出された距離がクラッタの場合に算出される距離よりも大きいことで、受光した反射光が確かにクラッタによるものでないことを確認するために実行される。
 ステップS125により算出された距離が閾値距離よりも小さくないと判定された場合(ステップS130:NO)、上述のステップS115が実行され、通常の方法により対象物距離が算出される。これに対して、ステップS120により算出された距離が閾値距離よりも小さいと判定された場合(ステップS130:YES)、距離算出部270は、ステップS125により算出された距離を、対象物距離として特定する(ステップS135)。上述のステップS115またはステップS135の完了後、該当画素位置での測距処理は終了し、次の画素位置についての測距処理が開始されることとなる。
 以上説明した第1実施形態の測距装置10によれば、窓92や筐体内壁面等の筐体内部におけるパルス光の反射光であるクラッタを含む複合山部が存在すると推定された場合に、立ち下がり時間に基づき基飛行時間を決定するので、クラッタが生じる場合において対象物距離を精度良く測定できる。一般に、窓92や筐体内部におけるパルス光の反射光は、測距装置の外部に存在する対象物よりも飛行時間が短い。換言すると、外部に存在する対象物によるパルス光の反射光は、筐体内部におけるパルス光の反射光よりも飛行時間が長い。このため、複合山部における立ち下がり位置(時間)と、外部に存在する対象物によるパルス光の反射光における立ち下がり位置(時間)との誤差は、小さい。したがって、本実施形態の測距装置10によれば、複合山部が存在すると推定される場合における測距精度を向上できる。
 また、クラッタが存在しない場合には、立ち上がり時間と立ち下がり時間との間の中間時間を、基飛行時間として決定されるので、このような場合においても立ち上がり時間と立ち下がり時間とのうちのいずれかに基づき基飛行時間が決定される構成に比べて、対象物距離をより精度良く測定できる。その一因として、立ち上がり時間と立ち下がり時間とをそれぞれ特定する際に誤差が生じたとしても、その中間時間とすることにより、立ち上がり時間または立ち下がり時間のいずれか一方のみを基準として基飛行時間を決定する構成に比べて誤差の影響を抑えることができることが挙げられる。
B.第2実施形態:
 図8に示す第2実施形態の測距装置10aの構成は、演算部20が中間時間特定部235を追加して備える点において、第1実施形態の測距装置10の構成と異なる。測距装置10aのその他の構成は、第1実施形態の測距装置10と同じであるので、同一の構成には同一の符号を付し、その詳細な説明を省略する。
 中間時間特定部235は、複合山部において、第1閾値受光強度と一致する2つの飛行時間の中間時間(以下、「第1中間時間」と呼ぶ)と、第2閾値受光強度と一致する2つの飛行時間の中間時間(以下、「第2中間時間」と呼ぶ)とを特定する。第1閾値受光強度は、第1実施形態の第1閾値受光強度と同じである。第2閾値受光強度とは、第1閾値受光強度よりも小さな受光強度であって、ノイズ強度もより大きな強度である。具体的には、本実施形態では、第2閾値受光強度は、下記式(2)に示す値として設定される。
 第2閾値受光強度=(ピーク強度-ノイズ強度)×0.1+ノイズ強度・・・(2)
 図9に示す第2実施形態の測距処理は、ステップS110が省略されている点と、ステップS106、S107、S108、S123を追加して実行する点とにおいて、図6に示す第1実施形態の測距処理と異なる。第2実施形態の測距処理におけるその他の手順は、第1実施形態の測距処理と同じであるので、同一の手順には同一の符号を付し、その詳細な説明を両略する。
 ステップS105の完了後、中間時間特定部235は、第1中間時間と第2中間時間とを算出する(ステップS106)。
 図10に示すように、或る画素位置81に、2つの対象物OBJ1、OBJ2が存在し、且つ、各対象物OBJ1、OBJ2が測距装置10aから見て露出した領域Ar1、Ar2を有する場合、これら2つの領域Ar1、Ar2からそれぞれ出力される反射光が測距装置10aの受光部60に入光することとなる。図10の例では、対象物OBJ1が測距装置10aから見て手前側に位置し、対象物OBJ2が測距装置10aから見て奥側に位置している。また、画素位置81において、測距装置10aから見た領域Ar1の露出面積は、領域Ar2の露出面積よりも大きい。このような状況において生成されるヒストグラムについて図11を用いて説明する。
 図11における縦軸および横軸は、図5と同じであるので、その詳細な説明を省略する。図10の状況において生成されるヒストグラムhr3は、大きな山部と、それよりも時間的に後にピークのある小さな山部とが合成されて生じた複合山部mp3を有する。具体的には、複合山部mp3は、細い一点鎖線により示される領域Ar1からの反射光によるヒストグラムhr31と、細い実線により示される領域Ar2からの反射光によるヒストグラムhr32とが合成されて得られるヒストグラムhr3における山部に相当する。上述のように、領域Ar1はより手前に位置して露出面積がより大きい。これに対して領域Ar2はより奥側に位置して露出面積が小さい。したがって、ヒストグラムhr31の山部は、ヒストグラムhr32の山部に比べて時間的に前に位置し、且つ、ピーク受光強度が大きい。このため、複合山部mp3は、時間的に前の位置に比較的大きなピークを有し、時間的に後の位置に比較的小さなピークを有するいびつな山型の形状を有する。
 このような形状の複合山部mp3に対し、上述のステップS106において中間時間特定部235は、複合山部mp3において、第1閾値受光強度I11と一致する2つの時間(立ち上がり時間Tu3および立ち下がり時間Td3)を特定し、これら2つの時間の中間時間Tc31を特定する。また、中間時間特定部235は、複合山部mp3において、第2閾値受光強度I12と一致する2つの時間(立ち上がり時間Tu4および立ち下がり時間Td4)とを特定し、これら3つの時間の中間時間Tc32を特定する。なお、本実施形態では、第1閾値受光強度と一致する2つの時間の中間時間を第1中間時間と呼び、第2閾値受光強度と一致する2つの時間の中間時間を第2中間時間と呼ぶ。したがって、中間時間Tc31は第1中間時間Tc31とも呼ばれ、また、中間時間Tc32は第2中間時間Tc32とも呼ばれる。
 図9に示すように、複合山部推定部240は、第2閾値受光強度と一致する2つの時間の間の幅(以下、「根元幅」と呼ぶ)が閾値幅よりも大きいか否かを判定する(ステップS107)。図11の例では、複合山部mp3において、第2閾値受光強度I12と一致する2つの時間Tu4、Td4の間の幅が根元幅dt31に相当する。複合山部mp3のように、複数の対象物からの反射光のヒストグラムに含まれる山部が合成して得られる複合山部では、根本幅は大きい。これに対して、単一の対象物からの反射光による山部の根本幅は小さい。そこで、複合山部であることを特定するための閾値幅が、予め実験等により特定されて設定されている。図9に示すように根元幅が閾値幅よりも大きくないと判定された場合(ステップS107:NO)、上述のステップS115が実行される。この場合、単一の対象物からの反射光が入射されているので、通常の方法により距離が算出されて対象物距離として特定される。
 根元幅が閾値幅よりも大きいと判定された場合(ステップS107:YES)、この場合、山部が複合山部である可能性が高い。そして、この場合、基飛行時間決定部260は、第1中間時間が第2中間時間以下であるか否かを判定する(ステップS108)。
 図11の例では、第1中間時間Tc31は、第2中間時間Tc32よりも小さい(早い)。また、図13の例では、第1中間時間と第2中間時間とが一致している。図13に示すヒストグラムhr5は、ヒストグラムhr51と、ヒストグラムhr52とを合成して得られたヒストグラムである。ヒストグラムhr51は、図10に示す対象物OBJ1からの反射光によるものであり、ヒストグラムhr52は、図10に示す対象物OBJ2からの反射光によるものである。但し、図13の例では、図10の例とは異なり、画素位置81において、測距装置10aから見た領域Ar1の露出面積は、領域Ar2の露出面積と等しい。このような状況においては、対象物OBJ1からの反射光によるヒストグラムhr51のピーク受光強度は、対象物OBJ2からの反射光によるヒストグラムhr52のピーク受光強度と等しくなる。このため、複合山部mp5は、ピークが一つの大きな山型の形状を有する。そして、この例では、第1中間時間と第2中間時間とはいずれも中間時間Tc51となる。
 図11および図13の例では、ステップS108において、第1中間時間Tc31は、第2中間時間Tc32よりも小さい(早い)と判定されることとなる。このように、第1中間時間Tc31は、第2中間時間Tc32よりも小さい(早い)と判定された場合(ステップS108:YES)、基飛行時間決定部260は、ステップS105で特定された立ち上がり時間を基飛行時間として決定する(ステップS123)。例えば、図11の例では、立ち上がり時間Tu3が基飛行時間として決定されることとなる。また、図13の例では、立ち上がり時間Tu7が基飛行時間として決定されることとなる。
 図11および図13の例とは異なり、第1中間時間が第2中間時間よりも大きい(遅い)場合について、図12を用いて説明する。図12に示すヒストグラムhr4は、ヒストグラムhr41と、ヒストグラムhr42とを合成して得られたヒストグラムである。ヒストグラムhr41は、図10に示す対象物OBJ1からの反射光によるものであり、ヒストグラムhr42は、図10に示す対象物OBJ2からの反射光によるものである。但し、図12の例では、図10の例とは異なり、画素位置81において、測距装置10aから見た領域Ar1の露出面積は、領域Ar2の露出面積よりも小さい。このような状況においては、対象物OBJ1からの反射光によるヒストグラムhr41のピーク受光強度は、対象物OBJ2からの反射光によるヒストグラムhr42のピーク受光強度に比べて小さくなる。このため、複合山部mp4は、時間的に前の位置に比較的小さなピークを有し、時間的に後の位置に比較的大きなピークを有するいびつな山型の形状を有する。
 このような図12の例では、複合山部mp4において、第1閾値受光強度I21と一致する2つの時間(立ち上がり時間Tu5および立ち下がり時間Td5)が特定され、また、これら2つの時間の中間時間である第2中間時間Tc41が特定される。さらに、複合山部mp4において、第2閾値受光強度I22と一致する2つの時間(立ち上がり時間Tu6および立ち下がり時間Td6)が特定され、また、これら2つの時間の中間時間である第2中間時間Tc42が特定される。そして、この例では、第1中間時間Tc41は、第2中間時間Tc42よりも大きい(遅い)。
 図9に示すように、第1中間時間が第2中間時間以下でないと判定された場合(ステップS108:NO)、上述のステップS120が実行され、ステップS105で特定された立ち下がり時間が基飛行時間として決定される。図12の例では、立ち下がり時間Td5が基飛行時間として決定される。
 図9に示すように、ステップS120またはステップS123の完了後、上述のステップS125、およびS135が実行される。
 図11に示すように、第1中間時間が第2中間時間以下(早いまたは同じ)場合、図10に示すように、測距装置10aから見たときに、より手前側の対象物OBJ1の露出面積が大きい。そして、この場合、対象物OBJ1の表面におけるパルス光の反射光の飛行時間を基時間として対象物距離が算出されるので、測距装置10aが搭載された車両にとってより影響の大きな対象物に基づき対象物距離が算出される。また、図12に示すように、第1中間時間が第2中間時間よりも大きい(遅い)場合、図10の例とは異なり、測距装置10aから見たときに、より奥側の対象物OBJ2の露出面積が大きい。そして、この場合、対象物OBJ2の表面におけるパルス光の反射光の飛行時間を基飛行時間として対象物距離が算出されるので、この場合もまた、測距装置10aが搭載された車両にとってより影響の大きな対象物に基づき対象物距離が算出される。
 以上説明した第2実施形態の測距装置10aは、第1実施形態の測距装置10と同様な効果を有する。加えて、第1中間時間が第2中間時間よりも小さい場合には、立ち上がり時間に基づき基飛行時間を決定するので、測距装置10aを基準として互いに同様な方角で距離が異なる複数の対象物が存在する場合であって、その方角におけるパルス光の照射面積が前方側の対象物の方が大きい場合において、前方側の対象物までの距離を対象物距離として求めることができる。また、第1中間時間が第2中間時間よりも大きい場合には、立ち下がり時間に基づき基飛行時間を決定するので、測距装置10aを基準として互いに同様な方角で距離が異なる複数の対象物が存在する場合であって、その方角におけるパルス光の照射面積が後方側の対象物の方が大きい場合において、後方側の対象物までの距離を対象物距離として求めることができる。このように、第2実施形態の測距装置10aによれば、測距装置10aを基準として互いに同様な方角で距離が異なる複数の対象物が存在する場合において、前方または後方のいずれかの対象物までの距離、すなわち、実際に対象物が存在する位置までの距離を、対象物距離として求めることができる。加えて、パルス光の照射面積がより大きな対象物までの距離を対象物距離として求めるので、測距装置10aから見て露出面積がより大きな対象物までの距離を対象物距離として求めることができる。このため、例えば、対象物距離を利用して車両の安全機能を発揮させるような構成において、かかる機能をより安全を確保するように発揮させることができる。また、第1中間時間と第2中間時間とが同じ場合には、立ち上がり時間に基づき基飛行時間を決定するので、測距装置10aを基準として互いに同様な方角で距離が異なる複数の対象物が存在し、かつ、パルス光の照射面積が互いに等しいような場合において、前後のいずれかの対象物までの距離を対象物距離として求めることができる。このため、対象物が存在しない位置までの距離が対象物距離として求められることを抑制できる。
 なお、第2実施形態および上述の第1実施形態からも理解できるように、基飛行時間決定部260は、複合山部が存在する場合に、複合山部について第1閾値受光強度に基づき特定された立ち上がり時間と立ち下がり時間とのうちのいずれかに基づき基飛行時間を決定することができる。
C.第3実施形態:
 第3実施形態の測距装置10aの構成は、図8に示す第2実施形態の測距装置10aと同じであるので、同一の構成要素には同一の符号を付し、その詳細な説明を省略する。第2実施形態の測距装置10aでは、立ち上がり時間Tu3および時間Td3と、中間時間(第1中間時間)Tc31とを特定するのために用いられる閾値受光強度は、第1閾値受光強度I11で共通していた。これに対して、第3実施形態では、立ち上がり時間Tu3および時間Td3を特定するために用いられる閾値受光強度と、中間時間(第1中間時間)を特定するために用いられる閾値受光強度とは、互いに相違する。以下、図14を用いて具体的に説明する。
 図14に示すように、第3実施形態では、第2実施形態と同様に、立ち上がり時間Tu3および立ち下がり時間Td3は、第1閾値受光強度I11を用いて特定される。これに対して、第3実施形態では、中間時間Tc33を特定するための閾値受光強度I14が、上述の第1閾値受光強度I11とは別に設定されている。そして、複合山部mp3においてかかる閾値受光強度I14と一致する2つの時間が特定され、これら2つの時間の中間時間が、中間時間Tc33として特定される。この中間時間Tc33は、第2実施形態の中間時間Tc31の代わりに第1中間時間として用いられる。したがって、例えば、図9に示す測距処理のステップS108では、第1中間時間である中間時間Tc33が、第2中間時間である第2中間時間Tc32以下であるか否かが判定される。第3実施形態において、閾値受光強度I14は、第1閾値受光強度I11よりも大きく、且つ、複合山部mp3のピーク受光強度である強度I13以下の値として、予め設定されている。なお、閾値受光強度I14を、第1閾値受光強度I11よりも小さく、かつ、第2閾値受光強度I12よりも大きな値として設定してもよい。
 以上説明した第3実施形態の測距装置10aは、第2実施形態の測距装置10aと同様な効果を有する。なお、第3実施形態における閾値受光強度I14は、本開示における第3閾値受光強度に相当する。
D.第4実施形態:
D1.装置構成:
 第4実施形態の測距装置10の構成は、第1実施形態の測距装置10と同じであるので、同一の構成には同一の符号を付し、その詳細な説明を省略する。第1実施形態において述べたとおり、窓92や筐体内壁面によるパルス光の反射光に起因して発生するクラッタについては、どの方角の画素位置において、どの程度の飛行時間で特定されるかは、予め実験やシミュレーションにより特定できる。他方、窓92の外側表面に砂や泥などの異物が付着した場合に、かかる位置におけるパルス光の反射光が非常に大きくなり、クラッタとして新たに生じる場合がある。かかる場合には、異物の付着位置に応じてクラッタの生じる画素位置が異なるため、予め特定しておくことができない。第4実施形態の測距装置10では、後述のクラッタ位置設定処理を実行することにより、上述のように異物付着等に起因して生じるクラッタの位置(どの方角の画素位置か)を特定し設定する。
D2.クラッタ位置設定処理:
 図15に示すクラッタ位置設定処理は、測距装置10の電源がオンすると実行される。ヒストグラム生成部230は、図6に示す測距処理によって対象物距離の特定が完了したか否かを判定する(ステップS205)。対象物距離の特定が完了した時点で、各画素についての測距装置10から対象物までの距離が特定されている。
 ヒストグラム生成部230は、今回の測距処理によって特定された距離を含めて、各画素についての距離のヒストグラムを生成する(ステップS210)。本実施形態では、メモリ290には、各画素について、所定回数分の特定された対象物距離が記憶される。上記所定回数は、メモリ290の大きさによるが、例えば、最大100回分を記憶するようにしてもよい。そして、100回分の対象物距離が記憶された状態で新たに対象物距離が特定された場合には、最も古く記憶された対象物距離に上書きして新たに特定された対象物距離が記憶されてもよい。ステップS210において生成される各画素についての距離のヒストグラムは、本開示における「各画素について算出された対象物距離の統計値」に相当する。
 図16では、ステップS210において生成されるヒストグラムの例として、2つのヒストグラムhd1、hd2が表されている。太い実線で示すヒストグラムhd1は、予めクラッタが生じると把握されている位置とは異なる画素についてのヒストグラムであって、窓92に泥等の異物が付着した場合のヒストラムを示している。細い破線のヒストグラムhd2は、予めクラッタが生じると把握されている位置とは異なる画素についてのヒストグラムであって、窓92に泥等の異物が付着していない場合のヒストグラムを示している。
 窓92に異物が付着した場合には、かかる窓92の位置に相当する距離において頻度のピークが生じる。他方、他の距離については、異物が付着しているため、仮にその方角から反射光が向かってきても異物によって遮られてかかる反射光に基づく距離は特定されず、かかる距離における頻度は低いままとなる。したがって、ヒストグラムhd1に示すように、比較的小さな距離d1においてピークが生じることとなる。
 これに対して、窓92に異物が付着していない場合には、車両が動いていれば、当該画素に相当する方角にさまざまな距離で対象物が存在し得る。このため、ヒストグラムhd2に示すように、広い範囲の距離において頻度が0よりも大きくなり、例えば、なだらかな丘状のヒストグラムとなる。
 図15に示すように、ヒストグラム生成部230は、閾値距離以内の範囲に閾値頻度を超える頻度の画素があるか否かを判定する(ステップS215)。図16のヒストグラムhd1が得られた場合には、閾値距離THd以内の範囲に閾値頻度THnを超える頻度が記録されているので、ヒストグラムhd1に対応する画素が、閾値距離以内の範囲に閾値頻度を超える頻度の画素に相当する。他方、図16のヒストグラムhd2が得られた場合には、閾値距離THd以内の範囲に閾値頻度THnを超える頻度が記録されていないので、ヒストグラムhd2に対応する画素は、閾値距離以内の範囲に閾値頻度を超える頻度の画素に相当しない。
 図15に示すように、閾値距離以内の範囲に閾値頻度を超える頻度の画素があると判定された場合(ステップS215:YES)、ヒストグラム生成部230は、当該画素、すなわち、閾値距離以内の範囲に閾値頻度を超える頻度の画素の位置を、クラッタ位置として追加設定する(ステップS220)。上述のように、クラッタ位置はメモリ290に記憶させているので、ヒストグラム生成部230は、かかるクラッタ位置の情報を読みだして、ステップS215で特定された画素位置を追加してメモリ290に記憶させる。ステップS220の完了後、処理はステップS205に戻る。
 上記ステップS215において、閾値距離以内の範囲に閾値頻度を超える頻度の画素が無いと判定された場合(ステップS215:NO)、ヒストグラム生成部230は、すべての画素の位置は、クラッタ位置ではないと追加設定する(ステップS225)。ステップS225の完了後、処理はステップS205に戻る。本実施形態では、窓92を拭く動作を行うワイパーを備え、ステップS215において、閾値距離以内の範囲に閾値頻度を超える頻度の画素があると判定された場合に、かかるワイパーにより窓92を拭く構成としてもよい。窓92に付着した泥等の異物に起因して閾値距離以内の範囲に閾値頻度を超える頻度の画素が存在することとなった場合に、ワイパーの拭き動作により異物が除去された後には、ステップS215において、閾値距離以内の範囲に閾値頻度を超える頻度の画素が無いと判定され、ステップS225の実行により、異物に由来するクラッタ位置はクラッタ位置ではないと更新して設定されることとなる。
 なお、ヒストグラム生成部230に代えて、距離算出部270や他の機能部が上述のステップS205~S220を実行するようにしてもよい。また、図4に示す機能部とは別の新たな機能部が上述のステップS205~S220を実行するようにしてもよい。このように、上述のステップS205~S220を実行する機能部は、本開示の「クラッタ位置設定部」に相当する。
 以上説明した第4実施形態の測距装置10は、第1実施形態の測距装置10と同様な効果を有する。加えて、第4実施形態の測距装置10は、各画素の距離のヒストグラムを生成し、閾値距離THd以内の範囲に閾値頻度THnを超える頻度の画素があると判定された場合に、かかる画素の位置をクラッタ位置として設定するので、泥等の異物が窓92に付着したことに起因して生じたクラッタの位置を特定することができる。このため、図6に示す測距処理のステップS110において、反射光の位置とクラッタ位置とが一致するか否かの判定を精度よく実行できる。
E.第5実施形態:
 第5実施形態の測距装置10の構成は、第4実施形態の測距装置10と同じであるので、同一の構成には同一の符号を付し、その詳細な説明を省略する。第5実施形態の測距装置10は、クラッタ位置設定処理の手順において、第4実施形態の測距装置10と異なる。
 図17に示す第5実施形態のクラッタ位置設定処理は、ステップS205に代えてステップS205aを実行する点と、ステップS210に代えてステップS210aを実行する点と、ステップS215に代えてステップS215aを実行する点とにおいて、図15に示す第4実施形態のクラッタ位置設定処理と異なる。第5実施形態のクラッタ位置設定処理におけるその他の手順は、第4実施形態のクラッタ位置設定処理と同じであるので、同一の手順には同一の符号を付し、その詳細な説明を省略する。
 図17に示すように、ヒストグラム生成部230は、ヒストグラムが生成済みであるか否かを判定する(ステップS205a)。ヒストグラムが生成済みでないと判定された場合(ステップS205a:NO)、処理はステップS205aに戻る。すなわち、ヒストグラム生成部230は、ヒストグラムが生成済みとなるまで待機する。
 ヒストグラムが生成済みであると判定された場合(ステップS205a:YES)、ヒストグラム生成部230は、今回算出されたヒストグラムを含めて、各画素についての飛行時間と平均強度のヒストグラムを生成する(ステップS210a)。本実施形態では、ヒストグラム生成部230により生成された各画素についてのヒストグラムは、メモリ290に記憶されている。そして、ステップS210aにおいて、ヒストグラム生成部230は、メモリ290に記憶されている複数のヒストグラムについて、飛行時間Tgごとの強度の平均値を求め、それにより平均強度のヒストグラムを生成する。各画素について生成される平均強度のヒストグラムは、本開示における「各画素の受光強度の統計値」に相当する。
 図18では、ステップS210aにおいて生成されるヒストグラムの例として、2つのヒストグラムhr61、hr62が表されている。太い実線で示すヒストグラムhr61は、予めクラッタが生じると把握されている位置とは異なる画素についてのヒストグラムであって、窓92に泥等の異物が付着した場合のヒストラムを示している。細い破線のヒストグラムhr62は、予めクラッタが生じると把握されている位置とは異なる画素についてのヒストグラムであって、窓92に泥等の異物が付着していない場合のヒストグラムを示している。
 窓92に異物が付着した場合には、かかる窓92の位置に相当する飛行時間Tf(より正確には、かかる位置までの往復の飛行時間)において受光強度Iのピークが生じる。他方、他の飛行時間Tfについては、異物が付着しているため、仮にその方角から反射光が向かってきても異物によって遮られるため、かかる反射光に基づく飛行時間Tfは特定されず、かかる飛行時間Tfにおける受光強度Iは低いままとなる。したがって、ヒストグラムhr61に示すように、比較的短い時間T6においてピークが生じることとなる。
 これに対して、窓92に異物が付着していない場合には、車両が動いていれば、当該画素に相当する方角にさまざまな距離で対象物が存在し得る。このため、ヒストグラムhr62に示すように、広い範囲の飛行時間Tfにおいて受光強度Iが0よりも大きくなり、例えば、なだらかな丘状のヒストグラムとなる。
 図17に示すように、ヒストグラム生成部230は、閾値飛行時間以内の範囲に平均強度が閾値強度を超えるピークを有する画素があるか否かを判定する(ステップS215a)。図18のヒストグラムhr61が得られた場合には、閾値飛行時間THt以内の範囲に閾値強度It1を超える飛行時間T6が記録されているので、ヒストグラムhr61に対応する画素が、閾値飛行時間THt以内の範囲に閾値強度It1を超える頻度の画素に相当する。他方、図18のヒストグラムhr62が得られた場合には、閾値飛行時間THt以内の範囲に閾値強度It1を超える頻度が記録されていないので、ヒストグラムhr62に対応する画素は、閾値飛行時間THt以内の範囲に閾値強度It1を超える強度の画素に相当しない。
 閾値飛行時間THt以内の範囲に平均強度が閾値強度It1を超えるピークを有する画素があると判定されると(ステップS215a:YES)、ヒストグラム生成部230は、当該画素、すなわち、閾値飛行時間THt以内の範囲に平均強度が閾値強度It1を超えるピークを有する画素の位置を、クラッタ位置として追加設定する(ステップS220)。
 他方、上記ステップS215aにおいて、閾値飛行時間THt以内の範囲に平均強度が閾値強度It1を超えるピークを有する画素が無いと判定された場合(ステップS215a:NO)、処理はステップS205aに戻る。したがって、この場合、メモリ290に記憶されているクラッタ位置は更新されないこととなる。
 以上説明した第5実施形態の測距装置10は、第4実施形態の測距装置10と同様な効果を有する。加えて、第5実施形態の測距装置10は、各画素についての飛行時間と平均強度のヒストグラムを生成し、閾値飛行時間THt以内の範囲に閾値強度It1を超える強度の画素があると判定された場合に、かかる画素の位置をクラッタ位置として設定するので、泥等の異物が窓92に付着したことに起因して生じたクラッタの位置を特定することができる。このため、図6に示す測距処理のステップS110において、反射光の位置とクラッタ位置とが一致するか否かの判定を精度よく実行できる。
F.第6実施形態:
 第6実施形態の測距装置10の構成は、測距処理の詳細手順を除き、第1実施形態の測距装置10と同じであるので、同一の構成には同一の符号を付し、その詳細な説明を省略する。図19に示す第6実施形態の測距処理は、ステップS132を追加して実行する点において、図6に示す第1実施形態の測距処理と異なる。第6実施形態の測距処理のその他の手順は、第1実施形態と同じであるので、同一の手順には同一の符号を付し、その詳細な説明を省略する。
 図19に示すように、ステップS130において、算出された距離が閾値距離よりも小さいと判定された場合(ステップS130:YES)、距離算出部270は、複合山部のパルス幅が予め設定されている閾値幅よりも小さいか否かを判定する(ステップS132)。本実施形態において、ステップS132のパルス幅は、第1閾値受光強度、すなわち、山部における立ち上がり時間と立ち下り時間とを決定するために用いられる飛行時間の幅を意味する。また、本実施形態において、ステップS132の閾値幅は、複合山部が、クラッタと対象物体からの反射光とから形成された山部と、後述する多重反射により形成された山部とのいずれであるかを弁別するために設定される飛行時間の幅を意味する。一般には、多重反射により形成された山部の幅は、クラッタと対象物体からの反射光とから形成された山部の幅よりも大きい。そして、弁別可能な閾値幅は、例えば、実験やシミュレーションにより決定されて設定される。上述の「多重反射」とは、車両の後方に設けられたリフレクタ等の高い反射率を有する物体(以下、「高反射率物体」と呼ぶ)と測距装置10との間においてパルス光およびかかる反射光が複数回往復して反射し合うことを意味する。高反射率物体にパルス光が照射された場合に、その反射光の強度は大きい。このため、かかる反射光が窓92や筐体90内部において反射し、その反射光が再び高反射率物体を照射して、再び測距装置10に向かう反射光が生じる場合がある。このような反射を繰り返して多重反射が起こると、高反射率物体と測距装置10との間で往復するに従って、飛行時間Tfが長くなり、且つ、受光強度が小さくなるような複数のヒストグラムを統合したヒストグラムが生成されることとなる。
 例えば、図20の例では、1回目の反射光により得られるヒストグラムhr71と、2回目の反射光により得られるヒストグラムhr72と、3回目の反射光により得られるヒストグラムhr73とが合成されて、複合山部mp7を示す幅の広いヒストグラムhr7が得られている。図20の例では、第1閾値受光強度I3と一致する受光強度が得られた飛行時間Tu7、Td7のうち、時間的に小さな飛行時間Tu7が立ち上がり時間として特定され、時間的な大きな飛行時間Td7が立ち下り時間として特定される。
 図19に示すステップS132において、複合山部のパルス幅が予め設定されている閾値幅よりも小さいと判定された場合(ステップS132:YES)、上述のステップS135が実行される。複合山部のパルス幅が予め設定されている閾値幅よりも小さい場合には、複合山部は、クラッタと対象物体からの反射光とで構成されている可能性が高い。したがって、この場合、第1実施形態と同様に、ステップS125により算出された距離を、対象物距離として特定するようにしている。
 これに対して、複合山部のパルス幅が予め設定されている閾値幅よりも小さくないと判定された場合(ステップS132:NO)、上述のステップS115が実行され、通常の方法により対象物距離が算出される。すなわち、立ち上がり時間Tu7と立ち下り時間Td7との中間時間がピーク時間として特定され、かかるピーク時間を利用して対象物距離が特定される。複合山部のパルス幅が予め設定されている閾値幅よりも小さくないと判定された場合、すなわち、複合山部のパルス幅が大きく多重反射が生じている可能性が高い場合には、最初の反射光に基づき距離を求めることが望ましい。しかし、上述のステップS135では、立ち下り位置を基飛行時間として算出された距離が対象物距離として特定されるため、特定精度が低くなるおそれがある。これは、多重反射を繰り返したために飛行時間が長くなったことに起因して大きくなった立ち下り時間が基飛行時間となるためである。そこで、本実施形態では、このような場合にはステップS135を実行せず、通常の方法により距離を算出して対象物距離として扱うことにしている。
 以上説明した第6実施形態の測距装置10は、第1実施形態の測距装置10と同様な効果を有する。加えて、パルス幅が閾値幅よりも小さいか否かを判定し、小さくない場合には、通常の方法により算出された距離を対象物距離として特定するので、多重反射により得られた複合山部について、立ち下り時間を基飛行時間として対象物距離を算出して特定する構成に比べて、精度よく対象物距離を特定できる。
G.第7実施形態:
 第6実施形態の測距装置10aの構成は、測距処理の詳細手順を除き、第2実施形態の測距装置10aと同じであるので、同一の構成には同一の符号を付し、その詳細な説明を省略する。図21に示す第7実施形態の測距処理は、ステップS108、S120、S123に代えて、ステップS117を実行する点において、図9に示す第2実施形態の測距処理と異なる。第7実施形態の測距処理のその他の手順は、第2実施形態と同じであるので、同一の手順には同一の符号を付し、その詳細な説明を省略する。
 図21に示すように、ステップS107において、根元幅が閾値幅よりも大きいと判定された場合(ステップS107:YES)、基飛行時間決定部260は、立ち上がり時間と、立ち下がり時間とをそれぞれ基飛行時間に決定する(ステップS117)。ステップS117の完了後、上述のステップS125、S135が実行される。ステップS117において、立ち上がり時間と、立ち下がり時間とがそれぞれ基飛行時間に決定されるので、ステップS125では、それぞれの基飛行時間に基づき距離が算出され、ステップS135では、算出された2つの距離がいずれも対象物距離として特定されることとなる。図10に示すように、2つの対象物OBJ1、OBJ2が存在する場合、2つの対象物OBJ1、OBJ2のそれぞれについて距離を求めることが求められる場合がある。これは、例えば、前後に存在する2つの物体までの距離をそれぞれ正確に特定するという要請による。本実施形態によれば、立ち上がり時間、すなわち、手前側に位置する対象物に由来する時間を基時間として距離が求められるので、手前側に位置する対象物についての対象物距離が精度よく求められる。また、立ち下がり時間、すなわち、奥側に位置する対象物に由来する時間を基時間として距離が求められるので、奥側に位置する対象物についての対象物距離が精度よく求められる。
 以上説明した第7実施形態の測距装置10aは、第2実施形態の測距装置10aと同様な効果を有する。加えて、根元幅が閾値幅よりも大きいと判定された場合、すなわち、複数の対象物からの反射光のヒストグラムに含まれる山部が合成して複合山部が得られている可能性が高い場合には、立ち上がり時間と、立ち下がり時間とがそれぞれ基飛行時間に決定されるので、前後に存在する2つの物体までの距離をそれぞれ正確に特定できる。
H.その他の実施形態:
 (H1)複合山部が存在しないと推定される場合、すなわち、第1実施形態においては反射光の受光位置がクラッタの位置と一致しないと判定された場合(ステップS110:NO)、また、第2および第3実施形態においては根元幅が閾値幅よりも大きくないと判定された場合(ステップS107:NO)、ステップS115が実行されて、山部の立ち上がり時間と立ち下がり時間の中間時間を基飛行時間として対象物距離が算出されていたが、本開示はこれに限定されない。このような場合において、山部の立ち上がり時間と立ち下がり時間とのうちのいずれかを基飛行時間として対象物距離を算出してもよい。
 (H2)第2および第3実施形態では、複合山部の存在が推定される場合において、第1中間時間と第2中間時間との比較結果に応じて、基飛行時間を、複合山部の立ち上がり時間または立ち下がり時間に決定していたが、本開示はこれに限定されない。例えば、複合山部の存在が推定される場合においては、第1中間時間と第2中間時間との比較を実行することなく立ち上がり時間を基飛行時間として決定してもよい。かかる構成によれば、測距装置10aから見てより前方側に位置する対象物までの距離を対象物距離として求められ易くでき、より近くに位置して車両にとってより影響の大きな対象物までの距離を求めることができる。
 (H3)第2および第3実施形態では、基飛行時間決定部260は、第1中間時間と第2中間時間とが等しい場合には、立ち上がり時間を基飛行時間として決定していたが、立ち上がり時間に代えて、立ち下がり時間を基飛行時間として決定してもよい。
 (H4)各実施形態における測距装置10、10aの構成はあくまでも一例であり、様々に変更可能である。例えば、第1実施形態の測距処理において、ステップS130を省略してもよい。また、ステップS135において、改めて基飛行時間に基づき距離を算出し、対象物距離として特定してもよい。また、測距装置10、10aは、車両に限らず、飛行機や船舶など任意の移動体に搭載されてもよい。或いは、セキュリティ等の用途のため、固定設置されて用いられてもよい。また、測距装置10、10aの筐体が窓部材を有さず、単に開口が設けられている構成であってもよい。
 (H5)本開示に記載の演算部20及びこれらの手法は、コンピュータプログラムにより具体化された一つ乃至は複数の機能を実行するようにプログラムされたプロセッサ及びメモリを構成することによって提供された専用コンピュータにより、実現されてもよい。あるいは、本開示に記載の演算部20及びこれらの手法は、一つ以上の専用ハードウエア論理回路によってプロセッサを構成することによって提供された専用コンピュータにより、実現されてもよい。もしくは、本開示に記載の演算部20及びこれらの手法は、一つ乃至は複数の機能を実行するようにプログラムされたプロセッサ及びメモリと一つ以上のハードウエア論理回路によって構成されたプロセッサとの組み合わせにより構成された一つ以上の専用コンピュータにより、実現されてもよい。また、コンピュータプログラムは、コンピュータにより実行されるインストラクションとして、コンピュータ読み取り可能な非遷移有形記録媒体に記憶されていてもよい。
 (H6)第4実施形態のクラッタ位置設定処理のステップS210においてヒストグラムを生成する際に、基となる対象物距離が特定されたのが、車両が走行中であるか否かに応じて、積算する度数に重みづけを行ってもよい。具体的には、車両が走行中である場合には、走行中でない場合(つまり、停車中の場合)に比べて、より大きな度数を積算するようにしてもよい。車両が停車中の場合には、周囲に存在する対象物までの距離は経時的な変化は小さく、近距離の位置に対象物が存在する場合には、近距離の範囲において高い頻度となるヒストグラムが作成される可能性が高い。このため、実際には窓92に異物が付着していないにも関わらずに近距離の範囲に高い頻度となるヒストグラムが得られ、誤ったクラッタ位置を特定するおそれがある。他方、走行中である場合には、周囲に存在する対象物までの距離は経時的な変化は大きい。このため、近距離の位置に対象物が存在したとしても短時間であるため、近距離の範囲において高い頻度となるヒストグラムが生成される可能性は低い。しかし、このような走行中において近距離の範囲において高い頻度となるのは、窓92に異物が付着してクラッタが発生している可能性が高い。そこで、そのような場合に高い頻度となるように、走行中の場合には、走行中でない場合(つまり、停車中の場合)に比べて、より大きな度数を積算するようにしてもよい。なお、上述のように走行中であるか否かに応じて重みづけを行うことに代えて、走行中に対象物距離が特定された場合にのみ、クラッタ位置設定処理を行い、停車中にはクラッタ位置設定処理を行わないようにしてもよい。
 (H7)第6実施形態のステップS132において用いる「閾値幅」は、予め実験等により特定された1つの固定値であったが、本開示はこれに限定されない。閾値幅を、複合山部の立ち下がり時間に応じて変化させて用いても良い。具体的には、複合山部の立ち下がり時間が大きい(長い)ほど、閾値幅を大きく設定するようにしてもよい。複合山部の立ち下がり時間が大きい(長い)場合には、多重反射の回数が多いことが想定される。そこで、このような場合には、より大きな閾値幅を設定することにより、互いに同様な方角で距離が異なる複数の対象物が奥行方向の比較的長い範囲に存在する場合と、多重反射が生じている場合とを、精度よく弁別できる。
 (H8)第6実施形態では、ステップS132において、複合山部のパルス幅が予め設定されている閾値幅よりも小さくないと判定された場合(ステップS132:NO)、ステップS115が実行され、通常の方法、つまり、中間時間を基飛行時間として距離が算出されていたが、本開示はこれに限定されない。この場合には、通常の方法とは異なり、立ち上がり時間を基時間として距離を算出してもよい。上述のように、多重反射が生じている場合には、第1回目の反射光に基づき算出される距離が、対象物距離により近い値として求められる。これは、第2回、第3回と、回数が増えるにしたがって、反射光の往復時間が加算されてしまい、正確な飛行時間(反射物までの往復時間)からずれてしまうからである。そこで、上述のような構成とすることにより、多重反射が生じている状況において、対象物距離を精度良く特定できる。
 本開示は、種々の形態で実現することも可能である。例えば、測距システム、測距装置を搭載する移動体、測距方法、これらの装置や方法を実現するためのコンピュータプログラム、かかるコンピュータプログラムを記録した一時的でない記録媒体等の形態で実現することができる。
 本開示は、上述の実施形態に限られるものではなく、その趣旨を逸脱しない範囲において種々の構成で実現することができる。例えば、発明の概要の欄に記載した形態中の技術的特徴に対応する各実施形態中の技術的特徴は、上述の課題の一部又は全部を解決するために、あるいは、上述の効果の一部又は全部を達成するために、適宜、差し替えや、組み合わせを行うことが可能である。また、その技術的特徴が本明細書中に必須なものとして説明されていなければ、適宜、削除することが可能である。

Claims (8)

  1.  測距装置(10、10a)であって、
     パルス光を射出する発光部(40)と、
     前記パルス光の反射光を含む光を受光する受光部(60)と、
     前記受光部で受光された光の飛行時間を利用して、前記パルス光を反射して前記反射光を出力する対象物までの距離である対象物距離を演算する演算部(20)と、
     を備え、
     前記演算部は、
       複数の前記飛行時間における各々の前記受光部における受光強度を表すヒストグラムを生成するヒストグラム生成部(230)と、
       前記ヒストグラムにおいて、ピーク飛行時間を境に前記飛行時間に対する前記受光強度の変化が増加から減少に変わる山部であって前記測距装置からの距離の互いの差が所定範囲内である複数の物体からの光の強度の分布が合成されて得られる山部である複合山部が存在するかを推定する複合山部推定部(240)と、
       前記山部において、背景光を含むノイズ光の受光強度よりも大きな第1閾値受光強度と一致する2つの飛行時間であって、より短い飛行時間である立ち上がり時間と、より長い飛行時間である立ち下がり時間と、を特定する飛行時間特定部(250)と、
       前記複合山部が存在すると推定された場合に、前記複合山部について特定された前記立ち上がり時間と前記立ち下がり時間とのうちのいずれかに基づき、前記対象物距離を算出するための基となる基飛行時間を決定する基飛行時間決定部(260)と、
       前記基飛行時間を利用して、前記複数の物体のうちの一の物体の距離を、前記対象物距離として算出する距離算出部(270)と、
     を備える、測距装置。
  2.  請求項1に記載の測距装置において、
     前記基飛行時間決定部は、前記複合山部が存在しないと推定された場合に、前記立ち上がり時間と前記立ち下がり時間との中間時間を、前記基飛行時間として決定する、測距装置。
  3.  請求項1または請求項2に記載の測距装置において、
     前記基飛行時間決定部は、前記複合山部が存在すると推定された場合に、前記立ち上がり時間を、前記基飛行時間として決定する、測距装置。
  4.  請求項1から請求項3までのいずれか一考に記載の測距装置において、
     射出する前記パルス光と前記受光部で受光する光とを透過する窓(92)を有する筐体(90)であって、前記発光部と前記受光部とを収容する筐体を、さらに備え、
     前記複合山部推定部は、前記複合山部として、前記筐体内部における前記パルス光の反射光であるクラッタと、前記測距装置の外部の物体である外部物体からの前記パルス光の反射光と、を含む第1複合山部が存在するか否かを推定し、
     前記基飛行時間決定部は、前記第1複合山部が存在すると推定された場合に、前記立ち下がり時間に基づき、前記基飛行時間を決定する、測距装置。
  5.  請求項1から請求項4までのいずれか一項に記載の測距装置において、
     前記複合山部において、前記第1閾値受光強度と一致する2つの飛行時間の中間時間である第1中間時間と、前記第1閾値受光強度よりも小さな第2閾値受光強度と一致する2つの飛行時間の中間時間である第2中間時間とを特定する中間時間特定部(235)を、さらに備え、
     前記複合山部推定部は、前記複合山部として、前記測距装置の外部の物体である複数の外部物体からの前記パルス光の反射光を含む第2複合山部が存在するか否かを推定し、
     前記基飛行時間決定部は、前記第2複合山部が存在すると推定された場合において、
       前記第1中間時間が前記第2中間時間よりも小さい場合には、前記立ち上がり時間に基づき、前記基飛行時間を決定し、
       前記第1中間時間が前記第2中間時間よりも大きい場合には、前記立ち下がり時間に基づき、前記基飛行時間を決定する、測距装置。
  6.  請求項1から請求項4までのいずれか一項に記載の測距装置において、
     前記複合山部において、前記第1閾値受光強度とは異なる受光強度である第3閾値強度と一致する2つの飛行時間の中間時間である第1中間時間と、前記第3閾値強度よりも小さな第2閾値受光強度と一致する2つの飛行時間の中間時間である第2中間時間とを特定する中間時間特定部(235)を、さらに備え、
     前記複合山部推定部は、前記複合山部として、前記測距装置の外部の物体である複数の外部物体からの前記パルス光の反射光を含む第2複合山部が存在するか否かを推定し、
     前記基飛行時間決定部は、前記第2複合山部が存在すると推定された場合において、
       前記第1中間時間が前記第2中間時間よりも小さい場合には、前記立ち上がり時間に基づき、前記基飛行時間を決定し、
       前記第1中間時間が前記第2中間時間よりも大きい場合には、前記立ち下がり時間に基づき、前記基飛行時間を決定する、測距装置。
  7.  請求項5または請求項6に記載の測距装置において、
     前記基飛行時間決定部は、前記第2複合山部が存在すると推定された場合において、前記第1中間時間と前記第2中間時間とが同じ場合には、前記立ち上がり時間と前記立ち下がり時間とのうちのいずれかに基づき、前記基飛行時間を決定する、測距装置。
  8.  請求項4に記載の測距装置において、
     前記基飛行時間決定部は、前記反射光を含む光の受光位置が、設定されているクラッタ位置と一致するか否かを判定し、前記受光位置が前記クラッタ位置と一致すると判定された場合に、前記基飛行時間を決定し、
     前記距離算出部は、前記受光位置が前記クラッタ位置と一致すると判定された場合に、前記対象物距離を算出し、
     前記ヒストグラム生成部は、前記受光部による受光単位である画素ごとに、前記ヒストグラムを生成し、
     前記距離算出部は、前記画素ごとに、前記対象物距離を算出し、
     各前記画素の受光強度の統計値と、各前記画素について算出された前記対象物距離の統計値とのうちの少なくとも一方に基づき、前記クラッタ位置を設定するクラッタ位置設定部を、さらに備える、測距装置。
PCT/JP2021/020214 2020-06-12 2021-05-27 測距装置 WO2021251161A1 (ja)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202180042031.1A CN115698756A (zh) 2020-06-12 2021-05-27 测距装置
US18/064,010 US20230108583A1 (en) 2020-06-12 2022-12-09 Distance measurement device

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
JP2020102143 2020-06-12
JP2020-102143 2020-06-12
JP2021-070064 2021-04-19
JP2021070064A JP7375789B2 (ja) 2020-06-12 2021-04-19 測距装置

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US18/064,010 Continuation US20230108583A1 (en) 2020-06-12 2022-12-09 Distance measurement device

Publications (1)

Publication Number Publication Date
WO2021251161A1 true WO2021251161A1 (ja) 2021-12-16

Family

ID=78845472

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/020214 WO2021251161A1 (ja) 2020-06-12 2021-05-27 測距装置

Country Status (3)

Country Link
US (1) US20230108583A1 (ja)
CN (1) CN115698756A (ja)
WO (1) WO2021251161A1 (ja)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2019229891A1 (ja) * 2018-05-30 2019-12-05 株式会社ニコンビジョン 光検出装置及び方法並びに測距装置及び方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101196562A (zh) * 2007-12-14 2008-06-11 武汉大学 一种基于改进的em算法的激光雷达波形数据分解的方法
JP2011215005A (ja) * 2010-03-31 2011-10-27 Hokuyo Automatic Co 信号処理装置、及び走査式測距装置
US20170176578A1 (en) * 2015-12-18 2017-06-22 Stmicroelectronics (Research & Development) Limited Ranging apparatus
JP2018017534A (ja) * 2016-07-26 2018-02-01 オムロン株式会社 測距センサおよび測距方法
US20180253404A1 (en) * 2017-03-01 2018-09-06 Stmicroelectronics (Research & Development) Limited Method and apparatus for processing a histogram output from a detector sensor

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101196562A (zh) * 2007-12-14 2008-06-11 武汉大学 一种基于改进的em算法的激光雷达波形数据分解的方法
JP2011215005A (ja) * 2010-03-31 2011-10-27 Hokuyo Automatic Co 信号処理装置、及び走査式測距装置
US20170176578A1 (en) * 2015-12-18 2017-06-22 Stmicroelectronics (Research & Development) Limited Ranging apparatus
JP2018017534A (ja) * 2016-07-26 2018-02-01 オムロン株式会社 測距センサおよび測距方法
US20180253404A1 (en) * 2017-03-01 2018-09-06 Stmicroelectronics (Research & Development) Limited Method and apparatus for processing a histogram output from a detector sensor

Also Published As

Publication number Publication date
US20230108583A1 (en) 2023-04-06
CN115698756A (zh) 2023-02-03

Similar Documents

Publication Publication Date Title
JP7375789B2 (ja) 測距装置
US11609329B2 (en) Camera-gated lidar system
US7158217B2 (en) Vehicle radar device
WO2021193289A1 (ja) 測距装置
US8810445B2 (en) Method and apparatus for recognizing presence of objects
JP6824236B2 (ja) レーザ距離測定装置
JP2003302468A (ja) 測距装置
US10139477B2 (en) Erroneous detection restraining circuit for laser range finder
JP6874592B2 (ja) 時間測定装置、距離測定装置、移動体、時間測定方法、及び距離測定方法
JP7095640B2 (ja) 物体検出装置
US20230065210A1 (en) Optical distance measuring device
WO2021251161A1 (ja) 測距装置
US20210055406A1 (en) Object sensing device
US11802970B2 (en) Alternating power-level scanning for time-of-flight lidar systems
JP7452069B2 (ja) 道路勾配推定装置、道路勾配推定システムおよび道路勾配推定方法
JP2019028039A (ja) 距離測定装置及び距離測定方法
CN113678021A (zh) 测距装置以及测距装置中的异常判定方法
WO2023149335A1 (ja) 測距装置及び測距方法
US20230266450A1 (en) System and Method for Solid-State LiDAR with Adaptive Blooming Correction
WO2021235317A1 (ja) 光学測距装置
US20240067094A1 (en) Gating camera, vehicle sensing system, and vehicle lamp
US20240094365A1 (en) Contamination detection apparatus
WO2024090116A1 (ja) データ処理装置、光学センサ、データ処理方法、データ処理プログラム
EP4203454A1 (en) Vehicle-mounted sensing system and gating camera
CN116615644A (zh) 用于检测激光雷达视窗上的污垢的方法和装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 21822635

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21822635

Country of ref document: EP

Kind code of ref document: A1