CN108921832A - Femoral head image analysis method, device, server and medium - Google Patents

Femoral head image analysis method, device, server and medium Download PDF

Info

Publication number
CN108921832A
CN108921832A CN201810668292.4A CN201810668292A CN108921832A CN 108921832 A CN108921832 A CN 108921832A CN 201810668292 A CN201810668292 A CN 201810668292A CN 108921832 A CN108921832 A CN 108921832A
Authority
CN
China
Prior art keywords
femoral head
threedimensional model
model
thighbone
analysis
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201810668292.4A
Other languages
Chinese (zh)
Inventor
陈卫衡
陈端端
何海军
李泰贤
李*
王荣田
闫天翼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Individual
Original Assignee
Individual
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN201810668292.4A priority Critical patent/CN108921832A/en
Publication of CN108921832A publication Critical patent/CN108921832A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • G06T2207/30008Bone

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Quality & Reliability (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Graphics (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Image Processing (AREA)

Abstract

The embodiment of the invention discloses a kind of femoral head image analysis method, device, server and media, wherein this method includes:Based on the target image for the femoral head for including necrotic area, the threedimensional model of the femoral head is created, wherein threedimensional model includes near end of thighbone threedimensional model and necrotic bone threedimensional model;According to threedimensional model, the analysis model of the femoral head is created;According to analysis model, extract the weight bearing area of the femoral head equivalent stress and total deformation quantity, and carry out data analysis.The embodiment of the present invention solves the problems, such as that the accuracy of femoral head image analysis result in the prior art is lower, improves the accuracy of femoral head image analysis result.

Description

Femoral head image analysis method, device, server and medium
Technical field
The present embodiments relate to image analysis technology field more particularly to a kind of femoral head image analysis method, device, Server and medium.
Background technique
Deep analysis should be carried out to treatment case in clinical research, illustrate its mechanism, the development to medicine technology field With very big directive significance, medical technology progress and the better therapeutic device of invention and material can be promoted, for example, much grinding It is because which kind of mechanism prevents collapse of femoral head that the person of studying carefully, which has illustrated from biology level and protected hip successful case, this is for orthopaedic technique There is very big directive significance.
Caput femoris necrosis seriously endangers patient health and quality of life, has been undisputable fact, has also become about its treatment Global problem, difficult point is nothing more than being:The Patients with Aseptic Necrosis of Femoral age is less than normal to be not easy to carry out joint replacement, if carrying out joint Displacement will face the risk repeatedly overhauled;Femoral head once collapses, and no matter which kind for the treatment of can not all restore femoral head initial configuration, Easy osteoarthritis secondary;Protecting hip success rate can not be completely satisfactory etc..Age patient less than normal is if protecting hip success It may not necessarily be in face of joint replacement and the risk repeatedly overhauled;If femoral head does not collapse, so that it may not have to face femur Head morphology changes the risk of the irreversible damage and its secondary osteoarthritis that cause, therefore, prevents and treats collapse of femoral head, makes hip Joint keeps initial form, is to maintain the cardinal principle of femoral head health.
Currently, the Prediction of Collapse with Magnetic Resonance Imaging analysis for caput femoris necrosis is mainly bad using caput femoris necrosis Area Prediction, femoral head The methods of dead position prediction, the prediction of caput femoris necrosis cartilage, caput femoris necrosis volume predictions and the prediction of caput femoris necrosis angle, but Be predictive analysis results accuracy it is barely satisfactory.However, at the analysis of the prediction result of collapse of femoral head and femoral head image Manage closely related, therefore, it is to need the critical issue that solves that femoral head image analysis, which how is effectively performed,.
Summary of the invention
The embodiment of the present invention provides a kind of femoral head image analysis method, device, server and medium, to solve existing skill The lower problem of the accuracy of femoral head image analysis result in art.
In a first aspect, the embodiment of the invention provides a kind of femoral head image analysis method, this method includes:
Based on the target image for the femoral head for including necrotic area, the threedimensional model of the femoral head is created, wherein described three Dimension module includes near end of thighbone threedimensional model and necrotic bone threedimensional model;
According to the threedimensional model, the analysis model of the femoral head is created;
According to the analysis model, extract the weight bearing area of the femoral head equivalent stress and total deformation quantity, line number of going forward side by side According to analysis.
Second aspect, the embodiment of the invention also provides a kind of femoral head image analysis apparatus, which includes:
Threedimensional model creation module creates the femoral head for the target image based on the femoral head for including necrotic area Threedimensional model, wherein the threedimensional model includes near end of thighbone threedimensional model and necrotic bone threedimensional model;
Analysis model creation module, for creating the analysis model of the femoral head according to the threedimensional model;
Data analysis module, for according to the analysis model, extract the femoral head weight bearing area equivalent stress with Total deformation quantity, and carry out data analysis.
The third aspect, the embodiment of the invention also provides a kind of servers, including:
One or more processors;
Storage device, for storing one or more programs,
When one or more of programs are executed by one or more of processors, so that one or more of processing Device realizes the femoral head image analysis method as described in any embodiment of the present invention.
Fourth aspect, the embodiment of the invention also provides a kind of computer readable storage mediums, are stored thereon with computer Program realizes the femoral head image analysis method as described in any embodiment of the present invention when the program is executed by processor.
The embodiment of the present invention is primarily based on the target image of the femoral head including necrotic area, creates the three-dimensional mould of the femoral head Type, wherein threedimensional model includes near end of thighbone threedimensional model and necrotic bone threedimensional model;Then, according to threedimensional model, creation should The analysis model of femoral head;Finally, according to analysis model, extract the weight bearing area of femoral head equivalent stress and total deformation quantity, and Data analysis is carried out, solves the problems, such as that the accuracy of femoral head image analysis result in the prior art is lower, from analysis femur Postnecrotic weight bearing area stress condition is started with, and the accuracy of femoral head image analysis result is improved.
Detailed description of the invention
Fig. 1 is the flow chart for the femoral head image analysis method that the embodiment of the present invention one provides;
Fig. 2 is the flow chart of femoral head image analysis method provided by Embodiment 2 of the present invention;
Fig. 3 is the structural schematic diagram for the femoral head image analysis apparatus that the embodiment of the present invention three provides;
Fig. 4 is a kind of structural schematic diagram for server that the embodiment of the present invention four provides.
Specific embodiment
The present invention is described in further detail with reference to the accompanying drawings and examples.It is understood that this place is retouched The specific embodiment stated is used only for explaining the present invention rather than limiting the invention.It also should be noted that in order to just Only the parts related to the present invention are shown in description, attached drawing rather than entire infrastructure.
Embodiment one
Fig. 1 is the flow chart for the femoral head image analysis method that the embodiment of the present invention one provides, and the present embodiment is applicable to The case where analyzing femoral head image, this method can be executed by femoral head image analysis apparatus, which can adopt It is realized, and can be integrated in the server with the mode of software and/or hardware.As shown in Figure 1, this method specifically includes:
S110, the target image based on the femoral head for including necrotic area, create the threedimensional model of femoral head, wherein three-dimensional Model includes near end of thighbone threedimensional model and necrotic bone threedimensional model.
The target image of femoral head including necrotic area includes CT image, such as the CT image of hip joint.User can lead to It crosses enabling image processing software and imports target image, be the format of image processing software identification by image format conversion, for example, figure As processing software includes Mimics software, the picture format of identification includes bin format;Then software operation is carried out, to target The different zones of image are identified and are divided, and near end of thighbone region and necrotic bone region are obtained, and it is close to respectively correspond creation femur Threedimensional model and necrotic bone threedimensional model are held, and is saved.
Image in the present embodiment for the femoral head for having existed necrotic area is analyzed, compared with the prior art middle base In the image of healthy femoral head is analyzed the case where, the three-dimensional reduction and simulation to necrotic femoral head may be implemented, be subsequent Accurate image analysis result is obtained to lay the foundation.
Optionally, target image includes the femoral head image scanned on a preset condition based, wherein preset condition includes It is at least one following:Scanning voltage, sweep current, Pixel Dimensions, scanning support gradient, Scan slice thickness and scanning slice interval.
Different sweep parameters is set, obtained image information can have differences, and then influence image analysis as a result, excellent Choosing, the preset condition of image scanning is as follows:Scanning voltage 120KV, sweep current 244.80mAs, Pixel Dimensions 0.840mm, 0.000 ° of scanning support gradient, bone tissue window scanning, thickness 1.25mm, interlayer is divided into 0.625mm.Based on scanning under this condition To femoral head image analyzed, the accuracy of image analysis result is higher.
S120, the threedimensional model according to femoral head, create the analysis model of femoral head.
User can import the near end of thighbone threedimensional model saved and necrotic bone is three-dimensional by enabling image analysis software Model carries out software operation, creates the analysis model of femoral head, wherein according to the analysis method of use, can be used different Image analysis software, image analysis method include finite element method.Finite element analysis can be by dividing human body difference group It knits and establishes threedimensional model, and the comparable load of the various physical activities of simulation loading, biomechanics Research is made not to be confined to entity Mechanical stress test, but the data that need are obtained in clinical research by computer and provide biomethanics support for clinic.
S130, according to the analysis model of creation, extract the weight bearing area of femoral head equivalent stress and total deformation quantity, and carry out Data analysis.
After analysis model creation, specific image analysis can be carried out.To the stock including necrotic area in the present embodiment Bone carries out reduction and simulation, analyzes its equivalent stress distribution and deformation under stress condition, specifically, extracting its area of bearing a heavy burden Equivalent stress and total deformation quantity.When the weight bearing area of femoral head is unable to get effective support of necrotic bone, will go out in weight bearing area Existing region of stress concentration, this region of stress concentration bear higher stress repeatedly and generate the deformation quantity more than normal range (NR), Femoral head will collapse, i.e. the stress and deformation quantity in weight bearing area are the immediate causes for leading to collapse of femoral head, therefore, to stock Bone image is effectively analyzed, and the equivalent stress of femoral head at weight loading region and total deformation quantity under load situation are accurately obtained, right There is important directive function in prediction collapse of femoral head, the accuracy of image analysis is higher, then the standard of collapse of femoral head prediction True property is higher.
Optionally, according to creation analysis model, extract femoral head weight bearing area equivalent stress and total deformation quantity, go forward side by side The analysis of row data, including:
Using pre-set angle, the weight bearing area of femoral head is determined in analysis model;
Based on analysis model, stress condition of the femoral head under walking states is simulated;
According to analog result, extract femoral head weight bearing area equivalent stress and total deformation quantity;
The equivalent stress of extraction and total deformation quantity are compared with corresponding preset threshold respectively, obtain the figure of femoral head As analysis result.
Wherein, pre-set angle is used to determine the weight bearing area of current femur head, and specific numerical value can be by user's root It is configured according to needs.The equivalent stress of extraction and total deformation quantity are compared with corresponding preset threshold respectively, can be obtained It is to maximum equivalent and maximum total deformation quantity of the weight bearing area for the femoral head for including necrotic area under load situation, i.e., required Image analysis result, wherein the preset threshold of equivalent stress and the preset threshold of total deformation quantity can by user according to demand into The setting of row adaptability.
Illustratively, true hip joint stress condition is simulated, in finite element analysis software, such as ANSYS software, into The following setting of row:Femur holding is leaned forward for 15 ° of neck of femur, the distal sections of femur model are set as staff cultivation and are fixed, if The weight bearing area (acetabular bone and femoral head contact area) of femoral head is set, this negative region is located at femur inward-outward direction opposite direction late-segmental collapse Radian be 40 °, the radian of front-rear direction opposite direction late-segmental collapse is 80 °.Hip joint will be applied with the 2.5 of weight when walking Corresponding pressure is simulated again, i.e., the pressure generated in the weight that the load of weight bearing area is equivalent to 2.5 times, and pressure calculation formula is 2.5*m*G/S, wherein m indicates weight, and G indicates that gravity coefficient, S indicate weight bearing area's area.Finally, utilizing finite element analysis software Included detecting function automatically extract the weight bearing area of femoral head equivalent stress and total deformation quantity, and then obtain the weight bearing area most Big equivalent stress and maximum total deformation quantity.
User can be using obtained maximum equivalent and maximum total deformation quantity as reference standard, i.e., the maximum equivalent is answered Power and maximum total deformation quantity can be used as the threshold value that caput femoris necrosis collapses, and carry out image analysis to using the present embodiment method Other femoral head images are made whether that there is a possibility that collapsing carries out prediction judgement.With regard to collapse of femoral head prediction angle and It says, carries out the simulation of stress condition in the present embodiment to the weight bearing area for the femoral head for including necrotic area, being equivalent to can direct base It is predicted in the reason of collapse of femoral head, therefore, the accuracy of prediction result, which is higher than, not to be related to femoral head in the prior art and collapse Fall into the result of the prediction technique of reason.
The technical solution of the present embodiment is primarily based on the target image of the femoral head including necrotic area, creates near end of thighbone three Dimension module and necrotic bone threedimensional model;Then, the threedimensional model according to creation, creates the analysis model of the femoral head;Finally, according to According to analysis model, extract the weight bearing area of femoral head equivalent stress and total deformation quantity, and carry out data analysis, solve existing skill The lower problem of the accuracy of femoral head image analysis result in art enters from the weight bearing area stress condition after analysis caput femoris necrosis Hand improves the accuracy of femoral head image analysis result, realizes the foundation that caput femoris necrosis collapses threshold value, can be femur Head Prediction of Collapse with Magnetic Resonance Imaging research provides crucial directive function, and then improves the accuracy of collapse of femoral head prediction.
Embodiment two
Fig. 2 is the flow chart of femoral head image analysis method provided by Embodiment 2 of the present invention, and the present embodiment is above-mentioned Further progress optimizes on the basis of embodiment.As shown in Fig. 2, this method specifically includes:
S210, according to the first image grayscale threshold value, the target image for the femoral head for including necrotic area is split, is determined Near end of thighbone region out, and create near end of thighbone threedimensional model.
Near end of thighbone region on target image is different with the grey scale pixel value in necrotic bone region, can be based on gray value differences It is different that target image is split.First image grayscale threshold value can be distributed according to the overall pixel near end of thighbone region to be carried out really It is fixed.
Optionally, according to the first image grayscale threshold value, the target image for the femoral head for including necrotic area is split, really Near end of thighbone region is made, and creates near end of thighbone threedimensional model, including:
According to the first image grayscale threshold value, and preset hatching is combined, the target figure to the femoral head for including necrotic area As being split, near end of thighbone region is determined, wherein hatching is the element marking line near end of thighbone region;
The soft tissue near end of thighbone region is removed, and creates near end of thighbone threedimensional model.
Wherein, while according to the first image grayscale threshold value and preset hatching, the accuracy of image segmentation can be improved; Remove soft tissue, can influence to avoid soft tissue to skeletal reconstruction, improve the accuracy of the near end of thighbone threedimensional model of creation.
Specific software operation is as follows:User utilizes trimming mask (Crop Mask) tool in image processing software by mesh Interested region is limited in the near end of thighbone region for needing to model in logo image, reduces the mask area for needing to edit;Then, It is filled using mask editor (Edit masks) tool non-in the near end of thighbone or erasing (Erase) mask of (Draw) each level The soft tissue of target area;After the completion of the image segmentation near end of thighbone region, calculating threedimensional model function (Calculate is used The near end of thighbone threedimensional model for 3D) rebuilding high quality, is finally exported and is saved with stl formatted file.
S220, according to the second image grayscale threshold value, the target image for the femoral head for including necrotic area is split, is determined Necrotic bone region out, and create necrotic bone threedimensional model.
Second image grayscale threshold value can be distributed according to the overall pixel in necrotic bone region to be determined.Specifically, user A new mask is created in the same project manager mask of image processing software, and the second image grayscale threshold value is set;So Afterwards, the necrotic zone of each level is filled using mask edit tool, wherein do not need to distinguish bone and soft tissue.Necrotic bone covers It is finally defeated with stl formatted file using calculating threedimensional model reconstruction high quality necrotic bone threedimensional model after the completion of plate editor Out and save.
S230, surface lubrication and grid dividing are carried out to the threedimensional model of creation respectively, the threedimensional model that obtains that treated.
The initial threedimensional model rough surface rebuild through image processing software Mimics, needs to handle through surface lubrication, The accurate grid division of energy, and then three dimensional analysis is carried out, exemplary, concrete operations are as follows:
Application image handles reverse engineering software, such as Geomagic Studio, imports with the stock of the stl format saved Proximal bone threedimensional model and necrotic bone threedimensional model, enabling grid doctor (Automatically Repairing a Mesh), Quick smooth (Using QuickSmooth), spike (Removing Spikes), relaxation (Minimizing are deleted Crease Angles Between Polygons) and the tools such as sand paper (Using Sandpaper) near end of thighbone three-dimensional mould Type and necrotic bone threedimensional model carry out surface lubrication operation.Surface lubrication operate when need to pay attention to lubrication dynamics be not easy it is excessive, It is not easy too small, dynamics crosses the original shape that senior general substantially change model, and it is distorted subsequent analysis result, and dynamics is too small It will be unable to carry out analytical calculation.Therefore, the three-dimensional mould that lubrication operation is in the form of obtaining keeping model original as far as possible and surface is smooth For the purpose of type.Then, analysis software, for example, ANSYS Mechanical APDL, is directed respectively at surface lubrication Threedimensional model after reason successively performs the following operation selection in preprocessing module:Single choice type (Element type) adds Add/editor/to delete (Add/Edit/Delete), solid unit (Solid), operation selecting unit type and surface mesh to divide (Meshing), the grid dividing of model is completed.
S240, the Young's modulus obtained according to the gray value based on target image, to treated, threedimensional model carries out material Expect assignment, obtain analysis model, wherein treated near end of thighbone threedimensional model and necrotic bone threedimensional model are assigned to difference respectively Young's modulus, analysis model includes finite element analysis model.
Application image processing software, such as the FEA/CFD function of Mimics, to being drawn by surface lubrication and grid for importing / after threedimensional model carry out Young's modulus assignment.It needs femoral head Accurate Segmentation to be necrosis using traditional assignment method The accuracy of bone, cancellous bone and cortex bone etc., this result for causing simulation to calculate and initial image segmentation is closely related, thus, Traditional assignment method universality is poor.In addition, Density Distribution is uneven in bone after the necrosis of femoral head, cause even same One organizes its biomechanics characteristic, and there is also notable differences, therefore, regard entire necrotic area as homogeneous material and give same Young Modulus value directly affects the accuracy of sunykatuib analysis result.In order to overcome the problems referred above, according to gray scale (the i.e. CT of target image Value), the numerical relation of bone density and Young's modulus, the Young's modulus of corresponding region can be directly obtained by target image gray value Value, to keep the difference of the biomechanics characteristic of different zones.Target image includes patient's femur CT image.The meter of Young's modulus Calculating formula is:ρ [kg/m3]=0.0464HU+1000, E [Pa]=3790 ρ3, wherein HU indicates image CT value, and ρ indicates bone Head density, E indicate Young's modulus.
Specifically, bone density is expressed as ρ near end of thighbone model1, correspondence image CT value is expressed as HU1, Young mould Amount is expressed as E1, then ρ1[kg/m3]=0.0464HU1+ 1000,For downright bad bone model, bone Density is expressed as ρ2, correspondence image CT value is expressed as HU2, Young's modulus is expressed as E2, then ρ2[kg/m3]=0.0464HU2+ 1000,Each image-region includes multiple pixel units, therefore, HU1And HU2It corresponds respectively to A series of grey scale pixel values of near end of thighbone model and downright bad bone model, ρ1And ρ2The respectively image of near end of thighbone and necrotic bone The corresponding bone density value of pixel unit.Material is carried out near end of thighbone model and downright bad bone model respectively according to above-mentioned calculation formula Expect assignment.
In addition, needing to select the Young's modulus number pair of different number respectively according to the property difference of necrotic bone and healthy bone Downright bad bone model and near end of thighbone model carry out initialization assignment.Optionally, initialization assignment includes:Entire necrosis bone model is assigned The Young's modulus of the first preset quantity is given, entire near end of thighbone model assigns the Young's modulus of the second preset quantity.Then, herein On the basis of according to above-mentioned formula continue material assignment.It is of course also possible to completing material assignment according to above-mentioned calculation formula And then initialization assignment is carried out, the present embodiment is not particularly limited this.Material assignment can accurately guarantee to analyze result Accuracy.The material assignment method of the present embodiment effectively improves conventional images cutting techniques to bone material assigned result accuracy Influence, therefore, this assignment method is more accurate.
Optionally, after carrying out surface lubrication and grid dividing to the threedimensional model of creation respectively, the present embodiment method is also Including:
Boolean calculation is carried out near end of thighbone threedimensional model and necrotic bone threedimensional model, is obtained three-dimensional without downright bad near end of thighbone Model, wherein no downright bad near end of thighbone threedimensional model is obtained later by removing necrotic bone near end of thighbone threedimensional model.Specifically , it can use image processing software, such as Boolean calculation (Boolean Command) tool in Mimics, by entire femur Proximal end and necrotic area subtraction are obtained without downright bad near end of thighbone.
Correspondingly, the Young's modulus obtained according to the gray value based on target image, to treated, threedimensional model is carried out Material assignment, including:
According to the Young's modulus that the gray value based on target image obtains, respectively to without downright bad near end of thighbone threedimensional model and Necrotic bone threedimensional model is assigned to different Young's modulus.
Illustratively, according to above-mentioned assignment method, 10 Young's modulus E are integrally assigned to necrotic bone threedimensional model first2, 30 Young's modulus E are integrally assigned to without downright bad near end of thighbone threedimensional model1, right then according to the calculation formula of Young's modulus Necrotic bone threedimensional model and the Young's modulus assignment that corresponding region is carried out without downright bad near end of thighbone threedimensional model.
S250, according to analysis model, extract the weight bearing area of femoral head equivalent stress and total deformation quantity, and carry out data point Analysis.
The technical solution of the present embodiment is split the femoral head image for including necrotic area based on image grayscale threshold value, obtains To near end of thighbone and necrotic bone region, and create threedimensional model;Secondly, the Young's modulus obtained using the gray value based on image Material assignment for the distinguishing property of threedimensional model divided by surface lubrication and grid dividing, obtains needed for data analysis Analysis model.The present embodiment solves the problems, such as that the accuracy of femoral head image analysis result in the prior art is lower, improves The accuracy of femoral head image analysis result can provide crucial directive function, Jin Erti for collapse of femoral head forecasting research The high accuracy of collapse of femoral head prediction.
Embodiment three
Fig. 3 is the structural schematic diagram for the femoral head image analysis apparatus that the embodiment of the present invention three provides, and the present embodiment can fit The case where for analyzing femoral head image.Femoral head image analysis apparatus provided by the embodiment of the present invention is executable originally Femoral head image analysis method provided by invention any embodiment has the corresponding functional module of execution method and beneficial to effect Fruit.As shown in figure 3, the device includes threedimensional model creation module 310, analysis model creation module 320 and data analysis module 330, wherein:
Threedimensional model creation module 310 creates femoral head for the target image based on the femoral head for including necrotic area Threedimensional model, wherein threedimensional model includes near end of thighbone threedimensional model and necrotic bone threedimensional model.
Optionally, the target image in threedimensional model creation module 310 includes the femur scanned on a preset condition based Head image, wherein preset condition includes at least one following:Scanning voltage, sweep current, Pixel Dimensions, scanning support gradient, Scan slice thickness and scanning slice interval.
Analysis model creation module 320 creates the analysis model of femoral head for the threedimensional model according to femoral head.
Data analysis module 330, for according to creation analysis model, extract femoral head weight bearing area equivalent stress with Total deformation quantity, and carry out data analysis.
Optionally, data analysis module 330 includes weight bearing area's determination unit, analogue unit, data extracting unit and result Analytical unit, wherein:
Area's determination unit of bearing a heavy burden determines the weight bearing of femoral head for utilizing pre-set angle in analysis model Area;
Analogue unit simulates stress condition of the femoral head under walking states for being based on analysis model;
Data extracting unit, for according to analog result, extract the weight bearing area of femoral head equivalent stress and total deformation quantity;
Interpretation of result unit, for respectively comparing the equivalent stress of extraction and total deformation quantity with corresponding preset threshold Compared with obtaining the image analysis result of femoral head.
Based on the above technical solution, optionally, threedimensional model creation module 310 includes near end of thighbone threedimensional model Creating unit and necrotic bone threedimensional model creating unit, wherein:
Near end of thighbone threedimensional model creating unit is used for according to the first image grayscale threshold value, to the femur including necrotic area The target image of head is split, and determines near end of thighbone region, and create near end of thighbone threedimensional model;
Necrotic bone threedimensional model creating unit is used for according to the second image grayscale threshold value, to the femoral head including necrotic area Target image be split, determine necrotic bone region, and create necrotic bone threedimensional model.
Optionally, near end of thighbone threedimensional model creating unit includes:
Near end of thighbone region determines subelement, is used for according to the first image grayscale threshold value, and combines preset hatching, right The target image of femoral head including necrotic area is split, and determines near end of thighbone region, wherein hatching is near end of thighbone Element marking line in region;
Soft tissue removes subelement, for removing the soft tissue near end of thighbone region, and creates near end of thighbone three-dimensional mould Type.
Optionally, analysis model creation module 320 includes surface treatment unit and material assignment unit, wherein:
Surface treatment unit, for obtaining to the threedimensional model progress surface lubrication of creation and grid dividing, treated Threedimensional model;
Material assignment unit, the Young's modulus for being obtained according to the gray value based on target image, to treated three Dimension module carries out material assignment, obtains analysis model, wherein treated near end of thighbone threedimensional model and necrotic bone threedimensional model It is assigned to different Young's modulus respectively, analysis model includes finite element analysis model.
Optionally, analysis model creation module 320 further includes Boolean calculation unit, is used for:
Boolean calculation is carried out near end of thighbone threedimensional model and necrotic bone threedimensional model, is obtained three-dimensional without downright bad near end of thighbone Model, wherein no downright bad near end of thighbone threedimensional model is obtained later by removing necrotic bone near end of thighbone threedimensional model;
Correspondingly, material assignment unit is specifically used for:
According to the Young's modulus that the gray value based on target image obtains, respectively to without downright bad near end of thighbone threedimensional model and Necrotic bone threedimensional model is assigned to different Young's modulus.
The present embodiment technical solution is primarily based on the target image of the femoral head including necrotic area, creates the three of the femoral head Dimension module, wherein threedimensional model includes near end of thighbone threedimensional model and necrotic bone threedimensional model;Then, according to threedimensional model, wound Build the analysis model of the femoral head;Finally, extracting the equivalent stress and total deformation in the weight bearing area of femoral head according to analysis model Amount, and data analysis is carried out, solve the problems, such as that the accuracy of femoral head image analysis result in the prior art is lower, from analysis Weight bearing area stress situation after caput femoris necrosis is started with, and the accuracy of femoral head image analysis result is improved.
Example IV
Fig. 4 is a kind of structural schematic diagram for server that the embodiment of the present invention four provides.Fig. 4, which is shown, to be suitable for being used to realizing The block diagram of the exemplary servers 412 of embodiment of the present invention.The server 412 that Fig. 4 is shown is only an example, should not be right The function and use scope of the embodiment of the present invention bring any restrictions.
As shown in figure 4, server 412 is showed in the form of generic server.The component of server 412 may include but not It is limited to:One or more processor 416, storage device 428 connect different system components (including storage device 428 and processing Device 416) bus 418.
Bus 418 indicates one of a few class bus structures or a variety of, including storage device bus or storage device control Device processed, peripheral bus, graphics acceleration port, processor or total using the local of any bus structures in a variety of bus structures Line.For example, these architectures include but is not limited to industry standard architecture (Industry Subversive Alliance, ISA) bus, microchannel architecture (Micro Channel Architecture, MAC) bus is enhanced Isa bus, Video Electronics Standards Association (Video Electronics Standards Association, VESA) local are total Line and peripheral component interconnection (Peripheral Component Interconnect, PCI) bus.
Server 412 typically comprises a variety of computer system readable media.These media can be it is any being capable of bedding and clothing The usable medium that business device 412 accesses, including volatile and non-volatile media, moveable and immovable medium.
Storage device 428 may include the computer system readable media of form of volatile memory, such as arbitrary access Memory (Random Access Memory, RAM) 430 and/or cache memory 432.Server 412 can be further Including other removable/nonremovable, volatile/non-volatile computer system storage mediums.Only as an example, storage system System 434 can be used for reading and writing immovable, non-volatile magnetic media (Fig. 4 do not show, commonly referred to as " hard disk drive ").To the greatest extent It is not shown in pipe Fig. 4, the disc driver for reading and writing to removable non-volatile magnetic disk (such as " floppy disk ") can be provided, with And to removable anonvolatile optical disk, such as CD-ROM (Compact Disc Read-Only Memory, CD-ROM), number Word optic disk (Digital Video Disc-Read Only Memory, DVD-ROM) or other optical mediums) read-write CD Driver.In these cases, each driver can be connected by one or more data media interfaces with bus 418. Storage device 428 may include at least one program product, which has one group of (for example, at least one) program module, These program modules are configured to perform the function of various embodiments of the present invention.
Program/utility 440 with one group of (at least one) program module 442 can store in such as storage dress It sets in 428, such program module 442 includes but is not limited to operating system, one or more application program, other program moulds It may include the realization of network environment in block and program data, each of these examples or certain combination.Program module 442 usually execute function and/or method in embodiment described in the invention.
Server 412 can also be with one or more external equipments 414 (such as keyboard, direction terminal, display 424 etc.) Communication, can also be enabled a user to one or more terminal interact with the server 412 communicate, and/or with make the clothes Any terminal (such as network interface card, modem etc.) that business device 412 can be communicated with one or more of the other computing terminal Communication.This communication can be carried out by input/output (I/O) interface 422.Also, server 412 can also be suitable by network Orchestration 420 and one or more network (such as local area network (Local Area Network, LAN), wide area network (Wide Area Network, WAN) and/or public network, such as internet) communication.As shown in figure 4, network adapter 420 passes through bus 418 It is communicated with other modules of server 412.It should be understood that although not shown in the drawings, can be used in conjunction with server 412 other hard Part and/or software module, including but not limited to:Microcode, terminal driver, redundant processor, external disk drive array, magnetic Disk array (Redundant Arrays of Independent Disks, RAID) system, tape drive and data backup Storage system etc..
The program that processor 416 is stored in storage device 428 by operation, thereby executing various function application and number According to processing, such as realize femoral head image analysis method provided by the embodiment of the present invention, this method includes:
Based on the target image for the femoral head for including necrotic area, the threedimensional model of the femoral head is created, wherein described three Dimension module includes near end of thighbone threedimensional model and necrotic bone threedimensional model;
According to the threedimensional model, the analysis model of the femoral head is created;
According to the analysis model, extract the weight bearing area of the femoral head equivalent stress and total deformation quantity, line number of going forward side by side According to analysis.
Embodiment five
The embodiment of the present invention five additionally provides a kind of computer readable storage medium, is stored thereon with computer program, should The femoral head image analysis method as provided by the embodiment of the present invention is realized when program is executed by processor, this method includes:
Based on the target image for the femoral head for including necrotic area, the threedimensional model of the femoral head is created, wherein described three Dimension module includes near end of thighbone threedimensional model and necrotic bone threedimensional model;
According to the threedimensional model, the analysis model of the femoral head is created;
According to the analysis model, extract the weight bearing area of the femoral head equivalent stress and total deformation quantity, line number of going forward side by side According to analysis.
The computer storage medium of the embodiment of the present invention, can be using any of one or more computer-readable media Combination.Computer-readable medium can be computer-readable signal media or computer readable storage medium.It is computer-readable Storage medium for example may be-but not limited to-the system of electricity, magnetic, optical, electromagnetic, infrared ray or semiconductor, device or Device, or any above combination.The more specific example (non exhaustive list) of computer readable storage medium includes:Tool There are electrical connection, the portable computer diskette, hard disk, random access memory (RAM), read-only memory of one or more conducting wires (ROM), erasable programmable read only memory (EPROM or flash memory), optical fiber, portable compact disc read-only memory (CD- ROM), light storage device, magnetic memory device or above-mentioned any appropriate combination.In this document, computer-readable storage Medium can be any tangible medium for including or store program, which can be commanded execution system, device or device Using or it is in connection.
Computer-readable signal media may include in a base band or as carrier wave a part propagate data-signal, Wherein carry computer-readable program code.The data-signal of this propagation can take various forms, including but unlimited In electromagnetic signal, optical signal or above-mentioned any appropriate combination.Computer-readable signal media can also be that computer can Any computer-readable medium other than storage medium is read, which can send, propagates or transmit and be used for By the use of instruction execution system, device or device or program in connection.
The program code for including on computer-readable medium can transmit with any suitable medium, including --- but it is unlimited In wireless, electric wire, optical cable, RF etc. or above-mentioned any appropriate combination.
The computer for executing operation of the present invention can be write with one or more programming languages or combinations thereof Program code, described program design language include object oriented program language-such as Java, Smalltalk, C++, It further include conventional procedural programming language-such as " C " language or similar programming language.Program code can be with It fully executes, partly execute on the user computer on the user computer, being executed as an independent software package, portion Divide and partially executes or executed on remote computer or terminal completely on the remote computer on the user computer.It is relating to And in the situation of remote computer, remote computer can pass through the network of any kind --- including local area network (LAN) or extensively Domain net (WAN)-be connected to subscriber computer, or, it may be connected to outer computer (such as provided using Internet service Quotient is connected by internet).
Note that the above is only a better embodiment of the present invention and the applied technical principle.It will be appreciated by those skilled in the art that The invention is not limited to the specific embodiments described herein, be able to carry out for a person skilled in the art it is various it is apparent variation, It readjusts and substitutes without departing from protection scope of the present invention.Therefore, although being carried out by above embodiments to the present invention It is described in further detail, but the present invention is not limited to the above embodiments only, without departing from the inventive concept, also It may include more other equivalent embodiments, and the scope of the invention is determined by the scope of the appended claims.

Claims (10)

1. a kind of femoral head image analysis method, which is characterized in that including:
Based on the target image for the femoral head for including necrotic area, the threedimensional model of the femoral head is created, wherein the three-dimensional mould Type includes near end of thighbone threedimensional model and necrotic bone threedimensional model;
According to the threedimensional model, the analysis model of the femoral head is created;
According to the analysis model, extract the weight bearing area of the femoral head equivalent stress and total deformation quantity, and carry out data point Analysis.
2. the method according to claim 1, wherein target image includes the stock scanned on a preset condition based Bone image, wherein the preset condition includes at least one following:Scanning voltage, sweep current, Pixel Dimensions, scanning support Gradient, Scan slice thickness and scanning slice interval.
3. the method according to claim 1, wherein creating point of the femoral head according to the threedimensional model Model is analysed, including:
Surface lubrication and grid dividing are carried out to the threedimensional model, the threedimensional model that obtains that treated;
According to the Young's modulus that the gray value based on target image obtains, treated that threedimensional model carries out material tax to described Value, obtains the analysis model, wherein treated near end of thighbone threedimensional model and necrotic bone threedimensional model are assigned to difference respectively Young's modulus, the analysis model includes finite element analysis model.
4. according to the method described in claim 3, it is characterized in that, carrying out surface lubrication and grid dividing to the threedimensional model Later, the method also includes:
Boolean calculation is carried out to the near end of thighbone threedimensional model and necrotic bone threedimensional model, is obtained three-dimensional without downright bad near end of thighbone Model, wherein it is described without downright bad near end of thighbone threedimensional model be by the near end of thighbone threedimensional model remove necrotic bone after It obtains;
Correspondingly, the Young's modulus that the basis is obtained based on the gray value of target image, to treated the threedimensional model Material assignment is carried out, including:
According to the Young's modulus that the gray value based on target image obtains, respectively to described without downright bad near end of thighbone three-dimensional mould Type and necrotic bone threedimensional model are assigned to different Young's modulus.
5. the method according to claim 1, wherein the target figure based on the femoral head for including necrotic area Picture creates the threedimensional model of the femoral head, including:
According to the first image grayscale threshold value, target image is split, determines near end of thighbone region, and create near end of thighbone Threedimensional model;
According to the second image grayscale threshold value, target image is split, determines necrotic bone region, and creates necrotic bone three-dimensional Model.
6. according to the method described in claim 5, it is characterized in that, described according to the first image grayscale threshold value, to target image It is split, determines near end of thighbone region, and create near end of thighbone threedimensional model, including:
According to the first image grayscale threshold value, and preset hatching is combined, target image is split, determines near end of thighbone Region, wherein the hatching is the element marking line in the near end of thighbone region;
The soft tissue in the near end of thighbone region is removed, and creates near end of thighbone threedimensional model.
7. the method according to claim 1, wherein extracting the negative of the femoral head according to the analysis model The equivalent stress in weight area and total deformation quantity, and data analysis is carried out, including:
Using pre-set angle, the weight bearing area of the femoral head is determined in the analysis model;
Based on the analysis model, stress condition of the femoral head under walking states is simulated;
According to analog result, the equivalent stress in the extraction weight bearing area and total deformation quantity;
The equivalent stress and total deformation quantity are compared with corresponding preset threshold respectively, obtain the image of the femoral head Analyze result.
8. a kind of femoral head image analysis apparatus, which is characterized in that including:
Threedimensional model creation module creates the three of the femoral head for the target image based on the femoral head for including necrotic area Dimension module, wherein the threedimensional model includes near end of thighbone threedimensional model and necrotic bone threedimensional model;
Analysis model creation module, for creating the analysis model of the femoral head according to the threedimensional model;
Data analysis module, for extracting the equivalent stress and total shape in the weight bearing area of the femoral head according to the analysis model Variable, and carry out data analysis.
9. a kind of server, which is characterized in that including:
One or more processors;
Storage device, for storing one or more programs,
When one or more of programs are executed by one or more of processors, so that one or more of processors are real The now femoral head image analysis method as described in any in claim 1~7.
10. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the program is by processor The femoral head image analysis method as described in any in claim 1~7 is realized when execution.
CN201810668292.4A 2018-06-26 2018-06-26 Femoral head image analysis method, device, server and medium Pending CN108921832A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810668292.4A CN108921832A (en) 2018-06-26 2018-06-26 Femoral head image analysis method, device, server and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810668292.4A CN108921832A (en) 2018-06-26 2018-06-26 Femoral head image analysis method, device, server and medium

Publications (1)

Publication Number Publication Date
CN108921832A true CN108921832A (en) 2018-11-30

Family

ID=64421278

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810668292.4A Pending CN108921832A (en) 2018-06-26 2018-06-26 Femoral head image analysis method, device, server and medium

Country Status (1)

Country Link
CN (1) CN108921832A (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110660137A (en) * 2019-09-02 2020-01-07 北京工业大学 Microscopic biomechanical research method for evaluating femoral fracture risk under different-angle lateral falling impact loads of human body
CN110680353A (en) * 2019-09-17 2020-01-14 中国人民解放军总医院第四医学中心 Analysis device and method for femoral head integrity rate
CN110728029A (en) * 2019-09-17 2020-01-24 中国人民解放军总医院第四医学中心 Femur integrity analysis system and femur integrity analysis model construction method
CN110811829A (en) * 2019-11-06 2020-02-21 中国人民解放军总医院第四医学中心 Construction method and system based on femoral rotation axis and varus analysis model
CN111839570A (en) * 2020-08-12 2020-10-30 北京中医药大学第三附属医院 Device for predicting risk of femoral head collapse and using method and application thereof
CN112712594A (en) * 2020-12-29 2021-04-27 深圳技术大学 Skeleton model construction method and device, computer equipment and storage medium
CN115131276A (en) * 2022-03-03 2022-09-30 中国人民解放军总医院第四医学中心 Method, device and equipment for acquiring stress distribution of skeleton and storage medium

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462636A (en) * 2014-09-30 2015-03-25 西安科技大学 Modeling method of necrosis caput femoris restoring model based on umbrella-shaped caput femoris supporter

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462636A (en) * 2014-09-30 2015-03-25 西安科技大学 Modeling method of necrosis caput femoris restoring model based on umbrella-shaped caput femoris supporter

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
何海军: "基于CT图像建立可供动态力学模拟分析的股骨头坏死有限元模型", 《中国中医骨伤科杂志》, vol. 24, no. 9, 30 September 2016 (2016-09-30), pages 5 - 10 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110660137A (en) * 2019-09-02 2020-01-07 北京工业大学 Microscopic biomechanical research method for evaluating femoral fracture risk under different-angle lateral falling impact loads of human body
CN110680353A (en) * 2019-09-17 2020-01-14 中国人民解放军总医院第四医学中心 Analysis device and method for femoral head integrity rate
CN110728029A (en) * 2019-09-17 2020-01-24 中国人民解放军总医院第四医学中心 Femur integrity analysis system and femur integrity analysis model construction method
CN110811829A (en) * 2019-11-06 2020-02-21 中国人民解放军总医院第四医学中心 Construction method and system based on femoral rotation axis and varus analysis model
CN112451090A (en) * 2019-11-06 2021-03-09 中国人民解放军总医院第四医学中心 Construction system of analysis model
CN112451090B (en) * 2019-11-06 2022-02-18 中国人民解放军总医院第四医学中心 Construction system of analysis model
CN111839570A (en) * 2020-08-12 2020-10-30 北京中医药大学第三附属医院 Device for predicting risk of femoral head collapse and using method and application thereof
CN111839570B (en) * 2020-08-12 2024-03-29 北京中医药大学第三附属医院 Device for predicting femoral head collapse risk and use method and application thereof
CN112712594A (en) * 2020-12-29 2021-04-27 深圳技术大学 Skeleton model construction method and device, computer equipment and storage medium
CN115131276A (en) * 2022-03-03 2022-09-30 中国人民解放军总医院第四医学中心 Method, device and equipment for acquiring stress distribution of skeleton and storage medium

Similar Documents

Publication Publication Date Title
CN108921832A (en) Femoral head image analysis method, device, server and medium
US11869202B2 (en) Method and system for processing multi-modality image
JP7346553B2 (en) Determining the growth rate of objects in a 3D dataset using deep learning
Veneziano et al. Surface smoothing, decimation, and their effects on 3D biological specimens
Park et al. Automatic reconstruction of the arterial and venous trees on volumetric chest CT
US20180012363A1 (en) Preview visualisation of tracked nerve fibers
US10013770B2 (en) Dynamic definition of a region of interest for tracking nerve fibers
Yosibash et al. Autonomous FEs (AFE)-A stride toward personalized medicine
Zhai et al. Automatic quantitative analysis of pulmonary vascular morphology in CT images
WO2016116136A1 (en) Atlas-based determination of tumour growth direction
Caprara et al. Automated pipeline to generate anatomically accurate patient-specific biomechanical models of healthy and pathological FSUs
CN112638262A (en) Similarity determination device, method, and program
US10769240B2 (en) Determining medical outcome quality
Schendel et al. A web-based, integrated simulation system for craniofacial surgical planning
Yu et al. Quantification of the bone lacunocanalicular network from 3D X‐ray phase nanotomography images
CN110728029B (en) Femur integrity analysis system and femur integrity analysis model construction method
KR20220106113A (en) Systems and methods for reconstruction and characterization of physiologically healthy and physiologically defective anatomy to facilitate preoperative surgical planning
US20230207106A1 (en) Image segmentation for sets of objects
JP5586953B2 (en) Access to medical image database using anatomical shape information
Chae et al. Medical image segmentation for mobile electronic patient charts using numerical modeling of IoT
EP3794550B1 (en) Comparison of a region of interest along a time series of images
Koutkalaki et al. Towards a foot bio-model for performing finite element analysis for footwear design optimization using a Cloud infrastructure
Plotkowiak et al. High performance computer simulations of cardiac electrical function based on high resolution MRI datasets
Ramme et al. Gaussian curvature analysis allows for automatic block placement in multi-block hexahedral meshing
Krol et al. Computer-aided osteotomy design for harvesting autologous bone grafts in reconstructive surgery

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination