CN102539532B - Ultrasonic C scanning imaging method based on two-dimensional neighborhood synthetic aperture focusing - Google Patents

Ultrasonic C scanning imaging method based on two-dimensional neighborhood synthetic aperture focusing Download PDF

Info

Publication number
CN102539532B
CN102539532B CN 201110458080 CN201110458080A CN102539532B CN 102539532 B CN102539532 B CN 102539532B CN 201110458080 CN201110458080 CN 201110458080 CN 201110458080 A CN201110458080 A CN 201110458080A CN 102539532 B CN102539532 B CN 102539532B
Authority
CN
China
Prior art keywords
ultrasonic
point
ripple signal
focusing
value
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.)
Active
Application number
CN 201110458080
Other languages
Chinese (zh)
Other versions
CN102539532A (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.)
Central South University
Original Assignee
Central South 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 Central South University filed Critical Central South University
Priority to CN 201110458080 priority Critical patent/CN102539532B/en
Publication of CN102539532A publication Critical patent/CN102539532A/en
Application granted granted Critical
Publication of CN102539532B publication Critical patent/CN102539532B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

The invention discloses an ultrasonic C scanning imaging method based on two-dimensional neighborhood synthetic aperture focusing, which includes five steps including performing visualization for an acoustic beam spread angle of a probe, building a mathematical model for a neighborhood of a focusing treatment point, collecting ultrasonic A wave signals, synthesizing focusing treatment ultrasonic A wave signals and performing ultrasonic C scanning imaging for a defect. The technical effect of the invention lies in that through focusing the two-dimensional neighborhood taking an objective point as a center and simultaneously taking the wafer diameter, the near field length and the acoustic beam spread angle of the probe into consideration, the neighborhood of the synthetic aperture focusing treatment point is confirmed, so that the scanning result is more accurate as compared with that of the conventional SAFT (synthetic aperture focusing technique) method, the adverse effects of the acoustic beam spread angle on the ultrasonic nondestructive detection can be effectively restrained, and the accuracy of the ultrasonic C scanning imaging is improved.

Description

