CN104049259A - Lidar three-dimensional imaging system based on virtual instrument - Google Patents

Lidar three-dimensional imaging system based on virtual instrument Download PDF

Info

Publication number
CN104049259A
CN104049259A CN201410308788.2A CN201410308788A CN104049259A CN 104049259 A CN104049259 A CN 104049259A CN 201410308788 A CN201410308788 A CN 201410308788A CN 104049259 A CN104049259 A CN 104049259A
Authority
CN
China
Prior art keywords
model
laser
dimensional
lidar
waveform
Prior art date
Legal status (The legal status 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 status listed.)
Granted
Application number
CN201410308788.2A
Other languages
Chinese (zh)
Other versions
CN104049259B (en
Inventor
王元庆
徐帆
曾真
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University
Original Assignee
Nanjing University
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 Nanjing University filed Critical Nanjing University
Priority to CN201410308788.2A priority Critical patent/CN104049259B/en
Publication of CN104049259A publication Critical patent/CN104049259A/en
Application granted granted Critical
Publication of CN104049259B publication Critical patent/CN104049259B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88Lidar systems specially adapted for specific applications
    • G01S17/89Lidar systems specially adapted for specific applications for mapping or imaging
    • 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/497Means for monitoring or calibrating

Abstract

A lidar three-dimensional imaging system based on a virtual instrument comprises a three-dimensional scene modeling module, a lidar testing system simulation environment modeling module, a full-waveform signal processing module and a three-dimensional reconstruction module. The three-dimensional scene modeling module comprises a three-dimensional model loading element used for setting background colors, visual angle control, illumination, a projection mode, a display model and other basic scene projects. The lidar testing system simulation environment modeling module comprises four modeling sub-modules which are a laser pulse model, an atmospheric transmission model, an object interaction model and a receiving unit model. The laser pulse model is used for simulating a laser source according to the wavelength, pulse width, energy and other characteristics of a laser. The atmospheric transmission model is used for simulating a tested atmospheric environment, so that a noise model is generated for acting on the laser. The object interaction model is used for simulating effects of the laser and a detected object, and the effects comprise mirror reflection, diffuse reflection, surface reflection and speckles. The receiving unit model is used for simulating detector noise and amplifier noise which are generated after sensing.

Description

Laser radar 3-D imaging system based on virtual instrument
One, technical field
The present invention is the laser radar 3-D imaging system based on virtual instrument, belongs to laser radar detection technical field.
Two, background technology
Laser radar detection and range measurement system (Light Detection and Ranging, be called for short Lidar) be a kind of active remote sensing technology, it has broken through traditional passive remote sensing imaging mechanism, be subject to that weather effect, illumination effect are little, automaticity high, there is unrivaled advantage.
Airborne imaging Lidar is by IMU aligning, DGPS acquisition device projection centre position, IMU/DGPS system refers to utilizes dress GPS receiver aboard and the GPS receiver being located on ground one or more base station synchronously and continuously to observe gps satellite signal, obtain the location parameter of aerial surveying camera by gps carrier phase measurement differential position, high precision Inertial Measurement Unit (the IMU that application and aerial surveying camera are closely connected, Inertial Measurement Uint) directly measure the attitude parameter of aerial surveying camera, pass through IMU, the associating post-processing technology of DGPS data obtains the photogrammetric measurement technology of every required photo high precision elements of exterior orientation of mapping, by Point Cloud Processing and filtering processing, three-dimensional reconstruction DEM.
In conjunction with multi-door ambits such as optics, radar, signal processing, computer visions, its range of application is across sea, land and sky for laser radar technique, is that recent domestic develops rapidly and the technical field of widespread use.Laser radar is applied to aerospace field the earliest, afterwards gradually for military affairs, fields such as laser guidance, battle reconnaissance, airplane anti-collision, land mine remote sensing.At present, laser radar technique is widely used in the various fields such as city digital Model Reconstruction, forest ecology monitoring, marine environment mapping, geology and geomorphology detection, space exploration.
The backscatter waveform of (Full-Waveform Lidar) the paired pulses transmitted waveform of all-wave laser radar and each pixel carries out data acquisition record, can decomposite flight time (TOF) and the strength information of each pixel on a time shaft.
The method of laser infrared radar imaging is mainly divided into three kinds:
(1) adopt single detector, only measure a pixel at every turn.This method is Detection Techniques the earliest, by the flight time of single measurement different pixels and strength information, recovers one by one distance and the reflectivity of each pixel.
(2) adopt planar array detector, each detectable multiple pixels.By utilizing emitted light is covered to whole target simultaneously, echoed signal contains the information of multiple pixels, and receiver can draw the distance of different pixels by demodulation, separation.This method is for large spot imaging, and more single of the pixel number of single measurement is measured obviously and improved.
(3) adopt APD (snow spring photodiode) detector array, the multiple pixels of standalone probe.Restraint on N the point of beating respectively in target by light beam being divided into N, receive by N unit of APD array, thus can the while information of a decomposing N point exactly.Scan all pixels of target by scanner, final imaging.The method has adopted highly sensitive APD sensor array, and two kinds of methods, can ensure high resolving power one-tenth figure above, can ensure again high hasty map.
Full wave shape echoed signal modeling and simulating is an important research content of Full wave shape Airborne Lidar survey technology with analyzing.The achievement of an actual measurement system often needs simulation modeling in earlier stage to carry out theoretical validation.About laser radar modeling and simulating aspect, there were a lot of successfully technological achievements abroad, and domestic to this rare systematic research achievement.
(1) the 3D Lidar system emulation of Lincoln laboratory based on OPENGL exploitation the hidden airborne 3D ladar test off the net of leafage, by interpretation of result, the result matching degree of simulation result and measured test is in 20%.
(2) Graham of Australia has carried out emulation based on Matlab to laser radar vegetation penetration test under carrying platform different motion track.
(3) Kim of the U.S., by setting up the Monte Carlo model of laser penetration trees, obtains Lidar full-wave simulation signal.This model, by Lidar system and natural scene, comprises atmosphere, and trees and landform finally obtain 3D model and the reflection characteristic of vegetation.
(4) German JUTZI, by the all-wave signal on the different gradients of emulation slope, utilizes and measures the gradient of determining plane with the correlativity of estimating ripple.
In existing technological achievement, not occurring can integrated control Lidar Simulation experiment scene, again can be for application platform provides hardware module interface, realize the Lidar 3-D imaging system of the expanded application of laser radar actual measurement experiment.Application of Virtual based on LabVIEW is a kind of brand-new exploitation implementation to laser infrared radar imaging field, by building of convenient system Visualization Platform to a great extent.LabVIEW is a graphical programming software based on G (Graphic) language.It adopts GRAPHICAL PROGRAMMING method, uses graphic symbolic expression program behavior, uses visualization technique to set up man-machine interface, more directly perceived, succinct.Abundant function library and subroutine that LabVIEW comprises as data acquisition, signal processing, I/O drive supervisor storehouse, provide multilingual development interface simultaneously, become the ideal tools of the application such as development and testing, measurement, data acquisition.
Three, summary of the invention
The object of the invention is the platform based on LabVIEW, develop the emulation test system of a set of complete laser radar three-dimensional imaging and the expansion of surveying 3-D imaging system based on the Lidar of virtual instrument is provided.
Laser radar detection based on virtual instrument and a range measurement system three-dimensional imaging analogue system, comprise three-dimensional scenic modeling, the modeling of Lidar pilot system simulated environment, the processing of all-wave signal, four functional modules of three-dimensional reconstruction; The feature of modules is described below:
(1) three-dimensional scenic MBM: read in three-dimension developing instrument 3DSMAX, the three-D elements of the forms such as the STL/WRL/ASE that AutoCAD produces, background colour is set, viewing angle control, illumination, projection mode, the basic scene project such as display model, to the landform that comprises of reading in, vegetation, the scope texture of aircraft, size, position, number is controlled, the wherein corresponding lasing light emitter in the position of aircraft position, dynamic change in test, finally present airborne Lidar imaging test dynamic scene, and export final digital elevation model DSM altitude figures, in order to follow-up test,
(2) Lidar pilot system simulated environment MBM: comprise laser pulse model, propagation in atmosphere model, target interaction models, four sub-module modelings of receiving element model; Laser pulse model is simulated lasing light emitter according to character such as sharp light wavelength, pulsewidth, energy; Propagation in atmosphere model is simulated the atmospheric environment of test, and generted noise model is to laser effect; Target interaction models, to laser and detection of a target effect simulation, comprises mirror-reflection, diffuse reflection, surface reflection, speckle; Noise of detector, amplifier noise that receiving element model produces after in order to analog sensed.By to the modeling of four model orders, the systematic parameters such as difference of elevation, temperature, flight angle are set, finally export the all-wave signal of each time period target respective pixel point;
(3) all-wave signal processing module: process the pixel degree of depth, the strength information information recovered for the echoed signal that each time period produces.This module mainly comprises signal enhancing, Wavelet Denoising Method, level and smooth, four parts of waveform fitting;
(4) three-dimensional reconstruction module: this module is tested final visual three-dimensional one-tenth figure is provided for airborne laser radar; The complete degree of depth and the strength information of target providing for all-wave signal processing module, generates Range Profile, intensity image; Be divided into non-ground point and ground point two category informations by segmentation and classification; Ground point information is set up grid, and ground point is carried out to triangle division, and aggregation strength information is finally rebuild generating digital elevation model DEM.
Further: in the corresponding different measuring cycle of aircraft position dynamic change in (1),, to contain pixel number in conjunction with hot spot and obtain pixel number and the echoed signal thereof that one-period records simultaneously laser light splitting according to APD number.
Further: (3) dynamically echoed signal search to each cycle, sampling unique point, specific aim processing, realizes the integrated function of signal enhancing, Wavelet Denoising Method, level and smooth, waveform fitting; In order to obtain higher measuring accuracy, waveform fitting module arranges matching template according to the actual waveform of laser pulse.
Further: based on the main development platform of LabVIEW, in conjunction with the mixed developing pattern of LabWindows/CVI, Matlab, C++, each function sub-modules (sub-VI) is called by father's module (father VI) according to inheritance, is finally integrated into the main gui interface of the superiors; Present the visual simulating test of complete laser radar three-dimensional imaging by carrying out total GUI.
Propose a kind of solution of the actual measurement of the Lidar based on virtual instrument imaging system, reserved TCP, serial communication interface, loaded measured data by collector, completed the function of Lidar eyeball cloud data three-dimensional imaging.
This platform utilizes the advantage of LabVIEW graphic language programming, uses dexterously three-dimensional value-added tax function control, builds adjustable laser radar experiment scene; To Lidar pilot system simulated environment, comprise that laser instrument, atmosphere, target are mutual, receiving element carries out comprehensive simulated, implement l-G simulation test and generate all-wave signal; Build the all-wave signal processing method of a cover system, all-wave signal is processed to resolve a picture dot point cloud information; Three-dimensional reconstruction function realizes module to cloud data classified filtering, final three-dimensional reconstruction DEM.
Beneficial effect of the present invention: virtual instrument framework: virtual instrument technique provides abundant high-performance modularized hardware interface, its Software tool LabVIEW provide pattern development platform, the simultaneously integrated abundant algorithms libraries such as data acquisition, signal processing, image processing that comprise, greatly shorten the construction cycle, there is very strong scalability simultaneously, realized software and hardware seamless integrated.Traditional upper imaging system of laser infrared radar imaging system is nearly all based on the platform such as Matlab, VC, realizing in experiments of measuring advantage not obvious.And it is very rare that virtual instrument is applied to laser infrared radar imaging field, and have Lidar visible, the system simultaneously with complete echo signal processing, imaging integrated treatment never occurs.
Three-dimensional visualization dynamic test scene construction: utilize LabVIEW tri-dimensional picture to load control, provide the three-dimensional model element form that professional cartography instrument is generated to load, by three-dimensional model deformation, Lidar l-G simulation test Visual Scene is built in the controls such as texture.Whole l-G simulation test is based on LabVIEW platform, can be in conjunction with multiple types of tools exploitations such as LabWindows/CVI, VC++, Matlab, in conjunction with three-dimension developing instruments such as 3DSMAX, AutoCAD, Pro/E, SolidWorks, VRMLPAD, by dissimilar 3D Unified Model and set up Lidar test scene, feature richness, extendability is strong, effect of visualization is good, complete and presented intuitively the test of all-wave laser radar three-dimensional imaging and functional verification is provided.Meanwhile, system provides the upper interface of actual measurement system, and to measured test, three-dimensional imaging provides the solution of expansion.
Four, brief description of the drawings
Fig. 1 is the structured flowchart of the imaging Lidar emulation platform based on virtual instrument;
Fig. 2 is that three-dimensional scenic module realizes three-dimensional modeling process flow diagram;
Fig. 3 is the process flow diagram that analogue system is set up;
Fig. 4 is the flow chart of data processing figure of all-wave signal processing module;
Fig. 5 is three-dimensional reconstruction module process flow diagram;
Fig. 6 is all-wave sampled signal waveform figure;
Fig. 7 is the imaging Lidar simulation software operational flowchart based on LabVIEW;
Fig. 8 is the main gui interface of the imaging Lidar simulation software based on LabVIEW;
Fig. 9 is three-dimensional scenic module modeling effect schematic diagram;
Figure 10 is the final three-dimensional reconstruction remodelling of Lidar l-G simulation test DEM;
Five, embodiment
Below in conjunction with accompanying drawing, implementer's case of the present invention is elaborated.
As shown in Figure 1, the present invention mainly comprises four functional modules: three-dimensional scenic modeling, the modeling of Lidar pilot system simulated environment, the processing of all-wave signal, three-dimensional reconstruction.Wherein three-dimensional scenic MBM can load the 3D model element of 3D mapping software output, and it is carried out to layout, adjusts, and makes it to produce required experiment scene.The establishment mode of 3D model element is to enrich flexibly.For example, the model element of landform can use the topography and geomorphology plug-in unit terrian plug-in unit in 3dsmax, generates the real figure landform of specifying longitude and latitude by connecting GIS; Vegetation model can utilize the tree plug-in unit structure of making of 3dsmax to specify vegetation group; Model aircraft can be at AutoCAD according to designing in dimensional requirement.
Shown in Fig. 2 is that three-dimensional scenic is built, and can load the 3D model file of 3D mapping software output, can obtain the implementation method of vertical array mould module by calling scene grid node.LabVIEW provides the information such as VRML, STL, tri-kinds of conventional 3-D geometric model loading methods of ASE, color array, texture coordinate, by setting up geometric model, add object, texture etc. is set calls nodal function and merge dissimilar 3D model buildings three-dimensional scenic.Meanwhile, LabVIEW provides abundant scene control method, comprises the basic setups such as background colour, light source, viewing angle control device, automatic projection mode, display mode, and the distortion control functionality such as scene convergent-divergent, translation, rotation.By merging these functionality controls in interface, foreground, can control flexibly three-dimensional scenic layout, finally present required emulation experiment scene, and export the digital elevation information of scene, use as target simulator experiment.In embodiment, setting aircraft flight height is 2km, and APD number is 8 array elements, and laser frequency is 10kHz, and atural object hot spot footprint is contained 8 pixel information.By the dynamic linear array scanning of for circulation control control aircraft, single sweep operation time corresponding real time 0.1ms.
Fig. 3 has provided Lidar pilot system simulated environment modeling embodiment, and this scheme is according to existing Laser beam energy distribution formula, the modeling of propagation in atmosphere formula.Lidar test simulation need for environment carries out modeling to laser pulse model, propagation in atmosphere model, target interaction models, receiving element model.The content of laser pulse modeling is the structure to Laser Time distribution, irradiance distribution, and input parameter comprises sharp light wavelength, half pulsewidth, the angle of divergence, laser energy, operating distance, frequency.The object of propagation in atmosphere model is the impact of simulation beam broadening, atmospheric turbulence, atmospheric attenuation, acts on laser pulse, and input parameter comprises spectral irradiance, turbulence intensity, visibility distance, luminance factor.The object of target interaction models is the impact that simulated laser pulse and the detection of a target are made mirror-reflection, diffuse reflection that the used time produces, surface reflection, speckle, and input parameter comprises specular reflection factor, the diffuse reflection factor, instantaneous ramp rate.Receiving element model is to want simulated laser effect target back echo signal to be subject to the noise effect of detector, amplifier through receiving trap, and the parameter of input comprises iris radius, instantaneous field of view, bandwidth, light transmission, responsiveness, focal length, dark face electric current, phaeodium electric current, excess noise fact or, detecting device internal gain, average background reflection, optical filter bandwidth, detector capacitors value, detecting device and amplifier total magnification.The correlation parameter of whole system model has difference of elevation, temperature, flight angle, target reflectivity, starts the each module of system simulation model, and order is carried out and produced Lidar Simulation echo waveform.In embodiment, (1) laser pulse time propagates irradiance distribution (2) propagation in atmosphere model comprises beam-broadening, atmospheric turbulence, energy attenuation: 1. beam-broadening rule: ω spread = ω · 1 + 1.624 · ( σ Rytov 2 ) 6 / 5 · Λ 1 , Wherein c nfor textural constant, the λ of atmospheric turbulence are that optical maser wavelength, R are target object distance; 2. atmospheric turbulence rule: atmospheric scinillation probability density function is P turb = 1 S · 2 πσ ln I 2 exp ( - ( ln ( S S av ) + 1 2 σ ln I 2 ) 2 2 σ ln I 2 ) , Wherein σ ln I 2 = ln ( 1 + σ I 2 ) = ln ( 1 + σ t arg et 2 ) , σ t arg et 2 = γ t arg er · σ sat , point 2 = γ t arg et · 1.3 σ point 2 1 + σ point 2 For the average rear effectively turbulent wave momentum of target, wherein r eff=min (r target, ω, υ fov/ 2R) be effective light spot radius (r targetfor the geometric radius of target, ω is the radius of laser at target place, υ fov/ 2R is the visual radius of target of receiver), ρ lfor (one dimension) effectively turbulent flow radius, determined by atmospheric turbulence constant C n, target range R, laser wavelength lambda and laser beam divergence φ; 3. energy attenuation rule: r is operating distance, σ aerit is decay factor; (3) target interaction models comprises mirror-reflection, irreflexive impact, is expressed as: ρ = ρ spec + ρ diff = A cos 6 ( θ ) e - tan 2 ( θ ) s 2 + B cos m ( θ ) , A and B are mirror-reflection and irreflexive association factor; (4) receiving element model is according to the noise equivalent power modeling of detector and amplifier: wherein e cfor electron charge, B is bandwidth (Hz), I dsfor surperficial dark current (A), I dbfor dark current (A), I bfor the electric current (A) that ground unrest produces, I sfor marking current (A), M is detector internal gain, and F is excess noise factor, for the response of detector.
Fig. 4 is the all-wave signal processing module of laser infrared radar imaging system.The function implementation procedure of this module is: all-wave signal, by specifying sampling rate to obtain all-wave sampled data, carries out signal enhancing to this signal, as methods such as threshold value enhancing, low-pass filtering, pulse accumulation.It is 1GSPS that embodiment sets sampling rate, and single-pulse laser pulsewidth is 10ns, and frequency is 10kHz.Echo data for model generation is sampled according to 10 sampled point values of taking out of monopulse, finally obtains actual emulation all-wave sample waveform, as shown in Figure 7.Carry out noise pre-service through low-pass filtering, then get the noise sample of 8 APD echoed signals, carry out region accumulation, y ' (i)=Σ na ny n(i), y n(i) represent input sample signal, a nrepresent weight, the signal to noise ratio (S/N ratio) of corresponding each group signal, n is sample number.
Ethernet host computer based on LabVIEW in embodiment is realized the input of actual measurement sampled data.Implementation method: call " opening TCP connects " function of labview, input server ip address, port address, by the transmission of three-way handshake command communication log-on data, completes the upper bit function of laser all-wave sampled data.The echo data position actual measurement laser all-wave sampled signal now obtaining.
All-wave sample waveform is carried out to noise pre-service, smoothing processing.Embodiment adopts wavelet de-noising disposal route.Wavelet de-noising has good time-frequency characteristic, selects base flexible, can process for ground the detailed information such as each edge, spike of non-stationary signal.In embodiment, call LabVIEW small echo kit, select " Advanced Signal Processing Toolkit → Wavelet Analysis → Feature Extraction → Wavelet Denoise " to call Wavelet Denoising Method function.Can select suitable de-noising level, threshold value selective rule and wavelet basis, input signal is processed.De-noising level is generally selected 4-5 level, threshold value is selected common several soft-thresholds, as without evaluation of risk threshold value (rigrsure) partially, to inspire three kinds of threshold value (heursure) and minimax threshold values (minimaxi), wavelet basis be mainly tetra-kinds of sym, db, bior, coif.According to actual effect, adopt 5 grades, rigrsure, tri-parameters of sym8 can reach comparatively desirable effect.Smoothing algorithm is in order to signal deburring, can select the method for moving average, 53 smoothing algorithms, vondrak algorithm in embodiment.Wherein vondrak algorithm smooth effect is best, effectively filtering burr, and the while can ensure again the integrality of signal well.Vondrak smoothing method is that Czech astronomer Vondrak proposes, and is applicable to uniformly-spaced be suitable for again the smoothing algorithm of unequal interval., its basic level and smooth criterion is: Q=F+ λ 2s=min, wherein F = Σ p i [ x ′ ( t i ) - x ( t i ) ] 2 , s = Σ i = 1 N - 3 [ Δ 3 x ′ ( t i ) ] 2 , X (t i) for treating smooth value.
Peak value, flex point detect for the number of definite waveform decomposition and initialized location and the peak value of nonlinear fitting, need before this to delete the wrong point of erroneous judgement, and calculated threshold are as matching constant.Implementation method: for peak value, search first order derivative=0; For rising edge flex point: first order derivative >0, second derivative=0, three order derivative <0; For negative edge flex point: first order derivative <0, second derivative=0, three order derivative >0.Nonlinear fitting is the core algorithm that waveform decomposes, and method is many, for example Levenberg-Marquardt nonlinear fitting, and LabVIEW provides this functionality controls directly to call.What existing approximating method used is Gauss curve fitting, there is certain error for laser pulse matching meeting.Actual Laser Time is distributed as: wherein t 1/2represent half-wavelength.Therefore, use as LM matching template, can obviously improve the goodness of fit, reduce error of fitting.Decompose the TOF (flight time), the strength information that obtain respective pixel by waveform.In embodiment, call " nonlinear fitting of belt restraining " function in LabVIEW, selected " LM boundary " pattern, sets matching template according to above-mentioned form, because adopting the pulse of 10ns pulse width laser, therefore set initialization condition a (3i+1)=peak-2a (3i+2), a (3i)=peak amplitude-4e -2, utilize LM iterative algorithm, output cycle index, best fit parameters.According to formula digital simulation goodness.Under suitable starting condition input, the goodness of fit can reach more than 0.999, can be used as final waveform decomposition result.To each wavelet of decomposing, its time delay is the elevation information of corresponding pixel points, and its amplitude is the strength information of corresponding pixel points.
Fig. 5 is the three-dimensional reconstruction module of Lidar imaging system.The implementation of this module is: utilize the elevation information of cloud data decomposition and the two-dimensional image of strength information generation target, as Range Profile, intensity image; Process is cut apart, is classified pixel, and filtering obtains non-ground point and ground point.Set up three-dimensional coordinate, generating mesh, to coordinate point cloud triangulation, merges degree of depth picture and intensity image, finally obtains digital elevation model.
Fig. 6 is the operating process of the laser radar 3-D imaging system based on LabVIEW, and its corresponding User Interface (GUI) as shown in Figure 7.
The first step: start after software platform, three-dimensional scenic basic parameter is carried out to initial configuration, load the three-dimensional model of landform, vegetation, aircraft, the situation elements generating is carried out to the controls such as size, position, texture, comprehensive layout, finally generates and presents Lidar three-dimensional scenic.
Second step: laser model parameter, Atmospheric models parameter, detector model parameter are set, and global parameter, for producing scenario triggered test, obtains and present all-wave signal.
The 3rd step: trigger the processing of all-wave signal, real-time resolving goes out the pixel depth information of each time cycle, and refreshes demonstration, can trigger data file.
The 4th step: trigger point cloud three-dimensional reconstruction, generate DEM and circle of equal altitudes, it can be set and present the demonstration of window shape formula, and can be to picture archiving.
The test effect that GUI platform based on Fig. 7 produces is as Fig. 8, shown in Fig. 9.Fig. 8 be three-dimensional scenic module construction in Fig. 7 Lidar l-G simulation test scene present window, Fig. 9 is three-dimensional DEM figure and the circle of equal altitudes that three-dimensional reconstruction module generates.
The above is embodiments of the invention, does not limit the present invention, all any amendments of making within the scope of thought of the present invention and principle, is equal to and replaces or invalid insertion etc., all should be included within protection scope of the present invention.

Claims (7)

1. the laser radar 3-D imaging system based on virtual instrument, is characterized in that comprising three-dimensional scenic modeling, the modeling of Lidar pilot system simulated environment, the processing of all-wave signal, four functional modules of three-dimensional reconstruction;
(1) three-dimensional scenic MBM: be written into three-dimensional model element, the basic scene project such as background colour, viewing angle control, illumination, projection mode, display model is set, the scope texture, size, position, the number that comprise landform, vegetation, aircraft that read in are controlled, the wherein corresponding lasing light emitter in the position of aircraft position, dynamic change in test, finally present airborne Lidar imaging test dynamic scene, and export final DSM surface altitude figures, in order to follow-up test;
(2) Lidar pilot system simulated environment MBM: comprise laser pulse model, propagation in atmosphere model, target interaction models, four sub-module modelings of receiving element model; Laser pulse model is simulated lasing light emitter according to character such as sharp light wavelength, pulsewidth, energy; Propagation in atmosphere model is simulated the atmospheric environment of test, and generted noise model is to laser effect; Target interaction models, to laser and detection of a target effect simulation, comprises mirror-reflection, diffuse reflection, surface reflection, speckle; Noise of detector, amplifier noise that receiving element model produces after in order to analog sensed.By to the modeling of four model orders, the systematic parameters such as difference of elevation, temperature, flight angle are set, finally export the all-wave signal of each time period target respective pixel point;
(3) all-wave signal processing module: process for the echoed signal that each time period produces, recover the pixel degree of depth, strength information.This module mainly comprises signal enhancing, denoising, level and smooth, four parts of decomposition;
(4) three-dimensional reconstruction module: this module is tested final visual three-dimensional one-tenth figure is provided for airborne laser radar; The complete degree of depth and the strength information of target providing for all-wave signal processing module, generates Range Profile, intensity image; Be divided into non-ground point and ground point two category informations by segmentation and classification; Ground point information is set up grid, and ground point is carried out to triangle division, and aggregation strength information is finally rebuild generating digital elevation model (DEM).
2. the laser radar 3-D imaging system based on virtual instrument according to claim 1, it is characterized in that the corresponding different measuring cycle of aircraft position dynamic change in (1),, cover pixel number in conjunction with hot spot and obtain pixel number and the echoed signal thereof that one-period records simultaneously laser light splitting according to APD number.
3. the laser radar 3-D imaging system based on virtual instrument according to claim 1, it is characterized in that in (3) dynamically the echoed signal search to each cycle, sampling unique point, specific aim processing, realizes the integrated function of signal enhancing, denoising, level and smooth, waveform fitting; In order to obtain higher measuring accuracy, waveform fitting module arranges matching template according to the actual waveform of laser pulse.
4. according to the laser radar detection based on virtual instrument described in claim 1 or 3 and range measurement system three-dimensional imaging analogue system, it is characterized in that reserved TCP, serial communication interface, load measured data by collector, complete the function of Lidar eyeball cloud data three-dimensional imaging.
5. according to the laser radar detection based on virtual instrument described in claim 1 or 3 and range measurement system three-dimensional imaging analogue system, the all-wave signal processing module implementation procedure that it is characterized in that imaging system is: all-wave signal is by specifying sampling rate to obtain all-wave sampled data, this signal is carried out to signal enhancing, threshold value enhancing, bandpass filtering, pulse accumulation method, noise reduction module is processed the detailed information such as each edge, spike for ground; Smoothing algorithm, to the level and smooth deburring of denoised signal, ensures the integrality of signal characteristic simultaneously.Peak value, flex point detect determines the number of waveform decomposition and initialized location and the peak value of nonlinear fitting, based on Levenberg-Marquardt Nonlinear Quasi hop algorithm, using practical laser pulse waveform as LM matching template, obviously improve the goodness of fit, reduce error of fitting.Decompose the depth information, the strength information that obtain respective pixel by waveform.
6. according to the laser radar 3-D imaging system based on virtual instrument described in claim 1 or 3, it is characterized in that the main development platform based on LabVIEW, utilize LabWindows/CVI, the Matlab of ActiveX technology, the mixed developing pattern of C++, build user interactions function of tonic chord interface.Can present the visual simulating test of complete laser radar three-dimensional imaging by carrying out main GUI, final three-dimensional one-tenth figure.Operating process is:
The first step: start after software platform, three-dimensional scenic basic parameter is carried out to initial configuration, load the three-dimensional model of landform, vegetation, aircraft, the situation elements generating is carried out to the controls such as size, position, texture, comprehensive layout, finally generates and presents Lidar three-dimensional scenic;
Second step: laser model parameter, Atmospheric models parameter, detector model parameter are set, and global parameter, for producing scenario triggered test, obtains and present all-wave signal;
The 3rd step: trigger the processing of all-wave signal, real-time resolving goes out the pixel depth information of each time cycle, and refreshes demonstration, can trigger data file;
The 4th step: trigger point cloud three-dimensional reconstruction, generate DEM and circle of equal altitudes, it is set and presents window display format, to picture archiving.
7. laser radar detection and the range measurement system three-dimensional imaging analogue system based on virtual instrument according to claim 5, the all-wave signal processing module implementation procedure that it is characterized in that imaging system is: noise reduction has good time-frequency characteristic, select base flexible, process the detailed information such as each edge, spike of non-stationary signal for ground.Peak value, flex point detect for the number of definite waveform decomposition and initialized location and the peak value of nonlinear fitting, need before this to delete the wrong point of erroneous judgement, and calculated threshold are as matching constant; Nonlinear fitting is that core algorithm that waveform decomposes adopts and comprises and provide this functionality controls directly to call by LabVIEW by Levenberg-Marquardt nonlinear fitting; Use Gauss curve fitting, there is certain error for laser pulse matching meeting; Actual Laser Time is distributed as: wherein t 1/2represent half-wave; Therefore, use as LM matching template, obviously improve the goodness of fit, reduce error of fitting.Decompose the TOF (flight time), the strength information that obtain respective pixel by waveform.
CN201410308788.2A 2014-07-01 2014-07-01 Laser radar 3-D imaging system based on virtual instrument Active CN104049259B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410308788.2A CN104049259B (en) 2014-07-01 2014-07-01 Laser radar 3-D imaging system based on virtual instrument

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410308788.2A CN104049259B (en) 2014-07-01 2014-07-01 Laser radar 3-D imaging system based on virtual instrument

Publications (2)

Publication Number Publication Date
CN104049259A true CN104049259A (en) 2014-09-17
CN104049259B CN104049259B (en) 2017-06-16

Family

ID=51502345

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410308788.2A Active CN104049259B (en) 2014-07-01 2014-07-01 Laser radar 3-D imaging system based on virtual instrument

Country Status (1)

Country Link
CN (1) CN104049259B (en)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105044731A (en) * 2015-08-31 2015-11-11 中国电子科技集团公司第十一研究所 Laser three-dimensional imaging system and imaging method
CN106143834A (en) * 2016-06-16 2016-11-23 上海电机学院 A kind of new energy marine unmanned intelligence pick-up boat
CN106530326A (en) * 2016-11-04 2017-03-22 中科宇图科技股份有限公司 Change detection method based on image texture features and DSM
CN107272447A (en) * 2017-08-07 2017-10-20 深圳市光速达机器人科技有限公司 A kind of emulation mode, simulator and robot emulation system
CN107644456A (en) * 2017-09-22 2018-01-30 南京理工大学 Sea-floor relief reconstructing system and method based on polarization characteristics of lasers
CN107784151A (en) * 2016-08-26 2018-03-09 福特全球技术公司 The physical modeling of radar and sonac
CN107817496A (en) * 2016-09-12 2018-03-20 德尔福技术有限公司 Laser radar object detection systems for automotive vehicle
CN107942315A (en) * 2017-11-24 2018-04-20 中船重工(武汉)凌久电子有限责任公司 A kind of radar background return based on satellite elevation data produces algorithm
CN109581408A (en) * 2018-12-10 2019-04-05 中国电子科技集团公司第十研究所 A kind of method and system carrying out target identification using laser complex imaging
CN109583049A (en) * 2018-11-14 2019-04-05 上海无线电设备研究所 A kind of multi-channel laser detection simulator
CN110488633A (en) * 2019-08-29 2019-11-22 中国航空工业集团公司沈阳飞机设计研究所 A kind of avionics weapon attacking test macro
CN110758925A (en) * 2019-10-29 2020-02-07 中国石油化工股份有限公司 Combine station vault storage tank sand removal system
CN110824443A (en) * 2019-04-29 2020-02-21 当家移动绿色互联网技术集团有限公司 Radar simulation method and device, storage medium and electronic equipment
CN111930161A (en) * 2020-06-12 2020-11-13 杭州奥创光子技术有限公司 Simulation control method and system for laser
CN112198495A (en) * 2020-09-04 2021-01-08 腾讯科技(深圳)有限公司 Laser radar simulation and parameter calibration method, device and equipment thereof
CN112307620A (en) * 2020-10-30 2021-02-02 幻境虚拟现实(广州)智能科技研究院有限公司 Radar virtual simulation method and device
CN112859099A (en) * 2021-02-25 2021-05-28 哈尔滨工业大学 High-precision distance measuring system based on real-time compensation of annular separator in cloud and mist environment
CN112904353A (en) * 2021-01-20 2021-06-04 南京理工大学 Laser radar distance signal simulation method and simulation signal generator
CN113357667A (en) * 2021-08-11 2021-09-07 四川大学 Laser absorption spectrum combustion diagnosis method for turbulence flame time-average characteristic
JP2022531627A (en) * 2019-07-09 2022-07-07 シーメンス インダストリー ソフトウェア ネザーランズ ベー.ヴェー. Simulation method of continuous wave rider sensor
CN115015888A (en) * 2022-05-17 2022-09-06 哈尔滨方聚科技发展有限公司 Laser three-dimensional dynamic scene simulation system
US11726190B2 (en) 2020-05-13 2023-08-15 National Instruments Corporation System for emulating an environment for testing a frequency modulated continuous wave (FMCW) detection and ranging (LiDAR) system
CN117098255A (en) * 2023-10-19 2023-11-21 南京波达电子科技有限公司 Edge calculation-based decentralization radar ad hoc network method

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111624583B (en) * 2020-07-30 2020-10-30 之江实验室 Rapid numerical simulation method for laser radar ranging considering speed factor

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201296224Y (en) * 2008-10-24 2009-08-26 张大彪 Automotive anti-collision alarm system based on a virtual instrument
CN103472443A (en) * 2013-09-24 2013-12-25 上海无线电设备研究所 Method for dynamic sea surface laser pulse wave beam echo simulation
US8681319B2 (en) * 2011-03-31 2014-03-25 Denso Wave Incorporated Laser radar for three-dimensional scanning

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN201296224Y (en) * 2008-10-24 2009-08-26 张大彪 Automotive anti-collision alarm system based on a virtual instrument
US8681319B2 (en) * 2011-03-31 2014-03-25 Denso Wave Incorporated Laser radar for three-dimensional scanning
CN103472443A (en) * 2013-09-24 2013-12-25 上海无线电设备研究所 Method for dynamic sea surface laser pulse wave beam echo simulation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨苏辉等: "基于LabVIEW的固体相干激光雷达多普勒测速的信号处理", 《北京理工大学学报》 *

Cited By (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105044731A (en) * 2015-08-31 2015-11-11 中国电子科技集团公司第十一研究所 Laser three-dimensional imaging system and imaging method
CN106143834A (en) * 2016-06-16 2016-11-23 上海电机学院 A kind of new energy marine unmanned intelligence pick-up boat
CN107784151A (en) * 2016-08-26 2018-03-09 福特全球技术公司 The physical modeling of radar and sonac
CN107784151B (en) * 2016-08-26 2023-06-30 福特全球技术公司 Physical modeling of radar and ultrasonic sensors
CN107817496B (en) * 2016-09-12 2021-05-07 安波福技术有限公司 Object detection system suitable for automatic vehicle
CN107817496A (en) * 2016-09-12 2018-03-20 德尔福技术有限公司 Laser radar object detection systems for automotive vehicle
CN106530326B (en) * 2016-11-04 2019-07-12 中科宇图科技股份有限公司 Change detecting method based on image texture feature and DSM
CN106530326A (en) * 2016-11-04 2017-03-22 中科宇图科技股份有限公司 Change detection method based on image texture features and DSM
CN107272447A (en) * 2017-08-07 2017-10-20 深圳市光速达机器人科技有限公司 A kind of emulation mode, simulator and robot emulation system
CN107644456A (en) * 2017-09-22 2018-01-30 南京理工大学 Sea-floor relief reconstructing system and method based on polarization characteristics of lasers
CN107644456B (en) * 2017-09-22 2020-12-18 南京理工大学 Laser polarization characteristic-based submarine topography reconstruction system and method
CN107942315B (en) * 2017-11-24 2019-08-27 中船重工(武汉)凌久电子有限责任公司 A kind of radar background return generation algorithm based on satellite elevation data
CN107942315A (en) * 2017-11-24 2018-04-20 中船重工(武汉)凌久电子有限责任公司 A kind of radar background return based on satellite elevation data produces algorithm
CN109583049A (en) * 2018-11-14 2019-04-05 上海无线电设备研究所 A kind of multi-channel laser detection simulator
CN109581408A (en) * 2018-12-10 2019-04-05 中国电子科技集团公司第十研究所 A kind of method and system carrying out target identification using laser complex imaging
CN109581408B (en) * 2018-12-10 2023-01-06 中国电子科技集团公司第十一研究所 Method and system for identifying target by using laser composite imaging
CN110824443A (en) * 2019-04-29 2020-02-21 当家移动绿色互联网技术集团有限公司 Radar simulation method and device, storage medium and electronic equipment
CN110824443B (en) * 2019-04-29 2020-07-31 当家移动绿色互联网技术集团有限公司 Radar simulation method and device, storage medium and electronic equipment
JP2022531627A (en) * 2019-07-09 2022-07-07 シーメンス インダストリー ソフトウェア ネザーランズ ベー.ヴェー. Simulation method of continuous wave rider sensor
US11460562B2 (en) 2019-07-09 2022-10-04 Siemens Industry Software Netherlands B.V. Methods to simulate continuous wave lidar sensors
JP7293488B2 (en) 2019-07-09 2023-06-19 シーメンス インダストリー ソフトウェア エヌヴェ How to simulate a continuous wave lidar sensor
CN110488633A (en) * 2019-08-29 2019-11-22 中国航空工业集团公司沈阳飞机设计研究所 A kind of avionics weapon attacking test macro
CN110758925A (en) * 2019-10-29 2020-02-07 中国石油化工股份有限公司 Combine station vault storage tank sand removal system
US11726190B2 (en) 2020-05-13 2023-08-15 National Instruments Corporation System for emulating an environment for testing a frequency modulated continuous wave (FMCW) detection and ranging (LiDAR) system
CN111930161B (en) * 2020-06-12 2021-10-15 杭州奥创光子技术有限公司 Simulation control method and system for laser
CN111930161A (en) * 2020-06-12 2020-11-13 杭州奥创光子技术有限公司 Simulation control method and system for laser
CN112198495A (en) * 2020-09-04 2021-01-08 腾讯科技(深圳)有限公司 Laser radar simulation and parameter calibration method, device and equipment thereof
CN112198495B (en) * 2020-09-04 2022-03-08 腾讯科技(深圳)有限公司 Laser radar simulation and parameter calibration method, device and equipment thereof
CN112307620A (en) * 2020-10-30 2021-02-02 幻境虚拟现实(广州)智能科技研究院有限公司 Radar virtual simulation method and device
CN112904353A (en) * 2021-01-20 2021-06-04 南京理工大学 Laser radar distance signal simulation method and simulation signal generator
CN112859099B (en) * 2021-02-25 2023-11-10 哈尔滨工业大学 High-precision ranging system based on real-time compensation of annular separator in cloud and fog environment
CN112859099A (en) * 2021-02-25 2021-05-28 哈尔滨工业大学 High-precision distance measuring system based on real-time compensation of annular separator in cloud and mist environment
CN113357667A (en) * 2021-08-11 2021-09-07 四川大学 Laser absorption spectrum combustion diagnosis method for turbulence flame time-average characteristic
CN113357667B (en) * 2021-08-11 2021-12-07 四川大学 Laser absorption spectrum combustion diagnosis method for turbulence flame time-average characteristic
CN115015888B (en) * 2022-05-17 2023-01-17 哈尔滨方聚科技发展有限公司 Laser three-dimensional dynamic scene simulation system
CN115015888A (en) * 2022-05-17 2022-09-06 哈尔滨方聚科技发展有限公司 Laser three-dimensional dynamic scene simulation system
CN117098255A (en) * 2023-10-19 2023-11-21 南京波达电子科技有限公司 Edge calculation-based decentralization radar ad hoc network method
CN117098255B (en) * 2023-10-19 2023-12-15 南京波达电子科技有限公司 Edge calculation-based decentralization radar ad hoc network method

Also Published As

Publication number Publication date
CN104049259B (en) 2017-06-16

Similar Documents

Publication Publication Date Title
CN104049259B (en) Laser radar 3-D imaging system based on virtual instrument
Giordan et al. The use of unmanned aerial vehicles (UAVs) for engineering geology applications
Auer et al. RaySAR-3D SAR simulator: Now open source
Yin et al. Simulation of satellite, airborne and terrestrial LiDAR with DART (II): ALS and TLS multi-pulse acquisitions, photon counting, and solar noise
Chandler et al. Structure from motion (SFM) photogrammetry vs terrestrial laser scanning
Guo et al. LiDAR principles, processing and applications in forest ecology
Kukko et al. Small-footprint laser scanning simulator for system validation, error assessment, and algorithm development
Gašparović et al. The evaluation of photogrammetry-based DSM from low-cost UAV by LiDAR-based DSM
Zhang et al. A combined deconvolution and Gaussian decomposition approach for overlapped peak position extraction from large-footprint satellite laser altimeter waveforms
Sun et al. Large-scale building height estimation from single VHR SAR image using fully convolutional network and GIS building footprints
McNicol et al. To what extent can UAV photogrammetry replicate UAV LiDAR to determine forest structure? A test in two contrasting tropical forests
Auer et al. 3D SAR simulation of urban areas based on detailed building models
Ajayi et al. Modelling 3D Topography by comparing airborne LiDAR data with Unmanned Aerial System (UAS) photogrammetry under multiple imaging conditions
Okojie et al. Relative canopy height modelling precision from UAV and ALS datasets for forest tree height estimation
Sittler et al. The potential of Lidar in assessing elements of cultural heritage hidden under forests
Castagnetti et al. 3D Reconstruction of rock paintings: a cost-effective approach based on modern photogrammetry for rapidly mapping archaeological findings
Zheng et al. Typical applications of airborne lidar technolagy in geological investigation
Steinle et al. Assessment of laser scanning technology for change detection in buildings
Devereux et al. Airborne LiDAR: instrumentation, data acquisition and handling
Filin Calibration of airborne and spaceborne laser altimeters using natural surfaces
Rodríguez-Bulnes et al. The documentation of archaeological heritage through aerial photogrammetry and UAS-based LiDAR: the case study of the Espique valley (La Peza, Granada, Spain)
Järvinen Airborne laser scanning data comparison based on roof features
Yin et al. Recent advances of modeling lidar data using dart and radiometric calibration coefficient from LVIS waveforms comparison
Novaković et al. New Approaches to Understanding the World Natural and Cultural Heritage by Using 3D Technology: UNESCO's Škocjan Caves, Slovenia
Grasso et al. The use of SLAM and UAV technology in geological field for monitoring: the case study of the Bossea Cave

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant