IMAGE DATA PROCESSING SYSTEMS

• US 20120224758A1
• Filed: 10/06/2010
• Published: 09/06/2012
• Est. Priority Date: 10/07/2009
• Status: Active Grant

• Current Assignee: Cambridge Enterprises Limited
• Original Assignee: Cambridge Enterprises Limited
• Inventors: Treece, Graham Michael, Poole, Kenneth
• Granted Patent: US 8,938,109 B2
• Time in Patent Office: Days
• US Class Current: 382/131
• CPC Class Codes: G06T 2200/24 involving graphical user interfaces [GUIs]
• CPC Class Codes: G06T 2207/10081 Computed x-ray tomography [CT]
• CPC Class Codes: G06T 2207/20108 Interactive selection of 2D slice in a 3D data set
• CPC Class Codes: G06T 2207/20168 Radial search
• CPC Class Codes: G06T 2207/30008 Bone
• CPC Class Codes: G06T 7/0012 Biomedical image inspection
• CPC Class Codes: G06T 7/12 Edge-based segmentation
• CPC Class Codes: G06T 7/149 involving deformable models, e.g. active contour models
• CPC Class Codes: G06T 7/62 of area, perimeter, diameter or volume

Abstract

We describe a method of estimating the thickness of a tissue structure, in particular cortical bone thickness from tomographic imaging data such as CT data. The method models the tissue structure as a variation of an imaging parameter along a line; models a variation of the tomographic imaging data along the line as a blurred version of the variation of the imaging parameter modelling the; and fits the blurred tomographic imaging model to data from the tomographic imaging by holding a tissue modelling parameter at a constraining value and allowing variation of a blurring parameter and at least one parameter defining the thickness of the tissue structure to determine an estimate of the thickness of the structure.

Claims

1. A method of determining data defining an estimate of the thickness of a tissue structure from tomographic imaging data, the method comprising:

inputting a tomographic imaging data set defining at least one slice through a region of tissue including said tissue structure;

modelling said tissue structure as a variation of an imaging parameter of said imaging data along a line crossing said tissue structure, a tissue structure model defined by said variation having at least a first tissue modelling parameter defining a value of said imaging parameter along said line and within said tissue structure and having at least one parameter defining said thickness of said tissue structure;

modelling a variation of said tomographic imaging data along said line as a blurred version of said variation of said imaging parameter modelling said tissue structure to provide a blurred tomographic imaging model, said blurred version of said variation of said imaging parameter being defined by a blurring function applied to said tissue structure model, wherein said blurred tomographic imaging model has a blurring parameter representing a degree of blurring applied by said blurring function;

determining a constraining value of said first tissue modelling parameter; and

fitting said blurred tomographic imaging model defining said blurred version of said variation of said imaging parameter along said line, to data derived from said tomographic imaging data set;

wherein said fitting comprises holding said first tissue modelling parameter substantially at said constraining value and allowing variation of said blurring parameter and said at least one parameter defining said thickness of said tissue structure to determine said at least one parameter defining said thickness of said tissue structure; and

wherein said at least one parameter defining said thickness of said tissue structure comprises said data defining said estimate of said thickness of said tissue structure.

2. A method as claimed in claim 1 wherein said tomographic imaging data comprises X-ray computed tomography data, wherein said tissue structure comprises cortical bone, and wherein said imaging parameter comprises a measure of bone density.

3. A method as claimed in claim 1 wherein said tissue structure model has second and third tissue modelling parameters defining values of said imaging parameter along said line outside said tissue structure respectively to either side of said tissue structure, wherein said at least one parameter defining said thickness of said tissue structure comprises two position values defining positions of boundaries along said line between regions where said imaging parameter has values defined by said first tissue modelling parameter and by said second and third tissue modelling parameters respectively, and wherein said fitting comprises determining values for said second and third modelling parameters.

4. A method as claimed in claim 1, wherein said line lies out of a plane defined by said at least one slice through said region of tissue, wherein said blurring parameter comprises an in-plane blurring parameter defining a degree of blurring in said plane of said slice, and wherein said blurred tomographic imaging model further comprises an out-of-plane blurring parameter defining a degree of blurring dependent on an angle between said line and said plane of said slice.

5. A method as claimed in claim 4 further comprising determining said angle between said line and a normal to said plane of said slice and constraining said out-of-place blurring parameter dependent on said determined angle during said fitting.

6. A method as claimed in claim 1 further comprising identifying a 3-D outer surface of said tissue structure from said tomographic imaging data set, identifying at least one surface normal to said outer surface, and aligning said line along a direction of said surface normal.

7. A method as claimed in claim 6 comprising identifying a plurality of said surface normals and, repeating said fitting of said blurred tomographic imaging model along a plurality of said lines each defined by one of said surface normals to determine a plurality of said estimates of thickness of said tissue structure, and storing or outputting thickness map data defining a map of said estimated thickness conforming to said outer surface of said tissue structure.

8. A method as claimed in claim 7 further comprising displaying said map to a user as a 3-D representation of said outer surface of said tissue structure painted with a colour or greyscale or patterned to represent variations of said estimated thickness over said outer surface.

9. A method as claimed in claim 1 as claimed in any preceding claim used to locate the position of one or both of an inner boundary and an outer boundary of said tissue structure, wherein said at least one parameter defining said thickness of said tissue structure comprises positions of said inner and outer boundaries.

10. A method as claimed in claim 1 as claimed in any preceding claim wherein said tissue structure comprises a layer of tissue, and wherein said tissue structure model defines a substantially constant value of said first tissue modelling parameter within said layer of tissue along said line, such that said model assumes a substantially constant value of said first tissue modelling parameter within a thickness of said layer of tissue.

11. A carrier carrying processor control code to, when running, implement the method of any preceding claim 1.

12. Tissue structure thickness estimation apparatus for determining data defining an estimate of the thickness of a tissue structure from tomographic imaging data, the apparatus comprising:

an input for inputting a tomographic imaging data set defining at least one slice through a region of tissue including said tissue structure;

working memory for storing data;

program memory storing processor control code; and

a processor coupled to said input, to said working memory, and to said program memory; and

wherein said stored processor control code comprises code to;

model said tissue structure as a variation of an imaging parameter of said imaging data along a line crossing said tissue structure, a tissue structure model defined by said variation having at least a first tissue modelling parameter defining a value of said imaging parameter along said line and within said tissue structure and having at least one parameter defining said thickness of said tissue structure;

model a variation of said tomographic imaging data along said line as a blurred version of said variation of said imaging parameter modelling said tissue structure to provide a blurred tomographic imaging model, said blurred version of said variation of said imaging parameter being defined by a blurring function applied to said tissue structure model, wherein said blurred tomographic imaging model has a blurring parameter representing a degree of blurring applied by said blurring function;

determine a constraining value of said first tissue modelling parameter; and

said blurred tomographic imaging model defining said blurred version of said variation of said imaging parameter along said line, to data derived from said tomographic imaging data set;

wherein said fitting comprises holding said first tissue modelling parameter substantially at said constraining value and allowing variation of said blurring parameter and said at least one parameter defining said thickness of said tissue structure to determine said at least one parameter defining said thickness of said tissue structure; and

wherein said at least one parameter defining said thickness of said tissue structure comprises said data defining said estimate of said thickness of said tissue structure.

13. A method of estimating the position of an edge of or the thickness of cortical bone from CT data, the method comprising:

modelling a CT image of said cortical bone using a piecewise constant bone density model blurred by a CT imaging blurring function;

fitting said CT data to said modelled blurred CT image by allowing said blurring function to vary whilst constraining a value of cortical bone density represented by said bone density model to be substantially fixed; and

estimating said position of said edge or said thickness of said cortical bone from a parameter of said fitted model.

14. A method as claimed in claim 13 further comprising estimating a fixed value for said cortical bone density represented by said bone density model from said CT data prior to said fitting.

15. A method as claimed in claim 13 wherein said piecewise constant bone density model is blurred by an in-plane blurring function and an out-of-plane blurring function, said in-plane blurring function representing blurring in a planar slice of said CT data, said out-of-plane representing blurring resulting from an angle between said planar slice and a normal to a surface of said cortical bone at a point where said CT data is fitted being less than 90 degrees.

16. A method as claimed in claim 13, comprising estimating said position of said edge thickness of said cortical bone at a plurality of positions, around a 2-D perimeter or over a 3-D surface of said cortical bone by, at each said position;

determining a perimeter normal or surface normal and performing said modelling and fitting along a direction defined by said perimeter normal or surface normal at said position.

17. A carrier carrying processor control code to, when running, implement the method of claim 13.

18. Apparatus for estimating the position of an edge of or the thickness of cortical bone from CT data, the apparatus comprising:

working memory for storing data;

program memory storing processor control code; and

a processor coupled to said input, to said working memory, and to said program memory; and

wherein said stored processor control code comprises code to;

model a CT image of said cortical bone using a piecewise constant bone density model blurred by a CT imaging blurring function;

fit said CT data to said modelled blurred CT image by allowing said blurring function to vary whilst constraining a value of cortical bone density represented by said bone density model to be substantially fixed; and

estimate said position of said edge or said thickness of said cortical bone from a parameter of said fitted model.

19. A method of providing a representation of cortical bone thickness from CT data, the method comprising:

inputting CT image data of said cortical bone;

generating a 2-D or 3-D needle map of an outer perimeter or surface of said cortical bone, said needle map defining a plurality of perimeter normals or surface normals to said cortical bone at a plurality of respective portions on said outer perimeter or surface;

determining a thickness estimate of said cortical bone at each of said positions along a direction of a respective said perimeter normal or surface normal for the position; and

generating mapping data representing a 2-D or 3-D map defining said estimated cortical thickness at each of said plurality of positions overlaid on a respective 2-D representation of said perimeter or 3-D representation of said surface of said cortical bone.

20. A method as claimed in claim 19 comprising displaying said map to a user as a representation of said cortical bone painted with colour or greyscale or patterned to represent variations of said estimated thickness over said cortical bone.

21. Apparatus for providing a representation of cortical bone thickness from CT data, the apparatus comprising:

working memory for storing data;

program memory storing processor control code; and

a processor coupled to said input, to said working memory, and to said program memory; and

wherein said stored processor control code comprises code to;

input CT image data of said cortical bone;

generate a 2-D or 3-D needle map of an outer perimeter or surface of said cortical bone, said needle map defining a plurality of perimeter normals or surface normals to said cortical bone at a plurality of respective portions on said outer perimeter or surface;

determine a thickness estimate of said cortical bone at each of said positions along a direction of a respective said perimeter normal or surface normal for the position; and

generate mapping data representing a 2-D or 3-D map defining said estimated cortical thickness at each of said plurality of position overlaid on a respective 2-D representation of said perimeter or 3-D representation of said surface of said cortical bone.