Small part X ray image background suppression method and device thereof

A background suppression and line image technology, applied in image enhancement, image data processing, instruments, etc., can solve the problems of loss of soft tissue information, incomplete background suppression, and uneven background intensity, achieve good suppression effect, increase display proportion, Effects that improve image quality

Active Publication Date: 2015-11-25
LIAONING KAMPO MEDICAL SYST
4 Cites 2 Cited by

AI-Extracted Technical Summary

Problems solved by technology

Usually, in order to prevent the tissue from being mistakenly judged as the background, the threshold is set relatively low. Usually, the ...
View more

Abstract

The invention discloses a small part X ray image background suppression method and a device thereof. The method comprises steps: rough background information is extracted from an original image; according to the rough background information, fine background information is acquired; according to the fine background information, a background image is synthesized; and a full background image is removed from the original image. The device comprises a rough background extraction module, a fine background extraction module, a background image synthesis module and a background removal module. The invention adopts an X ray image background extraction method for small parts, background information in an effective image region can be effectively removed, a display proportion of soft tissues and detail information in the image is enhanced, and the image quality is enhanced.

Application Domain

Image enhancement

Technology Topic

X ray imageBackground information +6

Image

  • Small part X ray image background suppression method and device thereof
  • Small part X ray image background suppression method and device thereof
  • Small part X ray image background suppression method and device thereof

Examples

  • Experimental program(1)

Example Embodiment

[0073] The present invention will be further described in detail below with reference to the drawings and embodiments.
[0074] At present, X-ray equipment produces X-rays by bombarding the anode target (usually a tungsten target) with high-speed electrons emitted from the cathode filament in the tube. The X-rays are not generated on the surface of the anode target, but by its atoms inside the anode target. And the interaction of high-speed electrons. The path that X-rays travel through when exiting the anode target varies according to their exit angle, such as figure 1 As shown, this effect is the heel effect. According to Lambert-Beer theorem I=I 0 e -μl That is, when the attenuation coefficient of the material is the same, the intensity of the rays is inversely proportional to the exponent of the path length of the rays in the material. Therefore, the intensity of the X-rays emitted from the anode target gradually decreases along the direction of the high-speed electrons. Obviously, the X-ray intensity difference caused by the heel effect will increase as the ray coverage increases. At present, the size of medical flat-panel detectors is usually 14 inches to 17 inches. When scanning small parts, the distance from the focus to the flat-panel detector is usually no more than 1 meter. It can be concluded that the angle covered by the heel effect should exceed 20 degrees. The influence on the X-ray intensity distribution is not negligible. When the X-ray emerges from the tube, its intensity change follows the inverse square law I∝D without considering the attenuation -2 , That is, the intensity of the ray is inversely proportional to the square of the distance to the focal point. The length of the filtering path that the X-ray travels through increases with the increase of the ray angle, and the attenuation of the ray is proportional to the index of the filtering path length, which leads to the difference of the ray intensity under different filtering paths. The X-rays will scatter after interacting with the scanned object, and the scattered X-rays are more randomly distributed.
[0075] During X-ray scanning, the above-mentioned factors will eventually produce an X-ray background with inconsistent intensity distribution in the image. When scanning large parts, due to the larger attenuation of the part, the proportion of the background in the image is smaller, and the coverage area is also smaller. After conventional processing, the influence of the background on the image quality can basically be eliminated Lost. However, when shooting a small part, the attenuation of the small part is very small, and sometimes the attenuation change of the part is even smaller than the change of the background intensity. In this case, the influence of the background intensity difference on the image quality cannot be ignored.
[0076] The small part X-ray image background suppression method of the present invention is mainly divided into two parts: background extraction and background suppression. The method flow is as follows figure 2 Shown
[0077] A. Extract the rough background information image from the original image;
[0078] B. Obtain a fine background information image according to the rough background information image;
[0079] C. Synthesize a background image based on the fine background information image;
[0080] D. Remove the background image from the original image to obtain an image with background suppression.
[0081] The specific algorithm flow of the present invention is as follows image 3 Shown:
[0082] Step 301. Obtain an absolute gradient image, first in the row direction according to the preset column interval c gap Compute absolute gradient image img grad (r,c):
[0083] img grad (r,c)=|img(r,c)-img(r,c+c gap )|
[0084] Among them, img represents the original image, img(r,c) represents the pixel value of the rth row and cth column in the original image.
[0085] Then based on the row gradient image according to the preset row interval r gap Calculate the final absolute gradient image img grad :
[0086] img grad =|img grad (r,c)-img grad (r+r gap ,c)|
[0087] Step 302. img the absolute gradient image grad Pixels that exceed the preset strong edge threshold grad_gt are set to 0, which means that the pixels are cleared and the image img is obtained no_eedge.
[0088] Step 303. Set the pixels in the original image img that exceed the preset strong tissue threshold tissue_gt to 0, which means that the pixels are cleared to obtain the image img no_etissue.
[0089] Step 304. Pass the image img no_eedge And image img no_etissue The set of pixels contained in it is intersected to obtain a rough background information image img rough , The image effect is like Figure 5 Shown:
[0090] img rough =img no_eedge ∩img no_etissue
[0091] Step 305. Perform row direction filtering on the rough background information image: first img rough Assign value to img tmp , Retain the original rough background information image img rough. Statistics rough background image img tmp The pixel value of the current line in, remove isolated points, and extract the starting position of the continuous segment seg_st i , End position seg_ed i And continuous segment length seg_len i.
[0092] Clear the background information segment whose length is less than the preset threshold seg_gt, i.e. img tmp Seg_st in the corresponding line i To seg_ed i The corresponding pixel value is set to 0.
[0093] Calculate the mean value of each segment in the remaining segment seg_mean i , Take the segment where the minimum mean value is located, and record it as the valley information segment seg valley. From the first background information segment to the valley information segment seg valley Derivation of background information between:
[0094] dev i = 2 X ( seg _ mean i - seg _ mean i + 1 ) seg _ st i + seg _ ed i - seg _ st i + 1 - seg _ ed i + 1
[0095] Clear derivative dev i The background information segment corresponding to greater than 0. Use the above formula to seg the valley information segment valley Go to the last paragraph of background information for derivation, and then clear dev i The background information segment corresponding to less than 0.
[0096] Among the remaining background information segments, take the last point c_ed in the previous information segment in the adjacent information segment i And its background value val_ed i , And the first point c_st in the following information segment i+1 And its background value val_st i+1 , Calculate the rate of change of adjacent segments and the allowable upper threshold:
[0097] change _ ratio = val _ ed i - val _ st i + 1 c _ ed i - c _ st i + 1
[0098] change_ratio_gt=(c_ed i -c_st i+1 )×adj_coef
[0099] Among them, adj_coef is the effective expansion coefficient. If the change rate change_ratio is greater than the allowable upper threshold change_ratio_gt, then the i-th information segment and the previous information segments are cleared; otherwise, if the change rate change_ratio is less than the allowable lower threshold -change_ratio_gt, then the i+1th information segment and The information segment that follows it.
[0100] After completing the above operations on the current row, perform the above operations on the remaining rows in turn, until all rows are covered, and then get the updated image img tmp.
[0101] Step 306. Image img tmp Perform the processing described in step 305 for each column of, until all the columns have been processed once, and then get the updated image img again tmp.
[0102] Step 307. Clear image img tmp Residual organization information in the. First calculate img rough Hist rough , Then calculate img tmp Hist tmp , Calculate the difference ratio of the two histograms:
[0103] hist dif _ ratio = hist rough - hist tmp hist rough
[0104] In hist dif_ratio Find the position tissue_lt where the tissue_ratio_lt continuously exceeds the preset threshold and the number of points that exceeds the threshold tissue_ratio_lt also exceeds the preset threshold len_lt. Clear img tmp Pixels whose mid-gray value exceeds tissue_lt.
[0105] Step 308. Clear img tmp Invalid row information in. The background information that the behavior is considered invalid when one of the following three rules is met, the specific rules are as follows:
[0106] The number of information points contained in the background information of the current row is less than the preset threshold;
[0107] The position of the starting background information segment of the current line is greater than a preset threshold;
[0108] The position of the current line ending background information segment is less than the preset threshold.
[0109] The image after the above processing is the fine background image img fine ,Such as Image 6 Shown.
[0110] Step 309. Perform row direction fitting to fine background img fine Polynomial fitting is performed on each row of information existing in, and second-order polynomial fitting is used in this embodiment, and then the fitting coefficients are used to generate target row information.
[0111] Step 310. Perform column-direction fitting, and perform row-direction fitting on the fine background image img fine Polynomial fitting is performed on each column of information in this embodiment. Second-order polynomial fitting is also used in this embodiment, and then fitting coefficients are used to generate target column information. After the above processing, an ideal background image is obtained, such as Figure 7 Shown.
[0112] Step 311. The background image is subtracted from the original image to obtain a background suppressed image, such as Figure 8 As shown in the right picture.
[0113] Such as Figure 4 As shown, the working mode of each module and unit in the device embodiment of the present invention corresponds to the method operation steps in the method embodiment, and will not be repeated here.
[0114] The modules and units described in the embodiments of the present invention may or may not be physically separated, and some or all of the modules may be selected according to actual needs to achieve the objectives of the solutions of the embodiments. According to the above-mentioned steps of the present invention, those of ordinary skill in the art can understand and implement them without creative work.

PUM

no PUM

Description & Claims & Application Information

We can also present the details of the Description, Claims and Application information to help users get a comprehensive understanding of the technical details of the patent, such as background art, summary of invention, brief description of drawings, description of embodiments, and other original content. On the other hand, users can also determine the specific scope of protection of the technology through the list of claims; as well as understand the changes in the life cycle of the technology with the presentation of the patent timeline. Login to view more.

Similar technology patents

Acne-removing composition and skin-care product containing same

InactiveCN107496247AGood antibacterialgood suppression effect
Owner:GUANGZHOU SALIAI STEMCELL SCI & TECH CO LTD

Algal bloom inhibitor, and preparation method and application thereof

ActiveCN110422919Agood suppression effectprevent algal blooms
Owner:CHENGDU ACADEMY OF ENVIRONMENTAL SCI

Preparation method of hydrate inhibitor modified on basis of polyvinyl pyrrolidone

InactiveCN107828003Ahigh puritygood suppression effect
Owner:SHAANXI HUANKE BIOLOGICAL TECH CO LTD

Method for inhibiting methane production by anaerobic methanogenesis mixed bacteria

PendingCN111087148Agood suppression effectlong-lasting inhibition
Owner:FUJIAN AGRI & FORESTRY UNIV

Application of acetyl bufalin in preparation of antitumor drugs

ActiveCN108354939Agood suppression effect
Owner:WENZHOU MEDICAL UNIV +1

Classification and recommendation of technical efficacy words

  • good suppression effect

Device and method for inhibiting vortex-induced vibration of spiral strake with rotating blade

InactiveCN103485732Alow costgood suppression effect
Owner:SOUTHWEST PETROLEUM UNIV

No-load magnetizing inrush current suppression method of YNd11 type three-phase transformer

ActiveCN107302216Aeasy accessgood suppression effect
Owner:XI AN JIAOTONG UNIV +1

Intelligent electrical valve locator main board hardware circuit

InactiveCN103343831AHighly integratedgood suppression effect
Owner:TIANJIN JINGTONG CONTROL INSTR TECH

Bacillus subtilis for preventing and treating phytophthora blight of pepper

InactiveCN104232511Aimprove securitygood suppression effect
Owner:YANGTZE UNIVERSITY

Pesticide composition for effectively preventing cucumber downy mildew

InactiveCN103749487Agood suppression effectEasy to manufacture
Owner:TIANMU COLLEGE OF ZHEJIANG A&F UNIV
Who we serve
  • R&D Engineer
  • R&D Manager
  • IP Professional
Why Eureka
  • Industry Leading Data Capabilities
  • Powerful AI technology
  • Patent DNA Extraction
Social media
Try Eureka
PatSnap group products