Bone density measurement method, system and device

Abstract

This application relates to the field of medical technology and discloses a bone density measurement method, system, and device. The method comprises: analyzing a captured medical image to obtain coordinate values and scan values for each scan point of a target area; determining a target three-dimensional area on the medical image in response to acquired three-dimensional area definition data; determining a target scan point within the target three-dimensional area based on the coordinate values of each scan point; and correcting the scan values of each target scan point using phantom correction data to determine the equivalent bone density of the target three-dimensional area based on the corrected scan values. This application allows for customizing the target three-dimensional area, thereby flexibly measuring bone density in the desired area.

Description

Technical Field

[0001] The present application relates to the field of medical technology, and in particular to a bone density measurement method, system, and device. Background Art

[0002] Existing technologies typically use dual-energy X-ray absorptiometry (DXA) and quantitative computed tomography (QCT) to measure bone density. However, these methods can only measure bone density in a fixed area, making them difficult to meet complex clinical needs such as local bone assessment and surgical planning. Summary of the Invention

[0003] In view of this, embodiments of the present application provide a bone density measurement method, system, and device.

[0004] In a first aspect, an embodiment of the present application provides a method for measuring bone density, comprising:

[0005] In response to the acquired medical image, analyzing the medical image to obtain coordinate values and scanning values of each scanning point of the target part;

[0006] determining a target three-dimensional region on the medical image in response to the acquired three-dimensional region definition data;

[0007] determining a target scanning point within the target three-dimensional area based on the coordinate value of each of the scanning points;

[0008] The scan value of each target scan point is corrected using the phantom correction data, so as to determine the equivalent bone density of the target three-dimensional area according to the corrected scan value.

[0009] In an optional embodiment, the three-dimensional region definition data includes any one of geometric parameters, function conditions, or a preset three-dimensional measurement model.

[0010] In an optional embodiment, determining the target scanning point within the target three-dimensional area based on the coordinate value of each scanning point includes:

[0011] Determining whether the coordinate value of each scanning point is within the range of the target three-dimensional area based on a spatial geometry screening algorithm;

[0012] The scanning points within the target three-dimensional area are used as target scanning points.

[0013] In an optional embodiment, before correcting the scan value of each target scan point using the phantom correction data, the method further includes:

[0014] Scanning the target area of the patient and a phantom with a preset number of reference rods to obtain a phantom medical image of the phantom and the medical image; wherein the reference rods are embedded in the phantom, and the equivalent density of potassium hydrogen phosphate and the equivalent density of water in each reference rod are both known;

[0015] Selecting a region of interest of each reference rod in the phantom medical image, and calculating a scan value of each reference rod;

[0016] Fitting the relationship between the equivalent density of dipotassium hydrogen phosphate in the reference rod and the scan value of the reference rod by linear regression to obtain a rate of change of the scan value with the equivalent density of dipotassium hydrogen phosphate and an offset value of the scan value of the reference rod based on the equivalent density of water;

[0017] The rate of change of the scanning value with the equivalent density of dipotassium hydrogen phosphate and the offset value of the scanning value of the reference rod based on the equivalent density of water are calibrated using scanning device parameters to obtain the phantom correction data; wherein the phantom correction data includes the scanning value of the reference rod, the response slope of the scanning device to dipotassium hydrogen phosphate, and the offset value of the scanning value under the scanning device based on the equivalent density of water as the baseline.

[0018] In an optional embodiment, correcting the scan value of each target scan point using phantom correction data includes:

[0019] The equivalent bone density of each target scanning point is calculated based on the scan value of the reference rod, the response slope of the scanning device to potassium dihydrogen phosphate, and the offset value of the scan value of the scanning device with the equivalent density of water as the baseline, in combination with a correction formula;

[0020] Wherein, the correction formula is: ;

[0021] In the formula, is the equivalent bone density of the target scanning point; is the scan value of the reference rod; is the offset value of the scan value after calibration; is the response slope of the scanning device to potassium dihydrogen phosphate.

[0022] In an optional embodiment, determining the equivalent bone density of the target three-dimensional area according to the corrected scan value of the target scan point includes:

[0023] The equivalent bone density of the target three-dimensional area is determined according to the equivalent bone density of each target scanning point and the weight of each scanning point.

[0024] In a second aspect, an embodiment of the present application provides a bone density measurement system, comprising: a processing module, a display module, and an input module;

[0025] The input module is used to input three-dimensional area definition data and transmit it to the processing module so that the processing module determines the target three-dimensional area;

[0026] The processing module is configured to execute the bone density measurement method described in any one of the aforementioned embodiments to obtain bone density information of the target three-dimensional area;

[0027] The display module is used to display the bone density distribution map generated by the processing module based on the bone density information of the target three-dimensional area.

[0028] In an optional embodiment, the processing module encapsulates the pydicom sub-library, pre-processing sub-module and bone density distribution map generation module of the Python library;

[0029] The preprocessing submodule is used to preprocess the acquired medical image and transmit the preprocessed medical image to the pydicom sub-library, so that the pydicom sub-library processes the preprocessed medical image to obtain the coordinate value and scanning value of each scanning point;

[0030] The bone density distribution map generating module is used to generate the corresponding bone density distribution map after obtaining the bone density of the target three-dimensional area, and send it to the display module for display.

[0031] In an optional embodiment, the system further includes a three-dimensional measurement model setting module;

[0032] The three-dimensional measurement model setting module is used to set and save the three-dimensional measurement model.

[0033] In a third aspect, an embodiment of the present application provides a bone density measurement device, which includes a processor and a memory, wherein the memory stores a computer program, and the processor is used to execute the computer program to implement the bone density measurement method described in any one of the aforementioned embodiments or to implement the functions of each module in the bone density measurement system described in any one of the aforementioned embodiments.

[0034] The embodiments of the present application have the following beneficial effects: The present application can define the target three-dimensional region in a customized manner, thereby achieving flexible definition of the target three-dimensional region. This avoids the problem in the prior art that QCT and DXA methods can only measure bone density in a fixed measurement area and cannot measure bone density according to actual needs. Furthermore, the present application corrects the scan values using a correction phantom, and then obtains the equivalent bone density of the target three-dimensional region based on the corrected scan values. This can eliminate equipment errors and environmental interference, ensure the accuracy of the scan values, and thus the accuracy of the final bone density. BRIEF DESCRIPTION OF THE DRAWINGS

[0035] In order to more clearly illustrate the technical solutions of the embodiments of the present application, the following is a brief introduction to the drawings required for use in the embodiments. It should be understood that the following drawings only show certain embodiments of the present application and therefore should not be regarded as limiting the scope. For ordinary technicians in this field, other relevant drawings can be obtained based on these drawings without creative work.

[0036] Figure 1 A schematic structural diagram of a bone density measurement system according to an embodiment of the present application is shown;

[0037] Figure 2 A first flow chart of the bone density measurement method according to an embodiment of the present application is shown;

[0038] Figure 3 A second flow chart of the bone density measurement method according to an embodiment of the present application is shown. DETAILED DESCRIPTION

[0039] The technical solutions in the embodiments of the present application will be described clearly and completely below in conjunction with the drawings in the embodiments of the present application. Obviously, the described embodiments are only part of the embodiments of the present application, rather than all the embodiments.

[0040] The components of the embodiments of the present application generally described and illustrated in the drawings herein may be arranged and designed in a variety of different configurations. Therefore, the following detailed description of the embodiments of the present application provided in the drawings is not intended to limit the scope of the claimed application, but rather merely represents selected embodiments of the present application. All other embodiments obtained by those skilled in the art based on the embodiments of the present application without creative effort are within the scope of protection of the present application.

[0041] Hereinafter, the terms "including", "having" and their cognates used in various embodiments of the present application are intended only to indicate specific features, numbers, steps, operations, elements, components or combinations of the aforementioned items, and should not be understood as excluding the existence of one or more other features, numbers, steps, operations, elements, components or combinations of the aforementioned items or adding the possibility of one or more features, numbers, steps, operations, elements, components or combinations of the aforementioned items. In addition, the terms "first", "second", "third" and the like are only used to distinguish descriptions and should not be understood as indicating or implying relative importance.

[0042] Unless otherwise defined, all terms used herein (including technical and scientific terms) have the same meaning as commonly understood by those skilled in the art to which the various embodiments of the present application belong. The terms (such as those defined in generally used dictionaries) will be interpreted as having the same meaning as in the context of the relevant technical field and will not be interpreted as having an idealized meaning or an overly formal meaning unless clearly defined in the various embodiments of the present application.

[0043] The following describes some embodiments of the present application in detail with reference to the accompanying drawings. In the absence of conflict, the following embodiments and features in the embodiments may be combined with each other.

[0044] Bone mineral density (BMD) measurement is an important tool for assessing osteoporosis, fracture risk, and diagnosing orthopedic diseases. Currently, commonly used bone density measurement technologies include dual-energy X-ray absorptiometry (DXA) and quantitative computed tomography (QCT). However, existing technologies have the following limitations in practical applications:

[0045] First, the measurement area is limited. Traditional QCT and DXA methods typically measure only fixed anatomical regions (such as the lumbar spine or hip joint) and lack the flexibility to define and measure bone density within arbitrary three-dimensional regions. This makes it difficult to meet the flexibility and accuracy requirements for complex clinical needs, such as local bone assessment and surgical planning.

[0046] Second, data processing efficiency is low. Existing methods rely on specialized software, resulting in complex and time-consuming data processing processes, making it difficult to achieve efficient and automated bone density measurement. Existing technical methods often require manual intervention when processing large amounts of imaging data, resulting in low measurement efficiency.

[0047] 3. Insufficient precision. The existing technology is easily affected by equipment errors and environmental interference during the calibration process, resulting in a decrease in the accuracy of the measurement results.

[0048] Fourth, visualization capabilities are lacking. Existing technologies typically only output numerical bone density results, lacking a visual representation of the measurement area, making it difficult to intuitively reflect the spatial distribution of bone density. For clinicians, the lack of visualization makes it difficult to directly use the results for diagnosis and surgical planning.

[0049] Based on this, the present application proposes a bone density measurement method, system, and device. By defining the target three-dimensional area in a customized manner, the problem that the existing QCT and DXA methods can only measure bone density in a fixed measurement area and cannot measure bone density according to actual needs is avoided. In addition, the present application corrects the scan values through a correction phantom, and then obtains bone density based on the corrected scan values. This can eliminate equipment errors and environmental interference, ensure the accuracy of the scan values, and thus ensure the accuracy of the final equivalent bone density.

[0050] The bone density measurement system will be described below.

[0051] Figure 1 FIG2 is a schematic diagram showing a structure of a bone density measurement system according to an embodiment of the present application. Schematically, the bone density measurement system includes a processing module 100 , a display module 300 and an input module 200 .

[0052] In this embodiment, when bone density measurement is required for a user's bones or bone samples, the area where bone density is required is first scanned using a scanning device, which may be a QCT device. While scanning the user, a phantom is also scanned. To ensure the accuracy of subsequent corrections to the bone scan values, the bones and phantom must be scanned synchronously during scanning. For example, if scanning the user's bones, the phantom can be placed below the scanning bed. If scanning a bone sample, the bone sample and phantom can be fixed side by side on the scanning bed to ensure that the bone sample and phantom are at the same scanning level. Furthermore, before scanning the bones, corresponding scanning parameters must be set based on the bones being scanned to ensure the quality of the image data. Scanning parameters include, but are not limited to, current, voltage, and slice thickness. For example, for adult lumbar spine scanning, the voltage can be 80 kVp-120 kVp, the current can be 200 mA-300 mA, and the slice thickness can be 1-3 mm to ensure clear trabecular bone structure; for children or small animals, the voltage can be 40 kVp-900 kVp, the current can be 50 mA-100 mA, and the slice thickness can be 2-5 mm. For another example, for low-density bone (such as osteoporosis patients), the voltage can be 80 kVp to enhance the contrast between bone and soft tissue (bone tissue is more sensitive to X-ray attenuation at low energy); for high-density bone (such as athletes or large animal specimens), the voltage can be 120 kVp to avoid artifacts caused by insufficient X-ray penetration. It will be understood that the above settings for scanning equipment parameters are only exemplary and the specific situation can be determined according to the clinical situation and are not limited here.

[0053] After scanning the bones with a scanning device, medical images in DICOM format can be exported through the software that comes with the scanning device. The exported data can be stored locally or in the cloud for subsequent processing. If you want to directly measure bone density, the scanned medical images can be directly input into the processing module 100 of the bone density measurement system through the input module 200.

[0054] The input module 200 is used to input three-dimensional region definition data and transmit the data to the processing module 100 so that the processing module 100 can determine a target three-dimensional region.

[0055] Exemplarily, in this embodiment, the target three-dimensional area can be customized through the input module 200, that is, the user can input three-dimensional area definition data (such as geometric shape parameters, function equations, etc.) to determine the target three-dimensional area to be measured.

[0056] When determining a target three-dimensional region by geometric shape parameters, the parameters for determining the geometric shape can be input through the input module 200. For example, in a femoral head CT scan, a target three-dimensional region can be defined as a cube with a side length of 10 mm, with the center located at the middle of the femoral neck. The coordinates of the center point can be first determined in the medical image through the input module 200, and then the length (X-axis), width (Y-axis), height (Z-axis) and other parameters of the cube can be input through the input module 200 to further determine the target three-dimensional region. For another example, if a target three-dimensional region is defined as a sphere with a radius of 5 mm in the lumbar vertebral body, covering the cancellous bone region, the radius (or XYZ semi-axis length) and coordinates of the center of the sphere can be input through the input module 200 to further determine the target three-dimensional region.

[0057] When determining the target three-dimensional region using a functional equation, this is achieved by inputting the function for determining the target three-dimensional region into the input module 200. For example, when obtaining the bone density of the subchondral bone region on the surface of the femoral head, if the desired target three-dimensional region is the sphere, then the functional equation can be defined as z = sqrt(1 - x² - y²), and limited to the range of x² + y² ≤ 1. The center point of the target three-dimensional region is then determined, and the desired target three-dimensional region can then be obtained.

[0058] It should be noted that some prompt data may be pre-set in the input module 200. For example, if a user is unsure of the required data when entering 3D region definition data through the input module 200, the prompt data can be viewed. For example, if the user desires a spherical target 3D region but is unsure of the parameters required to define the sphere, the user can search the prompt data for information related to the sphere and determine the required parameters, such as the center point coordinates and radius. The user can then enter the radius and center point coordinates based on the desired sphere size.

[0059] The processing module 100 is used for measuring bone density to obtain bone density information of the target three-dimensional area. The processing module 100 encapsulates the pydicom sub-library of the Python library, the pre-processing sub-module and the bone density distribution map generation module.

[0060] The preprocessing submodule is used for preprocessing the medical image according to the acquired medical image and then transmitting it to the pydicom sub-library, so that the pydicom sub-library processes the preprocessed medical image to obtain the coordinate value and the scan value of each scanning point. It can be understood that in the present embodiment, the preprocessing submodule and the pydicom sub-library are encapsulated in the processing module 100, and the processing rules are set. After the processing module 100 receives the medical image, it can be directly output to the preprocessing submodule for automatic processing, and then the preprocessed medical image is automatically transmitted to the pydicom sub-library to obtain the coordinate value of each scanning point and the scan value. In the present embodiment, by encapsulating the preprocessing submodule and the pydicom sub-library in the processing module 100, the automatic processing of medical images can be achieved, thereby improving the efficiency and accuracy of data processing.

[0061] In this embodiment, a bone density distribution map generation module is also encapsulated in the processing module 100. After processing the medical image and measuring the bone density, the corresponding bone density distribution map can be directly generated by the bone density distribution map generation module and sent to the display module 300 for display. Exemplarily, the bone density distribution map generation module includes Python's matplotlib library or plotly library to generate a three-dimensional bone density distribution map based on the bone density corresponding to each scan point. This three-dimensional bone density distribution map can be displayed through the display module 300 for clinicians to view. In addition, the module can also export the bone density data as a report in Excel or PDF format for clinicians to refer to.

[0062] In some embodiments, the bone density measurement system further includes a 3D measurement model setting module, which is used to set and save a 3D measurement model. This embodiment allows for customizing the target 3D area for bone density measurement. For customizing the target 3D area, the user can define the target 3D area by inputting geometric parameters that determine the geometric area, or by inputting a function equation. The system can also pre-set 3D measurement models. For example, commonly used 3D measurement models such as spheres and cubes can be pre-set. When acquiring a target 3D area of one of these shapes, the user can directly select from a library of pre-set 3D measurement models, determine the center point of the target 3D area, and then adjust or rotate the pre-set 3D measurement model to determine the target 3D area on the medical image, thus eliminating the need for customization. For example, for a pre-set 3D measurement model of a cylinder, the user can manually select or drag the specific coordinates of the center points of the two ends of the cylinder on the medical image. For example, the user can select a specific area of the femur as the center point of one end of the cylinder, with the other end located elsewhere on the femur. When the positions of the centers of the two ends of the cylinder change, the direction and axis of the cylinder will also be redefined, and the cylinder will rotate or translate accordingly to adapt to the new positions of the centers of the two ends.

[0063] Based on the above-mentioned bone density measurement system, the bone density measurement method will be described below in conjunction with some specific embodiments.

[0064] Figure 2 A flow chart of a bone density measurement method according to an embodiment of the present application is shown. Exemplarily, the bone density measurement method includes the following steps:

[0065] Step S100 : In response to the acquired medical image, the medical image is analyzed to obtain the coordinate value and scanning value of each scanning point of the target part.

[0066] Exemplarily, after the processing module acquires medical images of the skeleton and the phantom scanned by a scanning device, it first preprocesses the acquired medical images. The module then automatically processes the preprocessed images in DICOM format using the encapsulated pydicom library to extract the coordinates and scan value (HU value) of each scan point. A scan point refers to the spatial location and associated density information corresponding to each pixel or voxel in the medical image. These points are the basic units of three-dimensional image data and the basis for bone density measurement. Each scan point has a specific three-dimensional coordinate (x, y, z) representing its specific location. The scan point also contains a density value corresponding to that location, typically expressed in Hounsfield Units (HU). The HU value reflects the X-ray absorption capacity of the tissue or material at that point, indirectly indicating the density characteristics of that area.

[0067] Step S200 : determining a target three-dimensional region on a medical image in response to the acquired three-dimensional region definition data.

[0068] Exemplarily, in this embodiment, there are three ways to obtain the target three-dimensional region: determining the target three-dimensional region on the medical image based on acquired geometric parameters; determining the target three-dimensional region on the medical image based on acquired function conditions; or determining the target three-dimensional region on the medical image based on acquired reference point coordinates and a selected preset three-dimensional measurement model. These three methods for obtaining the target three-dimensional region have been described above and will not be repeated here. This embodiment, through this customized measurement region method, allows for the flexible definition of any three-dimensional region for bone density measurement, breaking through the problem of limited measurement areas in existing technical methods, and can be better applied to complex clinical needs (such as local bone assessment, surgical planning, etc.).

[0069] Step S300 : determining a target scanning point within a target three-dimensional area based on the coordinate value of each scanning point.

[0070] In the above steps, after the target 3D region is determined, the boundary range of the target 3D region is also determined. This step mainly involves filtering out all scan points (voxels) within the user-defined target 3D region from each scan point of the medical image.

[0071] In some implementations, a spatial geometry screening algorithm can be used to first determine whether the coordinates of each scan point are within the target 3D region. Scan points within the target 3D region are then used as target scan points. For example, a mathematical inequality can be used to determine whether a scan point is within the target 3D region, or an equation constraint can be used to determine whether a scan point is within the target 3D region.

[0072] The method of judging whether a scanning point is inside a target 3D region by mathematical inequality is usually applied to a target 3D region with a regular shape that can be determined by geometric parameters. For example, the target 3D region is a cube, and the boundaries of the X, Y, and Z axes of the cube can be given, that is, (minimum value of the X axis), (maximum value of the X axis), (minimum value of the Y axis), (maximum value of the Y axis), (minimum value of Z axis), (Z axis maximum value). When determining whether each scanning point is in the target three-dimensional area, it is only necessary to determine whether the coordinate value (x, y, z) of each scanning point is within the above-mentioned cubic range, that is, it is necessary to satisfy 、 and , if it satisfies, then the scan point is determined to be inside the cube. For another example, if the target three-dimensional area is a sphere, then the coordinates of the sphere center (x0, y0, z0) and the radius r can be given to check whether the coordinates (x, y, z) of each scan point satisfy , if satisfied, it means that the scanning point is within the target three-dimensional area of the sphere.

[0073] For determining whether the scanning point is within the target three-dimensional area through equation constraints, this method is usually applied to target three-dimensional areas with complex shapes that can be defined by function equations (such as surfaces, ellipsoids, etc.). For example, for the target three-dimensional area of an ellipsoid, the equation F(x,y,z)≤0 can be given to check whether the scanning point coordinates (x,y,z) satisfy the inequality, such as the ellipsoid ; If it is satisfied, it means that the scanning point is within the target three-dimensional area within the ellipsoid, where (x0, y0, z0) are the coordinates of the center of the ellipsoid, and a, b, and c are the semi-axis lengths in different directions.

[0074] For a target 3D region determined using a preset 3D measurement model, once the position and shape of the preset 3D measurement model on the medical image have been determined, the boundary of the target 3D region corresponding to the preset 3D measurement model is also determined. When determining target scan points, the coordinates of each scan point can be determined to be within the boundary. If so, the scan point is considered a target scan point. In the example above, where the preset 3D measurement model is a cylinder, once the shape and position of the cylinder are determined, the system will filter out target scan points from all scan points in the medical image based on the cylinder's boundary. When determining target scan points, the system checks the coordinates of each scan point one by one to determine whether it falls within the current cylinder's boundary. Specifically, for any scan point (x, y, z), the system first checks whether the point is within the axial range of the cylinder (i.e., the length between the two center points), and then checks whether the point is within the radial range of the cylinder (i.e., whether the distance from the cylinder's axis is less than the radius). If the scan point meets both of these conditions, it is marked as a target scan point.

[0075] After the scan points within the target three-dimensional area are screened out, the scan values of the screened out points may be stored as a data structure (such as a list or a matrix) for subsequent processing.

[0076] Step S400 , correcting the scan value of each target scan point using the phantom correction data, so as to determine the equivalent bone density of the target three-dimensional area according to the corrected scan value.

[0077] In some embodiments, as Figure 3 As shown, before step S400, the following steps are also included:

[0078] Step S500 : When scanning the target part of the patient, a phantom with a preset number of reference rods is simultaneously scanned to obtain a phantom medical image of the phantom.

[0079] The reference rods are embedded in the phantom, and the equivalent density of potassium hydrogen phosphate and water in each reference rod is known. That is, the reference rods are usually composed of materials with known low atomic number and high atomic number, and these reference materials can be measured by different concentrations of (Dipotassium hydrogen phosphate) solution or equivalent material is calibrated to simulate the X-ray attenuation characteristics of bone mineral (hydroxyapatite). The number of reference rods embedded in the phantom needs to cover the density range of the bone to be measured. For example, for the measurement of lumbar cancellous bone, its density is usually between 50 -300 Therefore, 5 reference rods can be set up, and the equivalent of each reference rod is The density is known, such as 0 , 50 , 100 , 200 , and 300 Because scan values from different scanning devices may differ, reference rods are used in this embodiment to provide a standardized benchmark and reduce fitting errors. The arrangement of the reference rods in the phantom can be determined based on the bones being measured. For example, to measure spinal bone density, the reference rods can be arranged in parallel within the phantom; to measure hip joint bone density, the reference rods can be arranged in a circular pattern within the phantom.

[0080] Step S600 : selecting a region of interest of each reference rod in the phantom medical image, and calculating a scan value of each reference rod.

[0081] In this step, a region of interest of each reference rod is selected from the acquired medical image of the phantom. Usually, the center of each reference rod is used as the center point, and the region of interest is selected by spreading outward, and the scanning value of each reference rod is obtained through the region of interest.

[0082] In step S700 , linear regression is used to fit the relationship between the equivalent density of dipotassium hydrogen phosphate in the reference rod and the scan value of the reference rod to obtain the rate of change of the scan value with the equivalent density of dipotassium hydrogen phosphate and the offset value of the scan value of the reference rod with the equivalent density of water as the baseline.

[0083] The fitting formula is: ;In the formula, is the scan value of the reference rod, is the equivalent density of water, is the rate of change of the scanning value with the equivalent density of dipotassium hydrogen phosphate, is the equivalent density of dipotassium hydrogen phosphate in the reference rod, is the offset value of the reference rod scan value based on the equivalent density of water as the baseline.

[0084] In step S800 , the rate of change of the scan value with the equivalent density of dipotassium hydrogen phosphate and the offset value of the reference rod scan value based on the equivalent density of water are calibrated by scanning device parameters to obtain phantom calibration data.

[0085] The phantom calibration data includes the scan value of the reference rod, the response slope of the scanning device to potassium dihydrogen phosphate, and the offset value of the scan value under the calibrated scanning device with the equivalent density of water as the baseline.

[0086] The rate of change of the scanning value with the equivalent density of potassium hydrogen phosphate The calibration formula is: ; Offset value of the reference rod scan value based on the equivalent density of water The calibration formula is: ; in the formula For The response slope of the scanning device to potassium dihydrogen phosphate obtained after calibration, For The offset values, e and f, obtained from the calibration of the scanning device using the equivalent density of water as the baseline are determined based on the calibration data for the corresponding phantom. Note that this calibration data is obtained according to the manufacturer's specifications. For example, for the Mindways Model 3 phantom, the parameters given in its specifications indicate e to be 0.2174 and f to be 999.6.

[0087] After obtaining the phantom correction data, this embodiment corrects the scan value of the target scanning point using the phantom correction data to avoid deviation of the scan value.

[0088] Exemplarily, the calibrated scan value of each target scanning point is calculated based on the scan value of the reference rod, the response slope of the scanning device to potassium dihydrogen phosphate, and the offset value of the scan value based on the equivalent density of water under the scanning device, combined with the correction formula.

[0089] The correction formula is: ;In the formula, It is the scan value after correction of the target scanning point, that is, the equivalent bone density of the target scanning point.

[0090] After obtaining the equivalent bone density of each target scanning point, in this embodiment, the equivalent bone density of the target three-dimensional area can be determined based on the equivalent bone density of each target scanning point and the weight of each scanning point.

[0091] If all scan points within the target three-dimensional area have the same weight, the equivalent bone density values of all target scan points can be directly averaged, and the calculation formula is: ;In the formula, is the equivalent bone density of the target three-dimensional area, N is the total number of all target scanning points in the target three-dimensional area, Scan point for the i-th target.

[0092] If different scan points within the target 3D region have different weights (for example, voxels with larger volumes contribute more to the total bone density), a weighted average method is required, which is calculated as follows: ;In the formula, is the equivalent bone density of the target three-dimensional area, N is the total number of all target scanning points in the target three-dimensional area, is the i-th target scanning point, is the weight of the i-th scan point (the weight is usually proportional to the volume of the voxel it corresponds to).

[0093] This application can customize the target three-dimensional area through user-input geometric parameters, function conditions, or by selecting a preset three-dimensional measurement model, thereby achieving flexible definition of the target three-dimensional area, avoiding the problem that the QCT and DXA methods in the prior art can only measure bone density in a fixed measurement area and cannot measure bone density according to actual needs. This application can meet diverse clinical needs through this customized target three-dimensional area. Furthermore, in this application, the preprocessing submodule and the pydicom sublibrary of the Python script are encapsulated in the processing module. In this way, after the processing module obtains the medical image, it can realize automated preprocessing of the medical image and automatically extract the coordinate values and scan values of the scanning points, thereby reducing manual operation time. In addition, after obtaining the bone density data, the present application can automatically generate a bone density distribution map and a numerical report, simplifying the workflow of the clinician. Furthermore, in this application, the scan value is corrected by a correction phantom, and then the bone density is obtained based on the corrected scan value. This can eliminate equipment errors and environmental interference, ensure the accuracy of the scan value, and thus ensure the accuracy of the final bone density.

[0094] The present application also provides a bone density measurement device. Exemplarily, the terminal device includes a processor and a memory, wherein the memory stores a computer program, and the processor runs the computer program to enable the terminal device to execute the functions of the various modules in the above-mentioned bone density measurement method or the above-mentioned bone density measurement system.

[0095] The processor can be an integrated circuit chip with signal processing capabilities. The processor can be a general-purpose processor, including at least one of a central processing unit (CPU), a graphics processing unit (GPU), a network processor (NP), a digital signal processor (DSP), an application-specific integrated circuit (ASIC), a field-programmable gate array (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, or discrete hardware components. The general-purpose processor can be a microprocessor or any conventional processor, etc., and can implement or execute the various methods, steps, and logic block diagrams disclosed in the embodiments of this application.

[0096] The memory may be, but is not limited to, random access memory (RAM), read-only memory (ROM), programmable read-only memory (PROM), erasable programmable read-only memory (EPROM), and electrically erasable programmable read-only memory (EEPROM). The memory is used to store computer programs, and the processor can execute the computer programs accordingly after receiving an execution instruction.

[0097] This application also provides a computer-readable storage medium for storing the computer program used in the above-mentioned bone density measurement device. For example, the computer-readable storage medium may include, but is not limited to, various media capable of storing program code, such as a USB flash drive, a mobile hard drive, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0098] In the several embodiments provided in this application, it should be understood that the disclosed devices and methods can also be implemented in other ways. The device embodiments described above are merely schematic. For example, the flowcharts and structure diagrams in the accompanying drawings show the possible architectures, functions and operations of the devices, methods and computer program products according to the multiple embodiments of the present application. In this regard, each box in the flowchart or block diagram can represent a module, a program segment or a part of the code, and the module, program segment or a part of the code contains one or more executable instructions for implementing the specified logical functions. It should also be noted that in an alternative implementation, the functions marked in the box can also occur in an order different from that marked in the accompanying drawings. For example, two consecutive boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the structure diagram and / or flowchart, and the combination of boxes in the structure diagram and / or flowchart, can be implemented using a dedicated hardware-based system that performs the specified function or action, or can be implemented using a combination of dedicated hardware and computer instructions.

[0099] In addition, the functional modules or units in the various embodiments of the present application can be integrated together to form an independent part, or each module can exist independently, or two or more modules can be integrated to form an independent part.

[0100] If the functions are implemented in the form of software function modules and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present application, or the part that contributes to the prior art, or part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes a number of instructions for enabling a computer device (which can be a smart phone, personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present application.

[0101] The above is only a specific implementation method of the present application, but the scope of protection of the present application is not limited thereto. Any technician familiar with this technical field can easily think of changes or replacements within the technical scope disclosed in this application, which should be covered by the scope of protection of the present application.

Claims

1. A method for measuring bone density, characterized in that: include: In response to the acquired medical image, analyzing the medical image to obtain coordinate values and scanning values of each scanning point of the target part; determining a target three-dimensional region on the medical image in response to the acquired three-dimensional region definition data; The three-dimensional area definition data is user-defined information input, including a preset three-dimensional measurement model. When the target three-dimensional area is determined by using the preset three-dimensional measurement model, the preset three-dimensional measurement model can be adjusted. determining a target scanning point within the target three-dimensional area based on the coordinate value of each of the scanning points; Correcting the scan value of each target scan point using phantom correction data to determine the equivalent bone density of the target three-dimensional area based on the corrected scan value; Before correcting the scan value of each target scanning point by using the phantom correction data, the method further includes: scanning the target part of the patient, and scanning a phantom with a preset number of reference rods to obtain a phantom medical image of the phantom and the medical image; wherein the reference rods are embedded in the phantom, and the equivalent density of potassium hydrogen phosphate and the equivalent density of water in each reference rod are known; selecting a region of interest of each reference rod in the phantom medical image, and calculating the scan value of each reference rod; fitting the equivalent density of potassium hydrogen phosphate in the reference rod by linear regression. The method comprises the steps of: calculating a relationship between the scanning value of the reference rod and the equivalent density of the dipotassium hydrogen phosphate, and calculating a relationship between the scanning value of the reference rod and the equivalent density of the dipotassium hydrogen phosphate, so as to obtain a rate of change of the scanning value with the equivalent density of the dipotassium hydrogen phosphate and an offset value of the scanning value of the reference rod with the equivalent density of water as a baseline; calibrating the rate of change of the scanning value with the equivalent density of the dipotassium hydrogen phosphate and the offset value of the scanning value of the reference rod with the equivalent density of water as a baseline by using scanning device parameters to obtain the phantom calibration data; wherein the phantom calibration data includes the scanning value of the reference rod, the response slope of the scanning device to the dipotassium hydrogen phosphate, and the offset value of the scanning value of the scanning device with the equivalent density of water as a baseline; The correcting the scan value of each target scanning point using the phantom correction data includes calculating the equivalent bone density of each target scanning point based on the scan value of the reference rod, the response slope of the scanning device to dipotassium hydrogen phosphate, and the offset value of the scan value under the scanning device with the equivalent density of water as the baseline, in combination with a correction formula.

2. The bone density measurement method according to claim 1, characterized in that: The three-dimensional region definition data includes any one of geometric parameters, function conditions or a preset three-dimensional measurement model.

3. The bone density measurement method according to claim 1, characterized in that: Determining a target scanning point within the target three-dimensional area based on the coordinate value of each scanning point includes: Determining whether the coordinate value of each scanning point is within the range of the target three-dimensional area based on a spatial geometry screening algorithm; The scanning points within the target three-dimensional area are used as target scanning points.

4. The bone density measurement method according to claim 1, characterized in that The correction formula is: ; In the formula, is the equivalent bone density of the target scanning point; is the scan value of the reference rod; is the offset value of the scan value after calibration; is the response slope of the scanning device to potassium dihydrogen phosphate.

5. The bone density measurement method according to claim 1, characterized in that: Determining the equivalent bone density of the target three-dimensional area according to the corrected scan value of the target scan point includes: The equivalent bone density of the target three-dimensional area is determined according to the equivalent bone density of each target scanning point and the weight of each scanning point.

6. A bone density measurement system, characterized in that: include: processing module, display module and input module; The input module is used to input three-dimensional area definition data and transmit it to the processing module so that the processing module determines the target three-dimensional area; The processing module is used to execute the bone density measurement method according to any one of claims 1 to 5 to obtain bone density information of the target three-dimensional area; The display module is used to display the bone density distribution map generated by the processing module based on the bone density information of the target three-dimensional area.

7. The bone density measurement system according to claim 6, characterized in that: The processing module encapsulates the pydicom sub-library of the Python library, the pre-processing sub-module and the bone density distribution map generation module; The preprocessing submodule is used to preprocess the acquired medical image and transmit the preprocessed medical image to the pydicom sub-library, so that the pydicom sub-library processes the preprocessed medical image to obtain the coordinate value and scanning value of each scanning point; The bone density distribution map generating module is used to generate the corresponding bone density distribution map after obtaining the bone density of the target three-dimensional area, and send it to the display module for display.

8. The bone density measurement system according to claim 6, characterized in that: The system also includes a three-dimensional measurement model setting module; The three-dimensional measurement model setting module is used to set and save the three-dimensional measurement model.

9. A bone density measuring device, characterized in that: The bone density measurement device includes a processor and a memory, the memory stores a computer program, and the processor is used to execute the computer program to implement the bone density measurement method according to any one of claims 1 to 5 or implement the functions of each module in the bone density measurement system according to any one of claims 6 to 8.

Images