A kind of Ultrasonic C Scanning Image method based on two-dimentional neighborhood synthetic aperture focusing
Technical field
The present invention relates to a kind of Ultrasonic NDT formation method, particularly a kind of Ultrasonic C Scanning Image method based on two-dimentional neighborhood synthetic aperture focusing.
Background technology
Ultrasound examination can be in the internal defects situation that judges of the condition of damaged member usability and shape not, it is the important means that guarantees product quality and performance, steady production technique in the many fields of modern industry, A ripple signal is as the raw data collected in the Ultrasonic Detection process, effect in ultrasonic automatic detection is very important: it is the foundation of Ultrasonic C Scanning Image, is also the basis of scanning result subsequent analysis and processing.Visible, obtaining accurate A ripple signal is the assurance that accurately detects defect.
And there is certain acoustic beam spread angle during ultrasonic probe emission sound wave, may receive the ultrasonic echo of same impact point when the diverse location of probe in measurand, make the flaw indication collected accurate not, thereby affect the precision of defect location, quantitative analysis results.Acoustic beam spread angle and wafer diameter, frequency are inversely proportional to, but large aperture probe is made complicated and application inconvenience, and there is the acoustic energy beam serious problem that decays in the too high meeting of frequency.Therefore, only depend on the probe of using different parameters instead to suppress the impact of spread angle unrealistic.Synthetic aperture focusing technology (Synthetic Aperture Focusing Technique, SAFT) be a kind of advanced person's signal processing method, can be under the prerequisite of application small-bore probe, signal by gathering the impact point neighborhood analyzing and processing in addition, imitate a large aperture probe and carry out work, thereby improve the resolution of ultrasonic signal.
At present, SAFT mainly concentrates in the application in Ultrasonic Detection field aspect ultrasonic B scanning imagery, less for the research of member C scan image.Also there is following deficiency in existing relevant SAFT algorithm application in ultrasonic field: the process points that 1. participates in focusing on is limited to the one dimension neighborhood of impact point, belongs in fact the 1D-SAFT algorithm.As shown in Figure 1: centered by impact point, only on a line segment, consider the distribution of process points, this algorithm has only suppressed the diffuse transmission influence of spread angle in some directions, is difficult to the focusing effect that reaches desirable.Therefore be necessary the two-dimentional neighborhood centered by impact point is focused on to (2D-SAFT), as shown in Figure 2.Do not take into full account the impact of the sound field characteristics of probe when 2. synthetic focusing is processed, often by virtue of experience value determines that participation focuses on number and the distribution of the ultrasonic A ripple signal of processing, and wafer diameter, near field length, the acoustic beam spread angle angle of considering each probe decide synthetic aperture focusing processing neighborhood of a point more reasonable, effect is more desirable.
Summary of the invention
In order to suppress the adverse effect of acoustic beam spread angle to Ultrasonic C-Scan, obtain ultrasonic A wave datum accurately, the invention provides a kind ofly towards ultrasonic A wave datum, consider that process points is made rational planning for and the Ultrasonic C Scanning Image method based on two-dimentional neighborhood synthetic aperture focusing of data problem of pretreatment simultaneously.
In order to realize above-mentioned technical purpose, technical scheme of the present invention is that a kind of Ultrasonic C Scanning Image method based on two-dimentional neighborhood synthetic aperture focusing comprises the following steps:
Step 1, as shown in Figure 3, a support bar is fixed on the tank base, it is 6mm-10mm smooth metal bead that supporting bar top is fixed a diameter, adjusting is connected in the columniform ultrasonic probe position on five degree of freedom motion stand, make the axis of ultrasonic probe in bead top and ultrasonic probe vertical with surface level, then control ultrasonic probe in the situation that axis keeps vertical with surface level, do the bow font scanning in same perpendicular along the radial line of the same horizontal direction of bead above bead, gather sound wave and meet the echoed signal that the bead surface produces, and extract this signal maximum amplitude, again by amplitude being mapped to color RGB (r, g, the method of the gray-scale value that b) three parameters are equal is carried out visual to the probe sound field, obtain the ultrasonic C-scan Image of the sound field that can mean to pop one's head in, as shown in Figure 4,
Step 2, for the resulting ultrasonic C-scan Image of step 1, use half-wave method to choose the threshold value of ultrasonoscopy binaryzation, the sound field image is separated from the background of C scan image, as shown in Figure 5.Then adopt digital image processing method to extract edge and the acoustic beam center line of sound field image, calculate on this basis the acoustic beam spread angle
Figure BDA0000128030940000031
near-zone length z, suppose that the probe wafer diameter is d, and the underwater sound journey between center wafer point and impact point is s, and the disc diameter of the focusing neighborhood that impact point is required is D, sets up the mathematical model that focuses on the processing vertex neighborhood and be
Figure BDA0000128030940000032
Step 3, workpiece for measurement is carried out discrete, be divided into the capable h of g row, row and the intersection point of row are decided to be the collection position of A ripple signal, and the ultrasonic A ripple signal indication that the capable j row intersection point of i place gathers is s (i, j)={ x ij, y ij, q ij(k) | k=1,2 ..., r}, wherein suppose i, j represents respectively row and the sequence number be listed as, x ij, y ijrepresent the coordinate figure of this collection position, q ij(k) mean k sampled point amplitude in this A ripple signal sampling point sequence, r means the sampled point number;
Step 4, ultrasonic A ripple signal is carried out to two-dimentional neighborhood synthetic aperture focusing process;
Step 5, ultrasonic A ripple signal after processing for two-dimentional neighborhood synthetic focusing, the boundary wave that the waveform of first determining first peak value maximum is tested member, then finding out the waveform that wave space with respect to boundary wave is less than measured workpiece thickness is defect waves, calculate square also cumulative energy eigenvalue that obtains defect of each sampled point amplitude in defect waves, and by the palette technology, the defect energy eigenvalue is mapped to color RGB (r, g, b) the color gray-scale value that three parameters are equal, point centered by the coordinate figure of A ripple signal acquisition point, take sampling interval as a little pel of square of length of side drafting, by all A ripple signals after the defect visualization processing, obtained the ultrasonic C-scan Image of measured workpiece.
A kind of described Ultrasonic C Scanning Image method based on two-dimentional neighborhood synthetic aperture focusing, in described step 4, ultrasonic A ripple signal is carried out to the synthetic focusing processing and comprise the following steps:
Step 1, by i, j assignment is 0;
Step 2, by interim ranks m, the n assignment is 0;
Whether step 3, judgement signals collecting position s (i, j) are less than D with the distance between s (m, n), are less than D execution step 4, otherwise perform step 5;
Step 4, s (m, n) is added into to the A ripple signal list L that s (i, j) participate in to focus on processes ij;
Step 5, judge whether m equals the maximal value g of row sequence number, if unequal perform step 6, otherwise would perform step 7;
Step 3 step of step 6, line number value m returning to after adding 1 continues to carry out;
Step 7, judge whether n equals the maximal value h of row sequence number, if unequal perform step 8, otherwise would perform step 9;
Step 8, column number value m add 1, and line number value n returns to step 3 step after composing 0 continues to carry out;
Step 9, to the list L of A ripple signal ijcarry out synthetic focusing processing assignment to the new ultrasonic A ripple signal S (i, j) in the capable j row intersection point of i place;
Step 10, empty list L ij;
Step 11, judge whether line number i equals the maximal value g of row sequence number, if unequal perform step 12, otherwise would perform step 13;
Step 12, line number value i return to step 2 and carry out after adding 1;
Step 13, judge whether row sequence number j equals the maximal value h of row sequence number, if equate whole flow process finishes, otherwise perform step 14;
Step 14, row sequence number j add 1, and line number value i composes 0 rear execution 002 step.
A kind of described Ultrasonic C Scanning Image method based on two-dimentional neighborhood synthetic aperture focusing, the execution flow process of described step 9 is:
From A ripple signal list L ijin choose the benchmark that s (i, j) processes as synthetic focusing, set the sampled point sequence to this benchmark A ripple signal after the sampled point threshold value and traveled through, can, referring to Fig. 7, extract the effective analystal section [k for defect b, k e], list of hypotheses L ijin the distance of certain A ripple signals collecting position and benchmark be l, ultrasonic card sample frequency is f, the velocity of propagation of ultrasound wave in measured workpiece is c, effective analystal section of this A ripple signal is
Figure BDA0000128030940000051
by the valid interval internal sort of each signal, identical sampled point amplitude is averaged after being sued for peace again, then using mean value assignment again to the corresponding sampled point of the effective analystal section of reference point signal as amplitude, obtain this position, collection point and focus on the A ripple signal S (i, j) after processing.
Technique effect of the present invention is, carry out synthetic focusing by the two-dimentional neighborhood centered by impact point, improved the signal focus effect, wafer diameter, near field length, the acoustic beam spread angle angle of simultaneously considering each probe decides synthetic aperture focusing to process neighborhood of a point, make scanning result compare existing SAFT method more accurate, effectively improved the precision of Ultrasonic C imaging.
The accompanying drawing explanation
Fig. 1 is target neighborhood of a point schematic diagram in existing SAFT disposal route;
Fig. 2 is target neighborhood of a point schematic diagram in SAFT disposal route of the present invention;
Fig. 3 is apparatus of the present invention schematic diagram;
Fig. 4 is that the present invention utilizes bead to scan drawn probe sound field ultrasonic C-scan Image;
Fig. 5 is the probe sound field image that the present invention separates from the background of C scan image;
Fig. 6 is ultrasonic probe sound field characteristic parameter schematic diagram of the present invention;
The collection schematic diagram that Fig. 7 is the ultrasonic A ripple of the present invention signal;
Fig. 8 is the process flow diagram that the present invention's two dimension neighborhood synthetic focusing is processed;
Fig. 9 is the schematic diagram figure that synthetic focusing of the present invention is processed;
Figure 10 is test block pictorial diagram of the present invention;
Figure 11 is the flaw size figure in test block of the present invention;
Figure 12 is the handled C scan image of the present invention;
Figure 13 is the original C scan image of processing without synthetic focusing;
The C scan image of Figure 14 for adopting traditional SAFT to process.
Embodiment
Referring to Fig. 3-Figure 14, this embodiment take the test block that there is an irregularly shaped groove on the surface sweeping surface as example illustrates detection method of the present invention, the present invention controls the ultrasonic probe be installed on five degree of freedom motion stand by the motion control card of installing on control computer, ultrasonic probe is movable in tank, and ultrasonic signal being transferred to the ultrasonic card of control computer through prime amplifier, ultrasonic card transfers to control computer by the ultrasonic signal after processing and is further processed.Detecting step is as follows:
Step 1: as shown in Figure 3, by diameter on a support bar that is fixed on base, be that 6mm-10mm smooth metal bead is placed in tank, the length of support bar is generally the 25cm left and right for good, ultrasonic probe is connected on five degree of freedom motion stand, make the axis of ultrasonic probe in bead top and ultrasonic probe vertical with surface level, then control ultrasonic probe in the situation that axis keeps vertical with surface level, do the bow font scanning in same perpendicular along the radial line of the same horizontal direction of bead above bead, be after ultrasonic probe keeps moving horizontally and carries out single pass, upwards vertically move a segment distance, and then again scan along the former opposite direction moved horizontally, gather sound wave and meet the echoed signal that the bead surface produces, and extract this signal maximum amplitude, by this maximum amplitude being mapped to color RGB (r, g, the method of the gray-scale value that b) three parameters are equal is carried out visual to the probe sound field, if the maximum amplitude gathered is V max, and the span of color RGB (r, g, b) three parameters is (0,255), therefore by V maxthe linear mapping of value to color RGB (r, g, b) three parameters, i.e. V max=0 o'clock, color RGB (r, g, b) three parameters assignment simultaneously were 0, V maxduring the maximal value that can gather for ultrasonic acquisition card, color RGB (r, g, b) three parameters assignment simultaneously are 255, V maxduring for other intermediate values, compose identical value in proportion color RGB (r, g, b) three parameters, amplitude can be converted to the color gray-scale value, probe acoustic beam spread angle is contained in the probe sound field, abstract sound field is carried out to the extraction of the visual acoustic beam spread angle of being convenient to pop one's head in, and use smooth ball as the target reflecting body, the acoustic beam of sound field optional position can the alignment with the ball target normal direction, can obtain to greatest extent reflection echo and guarantee the accuracy of measuring, and cost is low, simple to operate.
Step 2: for the resulting ultrasonic C-scan Image of step 1, use half-wave method, at first travel through ultrasonic signal value maximum in the sound field image, half of this value carried out to binary conversion treatment as threshold value to it, sound field is separated from background; Then to the C image after binaryzation by from left to right, order from top to bottom scanned, then adopt digital image processing method to extract edge and the acoustic beam center line of sound field image, on the acoustic beam center line last very big ultrasonic signal put to sound source the distance be called near field length z, as shown in Figure 6.The far-field region of sound field image perpendicular to the acoustic beam centerline direction on two diverse locations intercept respectively image outline, if s 1, s 2, d 1, d 2be respectively secondary when intercepting apart from length and the intercept of sound source, pop one's head in spread angle according to formula
Figure BDA0000128030940000071
calculate.Suppose that the probe wafer diameter is d, the underwater sound journey between center wafer point and impact point is s, and the disc diameter of impact point neighborhood is D, sets up the mathematical model that focuses on the processing vertex neighborhood to be
Figure BDA0000128030940000072
Step 3, as shown in Figure 7, carry out measured workpiece discrete, is divided into the capable h row of g, and row and the intersection point of row are decided to be the collection position of A ripple signal, and the ultrasonic A ripple signal indication that the capable j row intersection point of i place gathers is s (i, j)={ x ij, y ij, q ij(k) | k=1,2 ..., r}, wherein suppose i, j represents respectively row and the sequence number be listed as, x ij, y ijrepresent the coordinate figure of this collection position, q ij(k) mean k sampled point amplitude in this A ripple signal sampling point sequence, r means sampled point length;
Step 4, according to flow process shown in Fig. 8, ultrasonic A ripple signal is carried out to the synthetic focusing processing, comprises the following steps:
Step 1, by i, j assignment is 0;
Step 2, by interim ranks m, the n assignment is 0;
Whether step 3, judgement signals collecting position s (i, j) are less than D with the distance between s (m, n), are less than D execution step 4, otherwise perform step 5;
Step 4, s (m, n) is added into to the A ripple signal list L that s (i, j) participate in to focus on processes ij;
Step 5, judge whether m equals the maximal value g of row sequence number, if unequal perform step 6, otherwise would perform step 7;
The step 3 of step 6, line number value m returning to after adding 1 continues to carry out;
Step 7, judge whether n equals the maximal value h of row sequence number, if unequal perform step 8, otherwise would perform step 9;
Step 8, column number value m add 1, and line number value n returns to step 3 after composing 0 continues to carry out;
Step 9, to the list L of A ripple signal ijcarry out synthetic focusing processing assignment to the new ultrasonic A ripple signal S (i, j) in the capable j row intersection point of i place;
Step 10, empty list L ij;
Step 11, judge whether line number i equals the maximal value g of row sequence number, if unequal perform step 12, otherwise would perform step 13;
Step 12, line number value i return to step 2 and carry out after adding 1;
Step 13, judge whether row sequence number j equals the maximal value h of row sequence number, if equate that whole flow process finishes, otherwise perform step 14;
Step 14, row sequence number j add 1, and line number value i composes 0 rear execution step 2.
Wherein the execution flow process of step 9 is:
From A ripple signal list L ijin choose the reference point signal that s (i, j) processes as synthetic focusing, set the sampled point sequence to this benchmark A ripple signal after the sampled point threshold value and traveled through, referring to Fig. 9, extract the effective analystal section [k for defect b, k e], list of hypotheses L ijin the distance of certain A ripple signals collecting position and benchmark be l, ultrasonic card sample frequency is f, the velocity of propagation of ultrasound wave in measured workpiece is c, effective analystal section of this A ripple signal is to list L ijin after all participations focus on processing signals and carry out the valid interval analysis, by the valid interval internal sort of each signal, the amplitude of identical sampled point is averaged after being sued for peace again, then using mean value assignment again to corresponding sampled point in the reference point valid interval as amplitude, obtain this position, collection point and focus on the A ripple signal S (i, j) after processing.
Step 5, ultrasonic A ripple signal after processing for two-dimentional neighborhood synthetic focusing, the waveform of first determining first peak value maximum is boundary wave, then finding out the waveform that wave space with respect to boundary wave is less than measured workpiece thickness is defect waves, calculate square also cumulative energy eigenvalue that obtains defect of each sampled point amplitude in defect waves, and energy eigenvalue is mapped to color RGB (r by the palette technology by the defect characteristic value, g, b) the equal gray-scale value of three parameters, point centered by the coordinate figure of A ripple signal acquisition point, take sampling interval as a little pel of square of length of side drafting, by all A ripple signals after the defect visualization processing, obtained the ultrasonic C-scan Image of measured workpiece.
Defect in test block is 5.6mm apart from the surperficial degree of depth, and the defect bottom is smooth, test block thickness 15mm, the major diameter 35.2mm of defect, minor axis 16.4mm, area 550.28mm 2.
The six degree of freedom CYS-1100 model Ultrasonic C-scanning System that the present embodiment adopts Shanghai Cytrix Electrical Technology Co., Ltd. to produce is tested, the ultrasonic card of selecting is UT-2001 type high-speed sampling card, probe is the water logging flat probe of the 2.5Z14SJ model of Shantou ultrasound instrument research institute's production, the wafer diameter of this probe is 14mm, near field length is 81.6mm, steel ball size in ball target method acoustic field is 4mm, by the C scanning imagery, probe acoustic beam spread angle is carried out visually, calculate on this basis the acoustic beam spread angle of this flat probe
Figure BDA0000128030940000092
be 6.8 °, therefore, setting up and focusing on the mathematical model of processing vertex neighborhood is tan3.4 °=14+0.119 of D=14+2 (s-81.6) (s-81.6).50mm*50mm data acquisition rectangular area is set, and step pitch is 0.3mm, and the underwater sound is apart from being 150mm, and sampled point length is 2000.Calculating the radius of neighbourhood is D=22.13mm, surveyed area is bent to font scanning the ultrasonic A ripple of Real-time Collection signal, after been scanned, adopt the ultrasonic A ripple of two-dimentional neighborhood of the present invention signal focus method to be processed, finally extract energy (squared magnitude and) eigenwert of its defect, by the palette technology, defect is carried out to the C scanning imagery.Wherein Figure 11 is the C scan image through the 2D-SAFT processing, and Figure 12 is the original C scan image of processing without synthetic focusing, the C scan image of Figure 14 for adopting traditional 1D-SAFT to process.Under three kinds of methods, the defect characteristic contrast is shown in Table 1.
Figure BDA0000128030940000101
Table 1

Claims (3)

1. the Ultrasonic C Scanning Image method based on two-dimentional neighborhood synthetic aperture focusing, is characterized in that, comprises the following steps:
Step 1, a support bar is fixed on the tank base, it is 6mm-10mm smooth metal bead that supporting bar top is fixed a diameter, adjusting is connected in the columniform ultrasonic probe position on five degree of freedom motion stand, make the axis of ultrasonic probe in bead top and ultrasonic probe vertical with surface level, then control ultrasonic probe in the situation that axis keeps vertical with surface level, do the bow font scanning in same perpendicular along the radial line of the same horizontal direction of bead above bead, gather a plurality of sample point sound waves and meet the echoed signal that the bead surface produces, and extract this signal maximum amplitude, be mapped to color RGB(r by the maximum amplitude by extracted again, g, the method of the gray-scale value that b) three parameters are equal is carried out visual to the probe sound field, obtain the ultrasonic C-scan Image of the sound field that can mean to pop one's head in,
Step 2, for the resulting ultrasonic C-scan Image of step 1, use half-wave method to choose the threshold value of ultrasonic scanning image binaryzation, the sound field image of probe is separated from the background of C scan image, then adopt digital image processing method to extract edge and the acoustic beam center line of sound field image, calculate on this basis the acoustic beam spread angle
Figure FDA00002885944700012
, near-zone length z, suppose that the probe wafer diameter is d, the underwater sound journey between center wafer point and impact point is s, the disc diameter of the focusing neighborhood that impact point is required is D, the mathematical model that set up to focus on the processing vertex neighborhood is
Figure FDA00002885944700011
Step 3, workpiece for measurement is carried out discrete, be divided into the capable h of g row, row and the intersection point of row are decided to be the collection position of A ripple signal, and the ultrasonic A ripple signal indication that the capable j row intersection point of i place gathers is s (i, j)={ x ij, y ij, q ij(k) | k=1,2 ..., r}, wherein suppose i, j represents respectively row and the sequence number be listed as, x ij, y ijrepresent the coordinate figure of this collection position, q ij(k) mean k sampled point amplitude in this A ripple signal sampling point sequence, r means sampled point length;
Step 4, ultrasonic A ripple signal is carried out to the synthetic focusing processing;
Step 5, ultrasonic A ripple signal after processing for two-dimentional neighborhood synthetic focusing, the waveform of first determining first peak value maximum is boundary wave, then finding out the waveform that wave space with respect to boundary wave is less than measured workpiece thickness is defect waves, calculate square also cumulative energy eigenvalue that obtains defect of each sampled point amplitude in defect waves, and energy eigenvalue is mapped to color RGB(r by the palette technology by the defect characteristic value, g, b) the equal gray-scale value of three parameters, point centered by the coordinate figure of A ripple signal acquisition point, take sampling interval as a little pel of square of length of side drafting, by all A ripple signals after the defect visualization processing, obtained the ultrasonic C-scan Image of measured workpiece.
2. a kind of Ultrasonic C Scanning Image method based on two-dimentional neighborhood synthetic aperture focusing according to claim 1, is characterized in that, in described step 4, ultrasonic A ripple signal carried out to the synthetic focusing processing and comprise the following steps:
Step 1, by i, j assignment is 0;
Step 2, by interim ranks m, the n assignment is 0;
Whether step 3, judgement signals collecting position s (i, j) are less than D with the distance between s (m, n), are less than D execution step 4, otherwise perform step 5;
Step 4, s (m, n) is added into to the A ripple signal list L that s (i, j) participate in to focus on processes ij;
Step 5, judge whether m equals the maximal value g of row sequence number, if unequal perform step 6, otherwise would perform step 7;
The step 3 of step 6, line number value m returning to after adding 1 continues to carry out;
Step 7, judge whether n equals the maximal value h of row sequence number, if unequal perform step 8, otherwise would perform step 9;
Step 8, column number value m add 1, and line number value n returns to step 3 after composing 0 continues to carry out;
Step 9, to the list L of A ripple signal ijcarry out synthetic focusing processing assignment to the new ultrasonic A ripple signal S (i, j) in the capable j row intersection point of i place;
Step 10, empty list L ij;
Step 11, judge whether line number i equals the maximal value g of row sequence number, if unequal perform step 12, otherwise would perform step 13;
Step 12, line number value i return to step 2 and carry out after adding 1;
Step 13, judge whether row sequence number j equals the maximal value h of row sequence number, if equate whole flow process finishes, otherwise perform step 14;
Step 14, row sequence number j add 1, and line number value i composes 0 rear execution step 2.
3. a kind of Ultrasonic C Scanning Image method based on two-dimentional neighborhood synthetic aperture focusing according to claim 2, is characterized in that, the execution flow process of described step 9 is:
From A ripple signal list L ijin choose the reference point signal that s (i, j) processes as synthetic focusing, set the sampled point sequence to this benchmark A ripple signal after the sampled point threshold value and traveled through, obtain effective analystal section [k of defect waves b, k e], list of hypotheses L ijin the distance of certain A ripple signals collecting position and benchmark be l, ultrasonic card sample frequency is f, the velocity of propagation of ultrasound wave in measured workpiece is c, effective analystal section of this A ripple signal is [ 2 f ( k b c / 2 f ) 2 + l 2 / c , 2 f ( k b c / 2 f ) 2 + l 2 / c - k b + k e ] , To list L ijin after all participations focus on processing signals and carry out the valid interval analysis, by the valid interval internal sort of each signal, the amplitude of identical sampled point is averaged after being sued for peace again, then using mean value assignment again to corresponding sampled point in the effective analystal section of reference point signal as amplitude, obtain this position, collection point and focus on the A ripple signal S (i, j) after processing.
CN 201110458080 2011-12-31 2011-12-31 Ultrasonic C scanning imaging method based on two-dimensional neighborhood synthetic aperture focusing Active CN102539532B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201110458080 CN102539532B (en) 2011-12-31 2011-12-31 Ultrasonic C scanning imaging method based on two-dimensional neighborhood synthetic aperture focusing

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201110458080 CN102539532B (en) 2011-12-31 2011-12-31 Ultrasonic C scanning imaging method based on two-dimensional neighborhood synthetic aperture focusing

Publications (2)

Publication Number Publication Date
CN102539532A CN102539532A (en) 2012-07-04
CN102539532B true CN102539532B (en) 2013-06-26

Family

ID=46346941

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201110458080 Active CN102539532B (en) 2011-12-31 2011-12-31 Ultrasonic C scanning imaging method based on two-dimensional neighborhood synthetic aperture focusing

Country Status (1)

Country Link
CN (1) CN102539532B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107966497A (en) * 2017-10-30 2018-04-27 航天材料及工艺研究所 A kind of manual ultrasonic visualization non-destructive testing device and method based on camera positioning

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103018333B (en) * 2012-12-07 2014-10-29 清华大学 Synthetic aperture focused ultrasonic imaging method of layered object
CN103033816B (en) * 2012-12-07 2014-06-04 清华大学 Synthetic aperture focused ultrasonic imaging implementation method based on arc scanning transition
DE102013200974A1 (en) 2013-01-22 2014-07-24 Siemens Aktiengesellschaft Method and system for hand-held ultrasonic testing of a test object
CN107064303A (en) * 2017-04-24 2017-08-18 南通友联数码技术开发有限公司 A kind of solid wheel shaft non-pulling wheel ultrasonic testing system, detection method and imaging method
CN107576729B (en) * 2017-09-15 2021-02-12 南京中车浦镇城轨车辆有限责任公司 Ultrasonic phased array-based weld defect detection and rapid extraction system and method
CN109085237A (en) * 2018-06-20 2018-12-25 德淮半导体有限公司 A kind of ultrasonic scanning device and scan method
CN110361453B (en) * 2019-07-24 2021-09-14 四川幔壳科技有限公司 Array type ectopic receiving ultrasonic synthetic aperture focusing method based on nondestructive testing
CN112881537B (en) * 2021-01-22 2023-06-20 西安增材制造国家研究院有限公司 Elliptic synthetic aperture focusing laser ultrasonic signal imaging method
CN113533526B (en) * 2021-06-02 2022-10-11 中北大学 Ultrasonic phased array interface debonding defect full-focusing C scanning imaging method and system

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6128092A (en) * 1999-07-13 2000-10-03 National Research Council Of Canada Method and system for high resolution ultrasonic imaging of small defects or anomalies.
CN1804611A (en) * 2005-12-27 2006-07-19 哈尔滨工业大学 Ultrasonic time-of-flight diffraction detection method based on synthetic aperture focusing technique
EP1757928A1 (en) * 2004-06-14 2007-02-28 Kabushiki Kaisha Toshiba Three-dimensional ultrasonic imaging device
CN101930069A (en) * 2010-05-10 2010-12-29 中国人民解放军理工大学 Virtual instrument based system and method for detecting ultrasonic imaging by synthetic aperture focusing
CN102016566A (en) * 2008-03-27 2011-04-13 住友金属工业株式会社 Ultrasonic flaw detection method and device

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3525595B2 (en) * 1995-11-24 2004-05-10 Jfeスチール株式会社 Multi-channel automatic ultrasonic flaw detection method and apparatus for rolled metal sheet

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6128092A (en) * 1999-07-13 2000-10-03 National Research Council Of Canada Method and system for high resolution ultrasonic imaging of small defects or anomalies.
EP1757928A1 (en) * 2004-06-14 2007-02-28 Kabushiki Kaisha Toshiba Three-dimensional ultrasonic imaging device
CN1804611A (en) * 2005-12-27 2006-07-19 哈尔滨工业大学 Ultrasonic time-of-flight diffraction detection method based on synthetic aperture focusing technique
CN102016566A (en) * 2008-03-27 2011-04-13 住友金属工业株式会社 Ultrasonic flaw detection method and device
CN101930069A (en) * 2010-05-10 2010-12-29 中国人民解放军理工大学 Virtual instrument based system and method for detecting ultrasonic imaging by synthetic aperture focusing

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
JP特开平9-145686A 1997.06.06
合成孔径聚焦超声成像方法研究;杜英华 等;《海洋技术》;20100630;第29卷(第2期);全文 *
杜英华 等.合成孔径聚焦超声成像方法研究.《海洋技术》.2010,第29卷(第2期),94-96.

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107966497A (en) * 2017-10-30 2018-04-27 航天材料及工艺研究所 A kind of manual ultrasonic visualization non-destructive testing device and method based on camera positioning

Also Published As

Publication number Publication date
CN102539532A (en) 2012-07-04

Similar Documents

Publication Publication Date Title
CN102539532B (en) Ultrasonic C scanning imaging method based on two-dimensional neighborhood synthetic aperture focusing
CN101871895B (en) Laser scanning imaging nondestructive inspection method for hot continuous casting blank surface defects
CN108956761B (en) Steel plate full-coverage ultrasonic detection device and method
CN101101277B (en) High-resolution welding seam supersonic image-forming damage-free detection method
CN101806777B (en) Near surface flaw quantification detection method based on ultrasonic TOFD method
CN105973988B (en) A kind of ultrasonic three-dimensional imaging detection method of lamellar defect distribution
CN107576729B (en) Ultrasonic phased array-based weld defect detection and rapid extraction system and method
CN101655464A (en) Surface defect form and depth on-line lossless detection method of high-temperature continuous casting blank
CN102393422A (en) Ultrasonic time of flight diffraction (TOFD)-based offline defect judgment method
CN104634876A (en) Method for detecting inclusions in metal material by virtue of ultrasonic scanning microscope
CN112098526A (en) Near-surface defect feature extraction method for additive product based on laser ultrasonic technology
CN104698089A (en) Ultrasonic relative time propagation technology suitable for inclined crack quantifying and imaging
CN103033153B (en) Method for scanning ultrasonic microscope and meanwhile measuring mechanical property parameter of lamina material
CN111983030A (en) Friction welding seam defect quantitative detection method and system based on ultrasonic phased array
CN111398426A (en) Measurement and compensation calibration method for full-focusing phased array three-dimensional ultrasonic field
CN105403627A (en) Method for enhancing lateral resolution of ultrasonic testing images
CN111735774A (en) Method for quantifying size of crack defect based on time-distance curve of laser ultrasonic surface wave
CN114755298A (en) Method for detecting internal cracks of action rod of turnout switch machine based on ultrasonic technology
CN109030626A (en) A kind of detection method of stainless steel composite weld
CN104040329B (en) Method and apparatus for detecting the defect inside check object
KR101131994B1 (en) Real-time visualization system for automatically estimating ultrasonic signal in npp
CN106023293A (en) Three-dimensional reconstruction method based on C-scan ultrasonic image
CN113075297B (en) Titanium alloy phased array linear array ultrasonic detection sound field model construction method
CN114689584A (en) Method suitable for calculating grain size of silicon steel ferrite complete recrystallization structure
CN110823926B (en) Method for detecting quality of complex steel plate shear wall node based on ground penetrating radar scanning

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant