CN106023293A - Three-dimensional reconstruction method based on C-scan ultrasonic image - Google Patents
Three-dimensional reconstruction method based on C-scan ultrasonic image Download PDFInfo
- Publication number
- CN106023293A CN106023293A CN201610355273.7A CN201610355273A CN106023293A CN 106023293 A CN106023293 A CN 106023293A CN 201610355273 A CN201610355273 A CN 201610355273A CN 106023293 A CN106023293 A CN 106023293A
- Authority
- CN
- China
- Prior art keywords
- scan
- sample
- ultrasonoscopy
- image
- signal
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000012545 processing Methods 0.000 claims abstract description 6
- 238000002604 ultrasonography Methods 0.000 claims description 10
- 238000005070 sampling Methods 0.000 claims description 5
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 5
- 239000000284 extract Substances 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 3
- 240000007594 Oryza sativa Species 0.000 claims description 2
- 235000007164 Oryza sativa Nutrition 0.000 claims description 2
- 235000009566 rice Nutrition 0.000 claims description 2
- 241001269238 Data Species 0.000 claims 1
- 230000000877 morphologic effect Effects 0.000 claims 1
- 238000005516 engineering process Methods 0.000 abstract description 5
- 238000001514 detection method Methods 0.000 abstract description 3
- 238000009659 non-destructive testing Methods 0.000 abstract description 2
- 239000000523 sample Substances 0.000 description 20
- 239000000853 adhesive Substances 0.000 description 4
- 230000001070 adhesive effect Effects 0.000 description 4
- 230000007547 defect Effects 0.000 description 4
- 238000003384 imaging method Methods 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000009877 rendering Methods 0.000 description 2
- 241000406668 Loxodonta cyclotis Species 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 238000013332 literature search Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004449 solid propellant Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012800 visualization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/005—General purpose rendering architectures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T5/00—Image enhancement or restoration
- G06T5/77—Retouching; Inpainting; Scratch removal
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10132—Ultrasound image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20024—Filtering details
- G06T2207/20032—Median filtering
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Graphics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
A three-dimensional reconstruction method based on a C-scan ultrasonic image belongs to the field of image processing technology application study. The reconstruction method adopts an adjustable window threshold method to binarize a C-scan ultrasonic image and realizes three-dimensional reconstruction based on contour data. The reconstruction method is characterized by using the adjustable window threshold method to binarize the C-scan ultrasonic image, using a mathematical morphology method to process the binarized image, extracting a corresponding contour, so that realizing three-dimensional reconstruction based on the contour data. The method can highlight detailed information of three-dimensional reconstruction data, improves the display speed of three-dimensional reconstruction results, and has wide application in the fields of industrial non-destructive testing, medical ultrasonic image processing and medical ultrasonic equipment detection, etc.
Description
Technical field
A kind of method that the present invention relates to technical field of image processing, specifically a kind of based on C-scan ultrasonoscopy
Three-dimensional rebuilding method.
Background technology
Along with computer digital image technology is in the application of field of industry detection and popularization, ultrasonic imaging technique is increasingly becoming
In Dynamic Non-Destruction Measurement one evaluation means effectively.And combine Digital Image Processing, ultra sonic imaging, signal processing and
The ultrasound detection engineering that the ultrasonic C-scan Image processing method of computer technology has been successfully applied in Dynamic Non-Destruction Measurement
In.Owing to directly perceived and substantial amounts of material internal structure image can be provided, defect quantitative, qualitative, location is more accurate etc., super
Sound C-scan imaging is paid close attention to the most widely.Therefore, the three-dimensional reconstruction based on C-scan ultrasonoscopy of present invention research
Method has important practical value.
Through finding the literature search of prior art, at document, " adhesive structure ultrasonic C-scanning detects Ai Chunan et al.
3-D view reconstructs " (lathe and hydraulic pressure, volume 2011,39, the 16th phase, the page number: 80-83) on, elaborate to utilize ultrasound wave to examine
Examining system carries out, to solid propellant rocket adhesive structure, the image that C-scan obtains, and obtains figure by trilinear interpolation method
The volume data of pixel in Xiang, and carry out three-dimensional reconstruction based on volume data, it is thus achieved that the 3-D view of adhesive structure, more truly
Reflect the defect situation within adhesive structure.Sun Changli et al. is at document " Three dimensional
imaging based on ultrasonic linear phased array probe》(Proceedings,2014IEEE
Far East Forum on Nondestructive Evaluation/Testing:New Technology and
Application, Increasingly Perfect NDT/E, volume 2014,1, the page number: 98-101) on, elaborate to utilize line
Property phased array probe obtain image carry out artificial defect three-dimensional visualization.The method is first with the position of two dimensional image and side
To structure three-dimensional data structure, and obtain three-dimensional data point value according to the quadratic linear interpolation of adjacent two-dimensional images point, thus
Rebuild the three dimensional structure of artificial defect.
The three-dimensional rebuilding method based on C-scan image of the above field of non destructive testing all is intended to utilize interpolation method,
As trilinear interpolation, quadratic linear interpolation etc. obtain more dense data, and by the technique study described based on volume data
The three-dimensional reconstruction of ultrasonic C-scan Image.These use method based on volume data three-dimensional reconstruction to there is subject matter is that needs carry
The position of front known accurate two dimensional image and direction, and it is right to require to carry out the reconstruction of more directviewing description by interpolation increase data point
The three dimensional structure of elephant.These methods have impact on convenience based on view data three-dimensional reconstruction and flexibility ratio, and passes through interpolation
The dense data increased reduces speed during three dimensional display.
Summary of the invention
The present invention in view of the above problems, its object is to: studies the three-dimensional knot of a kind of reconstructed object that can be convenient, flexible
Structure, can comparatively fast show again three-dimensional reconstruction data based on C-scan ultrasonoscopy three-dimensional rebuilding method.
In accordance with the invention it is possible to provide convenient three-dimensional rebuilding method based on C-scan ultrasonoscopy flexibly；Can be prominent anti-
Mirror the detailed information of three-dimensional reconstruction data；Can improve and rebuild three-dimensional data display speed etc..
Described in application based on C-scan ultrasonoscopy three-dimensional rebuilding method, it is characterised in that step is as follows:
(1) the upper surface ultrasonic scanning signal data X=[x of sample is inputted1,x2,…xn-1,xn], wherein n represents all and sweeps
Retouch the sum of data.
(2) according to ultrasonic velocity v meter per second and the signal sampling time Δ t under frequency is f hertz of sample, calculate
The height (unit: rice) that sample is corresponding
Wherein vWaterRepresent ultrasonic spread speed in water under sample measuring condition.And utilize C-scan ultrasonoscopy principle
Obtaining the ultrasound image data of sample correspondence scan depths, C-scan ultrasound image data is described as
Xi=[xk,xk+1,…xk+m],(1<k<n；k+m<n)
Wherein i represents that amount of images, k represent the start point data of the i-th width image, and m represents the data total amount of the i-th width image.
(3) to C-scan ultrasound image data XiFirst carry out gaussian filtering, and take adjustable window threshold value method to be converted by image
For binary image；First, window function W is calculated1Interior signal XiMaximum
And minima
And to signal XiIt is ranked up；Then, med signal A is selected according to ranking results1, as signal XiInterior less than intermediate value
Signal A1Number reach the 50% of resultant signal and above time, select the window function W that size is little2, in re-starting sequence and selecting
Value signal A2As threshold value T；Otherwise, window function W is directly selected1Interior med signal A1As threshold value T.Finally, sliding window function
W1And W2Each width C-scan ultrasonoscopy is processed, selects to set threshold value and image is carried out binary conversion treatment, i.e. binaryzation
(u v) is result B
Wherein (u v) represents corresponding pixel points (u, gray value v) in C-scan ultrasonoscopy to I.
(4) to binary image B, (u v) carries out the opening operation in mathematical morphology and closed operation, deletes the speckle of little area
Point, and to the holes filling in closed figure, and extract corresponding profile.
(5) height that sample is corresponding, and the two-dimensional silhouette data that step (4) is extracted are calculated according to step (2), thus
(p, q, d), wherein (p, q) is two-dimensional silhouette location of pixels on image to the three-dimensional data points of acquisition sample, and d ∈ [0, h] is two
The actual height position of dimension image correspondence sample, the surface rendering method in then utilizing three-dimensional data field visualized obtains sample
Three-dimensional reconstruction result.
Detailed description of the invention
Below embodiments of the invention are elaborated: the present embodiment is carried out under premised on technical solution of the present invention
Implement, give detailed embodiment and concrete operating process, but protection scope of the present invention is not limited to following enforcement
Example.
Below in conjunction with instantiation a kind of based on C-scan ultrasonoscopy the three-dimensional rebuilding method content of the present invention done into
The detailed description of one step:
The first step: read C-scan ultrasonic signal data；
To certain imitative organization material sample, extracting its Ultrasonic C-Scan signal 900, wherein every bars has 800 samplings
Point.
Second step: calculate height and the C-scan ultrasound image data of correspondence thereof of sample；
Ultrasonic velocity 1528 meter per second of sample, water speed 1478 meter per second in ultrasonic, prolonging after sample insertion is ultrasonic
Time 0.4 microsecond late, calculating height corresponding to sample is 18 millimeters；Thus it is possible to it is ultrasonic to obtain representative cross sections survey C-scan
Image 800 width, image size is 30 × 30 pixels.
3rd step: use adjustable window threshold value method that image is converted into binary image；
First C-scan ultrasound image data being carried out gaussian filtering, arranging a width is W1(W1=21) window of pixel
It is processed by function, and maximum in calculation window function and minima are also ranked up；Then, according to ranking results,
Judge whether to become window function W2(W2=11) process；Finally, select corresponding threshold value T=0.32 to obtain C-scan to surpass
The binaryzation result of acoustic image.
4th step: binaryzation result is carried out opening operation and closed operation processes, and extract corresponding outline data, profile
Coordinate in data image be (p, q), 0≤p≤30,0≤q≤30.
5th step: according to the two-dimensional silhouette data of sample actual height and extraction, thus obtain the three-dimensional data points of sample
(p, q, h), wherein 0≤p≤30,0≤q≤30 are two-dimensional silhouette location of pixels on image, and 0≤h≤18 millimeter are X-Y scheme
Actual height position as corresponding sample；Then the field visualized middle surface rendering method of three-dimensional data is utilized to obtain the three-dimensional of sample
Reconstructed results.
Claims (4)
1. a three-dimensional rebuilding method based on C-scan ultrasonoscopy；It is characterized in that: according to ultrasonic velocity and the signal of sample
In the sampling time, calculate the height that sample is corresponding, and utilize C-scan ultrasonoscopy principle to obtain the super of sample correspondence scan depths
Acoustic image data；Utilize the adjustable window threshold value method binaryzation to C-scan ultrasonoscopy, and applied mathematics morphological method is to two
Value result processes, and extracts corresponding profile.
Method the most according to claim 1, it is characterised in that utilize adjustable window threshold value method two to C-scan ultrasonoscopy
Value, for the C-scan ultrasound image data X=[x that sample is total1,x2,…,xn-1,xn], select the data of the i-th width image
Xi=[xk,xk+1,…xk+m],(1<k<n；k+m<n)
Wherein i represents that amount of images, k represent the start point data of the i-th width image, and m represents the data total amount of the i-th width image, and n represents
The sum of all scan image datas；
Calculate window function W1Interior signal XiMaximum
And minima
And to window function W1Interior signal XiBeing ranked up, med signal is A1；
Work as XiAt window function W1Interior all signal amplitudes are than med signal A1The number that amplitude is little reach total number half and with
Time upper, reselected window function W2, W2Compare W1Little；And to window function W2Interior signal XiRe-start sequence, select window function W2
Interior med signal A2As threshold value T；Otherwise, window function W is directly selected1Interior med signal A1As threshold value T；
Sliding window function W1And W2Processing the i-th width C-scan ultrasonoscopy, threshold value T according to selecting carries out two-value to image
Change processes, and i.e. obtains the i-th width C-scan ultrasonoscopy binaryzation result Bi(u v) is
Wherein Ii(u v) represents corresponding pixel points (u, gray value v) in the i-th width C-scan ultrasonoscopy.
Method the most according to claim 1, it is characterised in that the ultrasonic velocity according to sample and the signal sampling time, meter
Calculate the height that sample is corresponding, and utilize C-scan ultrasonoscopy principle to obtain the ultrasound image data of sample correspondence scan depths；
Particularly as follows: first according to ultrasonic spread speed v in water under the conditions of sample measurement environmentWaterMeter per second, the ultrasonic velocity in sample
V meter per second and the signal Δ t second in sampling time under frequency is f hertz, calculating the height that sample is corresponding, unit is rice；
Obtain the scan depths information that sample is corresponding；Then, it is sample according to sample image in C-scan ultrasonoscopy principle
Internal cross section, it is vertical with ultrasonic wave acoustic beam；The coordinate in length and breadth of sample image, the respectively seat in length and breadth in representative sample cross section
Mark；The ultrasound image data of sample correspondence scan depths is:
Wherein h1,h2... h represents the respective heights of sample, I1(u,v),I2(u, v) ... (u v) represents the ultrasonic of respective heights to I
View data.
Method the most according to claim 1, it is characterised in that according to Mathematical Morphology Method to binaryzation result at
Reason, extracts corresponding profile, particularly as follows: first the i-th width C-scan ultrasonoscopy binaryzation result is utilized mathematical morphology side
Method processes, and obtains the outline data of the target coordinate (p in the i-th width imagei,qi)；Then, according to the i-th ultrasonic figure of width C-scan
It is (p as corresponding sample height value obtains target characteristic point three-dimensional coordinatei,qi,hi), wherein hi∈ [0, h] represents that the i-th width C sweeps
Retouch the scan depths information that ultrasonoscopy is corresponding；Obtaining, target characteristic point corresponding to sample all C-scan ultrasonoscopy is three-dimensional
After coordinate, use iso-surface patch algorithm in three-dimensional reconstruction, obtain the three-dimensional reconstruction result of sample.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610355273.7A CN106023293B (en) | 2016-05-26 | 2016-05-26 | A kind of three-dimensional rebuilding method based on C-scan ultrasound image |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610355273.7A CN106023293B (en) | 2016-05-26 | 2016-05-26 | A kind of three-dimensional rebuilding method based on C-scan ultrasound image |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106023293A true CN106023293A (en) | 2016-10-12 |
| CN106023293B CN106023293B (en) | 2018-11-30 |
Family
ID=57095070
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610355273.7A Active CN106023293B (en) | 2016-05-26 | 2016-05-26 | A kind of three-dimensional rebuilding method based on C-scan ultrasound image |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106023293B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106793447A (en) * | 2016-12-30 | 2017-05-31 | 中国科学技术大学 | A kind of fast two-dimensional Uniform Irradiation scan method |
| CN107993204A (en) * | 2017-11-28 | 2018-05-04 | 重庆大学 | A kind of MRI image reconstructing method based on image block enhancing rarefaction representation |
| CN116188459A (en) * | 2023-04-23 | 2023-05-30 | 合肥金星智控科技股份有限公司 | Line laser rapid identification method and system for belt tearing detection |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007113720A1 (en) * | 2006-03-30 | 2007-10-11 | Koninklijke Philips Electronics N.V. | Automatic cardiac band detection on breast mri |
| CN101327126A (en) * | 2008-07-23 | 2008-12-24 | 天津大学 | Method for extracting morphologic characteristic of human body bare footprint feature |
| CN104599270A (en) * | 2015-01-18 | 2015-05-06 | 北京工业大学 | Breast neoplasms ultrasonic image segmentation method based on improved level set algorithm |
| CN104657984A (en) * | 2015-01-28 | 2015-05-27 | 复旦大学 | Automatic extraction method of three-dimensional breast full-volume image regions of interest |
| CN104730148A (en) * | 2015-03-30 | 2015-06-24 | 北京科技大学 | Metal material inner inclusion three-dimensional reconstruction method based on ultrasonic testing technology |
-
2016
- 2016-05-26 CN CN201610355273.7A patent/CN106023293B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007113720A1 (en) * | 2006-03-30 | 2007-10-11 | Koninklijke Philips Electronics N.V. | Automatic cardiac band detection on breast mri |
| CN101327126A (en) * | 2008-07-23 | 2008-12-24 | 天津大学 | Method for extracting morphologic characteristic of human body bare footprint feature |
| CN104599270A (en) * | 2015-01-18 | 2015-05-06 | 北京工业大学 | Breast neoplasms ultrasonic image segmentation method based on improved level set algorithm |
| CN104657984A (en) * | 2015-01-28 | 2015-05-27 | 复旦大学 | Automatic extraction method of three-dimensional breast full-volume image regions of interest |
| CN104730148A (en) * | 2015-03-30 | 2015-06-24 | 北京科技大学 | Metal material inner inclusion three-dimensional reconstruction method based on ultrasonic testing technology |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106793447A (en) * | 2016-12-30 | 2017-05-31 | 中国科学技术大学 | A kind of fast two-dimensional Uniform Irradiation scan method |
| CN106793447B (en) * | 2016-12-30 | 2019-10-25 | 中国科学技术大学 | A kind of fast two-dimensional Uniform Irradiation scan method |
| CN107993204A (en) * | 2017-11-28 | 2018-05-04 | 重庆大学 | A kind of MRI image reconstructing method based on image block enhancing rarefaction representation |
| CN107993204B (en) * | 2017-11-28 | 2022-05-06 | 成都国一科技有限公司 | MRI image reconstruction method based on image block enhanced sparse representation |
| CN116188459A (en) * | 2023-04-23 | 2023-05-30 | 合肥金星智控科技股份有限公司 | Line laser rapid identification method and system for belt tearing detection |
| CN116188459B (en) * | 2023-04-23 | 2023-09-15 | 合肥金星智控科技股份有限公司 | Line laser rapid identification method and system for belt tearing detection |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106023293B (en) | 2018-11-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dinh et al. | Automated visualization of concrete bridge deck condition from GPR data | |
| US9207639B2 (en) | Transforming A-scan data samples into a three-dimensional space for facilitating visualization of flaws | |
| Fang et al. | Identification of two-phase flow regime using ultrasonic phased array | |
| CN101900695B (en) | System for recognizing pseudo dual-energy undersampling substance and method thereof | |
| CN102539532B (en) | Ultrasonic C scanning imaging method based on two-dimensional neighborhood synthetic aperture focusing | |
| CN101900696B (en) | Method for identifying material by means of dual-energy undersampling and system thereof | |
| CN104111285B (en) | Method for subjecting structure form of weld to imaging and device therefor | |
| McKee et al. | Volumetric imaging through a doubly-curved surface using a 2D phased array | |
| CN110363767A (en) | A kind of gridding ultrasound tomography detection method of shaft-like workpiece defect | |
| CN104535616A (en) | Window-scanning thermal imaging defect detecting and tomography method and system | |
| CN107255673A (en) | High temperature blade internal flaw three dimensional lossless detection method based on ultrasonic phase array | |
| US8805625B2 (en) | Three-dimensional visualization and analysis method and system for non-destructive examination of a rotor bore using ultrasound | |
| Posilović et al. | Generative adversarial network with object detector discriminator for enhanced defect detection on ultrasonic B-scans | |
| CN106023293A (en) | Three-dimensional reconstruction method based on C-scan ultrasonic image | |
| KR20150108872A (en) | System and method for identification, grouping and sizing of embedded flaws in rotor components using ultrasonic inspection | |
| CN107576729A (en) | Weld defect detection and quick extraction system and method based on ultrasonic phase array | |
| Guan et al. | Material damage diagnosis and characterization for turbine rotors using three-dimensional adaptive ultrasonic NDE data reconstruction techniques | |
| CN101190132B (en) | Pretreatment method and apparatus for ultrasonic imaging | |
| KR20010026857A (en) | Ultrasound imaging apparatus and method for a target separation from background | |
| Zou et al. | B-mode ultrasound imaging measurement and 3D reconstruction of submerged topography in sediment-laden flow | |
| CN114755298A (en) | Method for detecting internal cracks of action rod of turnout switch machine based on ultrasonic technology | |
| KR101131994B1 (en) | Real-time visualization system for automatically estimating ultrasonic signal in npp | |
| CN102435621B (en) | Dual-energy under-sampled matter identification method and system | |
| Misak et al. | Comparison of 2D and 3D ultrasonic pulse echo imaging techniques for structural assessment. | |
| Sun et al. | Three-dimensional reconstruction of ceramic membrane with internal defects based on ultrasound imaging technique applying triangular matrix-synthetic aperture focusing |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |