WO2019008670A1 - レーザレーダ装置 - Google Patents

レーザレーダ装置 Download PDF

Info

Publication number
WO2019008670A1
WO2019008670A1 PCT/JP2017/024496 JP2017024496W WO2019008670A1 WO 2019008670 A1 WO2019008670 A1 WO 2019008670A1 JP 2017024496 W JP2017024496 W JP 2017024496W WO 2019008670 A1 WO2019008670 A1 WO 2019008670A1
Authority
WO
WIPO (PCT)
Prior art keywords
light
signal
optical
laser radar
radar device
Prior art date
Application number
PCT/JP2017/024496
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
Application filed by 三菱電機株式会社 filed Critical 三菱電機株式会社
Priority to JP2017566166A priority Critical patent/JP6395958B1/ja
Priority to US16/618,303 priority patent/US11899112B2/en
Priority to EP17916881.0A priority patent/EP3637134B1/en
Priority to PCT/JP2017/024496 priority patent/WO2019008670A1/ja
Publication of WO2019008670A1 publication Critical patent/WO2019008670A1/ja

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
    • 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
    • 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/50Systems of measurement based on relative movement of target
    • G01S17/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • 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/95Lidar systems specially adapted for specific applications for meteorological use
    • 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/4812Constructional features, e.g. arrangements of optical elements common to transmitter and receiver transmitted and received beams following a coaxial path
    • 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/4814Constructional features, e.g. arrangements of optical elements of transmitters alone
    • G01S7/4815Constructional features, e.g. arrangements of optical elements of transmitters alone using multiple transmitters
    • 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/481Constructional features, e.g. arrangements of optical elements
    • G01S7/4817Constructional features, e.g. arrangements of optical elements relating to scanning
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Definitions

  • the present invention relates to radar technology, and more particularly to laser radar technology that uses laser light to observe the state of a target such as particulate matter (for example, an aerosol).
  • a target such as particulate matter (for example, an aerosol).
  • Laser light is emitted to the atmosphere with particulate matter such as aerosol in the air (fine particles consisting of liquid or solid suspended in the air, cloud particles or atmospheric molecules), and scattered light reflected by the target is
  • particulate matter such as aerosol in the air (fine particles consisting of liquid or solid suspended in the air, cloud particles or atmospheric molecules)
  • scattered light reflected by the target is
  • a laser radar device that can receive and measure information indicating atmospheric conditions such as wind speed based on the received scattered light.
  • Patent Document 1 Japanese Patent Laid-Open No. 2009-162678.
  • a beam scanning optical system that scans the atmosphere with the transmission light of a laser beam and receives scattered light reflected by the atmosphere as reception light, and two wedges An optical axis correction unit that refracts the received light with a prism to compensate for an angle difference (optical axis angular deviation) between the transmitted light and the received light, and a branch from the output light of the optical axis correction unit and the transmitted light And a signal processing unit for performing heterodyne detection on the output light of the optical coupler.
  • the optical axis correction unit can optically compensate for the angular difference generated between the transmission light and the reception light with beam scanning, so the signal reception unit It is possible to measure the wind speed without lowering it.
  • WO 2016/117159 (eg, FIGS. 1 and 3)
  • the optical axis correction unit compensates for the difference in angle by controlling the rotational states of the two wedge prisms.
  • the optical axis correction unit is configured by two wedge prisms, an optical component that guides the received light to the wedge prisms, and a mechanism that mechanically rotates the wedge prisms, the response is improved.
  • the distance to the target is 100 m
  • the light propagation time from the point when the transmission light is emitted toward the target to the point when the scattered light reflected by the target is received is about 0.66 ⁇ sec.
  • the response time required to compensate for the angular difference is on the order of tens of nanoseconds (nanoseconds), depending on the received field of view and the scanning speed.
  • nanoseconds nanoseconds
  • a radar device has a plurality of light emitting ends, and a light source array that simultaneously emits a plurality of laser beams from the plurality of light emitting ends, and the plurality of laser beams from the plurality of laser beams.
  • An optical splitter that separates a portion of the light as transmission light and separates a plurality of local light components from the plurality of laser lights, an optical modulator that modulates the transmission light to generate modulated transmission light, and A transmission / reception optical system for receiving the modulated transmission light reflected by a target present in the external space as reception light while scanning the outside space with the modulation transmission light, combining the reception light and the plurality of local light components
  • An optical multiplexer that generates a plurality of interference light components by wave and a plurality of interference light components that are disposed at positions that optically correspond to a plurality of reception fields different from one another of the transmission and reception optical system, respectively
  • An optical receiver array including a plurality of optical receivers that generate a plurality of detection signals by wave forming, and a detection signal from among the plurality of detection signals according to the scanning speed of the transmission / reception optical system with respect to the external space It is characterized by comprising: a switch circuit to be selected; and a signal processing unit that calculates an observation
  • the present invention it is possible to select a high SNR detected signal having a high signal strength from among a plurality of detected signals corresponding to a plurality of reception visual fields according to the scanning speed of the transmission and reception optical system with respect to the external space. Therefore, even when an angle difference occurs between the modulated transmission light and the reception light, the observation amount indicating the state of the target can be calculated with high accuracy based on the selected detection signal. In addition, since the phase difference is compensated without requiring a mechanical structure that compensates the angle difference optically, high responsiveness can be realized.
  • FIG. 2 is a view showing a specific configuration example of a main part of the laser radar device of the first embodiment.
  • FIG. 2 is a view schematically showing a configuration example of an optical scanner of the laser radar device of the first embodiment.
  • FIG. 7 is a diagram showing an example of a plurality of reception fields of view of the optical antenna in the first embodiment.
  • FIG. 2 schematically shows a configuration of a signal processing unit in the laser radar device of the first embodiment.
  • FIG. 2 is a block diagram showing an example of a hardware configuration of a signal processing unit in Embodiment 1.
  • FIG. 7 is a diagram for describing an example of a selection method of reception channels according to the first embodiment.
  • FIG. 5 is a diagram for explaining integration processing according to the first embodiment.
  • FIG. 7 is a block diagram schematically showing a configuration of a signal processing unit in a laser radar device of a second embodiment.
  • FIG. 14 is a diagram for describing an example of a method of selecting a reception channel according to Embodiment 2.
  • It is a block diagram which shows schematic structure of the laser radar apparatus of Embodiment 3 which concerns on this invention.
  • It is a block diagram which shows schematic structure of the laser radar apparatus of Embodiment 4 which concerns on this invention.
  • FIG. 1 is a block diagram showing a schematic configuration of a laser radar device 1 according to a first embodiment of the present invention.
  • FIG. 2 is a diagram showing a specific configuration example of the main part of the laser radar device 1.
  • the laser radar device 1 includes a light source array 10, an optical distributor 11, an optical modulator 12, an optical circulator 13, an optical antenna 14, an optical scanner 15, an optical multiplexer 16, an optical receiver array 17, A switch circuit 18, an A / D converter (analog-to-digital converter) 19, and a signal processing unit 20 are provided.
  • the transmission / reception optical system of the present embodiment is configured by the combination of the optical circulator 13, the optical antenna 14 and the optical scanner 15.
  • the light source array 10 has K light emitting ends (K is an integer of 3 or more) arranged spatially (two-dimensionally or three-dimensionally), and K lasers from these K light emitting ends
  • the lights L 1 to L K are emitted toward the light distributor 11.
  • the optical antenna 14 has K reception fields different from one another, and receives scattered light of K reception channels CH 1 , CH 2 ,..., CH K respectively through the K reception fields.
  • the number of laser beams L 1 to L K is the same as the number of reception channels CH 1 to CH K. Although the number of reception channels is three or more in this embodiment, the number of reception channels may be two without being limited thereto.
  • the light source array 10 includes a reference light source 101 that outputs a reference laser light, and one reference laser light input from the reference light source 101 via an optical waveguide as three laser lights L 1 to Yes and optical distribution device 102, and a light emitting end E1 ⁇ E3 respectively outputting toward the three laser beams L 1 ⁇ L 3 of input from the optical distribution device 102 to the optical divider 11 which distributes the L 3 doing.
  • a semiconductor laser or a solid state laser may be used as the reference light source 101.
  • the polarization state of the reference laser light is held so as to be a linear polarization state that is polarized in one direction.
  • a melting type optical fiber branch coupler or a filter type branch coupler an optical component that branches light using a dielectric multilayer film filter
  • An array laser having a plurality of light emitting surfaces may be employed as the light source array 10.
  • the light distributor 11 separates the transmission light components TL 1 to TL K from the input laser lights L 1 to L K , and the local light components LL from the laser lights L 1 to L K respectively. 1 to LL K are separated.
  • the optical distributor 11 outputs the transmission light TL including the transmission light components TL 1 to TL K to the optical modulator 12 and outputs the local light LL including the local light components LL 1 to LL K to the optical multiplexer 16.
  • a branch mirror or a beam splitter having a dielectric multilayer film filter may be used as the light distributor 11, for example.
  • the light distributor 11 is composed of a light guiding optical system 111 and a beam splitter 112.
  • the light guiding optical system 111 is an optical component for guiding the laser beams L 1 to L 3 input from the light source array 10 to the light incident surface of the beam splitter 112.
  • the beam splitter 112 transmits a part of the laser beams L 1 to L 3 input from the light guide optical system 111 as transmission light components TL 1 to TL 3 , and the laser beams L 1 to L K have a predetermined branching ratio.
  • the separated local light components LL 1 to LL 3 are reflected in the direction of the optical multiplexer 16.
  • the predetermined branching ratio may be determined by system design.
  • the optical modulator 12 performs frequency modulation and intensity modulation only on a part of transmission light components TL k of the transmission light components TL 1 to TL K input from the optical distributor 11. Accordingly, it shifts the frequency of the transmission light component TL k, and a modulator for pulsing the transmission light component TL k.
  • the optical modulator 12 generates modulated transmission light ML k (a series of pulse lights) having a predetermined pulse width T at a predetermined pulse repetition frequency (Pulse Repetitive Frequency, PRF) by frequency modulation and intensity modulation, and performs modulation transmission
  • the light ML k is output to the light circulator 13.
  • the optical modulator 12 also outputs a pulse trigger signal TS generated for pulsing the transmission light component TL k to the signal processing unit 20.
  • the pulse trigger signal TS is a signal representing the timing of pulsing of the transmission light component TL k .
  • the pulse width T corresponds to the distance resolution value. It is possible to set in advance a pulse width T corresponding to a desired distance resolution value. Also, in the light modulator 12, fixed pulse widths T and PRF set at the time of design may be used. Furthermore, when the output light intensity of the optical modulator 12 is insufficient, an optical amplifier for amplifying the output light of the optical modulator 12 may be provided downstream of the optical modulator 12. The optical modulator 12 may shift the frequency of transmission light using a phase modulator having a frequency shift function.
  • an optical modulator 12 for example, an AO frequency shifter (Acousto-Optical frequency shifter) or a known optical phase modulator may be used.
  • the pulse trigger signal TS for example, a TTL (Transistor-Transistor Logic) signal having a TTL level of 5 volts can be used.
  • the optical circulator 13 is a nonreciprocal optical device having three ports. That is, the optical circulator 13 has an optical input port through which light from the optical modulator 12 is incident, an optical input / output port through which light is incident from the optical antenna 14 and emits light to the optical antenna 14, and an optical multiplexer 16. And a light output port for emitting light.
  • the optical antenna 14 converts the modulated transmission light ML k input from the optical circulator 13 into parallel light and outputs the parallel light to the optical scanner 15.
  • Optical scanner 15 has a function of changing the irradiation direction of the outer space of the modulated transmission light ML k input from the optical antenna 14 (i.e., line-of-sight direction). Further, the optical scanner 15 can be repeatedly scanned a predetermined range of the external space modulation transmission light ML k input from the optical antenna 14.
  • Such an optical scanner 15 may be configured by, for example, a wedge prism, a galvano mirror or a polygon mirror, or a combination of these optical components.
  • FIG. 3 is a view schematically showing a configuration example of the light scanner 15 of the first embodiment.
  • the optical scanner 15 has a rotor 151 including a wedge prism 152, and a rotational drive unit 153 that rotationally drives the rotor 151 at an arbitrary speed about the optical axis OA.
  • the rotary drive unit 153 is configured to include a rotary motor such as a stepping motor and a rotary encoder (not shown).
  • the optical surface on the light antenna 14 side of the wedge prism 152 is perpendicular to the optical axis OA, and the optical surface on the outer space side of the wedge prism 152 is inclined to the optical axis OA.
  • the rotational drive unit 153 can scan the outside space with the laser beam by rotating the wedge prism 152 around the optical axis OA.
  • the rotation drive unit 153 outputs, to the signal processing unit 20, angle information AD indicating the current value of the scanning angle of the light scanner 15 detected by the built-in rotary encoder.
  • the scanning speed of the optical scanner 15 may be a value preset at the time of design.
  • the control signal for controlling the rotation of the rotor 151 is generated by the rotation drive unit 153, but is not limited to this.
  • the configuration of the signal processing unit 20 may be changed to generate a control signal that controls the rotation of the rotor 151.
  • the optical antenna 14 may be an aerosol in the external space (particulates suspended in the atmosphere, particles of liquid or solid), scattered light reflected back from targets such as cloud particles or atmospheric molecules, or structures or topography, etc.
  • the scattered light (diffuse light) reflected and returned by the hard target of (1) is received through the light scanner 15.
  • the optical circulator 13 shown in FIG. 1 outputs the scattered light, ie, the received light RL input from the optical antenna 14 to the optical multiplexer 16.
  • Such an optical antenna 14 can be configured, for example, using an optical telescope or a camera lens.
  • the optical antenna 14 may have a condensing adjustment function.
  • the optical circulator 13 is configured as a space propagation type optical device having a polarization beam splitter 131 and a wavelength plate 132.
  • the wave plate 132 converts linearly polarized light input from the light modulator 12 via the polarization beam splitter 131 into circularly polarized light, and converts circularly polarized light input from the optical antenna 14 into linearly polarized light (horizontally polarized light). It is a quarter wave plate to convert.
  • the polarization beam splitter 131 transmits linearly polarized light input from the light modulator 12 and reflects circularly polarized light input from the light antenna 14 in the direction of the optical multiplexer 16.
  • the optical antenna 14 shown in FIG. 1 has K reception fields different from each other for scattered light coming from a target.
  • Output to the FIG. 4 is a view showing an example of three reception visual fields (first to third reception visual fields) of the optical antenna 14.
  • the received viewing angle ⁇ per one reception channel is set equal to the divergence angle ⁇ of the transmitted laser beam (laser beam modulated transmission light ML k).
  • these reception fields of view are set so that different angular directions can be observed with an angular offset of ⁇ .
  • the value of the angular offset (for example, ⁇ / 2) may be set such that adjacent received fields of view overlap each other.
  • the optical multiplexer 16 multiplexes the local light LL input from the optical distributor 11 and the reception light RL input from the optical circulator 13 to form K interference light components.
  • the interference light CL including CL 1 ,..., CL K is generated, and the interference light CL, that is, the light beat signal is output to the light receiver array 17.
  • the optical multiplexer 16 can be configured, for example, using a branch mirror or a beam splitter having a dielectric multilayer film filter.
  • the light receiver array 17 has K light receivers arranged at positions respectively corresponding to the K reception fields of the light antenna 14. These K optical receivers perform heterodyne detection (frequency demodulation) on K interference light components CL 1 , CL 2 ,..., CL K to obtain K analog detection signals CS 1 , CS 2 ,. It may generate an analog detection signal group CS including CS K.
  • the k-th interference light component CL k is an interference light component generated by multiplexing the k-th local light component LL k and the k-th reception light component RL k .
  • the optical coupler 16 and the optical receiver array 17 may be configured to perform balanced detection.
  • Balanced detection makes it possible to reduce phase noise after heterodyne detection.
  • Switch circuit 18 in accordance with the scanning speed of the optical scanner 15, analog detection signals CS 1, CS 2, ..., and select the analog detection signal SS from the CS K, the selected analog detection signal SS A / D It can be output to the converter 19.
  • the signal processing unit 20 sequentially calculates the speed information SD indicating the scanning speed of the optical scanner 15 based on the angle information AD, and supplies the speed information SD to the switch circuit 18.
  • the angular misalignment theta d is the angle difference between the transmission direction of the laser beam (laser beam modulated transmission light ML k), and the arrival direction of the scattered light to the optical antenna 14.
  • the switch circuit 18 on the basis of the value of the angular misalignment theta d, and selects the reception channel CH p corresponding to the arrival direction of the scattered light, analog detection signals CS 1, CS 2, ..., received from the CS K
  • the analog detection signal SS corresponding to the channel CH p is selected.
  • the switch circuit 18 can select the output signal (analog detection signal) SS of the light receiver that optically corresponds to the reception field of view indicating the arrival direction of the scattered light.
  • analog detection signals CS 1, CS 2, ..., if CS K are added all for analog detection signal of the reception channel (reception field) to which the scattered light does not arrive is also added, SNR of the additional signal Will deteriorate.
  • the switch circuit 18 of this embodiment selects the analog detection signal SS of the reception channel (reception field of view) to which the scattered light arrives, the deterioration of the SNR caused by the angular deviation ⁇ d is avoided.
  • a delay time ⁇ (unit: second) from the time of transmission of the laser beam to the time of reception of scattered light
  • ⁇ ⁇ V scan unit: radian
  • the switch circuit 18 based on the angular deviation theta d, the delay time tau, whether to select which receiving channel can be predetermined.
  • is a reception viewing angle per reception channel.
  • the switch circuit 18 may add the analog detection signals of the two or more selected reception channels and output the addition signal to the A / D converter 19. Thereby, a high SNR signal can be generated while securing a wide reception field of view.
  • the optical multiplexer 16 is configured to include a beam splitter 161, a reflection mirror 162, and condensing optical systems 163 and 164.
  • the optical receiver array 17 of FIG. 1 is composed of optical receivers Da1 to Da3 and optical receivers Db1 to Db3.
  • the sum of the optical receiver Da1, Db1 output signal corresponds to the analog detection signal CL 1
  • the sum of the optical receiver Da2, Db2 of the output signal corresponds to the analog detection signal CL 2
  • the optical receiver Da3, Db3 output the sum of the signal corresponds to an analog detection signal CL 3.
  • the switch circuit 18 of FIG. 1 is composed of a first switch 181, a second switch 182, an adder 183, and a channel selection circuit 184.
  • the reflection mirror 162 reflects the reception light RL input from the polarization beam splitter 131 in the direction of the beam splitter 161.
  • the beam splitter 161 transmits a part of the received light RL incident from the polarization beam splitter 131 via the reflecting mirror 162 in the direction of one of the condensing optical systems 163, and the remaining part of the received light RL is collected by the other condensing optical system.
  • the beam splitter 161 reflects a part of the local light LL incident from the light distributor 11 in the direction of one of the focusing optical systems 163, and the remaining part of the local light LL is reflected by the other focusing optical system 164. Transmit in the direction.
  • the beam splitter 161 outputs the interference light generated by combining the part of the reception light RL and the part of the local light LL to one of the condensing optical systems 163, and at the same time, the remaining part of the reception light RL. And the remaining light of the local light LL are output to the other focusing optical system 164.
  • the condensing optical system 163 condenses the interference light incident from the beam splitter 161 on the light receivers Da1 to Da3, and the condensing optical system 164 is incident from the beam splitter 161.
  • the interference light is condensed on the light receivers Db1 to Db3.
  • the optical receivers Da1, Da2, Da3 are respectively disposed at positions that optically correspond to the first reception field, the second reception field, and the third reception field in FIG. For this reason, the optical receivers Da1, Da2, Da3 detect three interference light components including three scattered light components respectively propagated from the first reception visual field, the second reception visual field, and the third reception visual field, respectively, and detect three interference light components. Can output an analog detection signal of
  • the light receivers Db1, Db2 and Db3 are also arranged at positions optically corresponding to the first reception field, the second reception field and the third reception field in FIG. For this reason, the optical receivers Db1, Db2 and Db3 detect three interference light components including three scattered light components respectively propagated from the first reception visual field, the second reception visual field and the third reception visual field, and Can output an analog detection signal of
  • the channel selection circuit 184 calculates an angular deviation ⁇ d corresponding to the delay time ⁇ from the time of transmission of the laser beam to the time of reception of the scattered light based on the speed information SD sequentially supplied from the signal processing unit 20 based on the calculated angular deviation theta d, selects a reception channel CH p corresponding to the arrival direction of the scattered light. Then, channel selection circuit 184 supplies a selection control signal CS for selecting an output representative of the analog detection signal CS p reception channel CH p to the first switch 181 and second switch 182.
  • the first switch 181 and the second switch 182 select the outputs Ss1 and Ss2 designated by the selection control signal CS from the outputs of the optical receivers Da1 to Da3 and Db1 to Db3, and select the selected outputs Ss1 and Ss1.
  • the signal Ss2 is supplied to the adder 183.
  • the A / D converter 19 samples the analog detection signal SS at the sampling frequency fs to generate a digital detection signal DS, and outputs the digital detection signal DS to the signal processing unit 20.
  • the A / D converter 19 for example, a double integration type A / D converter, a successive approximation type A / D converter, or a parallel comparison type A / D converter can be used.
  • the A / D converter 19 may generate the digital detection signal DS by sampling the analog detection signal SS in synchronization with the pulse trigger signal TS using the pulse trigger signal TS as a trigger.
  • the signal processing unit 20 illustrated in FIG. 1 has a function of calculating an observation amount indicating the state of the target based on the digital detection signal DS. As described later, the signal processing unit 20 calculates the sight line velocity and velocity vector of the target as an observation quantity indicating the state of the target such as particulate matter (for example, aerosol, cloud particles or atmospheric molecules) suspended in the air. Although it can be configured as such, the calculated observation quantities are not limited to the eye velocity and the velocity vector.
  • the light propagation time between the laser radar device 1 and the target can be measured based on the digital detection signal DS, and the distance (ranging value) to the target can be calculated based on the measured value of the light propagation time.
  • the signal processing unit 20 may be configured.
  • the laser radar device 1 may be configured such that the optical scanner 15 scans the surface of the target with the modulated transmission light ML k and the signal processing unit 20 measures the distance distribution with the target. This makes it possible to measure the three-dimensional shape of a target such as a mobile (for example, a vehicle, a human body or an animal), a terrain object or a structure.
  • a target such as a mobile (for example, a vehicle, a human body or an animal), a terrain object or a structure.
  • FIG. 5 is a block diagram schematically showing a configuration example of the signal processing unit 20 in the present embodiment.
  • the signal processing unit 20 shown in FIG. 5 calculates the signal spectrum (amplitude spectrum or power spectrum) of the digital detection signal DS, and calculates the eye velocity and the velocity vector indicating the movement state of the target based on the signal spectrum. Can.
  • the signal processing unit 20 generates M range bin signals (time domain) which respectively represent signal waveforms of the continuous M gate time domains (M is a positive integer of 3 or more) in the digital detection signal DS. signal) B 1, B 2, ...
  • a signal splitter 21 which divides the B M, these range bin signals B 1, B 2, ..., the frequency domain signal B M F 1, F 2, ..., respectively F M .., F M , the integrator 23 for calculating the integrated spectrum IS as the signal spectrum of the digital detection signal DS based on the frequency domain signals F 1 , F 2 ,.
  • the angle information AD sequentially supplied from the light scanner 15 and the vector calculator 27 for calculating the velocity vector v of the target.
  • Te and a scanning speed calculator 28 for calculating the speed information SD indicating the scanning rate of the optical scanner 15.
  • the velocity calculator 26 and the vector calculator 27 constitute an observable calculator according to this embodiment.
  • the hardware configuration of such a signal processing unit 20 can be realized by, for example, a processor having a semiconductor integrated circuit such as a digital signal processor (DSP), an application specific integrated circuit (ASIC), or a field-programmable gate array (FPGA). Just do it.
  • the hardware configuration of the signal processing unit 20 is realized by a processor including an arithmetic device such as a central processing unit (CPU) or a graphics processing unit (GPU) that executes program code of software or firmware read from the memory. It may be done.
  • FIG. 6 is a block diagram schematically showing a signal processing device 40 which is an example of a hardware configuration for realizing the function of the signal processing unit 20.
  • the signal processing device 40 includes a processor 41, a memory 42, an input interface unit 43, an output interface unit 44, and a signal path 45.
  • the signal path 45 is a bus for mutually connecting the processor 41, the memory 42, the input interface unit 43 and the output interface unit 44.
  • the input interface unit 43 has a function of transferring the digital detection signal DS input from the outside and the angle information AD (not shown) to the processor 41 via the signal path 45.
  • the processor 41 performs digital signal processing on the transferred digital detection signal DS to calculate the velocity vector v of the target, and calculates the velocity information SD based on the transferred angle information AD.
  • the processor 41 can output data indicating the velocity vector v to an external device (for example, a display device) via the signal path 45 and the output interface unit 44.
  • the processor 41 can also output the speed information SD to the switch circuit 18 via the signal path 45 and the output interface unit 44.
  • the memory 42 is a data storage area used when the processor 41 executes digital signal processing.
  • the memory 42 may have a data storage area for storing a program code of software or firmware executed by the processor 41.
  • a semiconductor memory such as a ROM (Read Only Memory) and an SDRAM (Synchronous Dynamic Random Access Memory) can be used.
  • the scanning speed calculator 28 sequentially calculates the speed information SD indicating the scanning speed of the light scanner 15 based on the angle information AD sequentially supplied from the light scanner 15, and supplies the speed information SD to the switch circuit 18.
  • the signal divider 21 receives as input the digital detection signal DS selected by the switch circuit 18, and M range bin signals (time domain signals) representing the signal waveforms of M continuous gate time domains respectively. ) Divide into B 1 , B 2 , ..., B M.
  • FIG. 7 is a graph for explaining the concept of range bins. In this graph, the horizontal axis represents time, and the vertical axis represents signal strength.
  • the transmission pulse Tw indicating the modulated transmission light ML k is repeatedly transmitted with a period of 1 / PRF. For each transmission (shot) of the transmission pulse Tw, the reception signal Rw of the scattered light reflected by the target is observed. Time t s is a measurement start time for each transmission (shot) of the transmission pulse Tw.
  • the received signal Rw is divided into range bins Rb 1 , Rb 2 ,..., Rb M over M gated time domains.
  • Range bin signals B 1, B 2, ..., B M these range bins Rb 1, Rb 2, ..., are those corresponding respectively to the Rb M.
  • the gate time (time width) for determining each range bin can be set as 2R res / c (c: speed of light).
  • a distance value corresponding to a delay timing until the A / D converter 19 starts A / D conversion according to the pulse trigger signal TS is represented by R min. It shall be.
  • a value of 40 can be set as the distance value R min .
  • the switch circuit 18 as illustrated in Figure 8, the receiving channels CH 1, CH 2, ..., analog detection signals CS 1, CS 2 of CH K, ..., CS K are input.
  • the switch circuit 18 receives to compensate for the angular deviation corresponding to the delay time ⁇ 1 or ⁇ 2 select the channel CH 1, selects an analog detection signal CS 1 of the reception channel CH 1 as an analog detection signal SS.
  • the switch circuit 18 selects the reception channel CH 2 to compensate for the angular deviation corresponding to the delay time ⁇ 3 or ⁇ 4 , and the reception channel CH 2.
  • the analog detection signal CS 2 is selected as an analog detection signal SS of.
  • the signal divider 21, M-number of range bins signals analog detection signal SS B 1, B 2, ..., is divided into B M.
  • Domain transform section 22 shown in FIG. 5 range bin signal B 1 obtained for each shot, B 2, ..., respectively into the frequency domain signals B M F 1, F 2, ..., is converted to F M.
  • the domain transformation unit 22 performs discrete Fourier transform on each of the range bin signals B 1 , B 2 ,..., B M to generate frequency domain signals F 1 , F 2 ,. Can be calculated.
  • FFT Fast Fourier Transform
  • N FFT for example, 256 points
  • Integrator 23 for each shot, the frequency domain signal F 1, F 2, ..., range bin signals B 1, B 2, based on the F M, ..., B M each spectrum S (m, n) is calculated.
  • m 1 to M
  • N is the number of integrations designated by the user in advance.
  • the spectrum S (m, n) may be an amplitude spectrum that represents the relationship between frequency and signal amplitude, or a power spectrum that represents the relationship between frequency and power.
  • the integrator 23 integrates (incoherently integrates) high SNRs by integrating the spectra S (m, 1) to S (m, N) in the direction of the shot numbers for each range bin (range bin number m).
  • the signal spectrum IS (m) is calculated.
  • FIG. 9 is a diagram showing an outline of integration processing in the integrator 23.
  • integrator 23 can calculate the high SNR of the integration signal spectrum IS (1) by integrating the spectrum S (1,1) ⁇ S (1 , N). When this integration result is averaged by the number of integrations N, the SNR improves by N 1/2 times the number of integrations N.
  • incoherent integration is adopted, but instead of the incoherent integration, coherent integration may be adopted.
  • the integrator 23 integrates the integrated signal spectra IS (1) to IS (M) to calculate an integrated spectrum IS, and outputs the integrated spectrum IS to the velocity calculator 26.
  • Velocity calculator 26 detects the frequency of the peak position or gravity center position of the integrated spectrum IS as Doppler frequency f d. Then, the speed calculator 26 uses the Doppler frequency f d, the target velocity component, or the radial velocity Vr in the viewing direction corresponding to the angle information AD (Unit: m / s) is calculated.
  • the eye velocity Vr can be calculated according to the following equation (1).
  • is the wavelength of the laser beam.
  • Vr ⁇ ⁇ f d / 2 (1)
  • the vector calculator 27 can calculate the velocity vector of the target based on the eye velocity Vr using a known vector synthesis method or VAD (Velocity Azimuth Display) method.
  • VAD Vector Azimuth Display
  • the eye velocity Vr is the horizontal wind velocity in the east-west direction (U), the horizontal wind velocity in the north-south direction (V), the vertical wind velocity (W), the elevation angle ( ⁇ ), and It is expressed by the following equation (2) using an azimuth angle ( ⁇ ) based on the north.
  • Vr Usin ⁇ sin ⁇ + V cos ⁇ sin ⁇ + W cos ⁇ (2)
  • the optical receiver array 17 has K optical receivers disposed at positions respectively corresponding to the K reception fields of the optical antenna 14, and the optical receivers include K interference light components. Are respectively detected to generate analog detection signals CS 1 to CS K.
  • the switch circuit 18 selects a high SNR analog detection signal SS having a high signal strength from the analog detection signals CS 1 to CS K respectively corresponding to the K receiving fields of view according to the scanning speed of the optical scanner 15. be able to.
  • the signal processing unit 20 can achieve high SNR
  • the observation amount indicating the state of the target can be calculated with high accuracy.
  • the laser radar device 1 compensates for the angular deviation with a short response time without requiring a mechanical structure for compensating for the angular difference as disclosed in Patent Document 1 above. Can. Therefore, even if the distance to the target is short, a high SNR digital detection signal DS can be obtained, so it is possible to calculate the observation amount with high accuracy.
  • FIG. 10 is a block diagram showing a schematic configuration of a laser radar device 2 according to a second embodiment of the present invention.
  • the laser radar device 2 includes the light source array 10, the light distributor 11, the light modulator 12, the light circulator 13, the light antenna 14, the light scanner 15, and the light combining as in the first embodiment.
  • An optical receiver 16 and an optical receiver array 17 are provided.
  • the configuration of the laser radar device 2 of this embodiment is the A / D converter 19A and signal processing of FIG. 10 in place of the switch circuit 18, the A / D converter 19 and the signal processing unit 20 of the first embodiment.
  • the configuration is the same as that of the laser radar device 1 of the first embodiment except that the portion 20A is included.
  • a / D converter 19A is at a predetermined sampling frequency, analog detection signals CS 1, CS 2, ..., to sample each of the CS K, analog detection signals CS 1, CS 2, ..., CS K ( receiver channel CH 1, CH 2, ..., a digital detection signal DS 1, DS 2 respectively corresponding to the CH K), ..., and generates a DS K, these digital detection signals DS 1, DS 2, ..., signal processing DS K Output to section 20A.
  • the A / D converter 19A for example, a double integration type A / D converter, a successive approximation type A / D converter, or a parallel comparison type A / D converter can be used.
  • a / D converter 19A is the pulse trigger signal TS as a trigger, by sampling the analog detection signal group CS in synchronism with the pulse trigger signal TS, a digital detection signal DS 1, DS 2, ..., DS K May be generated.
  • the signal processing unit 20A a digital detection signal DS 1, DS 2, ..., DS K respective signal spectrum (integrated spectrum) IS 1, IS 2, ... , and calculates the IS K, these signals spectrum IS 1, IS 2, ..., digital detection signal DS 1, DS 2 based on the iS K, ..., DS K respective signal-to-noise power ratio R 1, R 2, ..., and calculates the R K.
  • FIG. 11 is a block diagram schematically showing a configuration example of the signal processing unit 20A in the second embodiment.
  • the signal processing unit 20A includes a signal divider 21A, a domain conversion unit 22A, an integrator 23A, an SNR calculator 24, a selector 25, a speed calculator 26, and a vector calculator 27.
  • the configuration of the signal processing unit 20A is the signal dividing unit 21A, the area converting unit 22A, the integrator 23A, and the SNR calculator in FIG. 11 instead of the signal dividing unit 21, the area converting unit 22 and the integrator 23 of the first embodiment.
  • the configuration is the same as the configuration of the signal processing unit 20 of the first embodiment except that it has a selector 24 and a selector 25.
  • Signal splitter 21A the digital detection signal DS 1, DS 2 obtained for each shot, ..., each digital detection signal DS k of DS K, M number respectively representing the signal waveform of the M gate time domain .., B k, and M.
  • the signal divider 21A outputs a range bin signal group RB k including M range bin signals B k, 1 , B k, 2 ,..., B k, M to the region conversion unit 22A for each reception channel CH k .
  • the range conversion unit 22A determines, for each shot, each of the range bin signal groups RF k of the range bin signal groups RF 1 , RF 2 ,..., RF K as M frequency domain signals F k, 1 ,. transformed into the frequency domain signal group RF k containing M, and outputs the frequency domain signal group RF k in the multiplier 23A. That is, region conversion unit 22A performs discrete Fourier transform on each of range bin signals B k, 1 , B k, 2 ,..., B k, M to obtain frequency domain signals F k, 1 , F k. , 2 ,..., F k, M can be calculated. As the discrete Fourier transform, Fast Fourier Transform (FFT) with FFT bin number N FFT (for example, 256 points) can be used.
  • FFT Fast Fourier Transform
  • the integrator 23A calculates an integrated spectrum IS k based on each frequency domain signal group RF k for each shot, and outputs the integrated spectrum IS k to the SNR calculator 24 and the selector 25.
  • the integrator 23A for each shot, the frequency domain signal F k, 1, F k, 2, ..., F k, the range bin signals B k based on M, 1, B k, 2, ..., Spectrum S k (m, n) of each of B k and M is calculated.
  • N is the number of integrations designated by the user in advance.
  • the spectrum S (k, m, n) may be an amplitude spectrum or a power spectrum.
  • the integrator 23A integrates the spectra S k (m, 1) to S k (m, N) in the direction of the shot numbers (incoherently integrates) for each range bin (range bin number m) to achieve high SNR Integral signal spectrum IS k (m) is calculated.
  • accumulator 23A is able to calculate the high SNR of the integration signal spectrum IS k (1) by integrating the spectrum S k (1,1) ⁇ S k (1, N) it can.
  • the SNR improves by N 1/2 times the number of integrations N.
  • incoherent integration is adopted, but instead of the incoherent integration, coherent integration may be adopted.
  • integrator 23A is integral signal spectrum IS k (1) to calculate an integrated spectrum IS k by integrating the ⁇ IS k (M), and outputs the integrated spectrum IS k in SNR calculator 24 and selector 25 Do.
  • the SNR calculator 24 the integrated spectrum IS 1, IS 2, ..., signal-to-noise power ratio from the IS K R 1, R 2, ..., and calculates the R K.
  • a known method may be used for the method of calculating the signal-to-noise power ratio R k , and it is not particularly limited.
  • the SNR calculator 24 can detect, for example, a peak value of the integrated spectrum IS k and calculate the ratio of the peak value to the noise level measured in advance as a signal-to-noise power ratio R k .
  • Selector 25 the signal-to-noise power ratio R 1, R 2, ..., integrated with R K spectrum IS 1, IS 2, ..., selects the integrated spectrum of the at least one high SNR from the IS K.
  • the selector 25 When one integrated spectrum is selected from the integrated spectra IS 1 , IS 2 ,..., IS K , the selector 25 outputs the selected integrated spectrum as it is as an integrated spectrum IS.
  • the selector 25 adds or averages the selected integrated spectra to calculate one integrated spectrum IS. And outputs the integrated spectrum IS to the velocity calculator 26.
  • selector 25 as illustrated in FIG.
  • reception channel CH 1, CH 2, ..., the integrated spectrum IS 1, IS 2 of CH K, ..., IS K are input.
  • the selector 25 uses the signal-to-noise power ratios R 1 , R 2 ,..., R K to compensate for the angular deviation corresponding to the delay time, and gives an integrated spectrum including frequency components exceeding the threshold TH. Select a channel. In the example of FIG.
  • the selector 25 selects the reception channel CH 1 to compensate for the angular deviation corresponding to the delay time tau 1 or tau 2, the integrated spectrum IS the integrated spectrum IS 1 of the reception channel CH 1 Can be output as Further, the selector 25 selects the reception channels CH 1 and CH 2 in order to compensate the angular deviation corresponding to the delay time ⁇ 3 , and adds the integrated spectra IS 1 and IS 2 of these reception channels CH 1 and CH 2. Alternatively, the noise reduction spectrum can be calculated by averaging, and this noise reduction spectrum can be output as an integrated spectrum IS. Further, selector 25 may select the reception channel CH 2 to compensate for the angular deviation corresponding to the delay time tau 4, it selects the integrated spectrum IS 2 of the reception channel CH 2 as the integrated spectrum IS.
  • the selector 25 calculates the noise reduction spectrum by adding or averaging the plurality of selected integrated spectra IS 1 and IS 2 , and this noise reduction spectrum is calculated.
  • the integrated spectrum IS is output to the speed calculator 26. This can improve the SNR.
  • Velocity calculator 26 the frequency of the peak position or gravity center position of the integrated spectrum IS detected as Doppler frequency f d, with reference to the Doppler frequency f d, it can be calculated radial velocity Vr corresponding to the angle information AD . Also, the vector calculator 27 can calculate the velocity vector v of the target based on the eye velocity Vr using a known vector synthesis method or VAD method.
  • the eye velocity Vr and the velocity vector v of the target are calculated as the observation amount indicating the state of the target, but the calculated observation amount is limited to the eye velocity Vr and the velocity vector v It is not a thing.
  • the signal processing unit 20A may be configured to have a functional block capable of calculating the distance (ranging value) between the laser radar device 2 and the target based on the integrated spectrum IS.
  • the hardware configuration of the signal processing unit 20A described above may be realized by, for example, a processor having a semiconductor integrated circuit such as a DSP, an ASIC, or an FPGA.
  • the hardware configuration of the signal processing unit 20A may be realized by a processor including an arithmetic device such as a CPU or a GPU that executes a program code of software or firmware read from a memory.
  • the function of the signal processing unit 20A may be realized using the signal processing device 40 shown in FIG.
  • selector 25 the signal-to-noise power ratio R 1, R 2, ..., the integrated spectrum IS 1, IS 2 with R K, ..., IS K And at least one high SNR integrated spectrum can be selected.
  • the selector 25 adds or averages the plurality of integrated spectra to reduce noise. Can be generated.
  • an angle difference (angle deviation) caused by the delay time occurs between the transmission direction of the modulated transmission light ML k and the arrival direction of the scattered light to the optical antenna 14, the integration output from the selector 25
  • An observable indicative of the state of the target can be calculated with high accuracy based on the spectrum or the noise reduction spectrum. In the present embodiment, it is possible to calculate the observation amount with high accuracy even under the situation where it is difficult to grasp the angle difference (angle deviation).
  • the laser radar device 2 of the present embodiment does not require a mechanical structure for compensating for the angle difference as disclosed in the above-mentioned patent document 1.
  • the angular deviation can be compensated with a short response time. Therefore, even if the distance to the target is short, the high SNR integrated spectrum IS can be obtained, and therefore, the observation amount can be calculated with high accuracy.
  • FIG. 13 is a block diagram showing a schematic configuration of a laser radar device 3 according to a third embodiment of the present invention.
  • the configuration of the laser radar device 3 shown in FIG. 13 is the laser radar device of the first embodiment except that the light modulator 12A of FIG. 13 is included instead of the light modulator 12 of the first embodiment. It is the same as the configuration of 1.
  • the optical modulator 12A performs frequency modulation and intensity modulation on the transmission light TL input from the optical splitter 11, thereby shifting the frequency of the transmission light TL, and transmitting the transmission light TL. It is a modulator that pulsates.
  • the optical modulator 12A generates modulated transmission light ML (a series of pulse lights) including a plurality of modulated light components ML 1 to ML K and having a predetermined pulse width T at a predetermined PRF by frequency modulation and intensity modulation.
  • the modulated transmission light ML is output to the optical circulator 13.
  • the optical modulator 12A also outputs a pulse trigger signal TS generated for pulsing the transmission light TL to the signal processing unit 20.
  • the pulse trigger signal TS is a signal that represents the timing of pulsing of the transmission light TL.
  • the optical divider 11 outputs a plurality of transmitting light components TL 1 ⁇ TL K to the optical modulator 12A, an optical modulator 12A, these transmission light component TL 1 ⁇ by performing frequency modulation and intensity modulation in TL K, to generate a modulated transmission light ML that includes a plurality of modulated light components ML 1 ⁇ ML K. Therefore, the optical antenna 14 can transmit the modulated transmission light ML having a wide beam diameter and a wide beam spread angle through the optical scanner 15. Therefore, as described below, the reception aperture D of the optical antenna 14 can be improved, and the SNR of the received signal can also be improved.
  • the SNR of the reception signal is as follows:
  • the larger the f the larger the receiving aperture D.
  • the optical receiver array 17 receives the reception bore D of the optical antenna 14 field ⁇ improvement in both the value of it can.
  • the optical receiver array 17 adopts a heterodyne detection method that can receive only single mode laser light, a wide receiving field and a wide receiving aperture can be secured, so that the SNR can be improved.
  • FIG. 14 is a block diagram showing a schematic configuration of a laser radar device 4 according to a fourth embodiment of the present invention.
  • the configuration of the laser radar device 4 shown in FIG. 14 is the laser radar device of the second embodiment except that the light modulator 12A of FIG. 14 is included instead of the light modulator 12 of the second embodiment. It is the same as the configuration of 2.
  • the optical modulator 12A performs a plurality of modulations by performing frequency modulation and intensity modulation on the transmission light TL input from the optical distributor 11.
  • the modulated transmission light ML (a series of pulse lights) including the light components ML 1 to ML K is output to the optical circulator 13. Further, the optical modulator 12A outputs the pulse trigger signal TS to the signal processing unit 20.
  • the reception aperture D of the optical antenna 14 can be improved, and the SNR of the reception signal (analog detection signal CS) can also be improved.
  • Embodiments 1 to 4 As mentioned above, although various embodiments according to the present invention have been described with reference to the drawings, these embodiments are merely examples of the present invention, and various embodiments other than these embodiments can be adopted.
  • the pulse type laser radar devices 1 to 4 use pulse waves as transmission light
  • the present invention is not limited to the pulse type.
  • a CW (Continuous Wave) method using continuous wave laser light as transmission light may be applied to the laser radar devices 1 to 4.
  • the space propagation structure is described as the optical connection means between the optical components, but the present invention is not limited to this.
  • an optical waveguide structure using an optical fiber may be employed.
  • the velocity vector v is measured as the observation amount, but the observation amount is not limited to the velocity vector v.
  • the configuration of the signal processing units 20 and 20A may be appropriately changed such that the moving speed or the distance measurement value of a hard target such as a structure or a terrain object can be calculated using the integrated spectrum IS.
  • the external space is scanned using the optical scanner 15.
  • the method of scanning the external space is not limited to the method using the optical scanner 15.
  • a method of scanning the outside space with a laser beam emitted from a moving body such as an aircraft or a car may be employed.
  • the configuration of the signal processing units 20 and 20A may be appropriately changed so that the observation amount is calculated using the relative moving velocity of the moving body with respect to the external space instead of the angle information AD.
  • the laser radar device is suitable for being used in a laser sensor system such as a three-dimensional measurement system for measuring the shape of an observation object, including an observation system for observing the atmospheric state using laser light. There is.
  • 1 to 4 laser radar device 10 light source array, 101 reference light source, 102 light distribution element, 11 light distributor, 111 light guiding optical system, 112 beam splitter, 12, 12A light modulator, 13 optical circulator, 131 beam splitter, 132 wavelength plate, 14 optical antennas, 15 optical scanners, 151 rotators, 152 wedge prisms, 153 rotary driving units, 16 optical multiplexers, 161 beam splitters, 162 reflecting mirrors, 163, 164 condensing optical systems, 17 optical receivers Array, 18 switch circuit, 181 first switch, 182 second switch, 183 adder, 184 channel selection circuit, 19, 19A A / D converter, 20, 20A signal processing unit, 21, 21A signal divider, 22, 22A area conversion unit, 23, Reference Signs List 3 A integrator, 24 SNR calculator, 25 selectors (selectors), 26 speed calculators, 27 vector calculators, 28 scan speed calculators, 40 signal processors, 41 processors, 42 memories, 43 input interface units, 44 outputs Interface unit 45 signal paths Da1 to Da

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

レーザレーダ装置(1)は、複数の光出射端から複数のレーザ光を同時に出射する光源アレイ(10)と、これら複数のレーザ光から分離された送信光を変調して変調送信光を生成する光変調器(12)と、変調送信光で外部空間を走査しつつ、ターゲットで反射された当該変調送信光を受信光として受信する送受信光学系(14,15)と、複数のレーザ光から分離された複数のローカル光成分と当該受信光とを合波することにより複数の干渉光成分を生成する光合波器(16)と、これら干渉光成分を検波することで複数の検波信号を生成する光受信器アレイ(17)と、外部空間に対する走査速度に応じて、複数の検波信号の中から検波信号を選択するスイッチ回路(18)と、当該選択された検波信号に基づいてターゲットの状態を示す観測量を算出する信号処理部(20)とを備える。

Description

レーザレーダ装置
 本発明は、レーダ技術に関し、特に、レーザ光を用いて粒子状物質(たとえば、エアロゾル)などのターゲットの状態を観測するレーザレーダ技術に関するものである。
 大気中のエアロゾル(当該大気に浮遊する、液体または固体からなる微粒子)、雲粒または大気分子などの粒子状物質をターゲットとしてレーザ光を大気中に照射し、当該ターゲットで反射された散乱光を受信し、その受信散乱光に基づいて風速などの大気状態を示す情報を計測することができるレーザレーダ装置が知られている。この種のレーザレーダ装置は、たとえば特許文献1(特開2009-162678号公報)に開示されている。
 特許文献1に開示されている従来のレーザレーダ装置は、レーザビームの送信光で大気を走査するとともに当該大気で反射された散乱光を受信光として受信するビーム走査光学系と、2枚のウェッジプリズムで当該受信光を屈折させて送信光と受信光との間の角度差(光軸角度ずれ)を補償する光軸補正部と、この光軸補正部の出力光と当該送信光から分岐されたローカル光とを合波する光カプラと、この光カプラの出力光に対してヘテロダイン検波を行う信号処理部とを備えている。この従来のレーザレーダ装置では、光軸補正部は、ビーム走査に伴い送信光と受信光との間で生じる角度差を光学的に補償することができるので、信号受信部は、受信信号強度を低下させることなく風速計測を行うことができる。
国際公開第2016/117159号(たとえば、図1及び図3)
 上記した従来のレーザレーダ装置では、光軸補正部は、2枚のウェッジプリズムの回転状態を制御することによって上記角度差を補償する。しかしながら、この光軸補正部は、2枚のウェッジプリズムと、これらウェッジプリズムに受信光を導光する光学部品と、それらウェッジプリズムを機械的に回転させる機構とで構成されるため、応答性が低いという課題がある。たとえば、ターゲットとの距離が100mの場合、当該ターゲットに向けて送信光が出射された時点から当該ターゲットで反射された散乱光が受信される時点までの光伝搬時間は、約0.66μ秒となるので、当該角度差の補償に必要な応答時間は、受信視野と走査速度とに依存するものの、数十n秒(ナノ秒)のオーダである。上記した従来のレーザレーダ装置では、そのような短い応答時間で当該角度差を正確に補償することが難しい。
 また、上記した従来のレーザレーダ装置では、2枚のウェッジプリズム及び光学部品の各々の加工精度あるいは組み込み精度が低いと、送信光と受信光との間の角度差を正確に補償することができないおそれがある。この場合、信号受信部における受信信号強度が低下し、信号対雑音電力比(Signal-to-Noise Ratio,SNR)が劣化する。
 上記に鑑みて本発明の目的は、送信光と受信光との間に角度差が生じた場合であっても、高い応答性で当該角度差を補償することができ、高SNRを得ることもできるレーザレーダ装置を提供することである。
 本発明の一態様によるレーダ装置は、複数の光出射端を有し、前記複数の光出射端からそれぞれ複数のレーザ光を同時に出射する光源アレイと、前記複数のレーザ光から当該複数のレーザ光の一部を送信光として分離し、かつ前記複数のレーザ光からそれぞれ複数のローカル光成分を分離する光分配器と、前記送信光を変調して変調送信光を生成する光変調器と、前記変調送信光で外部空間を走査しつつ、当該外部空間に存在するターゲットで反射された当該変調送信光を受信光として受信する送受信光学系と、前記受信光と前記複数のローカル光成分とを合波することにより複数の干渉光成分を生成する光合波器と、前記送受信光学系の互いに異なる複数の受信視野と光学的に対応する位置にそれぞれ配置され、前記複数の干渉光成分を検波することで複数の検波信号を生成する複数の光受信器を含む光受信器アレイと、前記外部空間に対する前記送受信光学系の走査速度に応じて、前記複数の検波信号の中から検波信号を選択するスイッチ回路と、当該選択された検波信号に基づいて前記ターゲットの状態を示す観測量を算出する信号処理部とを備えることを特徴とする。
 本発明によれば、外部空間に対する送受信光学系の走査速度に応じて、複数の受信視野に対応する複数の検波信号の中から、信号強度の高い高SNRの検波信号を選択することができる。よって、変調送信光と受信光との間に角度差が生じた場合であっても、当該選択された検波信号に基づいてターゲットの状態を示す観測量を高い精度で算出することができる。また、光学的に当該角度差を補償する機構的な構造を必要とせずに当該位相差が補償されることから、高い応答性を実現することができる。
本発明に係る実施の形態1のレーザレーダ装置の概略構成を示すブロック図である。 実施の形態1のレーザレーダ装置の主要部の具体的な構成例を示す図である。 実施の形態1のレーザレーダ装置の光スキャナの構成例を概略的に示す図である。 実施の形態1における光アンテナの複数の受信視野の例を示す図である。 実施の形態1のレーザレーダ装置における信号処理部の構成を概略的に示す図である。 実施の形態1における信号処理部のハードウェア構成例を示すブロック図である。 変調送信光及び受信光の信号強度分布の例を概略的に示すグラフである。 実施の形態1に係る受信チャンネルの選択方法の例を説明するための図である。 実施の形態1に係る積分処理を説明するための図である。 本発明に係る実施の形態2のレーザレーダ装置の概略構成を示すブロック図である。 実施の形態2のレーザレーダ装置における信号処理部の構成を概略的に示すブロック図である。 実施の形態2に係る受信チャンネルの選択方法の例を説明するための図である。 本発明に係る実施の形態3のレーザレーダ装置の概略構成を示すブロック図である。 本発明に係る実施の形態4のレーザレーダ装置の概略構成を示すブロック図である。
 以下、図面を参照しつつ、本発明に係る種々の実施の形態について詳細に説明する。なお、図面全体において同一符号を付された構成要素は、同一構成及び同一機能を有するものとする。
実施の形態1.
 図1は、本発明に係る実施の形態1であるレーザレーダ装置1の概略構成を示すブロック図である。図2は、レーザレーダ装置1の主要部の具体的な構成例を示す図である。図1に示されるようにレーザレーダ装置1は、光源アレイ10、光分配器11、光変調器12、光サーキュレータ13、光アンテナ14、光スキャナ15、光合波器16、光受信器アレイ17、スイッチ回路18、A/D変換器(Analog-to-Digital converter)19及び信号処理部20を備えている。光サーキュレータ13、光アンテナ14及び光スキャナ15の組合せによって本実施の形態の送受信光学系が構成される。
 光源アレイ10は、空間的(2次元的または3次元的)に配列されたK個の光出射端(Kは3以上の整数)を有し、これらK個の光出射端からK本のレーザ光L~Lを光分配器11に向けて出射する。後述するように光アンテナ14は、互いに異なるK個の受信視野を有し、これらK個の受信視野を通じてそれぞれK個の受信チャンネルCH,CH,…,CHの散乱光を受信することができる。レーザ光L~Lの本数は、受信チャンネルCH~CHの数と同じである。なお、本実施の形態では、受信チャンネル数は3以上であるが、これに限定されずに受信チャンネル数が2であってもよい。
 図2の例では、光源アレイ10は、基準レーザ光を出力する基準光源101と、この基準光源101から光導波路を介して入力された1本の基準レーザ光を3本のレーザ光L~Lに分配する光分配素子102と、この光分配素子102から入力された3本のレーザ光L~Lを光分配器11に向けてそれぞれ出力する光出射端E1~E3とを有している。基準光源101としては、半導体レーザまたは固体レーザが使用されればよい。基準レーザ光の偏光状態は、1つの方向に偏光する線偏光状態となるように保持されている。光分配素子102としては、たとえば、溶融型の光ファイバ分岐カプラまたはフィルタ型分岐カプラ(誘電体多層膜フィルタを使用して光を分岐する光学部品)が使用されればよい。なお、光源アレイ10として複数の発光面をもつアレイレーザが採用されてもよい。
 図1を参照すると、光分配器11は、入力されたレーザ光L~Lからそれぞれ送信光成分TL~TLを分離し、かつレーザ光L~Lからそれぞれローカル光成分LL~LLを分離する。光分配器11は、送信光成分TL~TLを含む送信光TLを光変調器12に出力し、ローカル光成分LL~LLを含むローカル光LLを光合波器16に出力する。光分配器11としては、たとえば、誘電体多層膜フィルタを有する分岐ミラーあるいはビームスプリッタが使用されればよい。
 図2の例では、光分配器11は、導光光学系111とビームスプリッタ112とで構成されている。導光光学系111は、光源アレイ10から入力されたレーザ光L~Lをビームスプリッタ112の光入射面に導光する光学部品である。ビームスプリッタ112は、導光光学系111から入力されたレーザ光L~Lの一部を送信光成分TL~TLとして透過させ、レーザ光L~Lから所定の分岐比で分離されたローカル光成分LL~LLを光合波器16の方向に反射させる。所定の分岐比は、システム設計により決定されればよい。
 図1を参照すると、光変調器12は、光分配器11から入力された送信光成分TL~TLのうちの一部の送信光成分TLにのみ周波数変調及び強度変調を実行することにより、送信光成分TLの周波数をシフトさせ、かつ送信光成分TLをパルス化する変調器である。光変調器12は、周波数変調及び強度変調により、所定のパルス繰り返し周波数(Pulse Repetitive Frequency,PRF)で所定のパルス幅Tを有する変調送信光ML(一連のパルス光)を生成し、変調送信光MLを光サーキュレータ13に出力する。また、光変調器12は、送信光成分TLのパルス化のために生成されたパルストリガ信号TSを信号処理部20に出力する。パルストリガ信号TSは、送信光成分TLのパルス化のタイミングを表す信号である。
 ここで、パルス幅Tは、距離分解能値に相当する。所望の距離分解能値に相当するパルス幅Tをあらかじめ設定しておくことが可能である。また、光変調器12においては、設計時に設定された固定のパルス幅T及びPRFが使用されてもよい。更に、光変調器12の出力光強度が不足している場合には、光変調器12の後段に、光変調器12の出力光を増幅する光増幅器が設けられてもよい。光変調器12は、周波数シフト機能を有する位相変調器を使用して送信光の周波数をシフトさせてもよい。
 このような光変調器12としては、たとえば、AO周波数シフタ(Acousto-Optical frequency shifter)または公知の光位相変調器が使用されればよい。パルストリガ信号TSとしては、たとえば、5ボルトのTTLレベルを有するTTL(Transistor-Transistor Logic)信号が使用可能である。
 光サーキュレータ13は、3ポートを有する非相反型の光デバイスである。すなわち、光サーキュレータ13は、光変調器12から光が入射される光入力ポートと、光アンテナ14から光が入射され、かつ光アンテナ14に光を出射する光入出力ポートと、光合波器16に光を出射する光出力ポートとを有する。光サーキュレータ13は、光変調器12から光入力ポートに入力された変調送信光MLと、光アンテナ14から光入出力ポートに入力された光とを光学的に互いに分離する。そして、光サーキュレータ13は、光変調器12から入力された変調送信光MLを光入出力ポートから光アンテナ14に出力し、同時に、光アンテナ14からの入力光を光出力ポートから光合波器16に出力する。
 光アンテナ14は、光サーキュレータ13から入力された変調送信光MLを平行光に変換して光スキャナ15に出力する。光スキャナ15は、光アンテナ14から入力された変調送信光MLの外部空間への照射方向(すなわち視線方向)を変化させる機能を有している。また、光スキャナ15は、光アンテナ14から入力された変調送信光MLで外部空間の所定範囲を繰り返し走査することができる。このような光スキャナ15は、たとえば、ウェッジプリズム、ガルバノミラーもしくはポリゴンミラー、またはこれら光学部品の組合せで構成されればよい。
 図3は、実施の形態1の光スキャナ15の構成例を概略的に示す図である。図3に示されるように、光スキャナ15は、ウェッジプリズム152を含む回転子151と、この回転子151を光軸OAの回りに任意速度で回転駆動させる回転駆動部153とを有する。回転駆動部153は、ステッピングモータなどの回転モータ及びロータリエンコーダ(図示せず。)を含んで構成されている。ウェッジプリズム152の光アンテナ14側の光学面は、光軸OAに対して垂直であり、ウェッジプリズム152の外部空間側の光学面は、光軸OAに対して傾斜している。このため、回転駆動部153は、ウェッジプリズム152を光軸OAの回りに回転させることで、レーザビームで外部空間を走査させることができる。回転駆動部153は、内蔵するロータリエンコーダで検出された、光スキャナ15の走査角度の現在値を示す角度情報ADを信号処理部20に出力する。
 なお、光スキャナ15の走査速度は、設計時にあらかじめ設定された値でよい。また、本実施の形態では、回転子151の回転を制御する制御信号は、回転駆動部153で生成されるが、これに限定されるものではない。この代わりに、回転子151の回転を制御する制御信号を生成するように信号処理部20の構成が変更されてもよい。
 光アンテナ14は、外部空間におけるエアロゾル(大気中に浮遊する、液体または固体からなる微粒子)、雲粒もしくは大気分子などのターゲットで反射されて戻ってきた散乱光、または、構造物もしくは地形物などのハードターゲットで反射されて戻ってきた散乱光(拡散光)を、光スキャナ15を介して受信する。図1に示される光サーキュレータ13は、光アンテナ14から入力された散乱光すなわち受信光RLを光合波器16に出力する。このような光アンテナ14は、たとえば、光学望遠鏡またはカメラレンズを用いて構成可能である。なお、光アンテナ14は、集光調整機能を有していてもよい。
 図2の例では、光サーキュレータ13は、偏光ビームスプリッタ131と波長板132とを有する空間伝搬型の光デバイスとして構成されている。この例では、波長板132は、光変調器12から偏光ビームスプリッタ131を介して入力された直線偏光を円偏光に変換し、光アンテナ14から入力された円偏光を直線偏光(水平偏光)に変換する1/4波長板である。偏光ビームスプリッタ131は、光変調器12から入力された直線偏光を通過させ、光アンテナ14から入力された円偏光を光合波器16の方向に反射させる。
 図1に示される光アンテナ14は、ターゲットから到来する散乱光に対して互いに異なるK個の受信視野を有している。光サーキュレータ13は、これら受信視野を通じて入力された、K個の受信チャンネルCH,CH,…,CHの受信光成分RL,RL,…,RLを含む受信光RLを光合波器16に出力する。図4は、光アンテナ14の3つの受信視野(第1~第3受信視野)の例を示す図である。図4の例では、1つの受信チャンネル当たりの受信視野角θは、送信レーザビーム(変調送信光MLのレーザビーム)の拡がり角θと同一となるように設定されている。また、これら受信視野は、θの角度オフセットで互いに異なる角度方向を観測できるように設定されている。ただし、この例の代わりに、隣接する受信視野が互いにオーバラップするように角度オフセットの値(たとえば、θ/2)が設定されてもよい。
 次に、図1を参照すると、光合波器16は、光分配器11から入力されたローカル光LLと光サーキュレータ13から入力された受信光RLとを合波して、K個の干渉光成分CL,…,CLを含む干渉光CLを生成し、この干渉光CLすなわち光ビート信号を光受信器アレイ17に出力する。光合波器16は、たとえば、誘電体多層膜フィルタを有する分岐ミラーまたはビームスプリッタを用いて構成可能である。
 光受信器アレイ17は、光アンテナ14のK個の受信視野と光学的に対応する位置にそれぞれ配置されたK個の光受信器を有する。これらK個の光受信器は、K個の干渉光成分CL,CL,…,CLをそれぞれヘテロダイン検波(周波数復調)することで、K個のアナログ検波信号CS,CS,…,CSを含むアナログ検波信号群CSを生成することができる。k番目の干渉光成分CLは、k番目のローカル光成分LLとk番目の受信光成分RLとの合波により生じた干渉光成分である。
 なお、バランスド検出ができるように光合波器16及び光受信器アレイ17が構成されてもよい。バランスド検出により、ヘテロダイン検波後の位相雑音の低減が可能となる。
 スイッチ回路18は、光スキャナ15の走査速度に応じて、アナログ検波信号CS,CS,…,CSの中からアナログ検波信号SSを選択し、選択されたアナログ検波信号SSをA/D変換器19に出力することができる。
 後述するように信号処理部20は、角度情報ADに基づいて光スキャナ15の走査速度を示す速度情報SDを逐次算出し、当該速度情報SDをスイッチ回路18に供給している。スイッチ回路18は、当該走査速度に基づいて、レーザビームの送信時から散乱光の受信時までの遅延時間τに対応する角度ずれθを算出する。ここで、角度ずれθとは、レーザビーム(変調送信光MLのレーザビーム)の送信方向と、光アンテナ14への散乱光の到来方向との間の角度差である。
 そして、スイッチ回路18は、角度ずれθの値に基づき、当該散乱光の到来方向に対応する受信チャンネルCHを選択し、アナログ検波信号CS,CS,…,CSの中から受信チャンネルCHに対応するアナログ検波信号SSを選択する。これにより、スイッチ回路18は、当該散乱光の到来方向を示す受信視野と光学的に対応する光受信器の出力信号(アナログ検波信号)SSを選択することができる。
 仮に、アナログ検波信号CS,CS,…,CSがすべて加算されるとすれば、当該散乱光が到来しない受信チャンネル(受信視野)のアナログ検波信号も加算されるため、加算信号のSNRは劣化してしまう。これに対し、本実施の形態のスイッチ回路18は、当該散乱光が到来する受信チャンネル(受信視野)のアナログ検波信号SSを選択するので、角度ずれθに起因するSNRの劣化を回避することができる。
 光スキャナ15によるビーム走査が所定の走査速度Vscan(単位:ラジアン/秒)で或る方向に行われた場合、レーザビームの送信時から散乱光の受信時までの遅延時間τ(単位:秒)に対応する角度ずれθは、τ×Vscan(単位:ラジアン)で与えられる。よって、スイッチ回路18は、角度ずれθに基づき、遅延時間τに対して、どの受信チャンネルを選択するのかを事前に決定することができる。たとえば、スイッチ回路18は、受信チャンネルを切り替えるタイミングを、tch=(θ/θ)の小数点切り上げ時の整数値で決めることができる。ここで、θは、1つの受信チャンネル当たりの受信視野角である。このとき、選択すべき受信チャンネルは1つのみでもよいが、これに限定されない。オーバーラップさせるために2つ以上の受信チャンネルが選択されてもよい。この場合、スイッチ回路18は、選択された2つ以上の受信チャンネルのアナログ検波信号を加算し、その加算信号をA/D変換器19に出力すればよい。これにより、広い受信視野を確保しつつ高SNRの信号を生成することができる。
 図2の例では、光合波器16は、ビームスプリッタ161、反射ミラー162及び集光光学系163,164を含んで構成されている。また、図1の光受信器アレイ17は、光受信器Da1~Da3と光受信器Db1~Db3とで構成されている。光受信器Da1,Db1の出力信号の和がアナログ検波信号CLに相当し、光受信器Da2,Db2の出力信号の和がアナログ検波信号CLに相当し、光受信器Da3,Db3の出力信号の和がアナログ検波信号CLに相当する。更に、図1のスイッチ回路18は、第1スイッチ181、第2スイッチ182、加算器183及びチャンネル選択回路184で構成されている。
 図2に示されるように、反射ミラー162は、偏光ビームスプリッタ131から入力された受信光RLをビームスプリッタ161の方向に反射させる。ビームスプリッタ161は、偏光ビームスプリッタ131から反射ミラー162を経て入射された受信光RLの一部を一方の集光光学系163の方向に透過させ、当該受信光RLの残部を他方の集光光学系164の方向に反射させる。同時に、ビームスプリッタ161は、光分配器11から入射されたローカル光LLの一部を一方の集光光学系163の方向に反射させ、当該ローカル光LLの残部を他方の集光光学系164の方向に透過させる。結果として、ビームスプリッタ161は、受信光RLの一部とローカル光LLの一部との合波により生成された干渉光を一方の集光光学系163に出力し、同時に、受信光RLの残部とローカル光LLの残部との合波により生成された干渉光を他方の集光光学系164に出力する。
 図2に示されるように、集光光学系163は、ビームスプリッタ161から入射された干渉光を光受信器Da1~Da3に集光させ、集光光学系164は、ビームスプリッタ161から入射された干渉光を光受信器Db1~Db3に集光させる。
 光受信器Da1,Da2,Da3は、図4の第1受信視野、第2受信視野及び第3受信視野と光学的に対応する位置にそれぞれ配置されている。このため、光受信器Da1,Da2,Da3は、第1受信視野、第2受信視野及び第3受信視野からそれぞれ伝搬した3つの散乱光成分を含む3つの干渉光成分をそれぞれ検出して3本のアナログ検波信号を出力することができる。一方、光受信器Db1,Db2,Db3も、図4の第1受信視野、第2受信視野及び第3受信視野と光学的に対応する位置にそれぞれ配置されている。このため、光受信器Db1,Db2,Db3は、第1受信視野、第2受信視野及び第3受信視野からそれぞれ伝搬した3つの散乱光成分を含む3つの干渉光成分をそれぞれ検出して3本のアナログ検波信号を出力することができる。
 チャンネル選択回路184は、信号処理部20から逐次供給される速度情報SDに基づいて、レーザビームの送信時から散乱光の受信時までの遅延時間τに対応する角度ずれθを算出し、当該算出された角度ずれθに基づき、当該散乱光の到来方向に対応する受信チャンネルCHを選択する。次に、チャンネル選択回路184は、受信チャンネルCHのアナログ検波信号CSを表す出力を選択させる選択制御信号CSを第1スイッチ181及び第2スイッチ182に供給する。第1スイッチ181及び第2スイッチ182は、光受信器Da1~Da3,Db1~Db3の出力の中から、選択制御信号CSで指定された出力Ss1,Ss2を選択し、当該選択された出力Ss1,Ss2を加算器183に供給する。加算器183は、供給された出力Ss1,Ss2を加算してアナログ検波信号SS(=CS)を生成し、アナログ検波信号SSをA/D変換器19に出力する。
 A/D変換器19は、サンプリング周波数fsでアナログ検波信号SSをサンプリングすることでディジタル検波信号DSを生成し、このディジタル検波信号DSを信号処理部20に出力する。A/D変換器19としては、たとえば、二重積分型A/D変換器、逐次比較形A/D変換器または並列比較型A/D変換器が使用可能である。なお、A/D変換器19は、パルストリガ信号TSをトリガとし、パルストリガ信号TSと同期してアナログ検波信号SSをサンプリングすることにより、ディジタル検波信号DSを生成してもよい。
 図1に示される信号処理部20は、ディジタル検波信号DSに基づいてターゲットの状態を示す観測量を算出する機能を有する。後述するように信号処理部20は、大気中に浮遊する粒子状物質(たとえば、エアロゾル、雲粒または大気分子)などのターゲットの状態を示す観測量として当該ターゲットの視線速度及び速度ベクトルを算出するように構成可能であるが、算出される観測量は、視線速度及び速度ベクトルに限定されるものではない。たとえば、ディジタル検波信号DSに基づいてレーザレーダ装置1とターゲットとの間の光伝搬時間を測定し、光伝搬時間の測定値に基づいて当該ターゲットとの距離(測距値)を算出できるように信号処理部20が構成されてもよい。光スキャナ15が変調送信光MLでターゲットの表面を走査し、信号処理部20が当該ターゲットとの距離分布を計測するようにレーザレーダ装置1を構成してもよい。これにより、移動体(たとえば、車両、人体もしくは動物)、地形物または構造物といったターゲットの3次元形状の計測が可能となる。
 図5は、本実施の形態における信号処理部20の構成例を概略的に示すブロック図である。図5に示される信号処理部20は、ディジタル検波信号DSの信号スペクトル(振幅スペクトルまたは電力スペクトル)を算出し、当該信号スペクトルに基づき、ターゲットの移動状態を示す視線速度及び速度ベクトルを算出することができる。図5に示されるように信号処理部20は、ディジタル検波信号DSを、連続するM個のゲート時間領域(Mは3以上の正整数)の信号波形をそれぞれ表すM個のレンジビン信号(時間領域信号)B,B,…,Bに分割する信号分割器21と、これらレンジビン信号B,B,…,Bを周波数領域信号F,F,…,Fにそれぞれ変換する領域変換部22と、周波数領域信号F,F,…,Fに基づき、ディジタル検波信号DSの信号スペクトルとして積算スペクトルISを算出する積算器23と、積算スペクトルISに基づいてターゲットの視線速度Vrを算出する速度算出器26と、ターゲットの速度ベクトルvを算出するベクトル算出器27と、光スキャナ15から逐次供給された角度情報ADに基づいて光スキャナ15の走査速度を示す速度情報SDを算出する走査速度算出器28とを含む。速度算出器26及びベクトル算出器27によって本実施の形態の観測量算出器が構成される。
 このような信号処理部20のハードウェア構成は、たとえば、DSP(Digital Signal Processor),ASIC(Application  Specific  Integrated  Circuit)またはFPGA(Field-Programmable Gate Array)などの半導体集積回路を有するプロセッサで実現されればよい。あるいは、信号処理部20のハードウェア構成は、メモリから読み出されたソフトウェアまたはファームウェアのプログラムコードを実行する、CPU(Central Processing Unit)またはGPU(Graphics Processing Unit)などの演算装置を含むプロセッサで実現されてもよい。
 図6は、信号処理部20の機能を実現するハードウェア構成例である信号処理装置40を概略的に示すブロック図である。信号処理装置40は、プロセッサ41、メモリ42、入力インタフェース部43、出力インタフェース部44及び信号路45を含んで構成されている。信号路45は、プロセッサ41、メモリ42、入力インタフェース部43及び出力インタフェース部44を相互に接続するためのバスである。入力インタフェース部43は、外部から入力されたディジタル検波信号DSと角度情報AD(図示せず。)とを信号路45を介してプロセッサ41に転送する機能を有する。プロセッサ41は、転送されたディジタル検波信号DSにディジタル信号処理を施してターゲットの速度ベクトルvを算出し、転送された角度情報ADに基づいて速度情報SDを算出する。プロセッサ41は、速度ベクトルvを示すデータを、信号路45及び出力インタフェース部44を介して外部機器(たとえば、表示装置)に出力することができる。また、プロセッサ41は、速度情報SDを信号路45及び出力インタフェース部44を介してスイッチ回路18に出力することもできる。
 メモリ42は、プロセッサ41がディジタル信号処理を実行する際に使用されるデータ記憶領域である。プロセッサ41がCPUなどの演算装置を内蔵する場合には、メモリ42は、プロセッサ41により実行されるソフトウェアまたはファームウェアのプログラムコードを記憶するデータ記憶領域を有していればよい。メモリ42としては、たとえば、ROM(Read Only Memory)及びSDRAM(Synchronous Dynamic Random Access Memory)などの半導体メモリを使用することが可能である。
 以下、図5に示される信号処理部20の構成及び動作を詳細に説明する。
 走査速度算出器28は、光スキャナ15から逐次供給された角度情報ADに基づいて光スキャナ15の走査速度を示す速度情報SDを逐次算出し、この速度情報SDをスイッチ回路18に供給する。
 信号分割器21は、スイッチ回路18で選択されたディジタル検波信号DSを入力とし、ディジタル検波信号DSを、連続するM個のゲート時間領域の信号波形をそれぞれ表すM個のレンジビン信号(時間領域信号)B,B,…,Bに分割する。図7は、レンジビンの概念を説明するためのグラフである。このグラフにおいて、横軸は時間を、縦軸は信号強度を示している。図7に示されるように、変調送信光MLを示す送信パルスTwが1/PRFの周期で繰り返し送信される。送信パルスTwの送信(ショット)ごとに、ターゲットで反射された散乱光の受信信号Rwが観測される。時刻tは、送信パルスTwの送信(ショット)ごとの計測開始時刻である。受信信号Rwは、M個のゲート時間領域上のレンジビンRb,Rb,…,Rbに分割される。レンジビン信号B,B,…,Bは、これらレンジビンRb,Rb,…,Rbにそれぞれ対応するものである。
 たとえば、30m程度に設定された距離分解能Rresについて、各レンジビンを定めるゲート時間(時間幅)を2Rres/c(c:光速)と設定することができる。また、観測開始距離を調整するための値として、A/D変換器19がパルストリガ信号TSに応じてA/D変換を開始するまでのディレイタイミングに相当する距離値がRminで表されるものとする。このとき、パルストリガ信号TSに対するA/D変換開始のための遅延量が0の場合であって、たとえば40mから測定開始する場合は、距離値Rminとして40の値を設定することができ、各ディジタル検波信号の時間波形に対する計測開始時間を、tstart=2Rmin/cと設定することができる。
 ところで、スイッチ回路18では、図8に例示されるように、受信チャンネルCH,CH,…,CHのアナログ検波信号CS,CS,…,CSが入力される。レーザビームの送信時から散乱光の受信時までの遅延時間τ=τまたはτの場合には、スイッチ回路18は、遅延時間τまたはτに相当する角度ずれを補償するために受信チャンネルCHを選択し、この受信チャンネルCHのアナログ検波信号CSをアナログ検波信号SSとして選択する。一方、遅延時間τ=τまたはτの場合には、スイッチ回路18は、遅延時間τまたはτに相当する角度ずれを補償するために受信チャンネルCHを選択し、受信チャンネルCHのアナログ検波信号CSをアナログ検波信号SSとして選択する。このようにスイッチ回路18で選択されたアナログ検波信号SSが入力されると、信号分割器21は、アナログ検波信号SSをM個のレンジビン信号B,B,…,Bに分割する。
 図5に示される領域変換部22は、各ショットについて得られたレンジビン信号B,B,…,Bをそれぞれ周波数領域信号F,F,…,Fに変換する。具体的には、領域変換部22は、レンジビン信号B,B,…,Bの各々に対して離散フーリエ変換を実行することにより、周波数領域信号F,F,…,Fを算出することができる。離散フーリエ変換としては、FFTビン数NFFT(たとえば、256点)での高速フーリエ変換(Fast Fourier Transform,FFT)が使用可能である。
 積算器23は、各ショットについて、周波数領域信号F,F,…,Fに基づいてレンジビン信号B,B,…,BそれぞれのスペクトルS(m,n)を算出する。ここで、mはレンジビン番号を示す正整数(m=1~M)、nはショット番号を示す正整数(n=1~N)である。Nは、あらかじめユーザによって指定された積分回数である。スペクトルS(m,n)は、周波数と信号振幅との関係を表す振幅スペクトル、または周波数と電力との関係を表す電力スペクトルであればよい。
 また、積算器23は、各レンジビン(レンジビン番号m)について、スペクトルS(m,1)~S(m,N)をショット番号の方向に積分して(インコヒーレント積分して)高SNRの積分信号スペクトルIS(m)を算出する。図9は、積算器23における積分処理の概要を示す図である。たとえば、レンジビンRbについては、積算器23は、スペクトルS(1,1)~S(1,N)を積分することで高SNRの積分信号スペクトルIS(1)を算出することができる。この積分結果を積分回数Nで平均化した場合、SNRは積分回数Nに対してN1/2倍で改善する。なお、積算器23では、インコヒーレント積分が採用されているが、当該インコヒーレント積分に代えてコヒーレント積分が採用されてもよい。
 更に、積算器23は、積分信号スペクトルIS(1)~IS(M)を積算することで積算スペクトルISを算出し、この積算スペクトルISを速度算出器26に出力する。
 速度算出器26は、積算スペクトルISのピーク位置または重心位置の周波数をドップラ周波数fとして検出する。そして、速度算出器26は、ドップラ周波数fを用いて、角度情報ADに対応する視線方向におけるターゲットの速度成分すなわち視線速度Vr(単位:m/s)を算出する。視線速度Vrは、次式(1)に従って算出可能である。ここで、λは、レーザビームの波長である。
 Vr=λ×f/2                  (1)
 ベクトル算出器27は、視線速度Vrに基づき、公知のベクトル合成法もしくはVAD(Velocity Azimuth Display)法を用いてターゲットの速度ベクトルを算出することができる。ベクトル合成法を採用する場合、たとえば、視線速度Vrは、東西方向の水平方向風速(U)、南北方向の水平方向風速(V)、鉛直方向の風速(W)、仰角(θ)、及び、北を基準とした方位角(φ)を用いて、以下の式(2)で表される。
 Vr=Usinφsinθ+Vcosφsinθ+Wcosθ  (2)
 ベクトル算出器27は、たとえば、少なくとも3方向の視線方向の視線速度に基づいて連立方程式を解くことで、ターゲットの速度ベクトルv=(U,V,W)を風速ベクトルとして算出することができる。
 以上に説明した実施の形態1のレーザレーダ装置1の効果は以下のとおりである。光受信器アレイ17は、光アンテナ14のK個の受信視野と光学的に対応する位置にそれぞれ配置されたK個の光受信器を有し、これら光受信器は、K個の干渉光成分をそれぞれ検波することでアナログ検波信号CS~CSを生成する。スイッチ回路18は、光スキャナ15の走査速度に応じて、K個の受信視野にそれぞれ対応するアナログ検波信号CS~CSの中から、信号強度の高い高SNRのアナログ検波信号SSを選択することができる。よって、変調送信光MLの送信方向と光アンテナ14への散乱光の到来方向との間に遅延時間に起因する角度差(角度ずれ)が生じたとしても、信号処理部20は、高SNRのアナログ検波信号SSから得られたディジタル検波信号DSに基づいてターゲットの状態を示す観測量を高い精度で算出することができる。
 また、本実施の形態のレーザレーダ装置1は、上記特許文献1に開示されているような当該角度差を補償する機構的な構造を必要とせずに、短い応答時間で角度ずれを補償することができる。よって、ターゲットとの距離が近くても、高SNRのディジタル検波信号DSが得られるため、高精度で観測量を算出することが可能である。
実施の形態2.
 図10は、本発明に係る実施の形態2であるレーザレーダ装置2の概略構成を示すブロック図である。図10に示されるようにレーザレーダ装置2は、上記実施の形態1と同様に、光源アレイ10、光分配器11、光変調器12、光サーキュレータ13、光アンテナ14、光スキャナ15、光合波器16及び光受信器アレイ17を備えている。本実施の形態のレーザレーダ装置2の構成は、上記実施の形態1のスイッチ回路18、A/D変換器19及び信号処理部20に代えて、図10のA/D変換器19A及び信号処理部20Aを有する点を除き、実施の形態1のレーザレーダ装置1の構成と同じである。
 A/D変換器19Aは、所定のサンプリング周波数で、アナログ検波信号CS,CS,…,CSの各々をサンプリングすることで、アナログ検波信号CS,CS,…,CS(受信チャンネルCH,CH,…,CH)にそれぞれ対応するディジタル検波信号DS,DS,…,DSを生成し、これらディジタル検波信号DS,DS,…,DSを信号処理部20Aに出力する。A/D変換器19Aとしては、たとえば、二重積分型A/D変換器、逐次比較形A/D変換器または並列比較型A/D変換器が使用可能である。なお、A/D変換器19Aは、パルストリガ信号TSをトリガとし、パルストリガ信号TSと同期してアナログ検波信号群CSをサンプリングすることにより、ディジタル検波信号DS,DS,…,DSを生成してもよい。
 信号処理部20Aは、ディジタル検波信号DS,DS,…,DSそれぞれの信号スペクトル(積算スペクトル)IS,IS,…,ISを算出し、これら信号スペクトルIS,IS,…,ISに基づいてディジタル検波信号DS,DS,…,DSそれぞれの信号対雑音電力比R,R,…,Rを算出する。そして、信号処理部20Aは、信号対雑音電力比R,R,…,Rを用いて、これら信号スペクトルIS,IS,…,ISの中から少なくとも1つの信号スペクトルを選択し、当該選択された少なくとも1つの信号スペクトルに基づいてターゲットの状態を示す観測量を算出する機能を有している。
 図11は、実施の形態2における信号処理部20Aの構成例を概略的に示すブロック図である。図11に示されるように信号処理部20Aは、信号分割器21A、領域変換部22A、積算器23A、SNR算出器24、選択器25、速度算出器26及びベクトル算出器27を備えて構成されている。信号処理部20Aの構成は、実施の形態1の信号分割器21、領域変換部22及び積算器23に代えて、図11の信号分割器21A、領域変換部22A、積算器23A、SNR算出器24及び選択器25を有する点を除き、実施の形態1の信号処理部20の構成と同じである。
 信号分割器21Aは、各ショットについて得られたディジタル検波信号DS,DS,…,DSのうちの各ディジタル検波信号DSを、M個のゲート時間領域の信号波形をそれぞれ表すM個のレンジビン信号(時間領域信号)Bk,1,Bk,2,…,Bk,Mに分割する。信号分割器21Aは、受信チャンネルCHごとに、M個のレンジビン信号Bk,1,Bk,2,…,Bk,Mを含むレンジビン信号群RBを領域変換部22Aに出力する。
 領域変換部22Aは、各ショットについて、レンジビン信号群RF,RF,…,RFのうちの各レンジビン信号群RFを、M個の周波数領域信号Fk,1,…,Fk,Mを含む周波数領域信号群RFに変換し、この周波数領域信号群RFを積算器23Aに出力する。すなわち、領域変換部22Aは、レンジビン信号Bk,1,Bk,2,…,Bk,Mの各々に対して離散フーリエ変換を実行することにより、周波数領域信号Fk,1,Fk,2,…,Fk,Mを算出することができる。離散フーリエ変換としては、FFTビン数NFFT(たとえば、256点)での高速フーリエ変換(Fast Fourier Transform,FFT)が使用可能である。
 積算器23Aは、各ショットについて、各周波数領域信号群RFに基づいて積算スペクトルISを算出し、この積算スペクトルISをSNR算出器24及び選択器25に出力する。
 具体的には、積算器23Aは、各ショットについて、周波数領域信号Fk,1,Fk,2,…,Fk,Mに基づいてレンジビン信号Bk,1,Bk,2,…,Bk,MそれぞれのスペクトルS(m,n)を算出する。ここで、kは受信チャンネル番号(k=1~K)、mはレンジビン番号(m=1~M)、nはショット番号(n=1~N)である。Nは、あらかじめユーザによって指定された積分回数である。スペクトルS(k,m,n)は、振幅スペクトルまたは電力スペクトルであればよい。
 また、積算器23Aは、各レンジビン(レンジビン番号m)について、スペクトルS(m,1)~S(m,N)をショット番号の方向に積分して(インコヒーレント積分して)高SNRの積分信号スペクトルIS(m)を算出する。たとえば、レンジビンRbについては、積算器23Aは、スペクトルS(1,1)~S(1,N)を積分することで高SNRの積分信号スペクトルIS(1)を算出することができる。この積分結果を積分回数Nで平均化した場合、SNRは積分回数Nに対してN1/2倍で改善する。なお、積算器23では、インコヒーレント積分が採用されているが、当該インコヒーレント積分に代えてコヒーレント積分が採用されてもよい。
 更に、積算器23Aは、積分信号スペクトルIS(1)~IS(M)を積算することで積算スペクトルISを算出し、この積算スペクトルISをSNR算出器24及び選択器25に出力する。
 SNR算出器24は、積算スペクトルIS,IS,…,ISから信号対雑音電力比R,R,…,Rを算出する。信号対雑音電力比Rの算出方法には公知の方法が使用されればよく、特に限定されるものではない。SNR算出器24は、たとえば、積算スペクトルISのピーク値を検出し、当該ピーク値とあらかじめ計測された雑音レベルとの比率を信号対雑音電力比Rとして算出可能である。
 選択器25は、信号対雑音電力比R,R,…,Rを用いて積算スペクトルIS,IS,…,ISの中から少なくとも1つの高SNRの積算スペクトルを選択する。積算スペクトルIS,IS,…,ISの中から1つの積算スペクトルを選択したとき、選択器25は、当該選択された積算スペクトルをそのまま積分スペクトルISとして出力する。一方、積算スペクトルIS,IS,…,ISの中から複数の積算スペクトルを選択したとき、選択器25は、当該選択された積分スペクトルを加算または平均化して1つの積算スペクトルISを算出し、当該積算スペクトルISを速度算出器26に出力する。選択器25においては、図12に例示されるように、受信チャンネルCH,CH,…,CHの積算スペクトルIS,IS,…,ISが入力される。選択器25は、遅延時間に相当する角度ずれを補償するために、信号対雑音電力比R,R,…,Rを用いて、閾値THを超える周波数成分を含む積算スペクトルを与える受信チャンネルを選択する。図12の例では、選択器25は、遅延時間τまたはτに相当する角度ずれを補償するために受信チャンネルCHを選択し、この受信チャンネルCHの積算スペクトルISを積算スペクトルISとして出力することができる。また、選択器25は、遅延時間τに相当する角度ずれを補償するために受信チャンネルCH,CHを選択し、これら受信チャンネルCH,CHの積算スペクトルIS,ISを加算または平均化してノイズ低減スペクトルを算出し、このノイズ低減スペクトルを積算スペクトルISとして出力することができる。更に、選択器25は、遅延時間τに相当する角度ずれを補償するために受信チャンネルCHを選択し、この受信チャンネルCHの積算スペクトルISを積算スペクトルISとして選択することができる。
 前述のとおり、遅延時間τの場合には、選択器25は、選択された複数の積算スペクトルIS,ISを加算または平均化することでノイズ低減スペクトルを算出し、このノイズ低減スペクトルを積算スペクトルISとして速度算出器26に出力する。これにより、SNRを改善することができる。
 速度算出器26は、積算スペクトルISのピーク位置または重心位置の周波数をドップラ周波数fとして検出し、このドップラ周波数fを用いて、角度情報ADに対応する視線速度Vrを算出することができる。また、ベクトル算出器27は、視線速度Vrに基づき、公知のベクトル合成法もしくはVAD法を用いてターゲットの速度ベクトルvを算出することができる。
 なお、本実施の形態では、ターゲットの状態を示す観測量として当該ターゲットの視線速度Vr及び速度ベクトルvが算出されているが、算出される観測量は、視線速度Vr及び速度ベクトルvに限定されるものではない。たとえば、積算スペクトルISに基づいてレーザレーダ装置2とターゲットとの間の距離(測距値)を算出できる機能ブロックを有するように信号処理部20Aが構成されてもよい。
 上記した信号処理部20Aのハードウェア構成は、たとえば、DSP,ASICまたはFPGAなどの半導体集積回路を有するプロセッサで実現されればよい。あるいは、信号処理部20Aのハードウェア構成は、メモリから読み出されたソフトウェアまたはファームウェアのプログラムコードを実行する、CPUまたはGPUなどの演算装置を含むプロセッサで実現されてもよい。実施の形態1の場合と同様に、図6に示した信号処理装置40を用いて信号処理部20Aの機能が実現されてもよい。
 以上に説明した実施の形態2のレーザレーダ装置2では、選択器25は、信号対雑音電力比R,R,…,Rを用いて積算スペクトルIS,IS,…,ISの中から少なくとも1つの高SNRの積算スペクトルを選択することができる。また、選択器25は、積算スペクトルIS,IS,…,ISの中から複数の積算スペクトルが選択された場合には、これら複数の積算スペクトルを加算または平均化することでノイズ低減スペクトルを生成することができる。よって、変調送信光MLの送信方向と光アンテナ14への散乱光の到来方向との間に遅延時間に起因する角度差(角度ずれ)が生じたとしても、選択器25から出力された積算スペクトルまたはノイズ低減スペクトルに基づいてターゲットの状態を示す観測量を高い精度で算出することができる。本実施の形態では、角度差(角度ずれ)を把握することが困難な状況下でも、観測量の高精度な算出が可能である。
 また、上記実施の形態1の場合と同様に、本実施の形態のレーザレーダ装置2は、上記特許文献1に開示されているような当該角度差を補償する機構的な構造を必要とせずに、短い応答時間で角度ずれを補償することができる。よって、ターゲットとの距離が近くても、高SNRの積算スペクトルISが得られるため、高精度で観測量を算出することができる。
実施の形態3.
 図13は、本発明に係る実施の形態3のレーザレーダ装置3の概略構成を示すブロック図である。図13に示されるレーザレーダ装置3の構成は、上記実施の形態1の光変調器12に代えて、図13の光変調器12Aを有する点を除いて、上記実施の形態1のレーザレーダ装置1の構成と同じである。
 本実施の形態の光変調器12Aは、光分配器11から入力された送信光TLに対して周波数変調及び強度変調を実行することにより、送信光TLの周波数をシフトさせ、かつ送信光TLをパルス化する変調器である。光変調器12Aは、周波数変調及び強度変調により、複数の変調光成分ML~MLを含み、所定のPRFで所定のパルス幅Tを有する変調送信光ML(一連のパルス光)を生成し、変調送信光MLを光サーキュレータ13に出力する。また、光変調器12Aは、送信光TLのパルス化のために生成されたパルストリガ信号TSを信号処理部20に出力する。パルストリガ信号TSは、送信光TLのパルス化のタイミングを表す信号である。
 以上に説明したように本実施の形態では、光分配器11は、複数の送信光成分TL~TLを光変調器12Aに出力し、光変調器12Aは、これら送信光成分TL~TLに周波数変調及び強度変調を施すことで、複数の変調光成分ML~MLを含む変調送信光MLを生成する。よって、光アンテナ14は、広いビーム径と広いビーム拡がり角を有する変調送信光MLを光スキャナ15を介して送信することができる。このため、以下に説明するように、光アンテナ14の受信口径Dを改善することができ、受信信号のSNRをも改善することができる。
 すなわち、受信信号(アナログ検波信号CS)のSNRは、非特許文献(S. Kameyama et al., Applied Optics, Vol. 46, No. 11, pp. 1953-1962, 2007.)に記載のとおり、受信口径Dに対して二乗(=D)で改善される。受信口径Dは、レンズの焦点距離をf、受光素子の開口数をNAとすると、D=NA×fで定義される。fが大きい方が受信口径Dは大きくなる。一方、受信視野αは、受信機の素子口径をdとすると、簡易的に、α=d/fで定義されるため、dが一定のままで、fを大きくすると受信視野αが小さくなる。
 本実施の形態では、光受信器アレイ17によって光アンテナ14の受信口径Dと受信視野α(アレイ数をNとすると、素子口径はd=d×N)の双方の値の改善ができる。本実施の形態では、光アンテナ14の受信視野とほぼ同等のビーム径とビーム拡がり角を有する変調送信光MLを送信することができる。したがって、光アンテナ14の受信口径Dを改善することができ、受信信号(アナログ検波信号CS)のSNRをも改善することができる。
 また、光受信器アレイ17が、シングルモードのレーザ光しか受信できないヘテロダイン検波法を採用する場合でも、広い受信視野と広い受信口径とを確保することができるため、SNRを改善することができる。
実施の形態4.
 図14は、本発明に係る実施の形態4であるレーザレーダ装置4の概略構成を示すブロック図である。図14に示されるレーザレーダ装置4の構成は、上記実施の形態2の光変調器12に代えて、図14の光変調器12Aを有する点を除いて、上記実施の形態2のレーザレーダ装置2の構成と同じである。
 本実施の形態の光変調器12Aは、上記実施の形態3の場合と同様に、光分配器11から入力された送信光TLに対して周波数変調及び強度変調を実行することにより、複数の変調光成分ML~MLを含む変調送信光ML(一連のパルス光)を光サーキュレータ13に出力する。また、光変調器12Aはパルストリガ信号TSを信号処理部20に出力する。
 したがって、上記実施の形態3の場合と同様に、光アンテナ14の受信視野とほぼ同等のビーム径とビーム拡がり角を有する変調送信光MLを送信することができる。よって、光アンテナ14の受信口径Dを改善することができ、受信信号(アナログ検波信号CS)のSNRをも改善することができる。
実施の形態1~4の変形例.
 以上、図面を参照して本発明に係る種々の実施の形態について述べたが、これら実施の形態は本発明の例示であり、これら実施の形態以外の様々な形態を採用することもできる。たとえば、上記実施の形態1~4では、パルス波を送信光として使用するパルス型のレーザレーダ装置1~4であったが、本発明はパルス型に限定されるものではない。連続波レーザ光を送信光として使用するCW(Continuous Wave)方式が上記レーザレーダ装置1~4に適用されてもよい。
 また、上記実施の形態1~4では、図2に示したように光学部品間の光接続手段として空間伝搬構造が記載されていたが、これに限定されるものではない。空間伝搬構造に代えて、光ファイバを使用した光導波路構造が採用されてもよい。
 また、上記実施の形態1~4は、速度ベクトルvを観測量として計測しているが、観測量は、速度ベクトルvに限定されるものではない。たとえば、積算スペクトルISを用いて、構造物または地形物などのハードターゲットの移動速度または測距値を算出することができるように信号処理部20,20Aの構成が適宜変更されてもよい。
 また、上記実施の形態1~4は、光スキャナ15を用いて外部空間を走査しているが、外部空間を走査する方法は、光スキャナ15を使用する方法に限定されるものではない。たとえば、航空機または自動車などの移動体の移動に伴い、この移動体から出射されるレーザビームで外部空間を走査する方法が採用されてもよい。この場合、角度情報ADに代えて、外部空間に対する移動体の相対移動速度を用いて観測量を算出するように信号処理部20,20Aの構成が適宜変更されればよい。
 なお、本発明の範囲内において、上記実施の形態1~4の自由な組み合わせ、各実施の形態の任意の構成要素の変形、または各実施の形態の任意の構成要素の省略が可能である。
 本発明に係るレーザレーダ装置は、レーザ光を用いて大気状態を観測する観測システムをはじめとする、観測対象物の形状を計測する3次元計測システムなどのレーザセンサシステムに用いられるのに適している。
 1~4 レーザレーダ装置、10 光源アレイ、101 基準光源、102 光分配素子、11 光分配器、111 導光光学系、112 ビームスプリッタ、12,12A 光変調器、13 光サーキュレータ、131 ビームスプリッタ、132 波長板、14 光アンテナ、15 光スキャナ、151 回転子、152 ウェッジプリズム、153 回転駆動部、16 光合波器、161 ビームスプリッタ、162 反射ミラー、163,164 集光光学系、17 光受信器アレイ、18 スイッチ回路、181 第1スイッチ、182 第2スイッチ、183 加算器、184 チャンネル選択回路、19,19A A/D変換器、20,20A 信号処理部、21,21A 信号分割器、22,22A 領域変換部、23,23A 積算器、24 SNR算出器、25 選択器(セレクタ)、26 速度算出器、27 ベクトル算出器、28 走査速度算出器、40 信号処理装置、41 プロセッサ、42 メモリ、43 入力インタフェース部、44 出力インタフェース部、45 信号路、Da1~Da3,Db1~Db3 受光素子、E1~E3 光出射端。

Claims (15)

  1.  複数の光出射端を有し、前記複数の光出射端からそれぞれ複数のレーザ光を同時に出射する光源アレイと、
     前記複数のレーザ光から当該複数のレーザ光の一部を送信光として分離し、かつ前記複数のレーザ光からそれぞれ複数のローカル光成分を分離する光分配器と、
     前記送信光を変調して変調送信光を生成する光変調器と、
     前記変調送信光で外部空間を走査しつつ、当該外部空間に存在するターゲットで反射された当該変調送信光を受信光として受信する送受信光学系と、
     前記受信光と前記複数のローカル光成分とを合波することにより複数の干渉光成分を生成する光合波器と、
     前記送受信光学系の互いに異なる複数の受信視野と光学的に対応する位置にそれぞれ配置され、前記複数の干渉光成分を検波することで複数の検波信号を生成する複数の光受信器を含む光受信器アレイと、
     前記外部空間に対する前記送受信光学系の走査速度に応じて、前記複数の検波信号の中から検波信号を選択するスイッチ回路と、
     当該選択された検波信号に基づいて前記ターゲットの状態を示す観測量を算出する信号処理部と
    を備えることを特徴とするレーザレーダ装置。
  2.  請求項1記載のレーザレーダ装置であって、前記光源アレイは、基準レーザ光を出力する単一の基準光源と、前記基準レーザ光を前記複数のレーザ光に分配する光分配素子とを含むことを特徴とするレーザレーダ装置。
  3.  請求項1記載のレーザレーダ装置であって、前記信号処理部は、前記スイッチ回路で選択された検波信号の信号スペクトルを算出し、前記信号スペクトルに基づいて前記観測量を算出することを特徴とするレーザレーダ装置。
  4.  請求項3記載のレーザレーダ装置であって、
     前記信号処理部は、
     前記スイッチ回路で選択された検波信号を、複数のゲート時間領域の信号波形をそれぞれ表す複数の時間領域信号に分割する信号分割器と、
     前記複数の時間領域信号を複数の周波数領域信号にそれぞれ変換する領域変換器と、
     前記複数の周波数領域信号から前記複数の時間領域信号それぞれの複数のスペクトルを算出し、前記複数のスペクトルを積算することにより前記信号スペクトルを算出する積算器と、
     前記信号スペクトルに基づいて前記観測量を算出する観測量算出器と
    を含むことを特徴とするレーザレーダ装置。
  5.  請求項1記載のレーザレーダ装置であって、前記信号処理部は、前記ターゲットの速度ベクトルを前記観測量として算出することを特徴とするレーザレーダ装置。
  6.  請求項1記載のレーザレーダ装置であって、
     前記光分配器は、前記複数のレーザ光から複数の送信光成分をそれぞれ分離し、当該複数の送信光成分を前記送信光として出力し、
     前記光変調器は、前記複数の送信光成分のうちの1つの送信光成分を変調して前記変調送信光を生成することを特徴とするレーザレーダ装置。
  7.  請求項1記載のレーザレーダ装置であって、
     前記光分配器は、前記複数のレーザ光から複数の送信光成分をそれぞれ分離し、当該複数の送信光成分を前記送信光として出力し、
     前記光変調器は、前記複数の送信光成分を変調して前記変調送信光を生成することを特徴とするレーザレーダ装置。
  8.  複数の光出射端を有し、前記複数の光出射端からそれぞれ複数のレーザ光を同時に出射する光源アレイと、
     前記複数のレーザ光から当該複数のレーザ光の一部を送信光として分離し、かつ前記複数のレーザ光からそれぞれ複数のローカル光成分を分離する光分配器と、
     前記送信光を変調して変調送信光を生成する光変調器と、
     前記変調送信光で外部空間を走査しつつ、当該外部空間に存在するターゲットで反射された当該変調送信光を受信光として受信する送受信光学系と、
     前記受信光と前記複数のローカル光成分とを合波することにより複数の干渉光成分を生成する光合波器と、
     前記送受信光学系の互いに異なる複数の受信視野と光学的に対応する位置にそれぞれ配置され、前記複数の干渉光成分を検波することで複数の検波信号を生成する複数の光受信器を含む光受信器アレイと、
     前記複数の検波信号それぞれの複数の信号スペクトルを算出し、前記複数の検波信号それぞれの複数の信号対雑音電力比を算出する信号処理部と
    を備え、
     前記信号処理部は、前記複数の信号対雑音電力比を用いて前記複数の信号スペクトルの中から少なくとも1つの信号スペクトルを選択し、当該選択された少なくとも1つの信号スペクトルに基づいて前記ターゲットの状態を示す観測量を算出する、
    ことを特徴とするレーザレーダ装置。
  9.  請求項8記載のレーザレーダ装置であって、前記光源アレイは、基準レーザ光を出力する単一の基準光源と、前記基準レーザ光を前記複数のレーザ光に分配する光分配素子とを含むことを特徴とするレーザレーダ装置。
  10.  請求項8記載のレーザレーダ装置であって、前記信号処理部は、前記複数の信号スペクトルの中から選択された少なくとも2つの信号スペクトルを加算または平均化することでノイズ低減スペクトルを算出し、前記ノイズ低減スペクトルに基づいて前記観測量を算出することを特徴とするレーザレーダ装置。
  11.  請求項8記載のレーザレーダ装置であって、
     前記信号処理部は、
     前記複数の検波信号の各々を、複数のゲート時間領域の信号波形をそれぞれ表す複数の時間領域信号に分割する信号分割器と、
     前記複数の検波信号の各々について、前記複数の時間領域信号を複数の周波数領域信号にそれぞれ変換する領域変換器と、
     前記複数の検波信号の各々について、前記複数の周波数領域信号から前記複数の時間領域信号それぞれの複数のスペクトルを算出し、当該算出された複数のスペクトルを積算することにより前記複数の検波信号にそれぞれ対応する複数の積算スペクトルを前記複数の信号スペクトルとして算出する積算器と、
     前記複数の積算スペクトルから前記複数の信号対雑音電力比をそれぞれ算出するSNR算出器と、
     前記複数の信号対雑音電力比を用いて前記複数の積算スペクトルの中から少なくとも1つの積算スペクトルを選択する選択器と、
     当該選択された少なくとも1つの積算スペクトルに基づいて前記観測量を算出する観測量算出器と
    を含むことを特徴とするレーザレーダ装置。
  12.  請求項11記載のレーザレーダ装置であって、
     前記選択器は、前記複数の積算スペクトルの中から選択された少なくとも2つの積算スペクトルを加算することでノイズ低減スペクトルを算出し、
     前記観測量算出器は、前記ノイズ低減スペクトルに基づいて前記観測量を算出する、
    ことを特徴とするレーザレーダ装置。
  13.  請求項8記載のレーザレーダ装置であって、前記信号処理部は、前記ターゲットの速度ベクトルを前記観測量として算出することを特徴とするレーザレーダ装置。
  14.  請求項8記載のレーザレーダ装置であって、
     前記光分配器は、前記複数のレーザ光から複数の送信光成分をそれぞれ分離し、当該複数の送信光成分を前記送信光として出力し、
     前記光変調器は、前記複数の送信光成分のうちの1つの送信光成分を変調して前記変調送信光を生成することを特徴とするレーザレーダ装置。
  15.  請求項8記載のレーザレーダ装置であって、
     前記光分配器は、前記複数のレーザ光から複数の送信光成分をそれぞれ分離し、当該複数の送信光成分を前記送信光として出力し、
     前記光変調器は、前記複数の送信光成分を変調して前記変調送信光を生成することを特徴とするレーザレーダ装置。
PCT/JP2017/024496 2017-07-04 2017-07-04 レーザレーダ装置 WO2019008670A1 (ja)

Priority Applications (4)

Application Number Priority Date Filing Date Title
JP2017566166A JP6395958B1 (ja) 2017-07-04 2017-07-04 レーザレーダ装置
US16/618,303 US11899112B2 (en) 2017-07-04 2017-07-04 Laser radar device
EP17916881.0A EP3637134B1 (en) 2017-07-04 2017-07-04 Laser radar device
PCT/JP2017/024496 WO2019008670A1 (ja) 2017-07-04 2017-07-04 レーザレーダ装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2017/024496 WO2019008670A1 (ja) 2017-07-04 2017-07-04 レーザレーダ装置

Publications (1)

Publication Number Publication Date
WO2019008670A1 true WO2019008670A1 (ja) 2019-01-10

Family

ID=63668514

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2017/024496 WO2019008670A1 (ja) 2017-07-04 2017-07-04 レーザレーダ装置

Country Status (4)

Country Link
US (1) US11899112B2 (ja)
EP (1) EP3637134B1 (ja)
JP (1) JP6395958B1 (ja)
WO (1) WO2019008670A1 (ja)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20210117907A (ko) * 2020-03-20 2021-09-29 한국과학기술원 셀프-헤테로다인 검파를 이용하여 거리를 측정하는 라이다 장치

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP6366880B2 (ja) * 2016-04-05 2018-08-01 三菱電機株式会社 レーザレーダ装置
US10761195B2 (en) 2016-04-22 2020-09-01 OPSYS Tech Ltd. Multi-wavelength LIDAR system
KR102592139B1 (ko) 2017-03-13 2023-10-23 옵시스 테크 엘티디 눈-안전 스캐닝 lidar 시스템
KR20220119769A (ko) 2017-07-28 2022-08-30 옵시스 테크 엘티디 작은 각도 발산을 갖는 vcsel 어레이 lidar 송신기
JP7117092B2 (ja) * 2017-09-25 2022-08-12 株式会社トプコン レーザ測定方法及びレーザ測定装置
KR102634870B1 (ko) 2017-11-15 2024-02-13 옵시스 테크 엘티디 잡음 적응형 솔리드-스테이트 lidar 시스템
KR102634880B1 (ko) 2018-04-01 2024-02-08 옵시스 테크 엘티디 잡음 적응형 솔리드-스테이트 lidar 시스템
US10838047B2 (en) 2018-04-17 2020-11-17 Santec Corporation Systems and methods for LIDAR scanning of an environment over a sweep of wavelengths
US20200041614A1 (en) * 2018-08-03 2020-02-06 OPSYS Tech Ltd. Distributed Modular Solid-State LIDAR System
WO2020062718A1 (en) * 2018-09-30 2020-04-02 SZ DJI Technology Co., Ltd. Optical scanning device with beam compression and expansion
CN113692540A (zh) 2019-04-09 2021-11-23 欧普赛斯技术有限公司 带激光控制的固态lidar发送器
WO2020242834A1 (en) 2019-05-30 2020-12-03 OPSYS Tech Ltd. Eye-safe long-range lidar system using actuator
US11513228B2 (en) * 2020-03-05 2022-11-29 Santec Corporation Lidar sensing arrangements
US11486792B2 (en) 2020-06-05 2022-11-01 Santec Corporation Tunable light source for optical fiber proximity and testing
JP7138828B2 (ja) 2020-08-05 2022-09-16 三菱電機株式会社 測距装置、測距方法及びレーダ装置
CN115685268B (zh) * 2023-01-03 2023-03-31 成都奇芯微电子有限公司 基于低轨mcsk调制的数据解调判决方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007248146A (ja) * 2006-03-14 2007-09-27 Omron Corp レーダ装置
JP2007316016A (ja) * 2006-05-29 2007-12-06 Mitsubishi Electric Corp レーダ装置
US20080278715A1 (en) * 2007-05-11 2008-11-13 Rosemount Aerospace Inc. Scanning ladar with adjustable operational parameters
JP2009162678A (ja) 2008-01-09 2009-07-23 Mitsubishi Electric Corp レーザレーダ装置
WO2014024508A1 (ja) * 2012-08-08 2014-02-13 三菱電機株式会社 レーダ装置
JP2016105082A (ja) * 2014-11-19 2016-06-09 株式会社豊田中央研究所 レーザレーダ装置およびレーザレーダ装置の光受信方法
WO2016092705A1 (ja) * 2014-12-12 2016-06-16 三菱電機株式会社 レーザレーダ装置
WO2016117159A1 (ja) 2015-01-21 2016-07-28 三菱電機株式会社 レーザレーダ装置
WO2017010176A1 (ja) * 2015-07-14 2017-01-19 コニカミノルタ株式会社 レーザレーダ装置
WO2017085799A1 (ja) * 2015-11-18 2017-05-26 三菱電機株式会社 レーザレーダ装置
WO2017098623A1 (ja) * 2015-12-10 2017-06-15 三菱電機株式会社 レーザレーダ装置

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6580497B1 (en) * 1999-05-28 2003-06-17 Mitsubishi Denki Kabushiki Kaisha Coherent laser radar apparatus and radar/optical communication system
FR2948459B1 (fr) * 2009-07-24 2011-08-19 Thales Sa Procede et systeme d'anemometrie mono-particule par lidar
EP2568314B1 (en) * 2010-05-07 2021-06-09 Mitsubishi Electric Corporation Laser radar device
US8294879B2 (en) * 2010-11-15 2012-10-23 Raytheon Company Multi-directional active sensor
US8659748B2 (en) * 2011-02-15 2014-02-25 Optical Air Data Systems, Llc Scanning non-scanning LIDAR
CA2857190A1 (en) * 2011-11-29 2013-06-06 Flidar Motion-stabilised lidar and method for wind speed measurement
WO2013094431A1 (ja) * 2011-12-21 2013-06-27 三菱電機株式会社 レーザレーダ装置
EP2896972B1 (en) * 2012-09-14 2019-10-30 Mitsubishi Electric Corporation Laser radar device and method for detecting speed of measured object
CN102890272B (zh) * 2012-11-05 2013-12-04 中国航天科工集团第二研究院二十三所 一种毫米波云雷达信号处理方法
KR102422784B1 (ko) * 2015-08-03 2022-07-19 엘지이노텍 주식회사 광파 탐지 및 거리 측정 장치

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007248146A (ja) * 2006-03-14 2007-09-27 Omron Corp レーダ装置
JP2007316016A (ja) * 2006-05-29 2007-12-06 Mitsubishi Electric Corp レーダ装置
US20080278715A1 (en) * 2007-05-11 2008-11-13 Rosemount Aerospace Inc. Scanning ladar with adjustable operational parameters
JP2009162678A (ja) 2008-01-09 2009-07-23 Mitsubishi Electric Corp レーザレーダ装置
WO2014024508A1 (ja) * 2012-08-08 2014-02-13 三菱電機株式会社 レーダ装置
JP2016105082A (ja) * 2014-11-19 2016-06-09 株式会社豊田中央研究所 レーザレーダ装置およびレーザレーダ装置の光受信方法
WO2016092705A1 (ja) * 2014-12-12 2016-06-16 三菱電機株式会社 レーザレーダ装置
WO2016117159A1 (ja) 2015-01-21 2016-07-28 三菱電機株式会社 レーザレーダ装置
WO2017010176A1 (ja) * 2015-07-14 2017-01-19 コニカミノルタ株式会社 レーザレーダ装置
WO2017085799A1 (ja) * 2015-11-18 2017-05-26 三菱電機株式会社 レーザレーダ装置
WO2017098623A1 (ja) * 2015-12-10 2017-06-15 三菱電機株式会社 レーザレーダ装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
S. KAMEYAMA ET AL., APPLIED OPTICS, vol. 46, no. 11, 2007, pages 1953 - 1962
See also references of EP3637134A4

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20210117907A (ko) * 2020-03-20 2021-09-29 한국과학기술원 셀프-헤테로다인 검파를 이용하여 거리를 측정하는 라이다 장치
KR102425721B1 (ko) * 2020-03-20 2022-07-29 한국과학기술원 셀프-헤테로다인 검파를 이용하여 거리를 측정하는 라이다 장치

Also Published As

Publication number Publication date
US11899112B2 (en) 2024-02-13
JPWO2019008670A1 (ja) 2019-07-04
US20210157000A1 (en) 2021-05-27
EP3637134A1 (en) 2020-04-15
EP3637134A4 (en) 2020-05-06
EP3637134B1 (en) 2021-04-28
JP6395958B1 (ja) 2018-09-26

Similar Documents

Publication Publication Date Title
JP6395958B1 (ja) レーザレーダ装置
US11940571B2 (en) Performing speckle reduction using polarization
US11762069B2 (en) Techniques for determining orientation of a target using light polarization
US11435453B1 (en) Techniques for simultaneous determination of range and velocity with active modulation
US11960032B2 (en) Techniques for increasing effective power in multi-beam LIDAR systems
US11054524B2 (en) Optimizing a lidar system using sub-sweep sampling
US11927699B1 (en) Techniques for range and velocity measurements in a lidar system
JP2023545775A (ja) コヒーレントlidarシステムにおけるゴースト低減技術
JP2021536586A (ja) 偏光エンコードされたビーム送出及び収集
US20220404498A1 (en) Distance and speed measuring apparatus
CN115407313B (zh) 多通道激光雷达
JP2023547877A (ja) コヒーレントlidarシステムにおける複数ターゲットのピーク関連付け技術
US11372105B1 (en) FMCW LIDAR using array waveguide receivers and optical frequency shifting
US11953627B2 (en) Techniques for multiplexing optical beams in coherent LiDAR systems
US12013466B2 (en) Techniques for using active amplitude modulation for simultaneous determination of range and velocity in an FMCW LIDAR system
US20240183947A1 (en) Techniques for providing a variety of lidar scan patterns
US20240036178A1 (en) Techniques for using a coherent receiver in a fmcw lidar system
JP7074889B2 (ja) センサ装置
JP2024512064A (ja) 同相・直交位相(iq)処理を用いたコヒーレントlidarシステムにおけるゴースト低減技術
CN116973926A (zh) 激光雷达
CN117337401A (zh) 利用同相/正交相位(iq)处理的用于减轻相干lidar系统中的重影的技术

Legal Events

Date Code Title Description
ENP Entry into the national phase

Ref document number: 2017566166

Country of ref document: JP

Kind code of ref document: A

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

Ref document number: 17916881

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 2017916881

Country of ref document: EP

Effective date: 20191211

NENP Non-entry into the national phase

Ref country code: DE