Ultrasonic imaging system with interpolated scan lines

Abstract: An ultrasonic diagnostic imaging system is provided which scans an image region to receive a plurality of spatially arranged lines of ultrasonic image information signals. These ultrasonic line information signals are used in an interpolater to interpolate one or more lines which are spatially interlineated between each pair of spatially adjoining received lines. These interpolated lines are produced using either the received RF or demodulated IF ultrasonic image information signals. The interpolated lines are produced prior to scan conversion, and preferably prior to nonlinear processing such as detection or log compression to reduce spatial aliasing artifacts. In one preferred embodiment the interpolater comprises a transversal filter of four taps which is responsive to received line information signals from a common range or depth. In a second embodiment the interpolater is adaptively responsive to motion to select samples for interpolation which would be spatially at the same range in the absence of motional effects.

Description: This invention relates to ultrasonic diagnostic imaging systems which scan the body by a plurality of scan lines received over an image plane or volume and, in particular, to the creation of additional image lines through interpolation prior to scan conversion. Ultrasonic diagnostic imaging systems create images of planes or three dimensional spaces within a body by transmitting ultrasonic waves into the body and receiving echoes returned from along various directional vectors or lines. As the lines of echoes are received the spatial locations from which they were reflected are either measured or computed and stored along with the echo information. The echo information is received at the frequencies of the transmitted ultrasonic waves (RF frequencies). The received RF echo information is conventionally amplified, demodulated to an intermediate frequency (IF) range, detected and further processed to provide a line of echo information at detected baseband frequencies. The lines of echo information may be received in a number of spatial formats. The lines may be arranged in a wedge-shaped sector format, in which the lines converge at a common point or region at the receiver and diverge as they extend to greater depths within the body. If the lines do not fully converge to a point at the receiver the sector format can have a trapezoidal appearance. If the receiving transducer is curved, the small side of the trapezoid can appear curved. The wedge shape can be extended over a 180° angle, giving the sector format the shape of a semicircle. Another common spatial format is the linear format, in which the lines are oriented parallel to each other over the image plane. These planar formats may be extended into a third dimension to project a volumetric image. Since the image is conventionally displayed on a video monitor it is necessary to convert the spatial format of the image plane, whatever its shape or line arrangement, to an arrangement of parallel lines which can be scanned onto a raster display. The lines of echo information are stored in a storage medium by using their spatial locational data to store the lines in the medium in their proper spatial orientations to each other. The fully assembled image plane can be read out of the storage medium in horizontal lines for scanning onto a video display monitor. Hence, the storage medium which assembles the image plane and reads it out in a video format is referred to as a scan converter, since the format of the ultrasonic scan lines in the body are converted into a format suitable for image display. As the image is assembled line by line in the scan converter storage medium, the interline spacing and its relation to the resolution of the display becomes readily apparent. Often the result will be that storage locations between adjacent lines will be empty due to the spatial spread of the lines of echo information. This is particularly the case with a sector format, in which the lines of echo information increasingly diverge in the far field. To improve the appearance of the displayed image, the storage areas between the lines are filled in with calculated interline display values. Usually the calculated values are computed through various interpolation techniques. This process of filling in the voids between (and sometimes even along) the lines of echo information is variously referred to as image smoothing or hole filling. However, by the time the lines of echo information are presented to the scan converter, a distortional effect, or aberration, has often been introduced into the image information. This aberration is due to the increase in frequency content resulting from nonlinear processing of the echo information signals, such as log compression or detection. This frequency content increase can often exceed the Nyquist limit of the spatial sampling rate, resulting in aliasing of the echo information signals. The aberration becomes especially noticeable when the ultrasonic probe is moved relative to the body being scanned. During such movement the artifact will give the image a scintillating appearance. In a paper entitled "Analysis of a Scan Conversion Algorithm for a Real-Time Sector Scanner" by Min Hwa Lee, Joo Han Kim and Song Bai Park, 8197 IEEE Transactions on Medical Imaging MI-S, No. 2, June 1996, there is described a digital scan converter in which the data points are so chosen that a scan conversion does not require any interpolation. The paper describes previous techniques wherein interpolation is required, and one such technique, wherein scanning lines are interpolated, is as described above, wherein interpolation takes place after nonlinear processing. In accordance with the present invention, there is provided an ultrasonic diagnostic imaging system as defined in Claim 1. In the drawings: FIGURE 1 is a block diagram of an ultrasonic imaging system constructed in accordance with the principles of the present invention; FIGURES 2a, 2b, 3a and 3b are graphical representations of the spectra of ultrasound signals before and after nonlinear processing; FIGURE 4 illustrates the line interpolater of FIGURE 1 in greater detail; FIGURE 5 illustrates the interpolation of scan lines between the received scan lines of a sector scanned image; FIGURE 6 illustrates an embodiment of the line interpolater of the present invention as a four tap filter; FIGURE 7 illustrates the axially adaptive interpolation of a scan line between two received scan lines of a sector scanned image to reduce motional effects; and FIGURE 8 illustrates in block diagram form the details of an adaptive line interpolater. Referring first to FIGURE 1, an ultrasonic imaging system constructed in accordance with the principles of the present invention is shown in block diagram form. An array 10 of transducer elements is activated by a beamformer 12 to transmit waves of ultrasonic energy into a body located adjacent to the transducer array 10. Ultrasonic echoes are reflected from particles and interfaces within the body and return along paths shown as lines l 1 , l 2 , l 3 , and so forth. Generally this transmission and reception is line sequential. A beam of ultrasonic energy is first sent out in the direction of line l 1 , which results in a return of echoes along the direction of line l 1 . Next, a beam of ultrasonic energy is transmitted along line l 2 resulting in a return of echoes from the direction of line l 2 . The body adjacent to the transducer 10 is scanned in this manner until the number of lines of coverage desired for the image plane area to be displayed have been gathered. The beamformer 12 receives the echo components returned following each transmission of ultrasonic energy and delays and combines the echo components as is known in the art to form a coherent sequence of echo information which returned from the direction of the scan line. As is well known, the beamformer 12 is used to focus and steer the transmitted beams of ultrasound and to combine and focus the received echo components of an array transducer to form a coherent line of echo information. The sequentially received and formed lines of echoes are then processed by a band pass filter and demodulator 14. Band pass filtering narrows and shapes the pass band of frequencies occupied by the echo signal and demodulation translates the received high frequency echo signals to an intermediate band or baseband of frequencies where they can be processed more efficiently and simply. In accordance with the principles of the present invention, the lines of received echo information are used to form spatially intermediate lines through interpolation by line interpolater 16. The line interpolater 16 can store all the lines of one complete sweep of lines (frame) over the area to be imaged, then the lines of echo information are used to calculate an interpolated line between adjacent line pair. However, to obviate the need to store a full frame of lines it is possible to compute interpolated lines as the line data is received. One way of doing this is shown in FIGURE 4. The signals from the currently received line from the bandpass filter and demodulator 14 are applied to the input of a delay line 30 with a delay of one line cycle (where a line cycle is the total of both the transmit and receive periods as well as any periods therebetween.) In this example the currently received line signal is identified as S 32 . The delay line 30 could be an analog delay line, a CCD delay line (for a sampled analog system), or a shift register (for a digital system). The currently received line signals are also applied to the input of a divide by two circuit 34. At the same time a signal from the same point (in time or distance) of the previously received line, S 31 , appears at the output of the one line delay 30 and is applied to the input of a second divide by two circuit 32. The weighted signals at the outputs of the divide by two circuits are applied to a summer 36 where they are combined to form an interpolated echo information signal X 3 , which is available for detection, as is one of the actually received lines of echoes signals, either at the output of the delay line as shown as S31, or the currently received line signal S 32 . The echo information signals S 31 and S 32 on each line are accompanied by data which provides their spatial coordinates, so that a spatial image can be properly assembled. In the same manner as the arrangement of FIGURE 4 interpolates to compute a signal value, the spatial coordinates of the two lines are used to compute an intermediate spatial location representing the location of the interpolated signal value X 3 in the image field. In the case of FIGURE 4, where the interpolated signal is computed by weighting each line value by one-half, the location of the interpolated line signal will be computed to be halfway between the two received lines. When the relative interpolation weighting is changed to other proportions, the location of the interpolated line signal will likewise be proportionally changed. The effect of the foregoing line interpolation technique is shown in FIGURE 5. In this figure an ultrasonic wave has been transmitted along line l 1 and echoes returned from along the line in the direction indicated by the arrowhead on the line. The echoes sequentially return from points S 11 , S 21 , S 31 and S 41 along the line l 1 . These points are at distances r 1 , r 2 , r 3 , and r 4 from the transducer, respectively. Since these distances are along the line they are referred to as indicators of axial direction. As echoes are received from along line 11 at points S 11 , S 21 , S 31 and S 41 they are applied to the one line delay 30. Echoes are then received from along the adjacent line, line l 2 . As these echoes are received the l 1 echoes from the same spatial distances "r" appear at the output of the delay line 30. Thus, echoes from point S 12 of the second line will be received as echoes from point S 11 appear at the output of the delay line. The two echoes are averaged (divided by two and summed) to compute the value of echoes at point X 1 of the interpolated line l i . Similarly, echoes are progressively received from along line 12 including sequential points S 22 , S 32 and S 42 . These echoes are used with the concurrently produced echoes from the same axial distances along line l 1 to compute a sequence of interpolated echoes including those at points X 2 , X 3 , and X 4 along the intermediate interpolated line l i . The interpolater 16 of FIGURE 4 is shown at the moment that signal S 32 is applied to the input of the one line delay 30 in synchronism with the production of signal S 31 at the output to produce interpolated signal X 3 . The interpolated line echo information is treated in the same manner as the echo information of the actually received lines in the formation of an image in the arrangement of FIGURE 1. The echo information is detected by a detector 18, then scaled in dynamic range by log compression circuit 20. (Two lines are shown in the signal path following the line interpolater 16 to denote the use of both received and interpolated line information.) The echo information is then locationally stored in the storage medium of a scan converter 22, enabling the sector scan to be displayed on the screen of a video display 24. The manner in which this interpolated line technique diminishes the disturbing aberration of scintillation is explained with reference to the graphs of FIGURES 2 and 3. Curve 40 represents the extent of the baseband spectrum of the received RF ultrasonic echo information with a line density in the lateral (azimuth) direction of an ultrasound image which spatially samples the subject at a spatial sampling frequency of f s . In FIGURE 2a this spatial sampling frequency f s has been chosen to be twice the highest frequency of the information spectrum, f o , in satisfaction of the Nyquist criterion. The choice of the spatial sampling frequency f s for the spectrum extending to frequency f o will result in a continuum of nonoverlapping spectra 40', 40'', and 40''' around multiples of the sampling frequency as shown in FIGURE 2a. As long as the signals of these spectra are not processed nonlinearly, no scintillation artifacts will develop. Suppose that the RF signals undergo detection before being displayed, a nonlinear operation. The detection process will create additional spectral components, extending the spectrum 40 to the greater spectrum illustrated by 40D in FIGURE 2b. With the same spatial sampling frequency f s the greater spectra resulting from sampling will now overlap, as indicted by the cross-hatched areas in FIGURE 2b. The spectral components of one spectrum will vectorially combine with the spectral components of the overlapping spectrum to form new, randomized vectorial spectral components. These randomized spectral components recur at each cross-hatched overlap area and manifest themselves in an image as scintillation effects as the object field moves in relation to the transducer 10. The creation of additional scan lines through interpolation increases the line density and effectively increases the spatial sampling frequency to a higher value f s ' as shown in FIGURE 3a. The same spectra 40, 40', 40'' and 40''' are obtained, but the higher spatial sampling frequency f s ' creates separations between the spectra as shown. Now if the signal in the spectra are subjected to detection, the spectra will be extended to 40D, 40D', 40D'', and 40D''' as shown in FIGURE 3b, but the inter-spectra spacing resulting from the higher spatial sampling frequency enables the Nyquist criterion for the spectra to continue to be satisfied. The signals are processed nonlinearly by detection, but since there is no spectra overlap and creation of randomized signals, no scintillation artifacts are produced. FIGURE 6 illustrates a preferred line interpolater for a sampled data or digital ultrasound embodiment. This line interpolater is a four tap linear filter. The linear filter can be constructed as either an infinite impulse response (IIR) filter or as a finite impulse response (FIR or transversal) filter, and in the embodiment of FIGURE 6 a transversal filter is shown. At least three lines of received ultrasonic signals are stored, enabling the transversal filter to operate upon samples of four lines at one time. The samples from the four lines, such as lines l 1 , l 2 , l 3 and l 4 of FIGURE 5 are applied sequentially to the input of the filter and shifted or transferred in until they are separated by respective line delay times Tau by delay elements 51, 52, and 53. The example of FIGURE 6 shows the samples from range r 1 of FIGURE 5 separated by delay elements 51-53. The samples are multiplied by filter coefficients C 1 -C 4 in multipliers 60-63 and the results are summed in summer 70 to produce the interpolated result Y 1 of a line which, in the preferred embodiment, is spatially centered in the middle of the four lines l 1 -l 4 . Other spatial locations for the "Y" line may also be obtained through choice of the filter coefficients. There are a variety of ways to operate the transversal filter. If the full field of ultrasonic signals is received and stored, the filter can traverse the entire lateral r 1 range by sequentially shifting in successive samples. For instance, after Y 1 is computed for line l i , the sample S 15 from the fifth line l 5 (not shown) is shifted in to compute the value for the interpolated line between lines 13 and l 4 at range r 1 . After all interpolated values at range r 1 have been computed the samples at the next range r 2 are applied to the transversal filter to compute the next value in each interpolated line, including value Y 2 in the "Y" line. Alternatively, if only four lines are available at one time, the filter would compute Y 1 , then Y 2 , then Y 3 , and so forth until all values for the "Y" line have been computed in axial succession. It is also possible to compute a number of interpolated lines between every pair of received lines. For instance, after computing value Y 1 the coefficients C 1 -C 4 of the filter could be changed, while retaining the same values at the delay element outputs, to compute a value between Y 1 and S 12 . Still keeping the sample values at the delay element outputs and changing the filter coefficients again, an interpolated value is computed between S 13 and Y 1 . Thus, three interpolated lines are computed in this manner between every pair of received lines, resulting in a four-fold increase in line density. A higher speed technique to achieve this result is to operate three filters in parallel. Each filter is structurally identical to the other two, differing only in filter coefficients used. Once the four samples from the four lines are available for processing, the three parallel filters would compute three interpolated line values simultaneously. For a system employing the transversal filter of FIGURE 6 for the line interpolater 16 of FIGURE 1, the bandpass filter of the bandpass filter and demodulator 14 could also be configured as a transversal filter, but operating on samples in the axial direction. For example, a twelve tap filter could be employed, through which the samples of the entire line are passed in sequence to filter the line data. While the line interpolater 16 can be located either before or after the bandpass filter and demodulator 14, in this instance it would be operationally simpler to have the line interpolater 16 follow the bandpass filter. This would obviate the need to run the samples of the interpolated lines through the bandpass filter, since the interpolater would be operating upon bandpass filtered data. The use of a line interpolater in accordance with the present invention can introduce erroneous results into the interpolated line values if the subject is moving relative to the transducer while the line data is acquired. Motion can result in the use of shifted adjacent samples in the computation of the interpolated lines. FIGURES 7 and 8 illustrate an adaptive system which compensates for the effects of motion in the axial direction. In FIGURE 8 the samples from consecutive received lines are transferred through a delay line 80 which is one line cycle plus one sample period Tau in length. For the sample sequence shown for lines l 1 and l 2 in FIGURE 7, sample S b2 of line l 2 is applied to the input of the delay line 80 at the same time as sample S a1 is produced at the output of the delay line. The delay line has two taps, separated from the output by one sample period and two sample periods, respective, at which samples S b1 and S c1 are produced at the time that sample S a1 appears at the output of the delay line. A number of samples from the two lines l 1 and 12, taken from the input and the first tap of the delay line, respectively, are applied to a correlator 82. The correlator performs a cross-correlation of range aligned data samples of the two lines l 1 and l 2 to detect the condition of relative motion between the two lines. This cross correlation is performed in the conventional manner by sequentially shifting sample sequences from the two lines relative to each other, multiplying aligned samples after each shift, and summing the products to produce a correlation factor. The value and direction of shift for which the correlation factor is at a maximum indicates the amount and direction of motion that has occurred in the period between the acquisition of the two lines l 1 and l 2 . The peak of the correlation factor is then used as the control input of a selector or multiplexer 84 to select the sample at the input of the selector which would, in the absence of motion, be range-aligned with the sample at the input of the delay line 80. Thus, when sample S b2 of line l 2 is at the input of the delay line and there is no motion when the line samples are acquired, there would be a high degree of correlation between the two lines when no relative shift occurred, indicating no motion. The selector 84 would then select sample S b1 to be used in calculating an interpolated value X b with sample S b2 at the output of summer 36. As shown in FIGURE 7, these two samples are at the same range r 4 in this example. But when there is motion away from the receiver (in a direction opposite that of the arrowheads in FIGURE 7), the correlator 82 will detect the motion and its direction. In the example the selector 84 would select sample S a1 to be used in interpolation with sample S b2 . An interpolated value X a is then computed at the output of summer 36 using these two sample values. By virtue of the motion, the value X a is seen to be at a half range increment between ranges r 3 and r 4 in FIGURE 7. In a similar manner motion toward the receiver would be detected by the correlator 82 and sample S c1 would be selected to be used with sample S b2 to interpolate value X c , again at a range half increment due to motion. The interpolated values may be at fractional range increments as the foregoing example illustrates. To bring the samples back into alignment with the same ranges as the received lines, the interpolated line values can be processed through an axial transversal filter with coefficients chosen to compute an interpolated value at each whole range increment along the interpolated line. This alignment would be useful if further line filtering were to be done using a multitap filter, for instance, or if the scan converter requires sample data points to be spatially organized in a uniform grid or pattern. While the embodiment of FIGURE 8 adapts the interpolater to account for axial motion, it will be apparent that an adaptive technique can be employed to account for motion in the lateral direction also. A correlation can be performed laterally across the aperture using the signal values at two adjacent r distances. If lateral motion is detected, the values of affected signal samples can be adjusted by a weighting or interpolative technique in consideration of the values of adjacent signal samples. Thus it is possible to employ the interpolative techniques of the present invention with a high degree of immunity to motional effects. One skilled in the art will recognize that techniques other than correlation may be used to detect relative motion to be used for interpolation adjustment. For instance, EP-A-0 487 339 describes a transducer assembly which includes accelerometers used for detecting motion of the transducer. U.S. Patent 5,127,409 describes a Doppler technique for sensing transducer motion. These and other techniques may be used to provide a selection signal for control of an interpolater in accordance with the principles of the embodiment of FIGURE 8. While the increase in line density achieved by line interpolation is employed in the foregoing examples to reduce scintillation artifact effects, advantage can be taken of the line density for other purposes. For instance, fewer received lines need to be acquired to create a frame of significantly greater line density when additional lines are interpolated. The time required to assemble all of the lines necessary to create a full image frame is reduced by the difference between the time that would have been required to transmit ultrasonic waves and receive echoes for the additional lines, and the time actually utilized to interpolate the lines. This results in a decrease in the amount of time required to form a complete frame and a consequent increase in the frame rate. Thus, line interpolation can be used where an increase in the image frame rate is useful or required. Stated another way, when the frame rate of realtime imaging drops to unacceptably low levels, the replacement of scanned lines with more quickly formed interpolated lines can boost the frame rate back to acceptable levels. Alternatively, a frame may be computed in a relatively short time through line interpolation, and the time savings used to perform other image enhancements, such as frame averaging, speckle reduction or elimination of other artifacts. The principles of the present invention may be used to improve Doppler scanning as well as B mode imaging. In two dimensional Doppler scanning for flow imaging, an ensemble of lines are acquired at different times in the same direction. The data from the lines of the ensemble are used to compute Doppler velocities of flow at points along the direction of the ensemble, and a plurality of spatially adjacent ensembles can be used to depict flow throughout an entire image field. The time required to compile flow data for a full image may be reduced by interpolating line ensembles intermediate to actually scanned ensembles. For example, assume that each spatial line location l 1 -l 4 in FIGURE 5 comprised an ensemble of eight lines acquired over time. Interlineated interpolated ensembles at the "X" and "Y" line locations can be formed by interpolating the first lines (in time) of adjacent scanned ensembles to compute the first scan line of the intermediate interpolated ensemble. The second lines (in time) of the scanned ensembles are then used to compute the second line of the interpolated ensemble, and so forth until an entire interpolated ensemble is formed. The interpolated ensembles are Doppler processed in the same manner as the scanned ensembles to depict flow in the image field.

Claims: An ultrasonic diagnostic imaging system comprising: means (10,12) for transmitting ultrasonic waves into a region of a subject to be imaged; means (10,12) for receiving ultrasonic information signals from along a plurality of spatial lines (l 1 ...l 4 ) of said region; interpolation means (16, Fig. 6) for utilizing said received ultrasonic line information signals to produce a plurality of interpolated lines spatially interlineated between said spatial lines (l 1 ...l 4 ); and scan conversion means (22) for utilizing said received signals and said interpolated lines to form a spatially arranged ultrasonic image of said region of said subject; characterised in that said interpolation means (16, Fig. 6) is arranged to produce said interpolated lines prior to nonlinear processing of said ultrasonic line information signals, thereby to provide a system with reduced susceptibility to spatial aliasing. The ultrasonic diagnostic imaging system of Claim 1, wherein said interpolation means (16, Fig. 6) comprises a filter responsive to a plurality of laterally aligned signals of said received ultrasonic line information signals. The ultrasonic diagnostic imaging system of Claim 2, wherein said filter is arranged to compute a plurality of interpolated lines spatially between each pair of received spatial lines (l 1 ...l 4 ). The ultrasonic diagnostic imaging system of Claim 2 or Claim 3, wherein said filter comprises a linear filter. The ultrasonic diagnostic imaging system of any one of Claims 2 to 4, wherein said filter comprises a transversal filter. The ultrasonic diagnostic imaging system of any one of Claims 2 to 5, wherein said filter has four or more taps. The ultrasonic diagnostic imaging system of any preceding claim, further comprising means (18,20) for processing said received signals and said interpolated lines nonlinearly, said scan conversion means (22) being arranged to form said image on the basis of the nonlinearly processed signals. The ultrasonic diagnostic imaging system of Claim 7, wherein said nonlinearly processing means (18,20) comprises a detector (18). The ultrasonic diagnostic imaging system of Claim 8, wherein said nonlinearly processing means (18,20) further comprises a compression circuit (20). The ultrasonic diagnostic imaging system of Claim 9, wherein said compression circuit (20) is arranged to rescale the amplitudes of time domain signals. The ultrasonic diagnostic imaging system of Claim 9 or Claim 10, wherein said compression circuit (20) comprises a log compression circuit. The ultrasonic diagnostic imaging system of any preceding claim, wherein said interpolating means (16, Fig. 6) is arranged to utilize received RF ultrasonic line information signals. The ultrasonic diagnostic imaging system of any preceding claim, further comprising demodulation means (14), responsive to said received ultrasonic line information signals, for producing intermediate frequency ultrasonic line information signals; and wherein said interpolating means (16, Fig. 6) is arranged to utilize said intermediate frequency ultrasonic line information signals to produce said interpolated lines. The ultrasonic diagnostic imaging system of any one of Claims 1 to 12, further comprising demodulation means (14), responsive to said received ultrasonic line information signals, for producing baseband ultrasonic line information signals; and wherein said interpolating means (16, Fig. 6) is arranged to utilize said baseband ultrasonic line information signals to produce said interpolated lines. The ultrasonic diagnostic imaging system of any preceding claim, wherein said receiving means (10,12) is arranged for receiving ensembles of said ultrasonic information signals at different times from a common direction; said interpolating means (16, Fig. 6) being arranged for utilizing the received ultrasonic line information signals from two or more ensembles to produce one or more spatially interlineated interpolated ensembles of corresponding lines; said system further comprising: means for producing information concerning said subject through Doppler processing of said received and interpolated ultrasonic line ensembles. The ultrasonic diagnostic imaging system of Claim 15, wherein said receiving means (10,12) further comprises means for receiving a plurality of ensembles of lines, each of which comprises lines received in a given time sequence; and wherein said interpolating means comprises means for producing an ensemble of interpolated lines of a time sequence corresponding to said given time sequence. The ultrasonic diagnostic imaging system of Claim 15, wherein said interpolating means (16, Fig. 6) further comprises means for producing one of said lines of said interpolated ensemble having a given position in said time sequence by using received lines having corresponding positions in the time sequences of said received line ensembles.

Inventors: COOLEY CLIFFORD R, ENTREKIN ROBERT R, HWANG JUIN-JET, ROBINSON BRENT STEPHEN, COOLEY, CLIFFORD R., ENTREKIN, ROBERT R., HWANG, JUIN-JET, ROBINSON, BRENT STEPHEN
Applicants: ADVANCED TECH LAB, ADVANCED TECHNOLOGY LABORATORIES, INC.
Application year: 1994
Main classification: G01S7
IPC codes: G01S7, G01S15