Collection, storage and calculation rapid imaging method for bladder B-ultrasonic mechanical sector scanning

Abstract

The invention discloses a collection, storage and calculation rapid imaging method for bladder B-ultrasonic mechanical sector scanning, and relates to the technical field of medical electronics. Comprising the steps that S1, bladder B-ultrasonic mechanical sector scanning information is configured, and a mapping table from rectangular coordinates to polar coordinates is generated according to the configuration information; s2, acquiring ultrasonic probe data based on a high-speed sampling clock, dividing a field of view according to focusing depths, and performing down-sampling with different multiplying power on different focusing depths to obtain sampling points; s3, calculating a weight value of R-theta under a polar coordinate; and S4, calculating the interpolation weight of the pixel point R-theta based on the weight value of the R-theta under the polar coordinate obtained in the S3, obtaining the value of the pixel point of the corresponding image, and storing the value into a cache region. According to the invention, the problems of fixed image size, high cache demand and large imaging delay of the traditional handheld bladder B ultrasonic equipment are effectively solved, and a more efficient and more accurate solution is provided.

Description

Technical Field

[0001] The invention relates to the field of medical electronic technology, and in particular to a rapid imaging method for acquisition, storage and calculation of bladder B-ultrasound mechanical fan scanning. Background Art

[0002] Bladder ultrasound examination is one of the commonly used non-invasive examination methods in clinical diagnosis and treatment, and is particularly important in assessing bladder capacity and detecting bladder lesions. Traditional handheld bladder ultrasound devices have some technical limitations when performing mechanical fan scanning imaging. The image size of traditional handheld bladder ultrasound devices is pre-fixed and cannot be flexibly adjusted according to specific application requirements. This fixed size limits the applicability and flexibility of the image and cannot provide the optimal image resolution and details under different diagnostic requirements.

[0003] Traditional handheld bladder ultrasound devices usually use an interpolation method in a rectangular coordinate system. When processing mechanical fan scan data in polar coordinates, this method must first cache all the data of a scan before performing interpolation calculations. This method requires a large amount of cache space. Since the data of the entire scan line needs to be cached before processing, it places high demands on the device's memory resources. At the same time, the delay is high, and the interpolation calculation is performed after the scan is completed, which increases the delay of the imaging process and affects the efficiency of real-time imaging and user experience.

[0004] In summary, traditional handheld bladder B-ultrasound devices have significant deficiencies in terms of flexible adjustment of image size, cache requirements and imaging delay. These technical limitations seriously affect the effectiveness of the device in practical applications and user experience.

[0005] Therefore, it is an urgent problem for technicians in this field to propose a rapid imaging method for mechanical fan scanning of bladder ultrasound to solve the difficulties existing in the prior art. Summary of the invention

[0006] In view of this, the present invention provides a rapid imaging method for bladder B-ultrasound mechanical fan scanning, which can dynamically adjust the sampling information and image size according to specific diagnostic needs, improve the applicability and flexibility of the image, and provide higher quality imaging effects in different diagnostic scenarios.

[0007] In order to achieve the above object, the present invention adopts the following technical solution:

[0008] A rapid imaging method for collecting, storing and calculating bladder B-ultrasound mechanical fan scanning, comprising the following steps:

[0009] S1: configuring bladder B-ultrasound mechanical sector scanning information, and generating a mapping table from rectangular coordinates to polar coordinates according to the configuration information, the configuration information including: ultrasound scanning depth, field of view angle and B-ultrasound image size;

[0010] S2: Acquire ultrasound probe data based on a high-speed sampling clock to obtain a field of view image, divide the field of view into a near field, a bladder focus area, and a far field according to the focus depth, and use different magnifications of downsampling for different focus depths to obtain sampling points;

[0011] S3: Using the improved R-theta interpolation method, the weight value of R-theta in polar coordinates is calculated for the sampling points obtained in S2 based on the image size set in S1;

[0012] S4: Based on the weight value of R-theta in polar coordinates obtained in S3, the interpolation weight of the pixel point R-theta is calculated, and the value of the corresponding image pixel point is obtained and stored in the cache area.

[0013] The above method optionally configures the scanning depth and field of view angle of the ultrasonic probe in S1:

[0014] According to the viewing angle Control the number of motor rotations:

[0015] Among them, the motor step angle is Δθ, and L is the total number of scanning lines of the B-ultrasound fan scan.

[0016] Number of high-speed sampling points on each scan line:

[0017] The sampling time of the ultrasonic probe is controlled according to the scanning depth H, and the sound speed is c, f s The sampling frequency is 10MHz-40MHz.

[0018] The above method optionally creates a mapping table from the rectangular coordinates of the image to the polar coordinates:

[0019] The rectangular coordinates are expressed as (x, y), that is, the pixel point at the yth row and the xth column. The upper left corner of the image is set as the rectangular coordinate origin, and the rectangular coordinates of the center position of the first row of the image are set as The point is the origin of polar coordinates;

[0020] Calculate the depth of the pixel point rectangular coordinates mapped to polar coordinates:

[0021]

[0022] Calculate the angle from the rectangular coordinates of the pixel to the polar coordinates:

[0023]

[0024] A mapping table of rectangular coordinates to polar coordinates of an image is generated according to the configured B-ultrasound scanning information and image information. The polar coordinate mapping of a pixel point (x, y) is represented as (h, θ), that is, the distance is h and the angle is θ.

[0025] In the above method, optionally, the downsampling strategy according to the depth in S2 includes:

[0026] A higher downsampling ratio is used in the near field, a moderate downsampling ratio is used in the bladder focus area, and a higher downsampling ratio is used in the far field.

[0027] For the near field within the focal depth of 5 cm, the downsampling factor is f s / 1000c, per f s / 1000c data are averaged to get a reduced sampling data, and each scanning line gets 100 reduced sampling data after reducing sampling in the near field area;

[0028] For the bladder region with a focal depth of 5cm-10cm, the downsampling factor is f s / 2000c, per f s / 2000c data are averaged to obtain a downsampled data, and each scanning line obtains 200 downsampled data after downsampled near the bladder focal area;

[0029] For the far field with a focal depth of more than 10 cm, the downsampling factor is f s (H-0.1) / 50c, per f s The (H-0.1) / 50c data are averaged to obtain a reduced sampling data. After each scanning line is reduced in the far field, 100 reduced sampling data are obtained.

[0030] The average value of the downsampled data is v a,b , which represents the value of the b-th downsampling point on the a-th scan line.

[0031] In the above method, optionally, calculating the weight value of R-theta in polar coordinates in S3 includes the following steps:

[0032] S301: Calculate the polar coordinates of the downsampled data (h b ,θ a ), represents the polar coordinates of the b-th downsampling point on the a-th scan line. According to the downsampling strategy, the depth h of the downsampling point is calculated. b , the downsampling point angle and the scan line l corresponding to the sampling point a The angle of the downsampling point is calculated as:

[0033]

[0034] S302: Obtain the sampling point (h) according to the polar coordinate mapping tablea ,θ a ) The interpolated pixel point (x, y) in the image, the rectangular coordinates of the pixel point meet the following conditions:

[0035] h b-1 <h (x,y) ≤h b ;

[0036] S303: Calculate the depth weight w of the R-theta interpolation algorithm h and the angle weight w θ ;

[0037] The depth value h of the polar coordinates mapped according to the pixel point (x, y) (x,y) , the weight on depth is calculated as:

[0038]

[0039] The depth value θ of the polar coordinates mapped according to the pixel point (x, y) (x,y) , the angle weight is calculated as:

[0040]

[0041] S304: Perform R-theta interpolation operation using the weights obtained in S303 to obtain the grayscale value of the pixel point (x, y):

[0042] image (x,y) =(1-w θ )·[(1-w h )·v a-1,b-1 +w h ·v a-1,b ]+w θ ·[(1-w h )·v a,b-1 +w h ·v a,b ];

[0043] Among them, v a-1,b-1 Indicates the b-1th sampling point on the a-1th scan line.

[0044] The above method, optionally, stores it in the cache area:

[0045] A storage area of ​​I×J is established to sequentially store images of rectangular coordinates I×J. Initially, the storage area is completely cleared to zero. After the grayscale value of the pixel point (x, y) is calculated, it is stored in the address space of x×y.

[0046] It can be seen from the above technical solution that, compared with the prior art, the fast imaging method for bladder B-ultrasound mechanical fan scanning of the present invention has the following beneficial effects:

[0047] Bilinear interpolation calculation is performed in polar coordinates, which can be dynamically adjusted according to the image size and adapt to different sampling strategies to achieve efficient B-ultrasound image restoration;

[0048] Real-time interpolation calculation without data caching: The weights of R-theta bilinear interpolation are calculated in real time, without the need to cache all scanned data for interpolation, which greatly reduces memory requirements and imaging delay, and improves imaging speed and real-time performance.

[0049] In summary, the present invention effectively solves the problems of fixed image size, high buffering requirements and large imaging delay of traditional handheld bladder B-ultrasound devices through an innovative real-time interpolation method, providing a more efficient and accurate solution for clinical diagnosis. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings required for use in the embodiments or the description of the prior art will be briefly introduced below. Obviously, the drawings described below are only embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying creative work.

[0051] Figure 1 A rapid imaging method for acquisition, storage and calculation of bladder B-ultrasound mechanical fan scanning provided by the present invention;

[0052] Figure 2 It is a schematic diagram of different depths of the present invention;

[0053] Figure 3 The rectangular coordinate and polar coordinate mapping diagrams used in the present invention;

[0054] Figure 4 Schematic diagram of the R-theta interpolation method of the present invention. DETAILED DESCRIPTION

[0055] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0056] In this application, relational terms such as first and second, etc. are used only to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply any such actual relationship or order between these entities or operations. The terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that a process, method, article or device including a series of elements includes not only those elements, but also other elements not explicitly listed, or also includes elements inherent to such process, method, article or device. In the absence of more restrictions, the elements defined by the sentence "comprise one..." do not exclude the existence of other identical elements in the process, method, article or device including the elements.

[0057] The present invention can be used in many general or special computing device environments or configurations, such as personal computers, server computers, handheld or portable devices, tablet devices, multi-processor devices, distributed computing environments including any of the above devices or devices, etc.

[0058] Reference Figure 1 As shown, the present invention discloses a rapid imaging method for bladder B-ultrasound mechanical fan scanning, comprising the following steps:

[0059] S1: configuring bladder B-ultrasound mechanical sector scanning information, and generating a mapping table from rectangular coordinates to polar coordinates according to the configuration information, the configuration information including: ultrasound scanning depth, field of view angle and B-ultrasound image size;

[0060] S2: Acquire ultrasound probe data based on a high-speed sampling clock to obtain a field of view image, divide the field of view into a near field, a bladder focus area, and a far field according to the focus depth, and use different magnifications of downsampling for different focus depths to obtain sampling points;

[0061] S3: Using the improved R-theta interpolation method, the weight value of R-theta in polar coordinates is calculated for the sampling points obtained in S2 based on the image size set in S1;

[0062] S4: Based on the weight value of R-theta in polar coordinates obtained in S3, the interpolation weight of the pixel point R-theta is calculated, and the value of the corresponding image pixel point is obtained and stored in the cache area.

[0063] Specifically, the size of the B-ultrasound image is I×J, which is represented by a B-ultrasound image with I rows and J columns;

[0064] The sampling frequency is f s , the high-speed sampling clock is above 10MHz. The depth of bladder B-ultrasound scanning is usually above 15cm, and the bladder is mainly located at 5-10cm.

