Classifications

• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01S—RADIO 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
  • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
  • G01S13/88—Radar or analogous systems specially adapted for specific applications
  • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
  • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
  • G01S13/904—SAR modes
  • G01S13/9054—Stripmap mode
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01S—RADIO 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
  • G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
  • G01S13/88—Radar or analogous systems specially adapted for specific applications
  • G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
  • G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
  • G01S13/9004—SAR image acquisition techniques
• • G—PHYSICS
  • G01—MEASURING; TESTING
  • G01S—RADIO 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
  • G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
  • G01S15/88—Sonar systems specially adapted for specific applications
  • G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
  • G01S15/8902—Side-looking sonar
  • G01S15/8904—Side-looking sonar using synthetic aperture techniques

Definitions

• the present invention relates to the field of sonar or radar imaging, and the detection of objects from an antenna, sonar or radar system, in particular a synthetic aperture system.
• Knowledge of an aquatic environment represents a major challenge for a large number of applications.
• the exploration of the underwater environment plays a very important role in many fields, such as the oil industry or the deployment of civil or military ships.
• Accurate imaging of submerged objects and/or the seabed makes it possible, in particular, to inspect structures and detect objects below the surface of the water, such as shipwrecks or underwater mines, underwater cables sailors or moving objects.
• acoustic waves are well suited to underwater detection and have quickly established themselves in this field.
• Sonar systems for “SOund Navigation And Ranging” have thus become privileged means for bathymetry and the detection of objects in aquatic environments.
• the angular resolution of a sonar system is linked to two main parameters: the frequency of the acoustic wave, and the length of the antenna, made up of one or more sensors.
• FIG. 1a An example of an SAS system is shown in Figure 1a, as part of a towed linear antenna, mounted on an autonomous underwater vehicle (or AUV, for "autonomous underwater vehicle” in English) moving along a rectilinear trajectory 101
• an autonomous underwater vehicle or AUV, for "autonomous underwater vehicle” in English
• incident waves in the form of sound pulses are emitted from the linear antenna at different positions 103a, 103b, 103c, 103d along the path 101.
• backscatter phenomenon or “backscatter” in English).
• the backscattered signals can be combined to produce an image of the seabed, and in particular of the object 102 (fixed or mobile) located in the field of the antenna, for example at the bottom of the sea. More precisely, by carrying out a coherent integration of the backscattered signals taking into account the delays (or path differences) of the signals relative to each other due to the successive positions of the antenna, it is possible to simulate a so-called "synthetic" large antenna (of LSA length much greater than the LR length of the actual antenna), which would be impossible to achieve for technical and cost reasons. The resolution of the SAS system is therefore greatly improved compared to a conventional sonar system.
• object 102 it is conceptually very convenient to treat object 102 as a fictitious source radiating towards the vehicle with a fictitious speed of sound equal to the real speed divided by two (this model is known in seismic as the "exploding reflector model”) , or “exploding reflector model” in English).
• This model is known in seismic as the "exploding reflector model”
• Exploding reflector model in English
• the factor two is introduced to compensate for the difference between the round trip path in the real case and the single path in the fictitious case.
• Imaging is then similar to the reconstruction of a field of such sources. When this field comprises only one point source, located for example at 102, the corresponding image is called the point spread function, or spatial impulse response, and is sufficient to fully characterize the imaging system, by the principle of superposition of linear systems.
• FIG. 1b represents an example of a theoretical trajectory and an example of a real trajectory of an autonomous marine or submarine vehicle.
• the autonomous underwater vehicle 104 is assumed to follow a theoretical rectilinear trajectory 105 represented at the top of FIG. 1b. Due to parasitic movements (for example, yaw, pitch, yaw, etc.) related to the aquatic environment, the real trajectory 106 actually followed by the vehicle 104 can be very different from the theoretical trajectory 105, as represented at the bottom of the Figure 1b.
• the vehicle 104 can be oriented along an axis making an angle ⁇ with respect to the direction of the nominal trajectory, as shown in Figure 1b. Such an angle can develop in particular in the presence of a current perpendicular to the trajectory.
• Several methods of the prior art make it possible to effectively measure the trajectory errors, that is to say the differences between the real trajectory 106 and the theoretical trajectory 105 followed by the vehicle 104, as well as the attitudes of the vehicle (for example the angle a represented at the bottom of Figure 1 b). These methods make it possible to determine the positions of the antenna and the successive path differences in a precise manner. However, knowledge of these elements is not sufficient, it is also necessary that the methods of forming images from the collected data can integrate them effectively.
• temporal methods i.e. methods in which data is processed only in the domain time domain
• frequency methods i.e. methods in which data is processed only in the frequency domain
• Temporal methods such as the backprojection algorithm (or “backprojection” in English), also called spatial filtering by delay and summation (or “delay and sum beamforming” in English), advantageously integrate all the determined path differences, leading thus to a very precise reconstruction of the seabed and the object.
• the back-projection algorithm is based on the calculation of the sum of the signals received by the different sensors of the antenna at different positions, to which respective delays have been applied to compensate for their round-trip propagation times from the sensor to to the considered image pixel.
• backpropagation such an application of a delay to a signal to compensate for its propagation time is called “backpropagation”, “backprojection” or even “migration” (terminology borrowed from the field of seismic).
• the received signals are modeled by cylindrical or spherical waves, which are migrated towards the place where the source is supposed to be located and that we sum, to obtain a new signal whose energy we can then calculate in order to deduce the scattering power of the source (and therefore the intensity of the pixel of the image associated with this source).
• Such methods are very computationally expensive. Indeed, because of the curvature of the wave surfaces, it is not possible to migrate the waves towards an entire group of pixels at the same time, for example an entire row or an entire column of the image. It is necessary to restart the migration, and the associated interpolations, for each pixel of the image, which requires a computer having a weight and a high energy consumption.
• frequency methods such as the Omega-k (or co-k) algorithm
• Omega-k or co-k
• these methods are mainly based on the Fourier diffraction theorem, which uses monochromatic plane waves, i.e. narrow-band plane waves (more precisely, this theorem provides the Fourier transform of the scattered or backscattered field associated with an object insonified by a wave monochromatic plane).
• Fourier's diffraction theorem presupposes an "ideal" antenna, i.e. an antenna following a uniform rectilinear motion. For this reason, the frequency methods take into account the antenna deformations with difficulty (i.e.
• the invention proposes to decompose broadband signals received by a synthetic antenna sonar or radar system from a source into a plurality of broadband plane waves in the domain temporal, thanks to a broadband and far-field temporal channel formation method. Then, for each plane wave formed, a migration of this wave towards a pixel of the image where the source is supposed to be located is implemented, also in the time domain. As detailed below, these two main steps can be implemented with a reasonable computational cost. Furthermore, the use of wideband signals is well suited to sonar signals (as well as to certain radar systems), and therefore allows more accurate results than the methods of the prior art based on narrowband signals. The proposed method is therefore both accurate and inexpensive in terms of calculations.
• Antenna system means a system comprising a physical antenna, or a synthetic aperture antenna system.
• scene it is understood a portion of the space of which one seeks to reconstruct an image.
• the scene can designate a limited portion of the seabed, of which one wishes to obtain an image.
• the present invention advantageously uses broadband signals, and makes it possible to obtain more precise results than methods which presuppose narrowband signals, this assumption often not being valid in many applications such as seabed imaging. For example.
• step b typically corresponds to a channel formation step, which is a technique well known in the field of antenna processing for locating sources, during which the signals received by the various sensors are combined. so as to form new signals associated with a set of given directions.
• Lane formation consists of compensating for the path differences between the sensors in the given directions.
• synthetic antenna processing which uses signals received from several successive transmissions, it is necessary to compensate for these path differences on the outward and return paths. It is possible to approximate the round-trip travel time to an object 102 by twice the "simple" travel time between a fictitious geometric point in space, called the phase center, and said object, as if the transmission and reception took place in the same place.
• the phase center corresponds to the midpoint between the position in space occupied by the transmitter at the time of emission and that occupied by the receiving sensor at the time of reception of the echo of the object 102.
• This approximation makes it possible to simplify the disclosure of the invention but it is by no means mandatory. It is also possible to calculate exact round-trip travel times by taking into account the spatial separation between transmission and reception.
• the new signals formed correspond to plane waves, as if the object 102 were located at infinity. This is a so-called far-field method of channel formation.
• the use of such a method for a synthetic antenna is part of the originality of the invention because it is well known that synthetic imaging is performed in the near field in the vast majority of cases.
• object 102 is a few hundred meters from the sonar while the far field of the synthetic antenna is several kilometers away. It is for this reason that all the temporal methods of the prior art use spherical or cylindrical waves, that is to say so-called near-field channel formation methods, which alone make it possible to obtain in a single stage a properly focused image.
• the formation of channels in the far field is only a second step b, the focused image being produced only in the next step c.
• the decomposition into plane waves does not aim to focus an image but to reduce the computational load.
• Stage c corresponds to a migration of the plane waves thus formed towards the pixels of the image where the sources which gave rise to them are supposed to be located.
• the migrated waves are the newly formed plane waves and not the spherical or cylindrical waves received, which proves to be much more efficient in computation time because the migration can be carried out simultaneously for all the pixels of a line or a column of the image by a simple expansion of the signal, without there being any need to start the migration again for each pixel.
• the luminous intensity of a pixel can correspond to its gray level.
• Such a method advantageously makes it possible to construct an image of a scene observed by a synthetic aperture system in a rapid, precise manner and with low energy consumption.
• the antenna system may be synthetic aperture radar or synthetic aperture sonar.
• the scene may correspond to a limited portion of a seabed.
• the method may further comprise receiving data relating to a variation in movement of the system relative to a theoretical trajectory.
• the new samples of the new time signals can then further be calculated from said data. Furthermore, since the processing is done in the time and non-frequency domain, these motion variations can themselves vary rapidly over time without the computational load being affected.
• the method advantageously integrates the motion errors and their variation over time, and provides much more precise results than the frequency methods which have great difficulty in dealing with the dynamic deformations of the synthetic antenna.
• the new samples of the new time signals may be computed using a broadband, far-field time channeling method.
• the channel formation method is based on the assumption that the sources are located in the far field, that is to say far enough from the sensors for the wave front of the signals received to be considered. as blueprint.
• broadband and far-field temporal channels can be based on an iterative algorithm of the kind described in "A fast beamforming algorithm” by Kenneth M. Houston, IEEE Oceans Conference Record 1 (1994) 211- 216 or similar algorithms that are known to those skilled in the art.
• this algorithm has O(Nlog2N) computational complexity, while time channeling methods typically have O(N 2 ) computational complexity. Such an algorithm therefore makes it possible to implement step b rapidly.
• step c can comprise an implementation of at least one Z-Chirp transform, which is a transform making it possible to expand or compress a signal in an effective manner.
• the Z-Chirp transform also has a computational complexity in O(Nlog2N), thus allowing rapid implementation of step c.
• the light intensity of said each pixel may be a function of a sum of values of the selected samples.
• Another aspect of the invention relates to a device for reconstructing an image of a scene from data coming from a synthetic aperture antenna system, said image comprising a plurality of pixels.
• This device may include: an input interface for: a. receiving broadband time signal samples backscattered from reflectors in the scene for one or more sensors of the antenna system.
• the device may further comprise a circuit for implementing, for each pixel among the plurality of pixels: b. for each direction among a plurality of directions, calculating, from the samples received by said one or more sensors, new samples of a new temporal signal respectively associated with said each direction, each new signal respectively corresponding to a broadband plane wave in said each direction; vs. for each direction among the plurality of directions, selecting, among the new samples of the respective new temporal signal, samples corresponding to a round trip time of said new signal from an origin of the image to said each pixel; and D. determining a light intensity of said each pixel from the selected samples.
• the device for reconstructing an image can be integrated into the synthetic aperture antenna system, or be a device separate from the synthetic aperture antenna system.
• a computer program, implementing all or part of the method described above, installed on pre-existing equipment, is in itself advantageous, since it makes it possible to obtain an image of the scene in a precise manner and with a relatively low computational cost.
• the present invention also relates to a computer program comprising instructions for implementing the method described above, when this program is executed by a processor.
• This program can use any programming language (for example, an object language or other), and be in the form of interpretable source code, partially compiled code or fully compiled code.
• Figure 5 described in detail below, can form the flowchart of the general algorithm of such a computer program, according to one or more embodiments.
• Another aspect relates to a non-transitory recording medium of a computer-executable program, comprising a set of data representing one or more programs, said one or more programs comprising instructions for, upon execution of said one or more programs by a computer comprising a processing unit operationally coupled to memory means and to an input/output interface module, to execute all or part of the method described above.
• Figure 1a shows an example of a synthetic aperture sonar system
• Figure 1b illustrates an example of a theoretical trajectory and an actual trajectory of an autonomous underwater vehicle on which an antenna is mounted as part of a synthetic aperture sonar system
• Figure 2 shows an example of a coordinate system that can be used to describe the motion of the sonar phase center C in embodiments of the invention
• Figure 3 is an illustration of Fourier's diffraction theorem
• Figure 4 represents an example of migration of a plane wave towards a pixel of the image
• Figure 5 is a flowchart of a method for reconstructing the image of an object according to one or more embodiments
• FIG. 6 represents the samples of the signal in the Fourier plane at different stages of the migration, according to one or more embodiments
• Figures 7a and 7b represent respectively the spread functions of the point obtained by a method of the prior art on the one hand, and of an image reconstruction method according to an embodiment of the invention d 'somewhere else ;
• Figure 8 represents a schematic block diagram of an information processing device for the implementation of one or more embodiments of the invention.
• the mobile device on which the antenna used to reconstruct an image is mounted follows a theoretical uniform rectilinear motion.
• the terms "antenna” and “sonar” are used interchangeably.
• the sonar phase center denoted C and defined as the point located in the middle of the position occupied by the sonar at the time of the transmission of the acoustic wave and the position occupied by the sonar at the time of the reception of the wave backscattered by an object or by the seabed in a given direction, theoretically moves at constant speed along an axis (Ox).
• the trajectory errors presented previously and illustrated in Figure 1b are assumed to be known, for example using a method of the prior art, such as the method described in patent application FR 1457951 filed on August 25, 2014 in the name of ECA ROBOTICS.
• Figure 2 shows an example of a coordinate system that can be used to describe the motion of the sonar phase center C in embodiments of the invention.
• point S represents the source of the signal received by one of the sensors of the antenna.
• source it is understood a point of the observed scene (i.e. space insonified by the synthetic antenna) from which one wishes to obtain an image, each pixel of the image corresponding for example to the power of the field diffused by a point of this scene.
• the source can correspond to a point of an object located at the bottom of the water.
• the phase center C of the sonar is supposed, in theory, to move along the axis (Ox) at constant speed.
• the distance x between the origin O of the coordinate system and the point C along the axis (Ox) of displacement of the sonar is commonly called “azimuth” in the SAR literature in English.
• the distance r corresponds to the orthogonal distance from S to the axis (Ox) and is commonly called “slant range” or simply “range” in the SAR literature in English.
• the coordinate system shown is a cylindrical coordinate system in which the axis of the cylinder corresponds to the axis of the path (Ox). Due to the symmetry of revolution around the axis (Ox), the study of the signals received from a source S towards the phase center C of the sonar can be carried out in the two-dimensional space Oxr.
• the signals used are generally considered to be narrowband. It is recalled that a signal can be considered to be narrowband if the time taken by a plane wave to insonify the entire antenna (ie the difference between the time taken by the wave to insonify the sensor of the antenna furthest from the source and the time taken by the wave to insonify the sensor of the antenna closest to the source) is very small in front of the duration of the sound pulse (ie the duration of the ping). In this case, it is considered that all the sensors of the antenna are insonified approximately at the same time.
• this approximation is valid when the ratio between the time taken by a plane wave to insonify the whole of the antenna and the duration of the sound pulse is less than 0.1.
• the above condition is often verified in the radar domain (although in some cases it is not applicable), so methods used in SAR imagery relying on narrowband signals generally provide good results. However, it is generally not applicable in the sonar field.
• an object 301 is shown insonified by a monochromatic plane wave 302.
• the scattered field 303 can then be measured. It is noted that in the context of sonar applications, one is interested in the backscattered field, but the principle remains the same as for a diffused field.
• the Fourier diffraction theorem states that the Fourier transform of the scattered field 303 along the axis 304 corresponds to the value of the two-dimensional (2D) Fourier transform of the object 301 along an arc of circle 305 in the frequency domain.
• 2D Fourier transform of the complete object 301 it is necessary to insonify this object 301 under different directions 304 and at different frequencies, in order to obtain a family of circular arcs 305.
• An image of the object 301 can then be obtained by applying an inverse Fourier transform to the family of circular arcs 305 obtained. It is noted that this theorem relies on monochromatic plane waves 302, and therefore applies to narrowband signals.
• the narrow-band SAS image I(x, r, k) can be obtained by weighting the angular spectrum W(u, k) by B(u, k), then performing a migration (or backpropagation) of narrow-band plane waves towards each point (or pixel) P of coordinates (x, r) of the image, by applying to W(u, k) a phase shift of 2k(xu + rv), that can be associated with the propagation delay 2(xu + rv)/c corresponding to the round trip travel time of the plane wave coming from the point O and heading in the direction 0 and the pixel P.
• the broadband SAS image i(x, r) can be obtained in an elementary way by linear superposition of the components I(x, r, k) of the narrowband SAS image, i.e. by back-propagation all the components W(u, k) of the wideband signal (i.e. all the narrowband plane waves resulting from the decomposition of the wideband signal) to which a weighting has been applied by a narrowband directivity diagram B(u , k).
• the broadband SAS image i(x, r) can be obtained by decomposing the broadband signal into a plurality of monochromatic (i.e. narrowband) plane waves, then performing a weighting of these narrow band waves followed by a migration towards each point P of coordinates (x, r) of the image.
• This migration can be performed as detailed above, i.e. by applying a phase shift corresponding to the round-trip travel time from O to P in the direction 0, which is equal to 2(xu + rv)/c .
• w(u, r) be the angular broadband plane wave spectrum, defined as follows: J where F/” 1 represents the one-dimensional inverse Fourier transform associated with the dual variables r and k.
• the broadband raw data d(x, r) is defined as follows:
• the convolution product of b(u, r) and w(u, r) can be written: which means that the weighting of the plane waves at narrow band by B(u, k) corresponds in wide band to a filtering of the angular spectrum by the function b(u,r).
• the broadband SAS image i(x, r) can be obtained directly by carrying out a back-propagation of the broadband plane waves, previously filtered by b(u,r), towards each point P of coordinates (x, r) picture:
• the processing of data from the synthetic aperture system can be performed directly in the time domain via far-field broadband channel formation (i.e., considering plane waves) followed by migration (or backpropagation) of the plane waves in the Cartesian grid (i.e., towards the points/pixels of the image).
• channels or formation of beams, or "beamforming” in English
• beamforming is a method well known to those skilled in the art for extracting from a signal the components propagating in a given direction 0.
• Channel formation and migration in the time domain can be very simply implemented for a plane wave (or, equivalently, cylindrical waves considered in the far field) which propagates in the direction perpendicular to the trajectory Ox. Indeed, in this case, the wave fronts are parallel to the direction axis (Ox), and are therefore aligned with the pixel grid of the image. It is therefore sufficient to choose the sampling period in oblique distance of the signal corresponding to the plane wave equal to the pitch between the lines of the image to carry out the migration without any calculation.
• broadband raw data is received at the level of the various sensors of the antenna (step 501 of FIG. 5).
• this raw data is received in the form of an NxP matrix, where P is the number of sensors of the antenna, and N the number of samples of the source signal received per sensor. These data are time data.
• the signal received on the sensor p can be written in the form: where tp is the propagation time taken by the wave between the source and the sensor p, and ⁇ p (t) models the ambient noise on the sensor p at time t.
• This signal r p (t) can be segmented into blocks (called “snapshots” in English) of the same duration with Typically, we can consider that the signal is stationary on during each snapshot k.
• the N-tuple of values for each sensor forms a vector r p (t k ) of dimension N*1.
• the raw data received at step 501 can then correspond to the concatenation of the vectors
• new wideband time signals are determined from the data received, by applying a method of forming broadband and far-field temporal channels.
• the wavefronts can be considered plane, so that the new temporal signals correspond to plane waves.
• time channels consists in delaying the signals received by the various sensors (ie the data received at step 501 ). The sum of all the signals delayed in a given direction 0 is then performed to obtain a new signal.
• the new signal can be of the form:
• each new signal represents a broadband plane wave received in a given 0-direction.
• the samples obtained during the formation of broadband temporal channels are back-propagate towards the points/pixels of the Cartesian grid of the image, thanks to a migration method (or back-propagation) illustrated in FIG. 4: a wave is emitted in a direction 0 from the point O.
• a migration method or back-propagation
• a set of delayed samples is therefore obtained, and each of these delayed samples corresponds to a contribution of a plane wave in a given direction to a pixel P(x,r) of the picture.
• the SAS image can then be reconstructed (step 505) from these back-propagated samples.
• the intensity (i.e. the gray level) of the image of a given pixel corresponds to the sum of all the samples received for this pixel, after formation of channels (502) and backpropagation ( 503).
• data relating to the actual movement of the antenna can be received at an optional step 502.
• the trajectory errors of the antenna with respect to a theoretical trajectory can be significant and must be taken into account for the processing of the data received, in particular at the level of the step 503 of formation of channels. Note that steps 501 and 502 can be performed in this order, in reverse order, or in parallel.
• raw data (resulting from step 501 of FIG. 5) is available, and possibly movement data (resulting from step 502 of FIG. 5).
• the raw data are for example data in (x,r).
• For broadband and far-field beamforming (step 502) it is possible to use the algorithm described in "A fast beamforming algorithm" by Kenneth M. Houston, IEEE Oceans Conference Record 1 (1994) 211- 216.
• This algorithm offers a fast method of forming beams from broadband signals, significantly less expensive in terms of calculations than the conventional beamforming method (called “sum-and-delay beamforming" in English, or DSBF), in which the signals are delayed according to the calculated and then summed rate differences.
• This algorithm is an iterative algorithm, in which a DSBF is first calculated for several subgroups of sensors and in a limited number of directions. Then, at each step, the pathways from pairs of subgroups are combined and the number of directions considered is increased.
• This algorithm provides very satisfactory results, for a computational complexity of O(Nlog2N), whereas the conventional DSBF algorithm has a computational complexity of O(N 2 ). Even if this algorithm has not been developed for large antenna arrays and has not been used in this context up to now, it proves to be particularly suitable for implementing step 503. By applying this algorithm, we obtain samples in (u,r).
• CZT Z-Chirp transform
• FFT Fast Fourier Transform
• FIG. 6 represents the output signal of the channel formation step in the Fourier plane.
• the data 601 resulting from the channel formation method 503 are represented by circles and are therefore located on a polar grid represented by solid lines.
• the circular arcs correspond to the different frequencies according to Fourier's diffraction theorem and the rays to the 0 directions of plane wave propagation.
• To perform the migration of these samples 601 towards the pixels of a Cartesian grid (represented by dashes in the figure), it is necessary to perform:
• the samples 603 represented by squares are obtained.
• the SAS image is then obtained by a two-dimensional Fourier transform. Contrary to the prior art, the transition to a Cartesian grid is obtained only by 1D expansions, which are rapid, and not by complicated interpolations (such as the Stolt interpolation) which are very costly in computation time.
• FIGS. 7a and 7b represent respectively the spread functions of the point obtained by the Omega-k method on the one hand, and by an image reconstruction method according to an embodiment of the present invention on the other hand. It appears from these figures that the point spread function (or PSF, for "point spread function" in English) is much more spread for the Omega-k algorithm than for the method according to the present invention, which means that the resolution of an imaging system according to the present invention is better than the resolution of an imaging system based on the Omega-k algorithm.
• PSF point spread function
• an advantage of the method described in Figure 5 is that the broadband plane waves can be calculated directly in the time domain for a deformed antenna array (which is the case for a towed sonar mounted on an underwater vehicle). autonomous sailor), if the deformations (i.e. the trajectory errors) are known, for example by using a motion estimation method of the prior art. Thus, the method is more accurate than prior art methods which use narrowband signals.
• Another advantage is that there are fast methods (in terms of computational cost) to perform steps 503 and 504 of Figure 5.
• the image reconstruction method shown in Figure 5 has linear computational complexity, i.e. in O(Nlog2N), ie of the same order of magnitude as the Omega-k algorithm.
• the back-projection algorithm which is a temporal method integrating the exact path differences, has a quadratic complexity, ie in O(N 2 ). The invention thus proposes a method that is both precise and rapid.
• Figure 8 shows a schematic block diagram of an information processing device 800 for implementing one or more embodiments of the invention.
• the device 800 can comprise a memory 805 for storing instructions allowing the implementation of the method, the data from the backscattered signal received, and temporary data for carrying out the different steps of a method as described previously.
• the device may further comprise a circuit 804.
• This circuit may be, for example:
• processors capable of interpreting instructions in the form of a computer program, or - an electronic card whose steps of the method of the invention are described in silicon, or else
• a programmable electronic chip such as an FPGA chip (for "Field-Programmable Gate Array” in English), a SOC (for “System On Chip” in English), a GPU (for “Graphics Processing Unit” in English) , or an ASIC (for "Application Specific Integrated Circuit”).
• FPGA chip for "Field-Programmable Gate Array” in English
• SOC for "System On Chip” in English
• GPU for "Graphics Processing Unit” in English
• ASIC Application Specific Integrated Circuit
• SOCs or system on a chip are embedded systems that integrate all the components of an electronic system into a single chip.
• An ASIC is a specialized electronic circuit that groups functionalities tailored to a given application. ASICs are usually configured when manufactured and can only be simulated by the user. Programmable logic circuits of the FPGA (Field-Programmable Gate Array) type are electronic circuits reconfigurable by the user.
• FPGA Field-Programmable Gate Array
• the device 800 may include an input interface 803 for receiving data from the backscattered signal, and an output interface 806 for supplying the image of an object located at the bottom of the sea and/or data from location of such an object.
• the computer can include, to allow easy interaction with a user, a screen 801 and a keyboard 802.
• the keyboard is optional, especially in the context of a computer in the form of a touch pad, For example.
• the device 800 can be a computer, a computer network, an electronic component, or another device comprising a processor operatively coupled to a memory, as well as, according to the chosen embodiment, a data storage unit, and other associated hardware such as a network interface and a media drive for reading removable storage media and writing to such media (not shown in the figure).
• the removable storage medium can be, for example, a compact disc (CD), a digital video/versatile disc (DVD), a flash disc, a USB key, etc.
• the memory, data storage unit, or removable storage medium contains instructions which, when executed by control circuit 804, cause this control circuit 804 to perform or control the input interface 803, output interface 806, data storage in memory 805 and/or data processing parts according to one or more embodiments of the proposed method.
• the functional diagram shown in Figure 5 is a typical example of a program of which certain instructions can be carried out at the device 800. As such, Figure 5 may correspond to the flowchart of the general algorithm of a computer program within the meaning of the invention.
• the present invention is not limited to the embodiments described above by way of examples; it extends to other variants.
• the present invention is applicable in the context of conventional fixed antenna (sonar or radar) systems, comprising a physical antenna comprising a transmitter and a plurality of sensors (receivers), in which a series of pulses are emitted and in which the return signals of these pulses are received by the various sensors.
• each phase center corresponds to the middle of the emitter and of each sensor.
• the present invention can be used in many other fields, such as ultrasonic medical imaging (echography, Doppler, etc.).

Landscapes

• Engineering & Computer Science (AREA)
• Remote Sensing (AREA)
• Radar, Positioning & Navigation (AREA)
• Physics & Mathematics (AREA)
• Computer Networks & Wireless Communication (AREA)
• General Physics & Mathematics (AREA)
• Electromagnetism (AREA)
• Acoustics & Sound (AREA)
• Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
• Radar Systems Or Details Thereof (AREA)

Applications Claiming Priority (2)

Publications (1)

Family

ID=81581242

Family Applications (1)

Country Status (4)

Families Citing this family (1)

Family Cites Families (1)

• 2022
  • 2022-02-09 FR FR2201132A patent/FR3132584B1/fr active Active
• 2023
  • 2023-01-27 WO PCT/FR2023/050108 patent/WO2023152434A1/fr active Application Filing
  • 2023-01-27 EP EP23705444.0A patent/EP4476567A1/de active Pending
  • 2023-01-27 AU AU2023218589A patent/AU2023218589A1/en active Pending

Also Published As

Similar Documents

Legal Events