[0065] Furthermore, the scanning depth and field of view angle of the ultrasonic probe in S1 are configured as follows:

[0066] According to the viewing angle Control the number of motor rotations:

[0067] Among them, the motor step angle is Δθ, and L is the total number of scanning lines of the B-ultrasound fan scan;

[0068] Number of high-speed sampling points on each scan line:

[0069] The sampling time of the ultrasonic probe is controlled according to the scanning depth H, and the sound speed is c, f s The sampling frequency is 10MHz-40MHz.

[0070] Furthermore, a mapping table from rectangular coordinates to polar coordinates of the image is established:

[0071] The rectangular coordinates are expressed as (x, y), that is, the pixel point at the yth row and the xth column. The upper left corner of the image is set as the rectangular coordinate origin, and the rectangular coordinates of the center position of the first row of the image are set as The point is the origin of polar coordinates;

[0072] Calculate the depth of the pixel point rectangular coordinates mapped to polar coordinates:

[0073]

[0074] Calculate the angle from the rectangular coordinates of the pixel to the polar coordinates:

[0075]

[0076] A mapping table of rectangular coordinates to polar coordinates of an image is generated according to the configured B-ultrasound scanning information and image information. The polar coordinate mapping of a pixel point (x, y) is represented as (h, θ), that is, the distance is h and the angle is θ.

[0077] Furthermore, the depth-based downsampling strategies in S2 include:

[0078] A higher downsampling ratio is used in the near field, a moderate downsampling ratio is used in the bladder focus area, and a higher downsampling ratio is used in the far field.

[0079] For the near field within the focal depth of 5 cm, the downsampling factor is f s / 1000c, per f s / 1000c data are averaged to get a reduced sampling data, and each scanning line gets 100 reduced sampling data after reducing sampling in the near field area;

[0080] For the bladder region with a focal depth of 5cm-10cm, the downsampling factor is f s / 2000c, per f s / 2000c data are averaged to obtain a downsampled data, and each scanning line obtains 200 downsampled data after downsampled near the bladder focal area;

[0081] For the far field with a focal depth of more than 10 cm, the downsampling factor is f s (H-0.1) / 50c, per f s The (H-0.1) / 50c data are averaged to obtain a reduced sampling data. After each scanning line is reduced in the far field, 100 reduced sampling data are obtained.

[0082] The average value of the downsampled data is v a,b , which represents the value of the b-th downsampling point on the a-th scan line.

[0083] Specifically, Figure 2 The depth-based downsampling strategy includes: using a higher downsampling ratio in the near field to reduce the amount of data and processing burden; using a moderate downsampling ratio in the bladder focus area to balance image quality and processing speed; and using a higher downsampling ratio in the far field to improve image quality.

[0084] Furthermore, the calculation of the weight value of R-theta in polar coordinates in S3 includes the following steps:

[0085] S301: Calculate the polar coordinates of the downsampled data (h b ,θ a ), represents the polar coordinates of the b-th downsampling point on the a-th scan line. According to the downsampling strategy, the depth h of the downsampling point is calculated. b , the downsampling point angle and the scan line l corresponding to the sampling point a The angle of the downsampling point is calculated as:

[0086]

[0087] S302: According to Figure 3 The polar coordinate mapping table shown in the figure is used to obtain the sampling point (h a ,θ a ) The interpolated pixel point (x, y) in the image, the rectangular coordinates of the pixel point meet the following conditions:

[0088] h b-1 <h (x,y) ≤h b ;

[0089] S303: Calculate the depth weight w of the R-theta interpolation algorithm h and the angle weight wθ , R-theta interpolation algorithm such as Figure 4 As shown;

[0090] The depth value h of the polar coordinates mapped according to the pixel point (x, y) (x,y) , the weight on depth is calculated as:

[0091]

[0092] The depth value θ of the polar coordinates mapped according to the pixel point (x, y) (x,y) , the angle weight is calculated as:

[0093]

[0094] S304: Perform R-theta interpolation operation using the weights obtained in S303 to obtain the grayscale value of the pixel point (x, y):

[0095] image (x,y) =(1-w θ )·[(1-w h )·v a-1,b-1 +w h ·v a-1,b ]+w θ ·[(1-w h )·v a,b-1 +w h ·v a,b ];

[0096] Among them, v a-1,b-1 Indicates the b-1th sampling point on the a-1th scan line.

[0097] Furthermore, it is stored in the cache area:

[0098] A storage area of ​​I×J is established to sequentially store images of rectangular coordinates I×J. Initially, the storage area is completely cleared to zero. After the grayscale value of the pixel point (x, y) is calculated, it is stored in the address space of x×y.

[0099] In a specific embodiment, a sampling, storage and computing architecture is proposed, including: front-end data acquisition, stepper motor control, ultrasonic probe transceiver control, and high-speed ultrasonic data acquisition; dynamic processing process, dynamic filtering of echo signals, envelope detection, and dynamic downsampling according to the focus depth; interpolation area calculation, when new data is collected, the image area of ​​the collected data interpolation is calculated, and the depth and angle of the pixel points mapped to the polar coordinates are calculated; weight calculation, the depth weight and angle weight of the pixel points are calculated; R-theta interpolation, the grayscale value of the image pixel points is calculated; storage process, the interpolated pixel points are stored in the storage space according to the rectangular coordinates.

[0100] Each embodiment in this specification is described in a progressive manner, and the same or similar parts between the embodiments can refer to each other, and each embodiment focuses on the differences from other embodiments. In particular, for the system or system embodiment, since it is basically similar to the method embodiment, the description is relatively simple, and the relevant parts can refer to the partial description of the method embodiment. The system and system embodiments described above are merely schematic, wherein the units described as separate components may or may not be physically separated, and the components displayed as units may or may not be physical units, that is, they may be located in one place, or they may be distributed on multiple network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Ordinary technicians in this field can understand and implement it without creative work.

[0101] Those skilled in the art may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein may be implemented by electronic hardware, computer software, or a combination of both.

[0102] In order to clearly illustrate the interchangeability of hardware and software, the above description has generally described the components and steps of each example according to function. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professional and technical personnel can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the present invention.

[0103] The above description of the disclosed embodiments enables one skilled in the art to implement or use the present invention. Various modifications to these embodiments will be apparent to one skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present invention. Therefore, the present invention will not be limited to the embodiments shown herein, but rather to the widest scope consistent with the principles and novel features disclosed herein.

Claims

1. A rapid imaging method for mechanical sector scanning of bladder ultrasound, characterized in that: The following steps are involved: S1: configuring bladder B-ultrasound mechanical sector scanning information, and generating a mapping table from rectangular coordinates to polar coordinates according to the configuration information, the configuration information including: ultrasound scanning depth, field of view angle and B-ultrasound image size; S2: Acquire ultrasound probe data based on a high-speed sampling clock to obtain a field of view image, divide the field of view into a near field, a bladder focus area, and a far field according to the focus depth, and use different magnifications of downsampling for different focus depths to obtain sampling points; S3: Using the improved R-theta interpolation method, the weight value of R-theta in polar coordinates is calculated for the sampling points obtained in S2 based on the image size set in S1; S4: Based on the weight value of R-theta in polar coordinates obtained in S3, the interpolation weight of the pixel point R-theta is calculated, and the value of the corresponding image pixel point is obtained and stored in the cache area.

2. The method for rapid imaging of bladder B-ultrasound mechanical sector scanning according to claim 1, characterized in that: The scanning depth and field of view angle configuration of the ultrasonic probe in S1: According to the viewing angle Control the number of motor rotations: Among them, the motor step angle is Δθ, and L is the total number of scanning lines of the B-ultrasound fan scan: Number of high-speed sampling points on each scan line: The sampling time of the ultrasonic probe is controlled according to the scanning depth H, c is the speed of sound, f s The sampling frequency is 10MHz-40MHz.

3. The method for rapid imaging of bladder B-ultrasound mechanical sector scanning according to claim 2, characterized in that: Create a mapping table from rectangular coordinates to polar coordinates of an image: The rectangular coordinates are expressed as (x, y), that is, the pixel point at the yth row and the xth column. The upper left corner of the image is set as the rectangular coordinate origin, and the rectangular coordinates of the center position of the first row of the image are set as The point is the origin of polar coordinates; Calculate the depth of the pixel point rectangular coordinates mapped to polar coordinates: Calculate the angle from the rectangular coordinates of the pixel to the polar coordinates: A mapping table of rectangular coordinates to polar coordinates of an image is generated according to the configured B-ultrasound scanning information and image information. The polar coordinate mapping of a pixel point (x, y) is represented as (h, θ), that is, the distance is h and the angle is θ.

4. The method for rapid imaging of bladder B-ultrasound mechanical sector scanning according to claim 1, characterized in that: The depth-based downsampling strategies in S2 include: A higher downsampling ratio is used in the near field, a moderate downsampling ratio is used in the bladder focus area, and a higher downsampling ratio is used in the far field. For the near field within the focal depth of 5 cm, the downsampling factor is f s / 1000c, per f s The average of 1000c data is used to obtain a reduced sampling data. After each scanning line is reduced in the near field, 100 reduced sampling data are obtained. For the bladder region with a focal depth of 5cm-10cm, the downsampling factor is f s / 2000c, per f s / 2000c data are averaged to obtain a downsampled data, and each scanning line obtains 200 downsampled data after downsampled near the bladder focal area; For the far field with a focal depth of more than 10 cm, the downsampling factor is f s (H-0.1) / 50c, per f s The (H-0.1) / 50c data are averaged to obtain a reduced sampling data. After each scanning line is reduced in the far field, 100 reduced sampling data are obtained. The average value of the downsampled data is v a,b , which represents the value of the b-th downsampling point on the a-th scan line.

5. The method for rapid imaging of bladder B-ultrasound mechanical sector scanning according to claim 2, characterized in that: Calculating the weight value of R-theta in polar coordinates in S3 includes the following steps: S301: Calculate the polar coordinates of the downsampled data (h b ,θ a ), represents the polar coordinates of the b-th downsampling point on the a-th scan line. According to the downsampling strategy, the depth h of the downsampling point is calculated. b , the downsampling point angle and the scan line l corresponding to the sampling point a The angle of the downsampling point is calculated as: S302: Obtain the sampling point (h) according to the polar coordinate mapping table a ,θ a ) The interpolated pixel point (x, y) in the image, the rectangular coordinates of the pixel point meet the following conditions: h b-1 <h (x,y) ≤h b ; S303: Calculate the depth weight w of the R-theta interpolation algorithm h and the angle weight w θ ; The depth value h of the polar coordinates mapped according to the pixel point (x, y) (x,y) , the weight on depth is calculated as: The depth value θ of the polar coordinates mapped according to the pixel point (x, y) (x,y) , the angle weight is calculated as: S304: Perform R-theta interpolation operation using the weights obtained in S303 to obtain the grayscale value of the pixel point (x, y): image (x,y) =(1-w θ )·[(1-w h )·v a-1,b-1 +w h ·v a-1,b ]+w θ ·[(1-w h )·v a,b-1 +w h ·v a,b ]; Among them, v a-1,b-1 Indicates the b-1th sampling point on the a-1th scan line.

6. The method for rapid imaging of bladder B-ultrasound mechanical sector scanning according to claim 2, characterized in that: Store in cache area: A storage area of ​​I×J is established to sequentially store images of rectangular coordinates I×J. Initially, the storage area is completely cleared to zero. After the grayscale value of the pixel point (x, y) is calculated, it is stored in the address space of x×y.

Images