WO2015138397A1 - Methods and systems for identifying anatomical landmarks in image data - Google Patents

Methods and systems for identifying anatomical landmarks in image data Download PDF

Info

Publication number
WO2015138397A1
WO2015138397A1 PCT/US2015/019615 US2015019615W WO2015138397A1 WO 2015138397 A1 WO2015138397 A1 WO 2015138397A1 US 2015019615 W US2015019615 W US 2015019615W WO 2015138397 A1 WO2015138397 A1 WO 2015138397A1
Authority
WO
WIPO (PCT)
Prior art keywords
reformatted
image data
oblique plane
tracing
soft tissue
Prior art date
Application number
PCT/US2015/019615
Other languages
French (fr)
Inventor
Yangqiu Hu
Angela B. MINES
Nathaniel M. Lenz
Original Assignee
Smith & Nephew, Inc.
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 Smith & Nephew, Inc. filed Critical Smith & Nephew, Inc.
Priority to US15/124,988 priority Critical patent/US10235757B2/en
Priority to EP15761756.4A priority patent/EP3116389B1/en
Publication of WO2015138397A1 publication Critical patent/WO2015138397A1/en
Priority to US16/152,211 priority patent/US10275884B2/en

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
    • G06T7/0014Biomedical image inspection using an image reference approach
    • G06T7/0016Biomedical image inspection using an image reference approach involving temporal comparison
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/05Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves 
    • A61B5/055Detecting, measuring or recording for diagnosis by means of electric currents or magnetic fields; Measuring using microwaves or radio waves  involving electronic [EMR] or nuclear [NMR] magnetic resonance, e.g. magnetic resonance imaging
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/45For evaluating or diagnosing the musculoskeletal system or teeth
    • A61B5/4504Bones
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/45For evaluating or diagnosing the musculoskeletal system or teeth
    • A61B5/4523Tendons
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/45For evaluating or diagnosing the musculoskeletal system or teeth
    • A61B5/4533Ligaments
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/48Other medical applications
    • A61B5/4887Locating particular structures in or on the body
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/742Details of notification to user or communication with user or patient ; user input means using visual displays
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/46Arrangements for interfacing with the operator or the patient
    • A61B6/461Displaying means of special interest
    • A61B6/466Displaying means of special interest adapted to display 3D data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/505Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of bone
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5217Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data extracting a diagnostic or physiological parameter from medical diagnostic data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis
    • A61B6/5211Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
    • A61B6/5223Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data generating planar views from image data, e.g. extracting a coronal view from a 3D image
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0875Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of bone
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/46Ultrasonic, sonic or infrasonic diagnostic devices with special arrangements for interfacing with the operator or the patient
    • A61B8/461Displaying means of special interest
    • A61B8/466Displaying means of special interest adapted to display 3D data
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/483Diagnostic techniques involving the acquisition of a 3D volume of data
    • 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
    • G06T19/00Manipulating 3D models or images for computer graphics
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T19/00Manipulating 3D models or images for computer graphics
    • G06T19/20Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/60Analysis of geometric attributes
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/70Determining position or orientation of objects or cameras
    • G06T7/73Determining position or orientation of objects or cameras using feature-based methods
    • G06T7/75Determining position or orientation of objects or cameras using feature-based methods involving models
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16HHEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
    • G16H50/00ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics
    • G16H50/30ICT specially adapted for medical diagnosis, medical simulation or medical data mining; ICT specially adapted for detecting, monitoring or modelling epidemics or pandemics for calculating health indices; for individual health risk assessment
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B2576/00Medical imaging apparatus involving image processing or analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/45For evaluating or diagnosing the musculoskeletal system or teeth
    • A61B5/4528Joints
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0833Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures
    • A61B8/085Detecting organic movements or changes, e.g. tumours, cysts, swellings involving detecting or locating foreign bodies or organic structures for locating body or organic structures, e.g. tumours, calculi, blood vessels, nodules
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5223Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for extracting a diagnostic or physiological parameter from medical diagnostic data
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • 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
    • 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/30204Marker
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2210/00Indexing scheme for image generation or computer graphics
    • G06T2210/32Image data format
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/004Annotating, labelling
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/008Cut plane or projection plane definition
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2219/00Indexing scheme for manipulating 3D models or images for computer graphics
    • G06T2219/028Multiple view windows (top-side-front-sagittal-orthogonal)

Definitions

  • Ligament attachment points can be used to facilitate consistent placement of implants on a patient's bone and may offer surgeons a reliable frame of reference for properly orienting implants.
  • the location of ligament attachment points may not be known by a surgeon before a surgical operation and may need to be assessed intraoperatively.
  • Noninvasive imaging techniques such as MRI scans and CT scans, allow image data representing a patient's joint to be collected in advance of a surgical procedure.
  • it can be difficult to use the raw image data collected from preoperative imaging to provide a surgeon with accurate preoperative information about the location of ligament attachment points.
  • the systems, devices, and methods include receiving image data, identifying a base oblique plane that is oblique to the orthogonal planes in the image data, generating an adjusted set of orthogonal planes based on the base oblique plane, and generating reformatted images sampled along the adjusted set of orthogonal planes.
  • These systems, devices, and methods may provide an oblique plane that more adequately contains an anatomical landmark such as a ligament attachment point to the tibia.
  • the method may be used to identify where the anterior cruciate ligament (ACL), medial collateral ligament (MCL) posterior cruciate ligament (PCL), or patellar tendon attaches to the tibia.
  • ACL anterior cruciate ligament
  • MCL medial collateral ligament
  • PCL posterior cruciate ligament
  • patellar tendon attaches to the tibia.
  • the systems, devices, and methods may allow tracing of an anatomical landmark in the reformatted image to generate a 3-D marking on a 3-D surface model of a patient's bone.
  • the attachment points may be useful landmarks for patient-matched instrumentation.
  • the systems, devices, and methods may also facilitate surgical planning.
  • the ligament locations may be used for preoperative planning and for designing patient-matched cutting blocks.
  • a method for locating anatomical landmarks in image data includes receiving, at a processor, image data representative of a patient's joint, identifying an oblique plane intersecting a soft tissue represented in the image data and normal to a longitudinal axis of the soft tissue, identifying, with the processor, a set of planes parallel or orthogonal to the oblique plane, generating, with the processor, a set of reformatted images, from the image data, sampled along the set of planes, tracing a feature of the soft tissue in the set of reformatted images, and generating, with the processor, a 3-D marking representing the traced feature on a 3-D model.
  • the image data is raw image data.
  • Tracing the feature may include tracing the feature in two or more reformatted images parallel to the oblique plane. Tracing the feature may include advancing through parallel reformatted images selected from the set of reformatted images. Tracing the feature may include displaying a first reformatted image on a screen and displaying a second reformatted image on the screen, wherein the second reformatted image is parallel to and offset from the first reformatted image. Tracing the feature may include outlining a cross-section of the feature in a reformatted image from the set of reformatted images. Tracing the feature may include storing data representing the tracing in a matrix.
  • generating the set of reformatted images includes rotating the image data.
  • Identifying the oblique plane may include identifying a first point and a second point along the longitudinal axis of the soft tissue, evaluating, with the processor, a line defined by the first and second points, and evaluating, with the processor, a plane normal to the line.
  • the first point may be located about where the soft tissue attaches to a patient's bone.
  • the second point may be located at about the center of the soft tissue.
  • the first point may be located at about the start of the soft tissue and the second point may be located at about the end of the soft tissue.
  • the longitudinal axis of the soft tissue may be curved.
  • identifying the oblique plane includes drawing a curve substantially parallel to the longitudinal axis of the tissue, and evaluating, with the processor, an oblique plane normal to the curve at the point where the oblique plane intersects the curve.
  • a plurality of reformatted images is generated, each reformatted image being normal to the curve at the point where each reformatted image intersects the curve.
  • Generating the 3-D marking may include generating a 3-D surface model from the matrix.
  • the method further includes smoothing the 3-D marking.
  • the 3-D model may be a model of a patient's bone.
  • the 3-D model may be a model of a patient's tibia.
  • the 3-D model may be displayed with the image data to orient a user.
  • generating the set of reformatted images includes identifying, with the processor, four corners of the oblique plane, generating, with the processor, texture coordinates using the corner locations, and interpolating, with the processor, pixel intensities for the texture coordinates using the image data.
  • the soft tissue may be an ACL, a MCL, a PCL, or a patellar tendon of a patient.
  • a system for locating anatomical landmarks in image data includes means for receiving image data representative of a patient's joint, means for identifying an oblique plane intersecting a soft tissue represented in the image data and normal to a longitudinal axis of the soft tissue, means for identifying a set of planes parallel or orthogonal to the oblique plane, means for generating a set of reformatted images, from the image data, sampled along the set of planes, means for tracing a feature of the soft tissue in the set of reformatted images, and means for generating a 3-D marking representing the traced feature on a 3-D model.
  • the image data is raw image data.
  • the means for tracing the feature may include means for tracing the feature in two or more reformatted images parallel to the oblique plane.
  • the means for tracing the feature may include means for advancing through parallel reformatted images selected from the set of reformatted images.
  • the means for tracing the feature includes means for outlining a cross-section of the feature in a reformatted image from the set of reformatted images.
  • the means for tracing the feature may include means for storing data representing the tracing in a matrix.
  • the means for generating the set of reformatted images may include means for rotating the image data.
  • the means for identifying the oblique plane may include means for identifying a first point and a second point along the longitudinal axis of the soft tissue, means for evaluating a line defined by the first and second points, and means for evaluating a plane normal to the line.
  • the first point may be located about where the soft tissue attaches to a patient's bone.
  • the second point may be located at about the center of the soft tissue.
  • the first point may be located at about the start of the soft tissue and the second point may be located at about the end of the soft tissue.
  • the longitudinal axis of the soft tissue may be curved.
  • the means for identifying the oblique plane includes means for drawing a curve substantially parallel to the longitudinal axis of the tissue, and means for evaluating an oblique plane normal to the curve at the point where the oblique plane intersects the curve.
  • a plurality of reformatted images is generated, each reformatted image being normal to the curve at the point where each reformatted image intersects the curve.
  • the means for generating the 3-D marking may include means for generating a 3-D surface model from the matrix.
  • system further includes means for smoothing the 3-D marking.
  • the 3-D model is a model of a patient's bone.
  • the 3-D model may be a model of a patient's tibia.
  • the 3-D model may be displayed with the image data to orient a user.
  • the means for generating the set of reformatted images includes means for computing four corners of the oblique plane, means for generating texture coordinates using the corner locations, and means for interpolating pixel intensities for the texture coordinates using the image data.
  • the soft tissue is an ACL, a MCL, a PCL, or a patellar tendon of a patient.
  • a method for locating anatomical landmarks in image data includes receiving, at a processor, image data representative of a patient's joint, identifying an oblique plane intersecting a soft tissue represented in the image data and normal to a longitudinal axis of the soft tissue, generating, with the processor, a reformatted image, from the image data, sampled along the oblique plane, tracing a feature of the soft tissue in the reformatted image, and generating, with the processor, a 3-D marking representing the traced feature on a 3-D model.
  • the image data is raw image data.
  • FIG. 1 shows an illustrative flowchart for identifying anatomical landmarks in image data
  • FIG. 2 shows an illustrative flowchart for identifying a base oblique plane
  • FIG. 3 shows an illustrative graphical user interface for identifying the base oblique plane of FIG. 2;
  • FIG. 4 A shows an illustrative graphical user interface for identifying a first point and a second point to identify a line parallel to a longitudinal axis of soft tissue
  • FIG. 4B shows an illustrative preview window that displays a candidate base oblique plane
  • FIG. 4C shows an illustrative graphical user interface for identifying a curve parallel to a longitudinal axis of soft tissue
  • FIG. 5 shows an illustrative flowchart for generating reformatted images sampled along an adjusted set of orthogonal planes
  • FIGS. 6A and 6B show screenshots of an illustrative graphical user interface displaying reformatted images sampled along the adjusted set of orthogonal planes;
  • FIG. 7 shows an illustrative graphical user interface for performing a tracing process using multiple views provided by the adjusted set of orthogonal planes.
  • FIG. 8 shows illustrative 3-D markings generated from tracing data.
  • anatomical landmarks such as ligament attachment points to the tibia
  • Various anatomical landmarks may be located in the image data using a set of reformatted images that provide a customized view of an anatomical landmark. It will be understood that the set of reformatted images may include two or more reformatted images.
  • a ligament attachment point for which this approach is well-suited is the attachment point between the anterior cruciate ligament (ACL) and the tibia.
  • any other suitable soft tissue, soft tissue attachment points, or other features of a patient's joint or bone may be identified including, for example, the medial sulcus of the tibia plateau or the lateral sulcus of the tibia plateau as well as attachment points between the medial collateral ligament (MCL) and the tibia, between the posterior cruciate ligament (PCL) and the tibia, between the patellar tendon and the tibia, or any other suitable feature or landmark.
  • MCL medial collateral ligament
  • PCL posterior cruciate ligament
  • FIG. 1 shows a flowchart 100 for identifying anatomical landmarks in image data according to certain embodiments.
  • image data representative of a patient's joint is received.
  • the representation of the patient's joint includes soft tissue, such as a ligament, and bone.
  • the soft tissue e.g. ligament
  • the image data received in step 102 can be raw image data.
  • the image data may be received from an imaging machine without significant post-processing.
  • the image data can be obtained using any suitable medical imaging technique including, for example, MRI, CT scan, ultrasound, X-ray, or any other suitable technique or combination thereof.
  • the image data may be topographic image data comprising multiple image slices taken in orthogonal planes. Orthogonal planes for tomographic imaging may include medial-lateral planes (i.e. sagittal planes), anterior- posterior planes (i.e.
  • a user can identify a base oblique plane, oblique to one or more of the orthogonal planes in step 104.
  • the customized set of views may include two or more views.
  • the base oblique plane identified in step 104 is normal to the longitudinal axis of the soft tissue and intersects the soft tissue.
  • Various methods may be used to identify a base oblique plane and are discussed in detail in relation to FIGS. 2-4.
  • an adjusted set of orthogonal planes that are orthogonal to the base oblique plane is generated in step 105.
  • the adjusted set of orthogonal planes may provide views that effectively rotate the coronal, sagittal, and axial views.
  • An adjusted set of orthogonal planes may be, for example, two or more planes, including the base plane.
  • reformatted images are sampled along the adjusted set of orthogonal planes generated in step 105. The reformatted images provide a customized set of views that facilitate identification of soft tissue attachment points on a patient's bone.
  • the reformatted images allow the landmark identification to be performed without relying on cross-sectional views of a 3-D model generated from segmentation of the image data.
  • errors associated with 3-D model generation e.g. segmentation errors and smoothing errors
  • accuracy may be increased by limiting the number of processing steps applied to the image data prior to landmark identification.
  • the generation of the reformatted images is further described in relation to FIG. 5 and an illustrative graphical user interface displaying a set of reformatted images is shown in FIGS 6B and 7.
  • the generation of the reformatted images may be performed by a graphics processing unit (GPU) to reduce processing time, thereby allowing dynamic, "on-the-fly," generation of views more conducive to landmark identification. While a central processing unit (CPU) may also be used, a GPU may be more computationally efficient due to the highly parallel architecture of GPUs. Additionally, generating the reformatted images along the adjusted set of orthogonal planes does not necessarily require the processing capabilities needed to form 3-D models of the joint.
  • GPU graphics processing unit
  • a user traces a feature of the soft tissue in one or more reformatted images in step 108.
  • the process for tracing the soft tissue feature is further described in relation to FIG. 7. Because the user traces this feature from reformatted image data without relying on cross-sectional views of 3-D models generated from segmentation data, segmentation errors and smoothing errors associated with 3-D model generation may be avoided.
  • a computer processor uses the tracings to generate a 3-D marking representing the traced feature on a 3-D model. If desired, the 3-D marking can be smoothed as further described in relation to FIG. 8.
  • FIG. 2 shows a flowchart 200 for identifying a base oblique plane according to certain embodiments.
  • the steps of flowchart 200 are related to step 104 of process 100.
  • the steps of process 200 may be performed with reference to the graphical user interface of FIGS. 3 and 4A.
  • the user identifies a first point and a second point along the longitudinal axis of the soft tissue.
  • the first and second points may be identified in a first orthogonal slice of raw image data.
  • the orthogonal plane used may be a medial-lateral plane.
  • a first point 410a and a second point 410b are selected in a medial-lateral plane.
  • first point 410a and the second point 410b is shown in FIG. 4A.
  • step 204 the processor evaluates a line defined by the first and second points. For example, as shown in FIG 4 A, line 412 is evaluated based on the location of point 410a and point 410b. The line is then used in step 206 to evaluate a base oblique plane normal to the line. For any given line there is a set of many planes normal to the line and parallel to each other along the length of the line. Therefore, in step 208 one of the planes from the set of planes is identified as the base oblique plane.
  • the selected base oblique plane is the plane from the set of planes that contains the first point. For example, in FIG. 4B, the selected base oblique plane 424 contains first point 410a. Any other suitable plane from the set of planes that adequately shows the soft tissue attachment point may be selected.
  • FIG. 4C shows a graphical user interface for identifying a curve parallel to a longitudinal axis of soft tissue according to certain embodiments. Because the longitudinal axis of the soft tissue may be curved, a user may be able to draw a curve parallel to the longitudinal axis of the soft tissue instead of a line. When a line is drawn to define the base oblique plane, planes offset from the oblique plane will also be normal to the line. However, when a curve is used to define the base oblique plane, planes offset from that base oblique plane may in some cases not be normal to the curve at the point that the offset plane intersects the curve. Therefore, to better align the adjusted set of orthogonal planes to the curve, multiple base oblique planes may be generated, each base oblique plane being normal to the curve at the point that each plane intersects the curve.
  • FIG. 3 shows a graphical user interface 300 for identifying the base oblique plane of FIG. 2 according to certain embodiments.
  • Graphical user interface 300 contains an image window 302, and buttons 306, 308, and 310.
  • the image window displays an orthogonal slice of raw image data.
  • a medial-lateral slice of a patient's knee joint is shown.
  • Located near the center of the knee joint is a soft tissue 304.
  • the soft tissue shown in the window is the anterior cruciate ligament (ACL).
  • the ACL is a narrow soft tissue structure which attaches to the tibia and the femur.
  • the ACL may not be adequately contained in the orthogonal planes.
  • an adjusted set of orthogonal views based on a base oblique plane that intersects the ACL and is normal to the longitudinal axis of the ACL is desired.
  • next case button 306 is provided so that a case can be skipped if the image quality is found to be inadequate for analysis.
  • An image can be inadequate for a number of reasons such as low image resolution, poor visibility of the structures to be identified, or absence of segmentation data.
  • the graphical user interface may provide the user with the ability to annotate an inadequate image to note the reason for skipping it. If, however, the image quality is found to be acceptable, the user can press "set rotation" button 308 to initiate the process of identifying an appropriate oblique plane as further described in FIG. 4A.
  • Save button 310 is provided so that users can save their progress at any point during the landmark identification process.
  • FIG. 4A shows a graphical user interface for identifying a first point and a second point to identify a line parallel to a longitudinal axis of soft tissue according to certain embodiments.
  • line 412 and points 410a and 410b appear on the graphical user interface.
  • the user can drag and drop the points 410a and 410b to align line 412 with the longitudinal axis of the soft tissue 304. It may be desirable for line 412 to be placed on the anterior side of the ACL.
  • the first point 410a may be selected to coincide with the attachment point of the ACL, and the second point 410b may be selected to coincide with the midpoint of the ACL although any suitable point along the ACL may be selected.
  • the first point 410a may be used later to select a plane from a set of planes that are normal to line 412.
  • the second point 410b may be used primarily for establishing the slope of line 412.
  • the exact position of point 410a may be less important than the position of point 410b or the slope of line 412. That is, the direction of the line may be more important than its length.
  • FIG. 4B shows a preview window 420 that displays a candidate base oblique plane 424 according to certain embodiments.
  • the candidate base oblique plane 424 provides a preview of the base oblique plane that would be generated from the user's current selection of points 410a and 410b.
  • the candidate base oblique plane 424 may be shown in the preview window 420 while the user positions the points 410a and 410b.
  • Such a preview window 420 may aid the user in the process of identifying the first and second points.
  • the candidate base oblique plane is the oblique plane selected from the set of oblique planes normal to line 412 that contains first point 410a.
  • the preview window 420 also displays a reformatted image on the candidate base oblique plane 424.
  • the user may accept or reject the candidate base oblique plane based on whether the reformatted image provides an adequate view of a ligament attachment point.
  • the reformatted image on the candidate base oblique plane 424 may be generated using a graphics processing unit (GPU) to reduce processing time.
  • the preview window 410 displays a 3-D view of both the candidate base oblique plane 424 and an orthogonal slice 422 together so that the user can verify the relative orientation of the candidate base oblique plane.
  • a 3-D model of the patient's bone 428 may also be displayed with the image data to orient the user.
  • the user may reposition points 410a and 410b to generate a new candidate base oblique plane 424. If the candidate base oblique plane 424 is found to be acceptable, the user may indicate the acceptance of the plane by clicking a button, and a separate window displaying the reformatted image may appear in the graphical user interface. Example reformatted image windows are shown in FIGS. 6B and 7.
  • FIG. 5 shows a flowchart 500 for generating reformatted images sampled along the adjusted set of orthogonal planes according to certain embodiments.
  • the steps of flowchart 500 are related to step 106 of process 100.
  • step 502 the location and orientation of the base oblique plane are received. This location and orientation information defines the base oblique plane.
  • An adjusted set of orthogonal planes may be generated based on the base oblique plane. Each plane in the adjusted set of orthogonal planes is either parallel or orthogonal to the base oblique plane.
  • the oblique planes in the set may be equally spaced and may span the entire volume of the raw image data.
  • the raw image data may be effectively rotated, resulting in a new stack of image slices oriented with the base oblique plane.
  • the oblique planes generated along a curve there may be multiple base oblique planes.
  • planes orthogonal to each base oblique plane may be generated to provide a set of adjusted orthogonal views that follows the curve. It is also envisioned that the steps here can be performed without identifying a set of orthogonal planes.
  • a reformatted image may be sampled along the base oblique plane only. [0040]
  • Four corners of the planes in the adjusted set of planes are computed in step 504.
  • the corners of the oblique planes are determined by the intersection of the infinite plane and the bounding volume of the raw image data.
  • candidate base oblique plane 424 has corner 426 at the intersection of the oblique plane and the bounding volume of the raw image data.
  • the corners are used to generate texture coordinates for the planes in the adjusted set of orthogonal planes.
  • the texture coordinates are given pixel intensity values for displaying the reformatted image on the oblique plane.
  • Pixel intensity values may represent, for example, radiodensity measurements, attenuation coefficients, or tissue relaxation time such as Tl, T2, and proton density.
  • the pixel intensity values for texture coordinates that do not coincide with points in the raw image data are inferred using interpolation.
  • interpolation methods include nearest-neighbor interpolation, bilinear interpolation, trilinear interpolation, cubic interpolation, and any other suitable technique for interpolation or any combination thereof.
  • One of the differences among the types of interpolation is the assumption of the continuity of the spatial distribution of the pixel values.
  • Step 508 may be performed using the graphical processing unit (GPU) or the central processing unit (CPU).
  • GPU graphical processing unit
  • CPU central processing unit
  • the output of the interpolation step may be a matrix of pixel intensity values, with each value in the matrix corresponding to a generated texture coordinate.
  • the GPU or the CPU sends the output to the display
  • FIGS. 6 A and 6B show screenshots of a graphical user interface displaying reformatted images sampled along the adjusted set of orthogonal planes.
  • Window 606 shows the reformatted image on the base oblique plane.
  • the image shown in window 606 is the same reformed image that was previewed on candidate base oblique plane 424.
  • 602a and 604 show the corresponding views orthogonal to the base oblique plane. Therefore, windows 602a, 604 and 606 together provide a set of adjusted orthogonal views.
  • the raw image data is effectively rotated to provide a new stack of image slices oriented with respect to the base oblique plane.
  • the user may be able to advance through the reformatted images in the same way that the user would normally advance through orthogonal slices of the raw image data.
  • the user may be provided with a tool bar having toggle buttons that allow different views to be turned on or off.
  • FIG. 7 shows a graphical user interface for performing a tracing process using multiple views provided by the adjusted set of orthogonal planes according to certain embodiments.
  • Reformatted image window 708 shows an enlarged reformatted image of a patient's knee joint and 2-D marking 710c.
  • 2-D markings may be created by tracing the ACL attachment point in image windows 702, 704, or 708.
  • 2-D marking 710c may be created by tracing the attachment point in image window 708.
  • the 2-D markings may be traced using a cursor controlled by computer mouse. For example, a user may click and hold a mouse button while tracing the ligament attachment point with the cursor to create the 2-D marking. Additionally, a user may advance through reformatted images during tracing.
  • the user may advance through reformatted images normal to the longitudinal axis of the soft tissue beginning just before the ligament attachment point and may stop after the end of the ligament. Advancing through images in this way may facilitate identification of the end of a ligament, and therefore the ligament attachment point, because there may be a sharp change in the grayscale of the cross-sections as the end of the ligament is passed.
  • a user may start just before the starting point of a ligament and advance through reformatted images until the end of the ligament is reached or just passed. Advancing in this way may allow user to identify attachment points on both ends of the ligament. For example, advancing from the start of the ACL to its end may facilitate identification of both the attachment point to the femur and the attachment point to the tibia.
  • the user may outline the cross- sections of a ligament in any of the adjusted orthogonal planes. Additionally, a user may also be provided with the ability to erase all or a part of a 2-D marking made in any view. It may be desirable to advance through 8 to 10 reformatted images during tracing to adequately mark an attachment point. As the user traces the attachment point in the cross-sections, a matrix representing the tracing data may be generated. For example, the data representing the points selected during tracing may be stored in a binary map. This tracing data may be used to generate a 3-D marking as discussed further in relation to FIG. 8.
  • the graphical user interface of FIG. 7 also contains surface model window 706 which displays a 3-D surface model 712 representing a patient's tibia.
  • the 3-D surface model may help orient the user to the patient's anatomy during tracing.
  • the 3-D surface model may be created prior to tracing based on segmentation of the raw image data.
  • the same 2-D marking traced in the adjusted orthogonal planes is displayed as marking 710b on the 3-D surface model 712.
  • the graphical user interface may facilitate creation of 2-D markings representing the ligament attachment point.
  • FIG. 8 shows 3-D markings generated from tracing data according to certain embodiments.
  • the 3-D markings may be generated using the 2-D markings created in the tracing process described in FIG. 7.
  • 3-D marking 802 which represents the ACL attachment point to the tibia
  • 3-D marking 802 may be generated from the tracing data using the marching cubes algorithm.
  • 3-D marking 804, which represents the patellar tendon attachment point to the anterior side of the tibia can be created by a similar process as that used to create 3-D marking 802.
  • point markings 806 and 808 may each be created from the identification of a single landmark point in one or more of the adjusted orthogonal planes.
  • the point marking 806 is selected to correspond to the medial sulcus of the tibia plateau and point marking 808 is selected to correspond to the lateral sulcus of the tibia plateau.
  • the point markings are selected to correspond to the lowest points on the tibia condyles.
  • smoothing may be performed on the surface model or on the tracing data before generation of the surface model. Such smoothing may be performed for processing efficiency because smoothed 3-D markings may not require subsequent smoothing when imported into a computer aided design (CAD) environment.
  • CAD computer aided design
  • a method for locating anatomical landmarks in image data comprising:
  • identifying, with the processor, a set of planes comprising the oblique plane and at least one other plane parallel or orthogonal to the oblique plane;
  • tracing the feature comprises tracing the feature in two or more reformatted images parallel to the oblique plane.
  • tracing the feature comprises: displaying a first reformatted image on a screen;
  • tracing the feature comprises outlining a cross-section of the feature in a reformatted image from the set of reformatted images.
  • tracing the feature comprises storing data representing the tracing in a matrix.
  • identifying the oblique plane comprises:
  • A10 The method of any of embodiments A1-A7, wherein the first point is located at about the start of the soft tissue and the second point is located at about the end of the soft tissue.
  • Al 1 The method of any of embodiments A1-A10, wherein the longitudinal axis of the soft tissue is curved.
  • A14 The method of any of embodiments A1-A13, wherein generating the 3-D marking comprises generating a 3-D surface model from the matrix.
  • A15 The method of any of embodiments A1-A14, further comprising smoothing the 3-D marking.
  • A16 The method of any of embodiments A1-A15, wherein the 3-D model is a model of a patient's bone.
  • A17 The method of any of embodiments A1-A16, wherein the 3-D model is a model of a patient's tibia.
  • A20 The method of any of embodiments A1-A19, wherein the soft tissue is an ACL, a MCL, a PCL, or a patellar tendon of a patient.
  • a system for locating anatomical landmarks in image data comprising:
  • Bl 1. The system of any of embodiments B1-B10, wherein the longitudinal axis of the soft tissue is curved.
  • B13 The system of embodiment Bl 1 or B12, wherein a plurality of reformatted images is generated, each reformatted image being normal to the curve at the point where each reformatted image intersects the curve.
  • B14 The system of any of embodiments B1-B13, wherein the means for generating the 3-D marking comprises means for generating a 3-D surface model from the matrix.
  • a method for locating anatomical landmarks in image data comprising:

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Radiology & Medical Imaging (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Optics & Photonics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Computer Graphics (AREA)
  • Rheumatology (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Human Computer Interaction (AREA)
  • Physiology (AREA)
  • Pulmonology (AREA)
  • Rehabilitation Therapy (AREA)
  • Quality & Reliability (AREA)
  • Architecture (AREA)
  • Primary Health Care (AREA)
  • Data Mining & Analysis (AREA)
  • Epidemiology (AREA)

Abstract

Systems, devices, and methods are described for locating and identifying anatomical landmarks, such as ligament attachment points, in image data. These systems, devices, and methods may provide an oblique plane that contains an anatomical landmark such as a ligament attachment point to the tibia. For example, the position at which the anterior cruciate ligament (ACL), medial collateral ligament (MCL) posterior cruciate ligament (PCL), or patellar tendon attaches to the tibia may be identified. The systems, devices, and methods allow for tracing of an anatomical landmark to generate a 3-D marking on a 3-D surface model of a patient's bone. The attachment points may be useful landmarks for patient-matched instrumentation.

Description

METHODS AND SYSTEMS FOR IDENTIFYING ANATOMICAL LANDMARKS IN
IMAGE DATA
Cross Reference to Related Applications
[0001] This application claims the benefit of U.S. Provisional Application No. 61/950,756, filed on March 10, 2014, which is hereby incorporated herein by reference in its entirety.
Background
[0002] Accurate location of anatomical landmarks, such as ligament attachment points, in patient image data is critical to the overall outcome of patient-matched instrumentation techniques. Ligament attachment points can be used to facilitate consistent placement of implants on a patient's bone and may offer surgeons a reliable frame of reference for properly orienting implants. However, the location of ligament attachment points may not be known by a surgeon before a surgical operation and may need to be assessed intraoperatively.
Noninvasive imaging techniques, such as MRI scans and CT scans, allow image data representing a patient's joint to be collected in advance of a surgical procedure. However, it can be difficult to use the raw image data collected from preoperative imaging to provide a surgeon with accurate preoperative information about the location of ligament attachment points.
[0003] The location of anatomical landmarks in such image data is evaluated using 2-D slices within the orthogonal planes of a 3-D volume (e.g., coronal, sagittal, and axial planes). This approach provides suboptimal visualization of ligament attachment points that are not adequately depicted in any of these planes. As a result, identifying such attachment points (e.g. attachments points on epicondyles) is a time-consuming process and suffers from high inter-operator variability. Therefore, to facilitate accurate preoperative location of ligament attachment points, there is a need for improved visualization of ligament attachment points. Summary
[0004] Disclosed herein are systems, devices, and methods for locating anatomical landmarks, such as ligament attachment points, in image data. In certain implementations, the systems, devices, and methods include receiving image data, identifying a base oblique plane that is oblique to the orthogonal planes in the image data, generating an adjusted set of orthogonal planes based on the base oblique plane, and generating reformatted images sampled along the adjusted set of orthogonal planes. These systems, devices, and methods may provide an oblique plane that more adequately contains an anatomical landmark such as a ligament attachment point to the tibia. In particular, the method may be used to identify where the anterior cruciate ligament (ACL), medial collateral ligament (MCL) posterior cruciate ligament (PCL), or patellar tendon attaches to the tibia. The systems, devices, and methods may allow tracing of an anatomical landmark in the reformatted image to generate a 3-D marking on a 3-D surface model of a patient's bone. The attachment points may be useful landmarks for patient-matched instrumentation. The systems, devices, and methods may also facilitate surgical planning. For example, the ligament locations may be used for preoperative planning and for designing patient-matched cutting blocks.
[0005] According to one aspect, a method for locating anatomical landmarks in image data includes receiving, at a processor, image data representative of a patient's joint, identifying an oblique plane intersecting a soft tissue represented in the image data and normal to a longitudinal axis of the soft tissue, identifying, with the processor, a set of planes parallel or orthogonal to the oblique plane, generating, with the processor, a set of reformatted images, from the image data, sampled along the set of planes, tracing a feature of the soft tissue in the set of reformatted images, and generating, with the processor, a 3-D marking representing the traced feature on a 3-D model. In some implementations, the image data is raw image data. Tracing the feature may include tracing the feature in two or more reformatted images parallel to the oblique plane. Tracing the feature may include advancing through parallel reformatted images selected from the set of reformatted images. Tracing the feature may include displaying a first reformatted image on a screen and displaying a second reformatted image on the screen, wherein the second reformatted image is parallel to and offset from the first reformatted image. Tracing the feature may include outlining a cross-section of the feature in a reformatted image from the set of reformatted images. Tracing the feature may include storing data representing the tracing in a matrix.
[0006] In certain implementations, generating the set of reformatted images includes rotating the image data. Identifying the oblique plane may include identifying a first point and a second point along the longitudinal axis of the soft tissue, evaluating, with the processor, a line defined by the first and second points, and evaluating, with the processor, a plane normal to the line. The first point may be located about where the soft tissue attaches to a patient's bone. The second point may be located at about the center of the soft tissue. The first point may be located at about the start of the soft tissue and the second point may be located at about the end of the soft tissue. The longitudinal axis of the soft tissue may be curved.
[0007] In certain implementations, identifying the oblique plane includes drawing a curve substantially parallel to the longitudinal axis of the tissue, and evaluating, with the processor, an oblique plane normal to the curve at the point where the oblique plane intersects the curve. In certain implementations, a plurality of reformatted images is generated, each reformatted image being normal to the curve at the point where each reformatted image intersects the curve. Generating the 3-D marking may include generating a 3-D surface model from the matrix. In certain implementations, the method further includes smoothing the 3-D marking. The 3-D model may be a model of a patient's bone. The 3-D model may be a model of a patient's tibia. The 3-D model may be displayed with the image data to orient a user.
[0008] In certain implementations, generating the set of reformatted images includes identifying, with the processor, four corners of the oblique plane, generating, with the processor, texture coordinates using the corner locations, and interpolating, with the processor, pixel intensities for the texture coordinates using the image data. The soft tissue may be an ACL, a MCL, a PCL, or a patellar tendon of a patient.
[0009] According to one aspect, a system for locating anatomical landmarks in image data includes means for receiving image data representative of a patient's joint, means for identifying an oblique plane intersecting a soft tissue represented in the image data and normal to a longitudinal axis of the soft tissue, means for identifying a set of planes parallel or orthogonal to the oblique plane, means for generating a set of reformatted images, from the image data, sampled along the set of planes, means for tracing a feature of the soft tissue in the set of reformatted images, and means for generating a 3-D marking representing the traced feature on a 3-D model. In some implementations, the image data is raw image data. The means for tracing the feature may include means for tracing the feature in two or more reformatted images parallel to the oblique plane. The means for tracing the feature may include means for advancing through parallel reformatted images selected from the set of reformatted images.
[0010] In certain implementations, the means for tracing the feature includes means for outlining a cross-section of the feature in a reformatted image from the set of reformatted images. The means for tracing the feature may include means for storing data representing the tracing in a matrix. The means for generating the set of reformatted images may include means for rotating the image data. The means for identifying the oblique plane may include means for identifying a first point and a second point along the longitudinal axis of the soft tissue, means for evaluating a line defined by the first and second points, and means for evaluating a plane normal to the line. The first point may be located about where the soft tissue attaches to a patient's bone. The second point may be located at about the center of the soft tissue. The first point may be located at about the start of the soft tissue and the second point may be located at about the end of the soft tissue. The longitudinal axis of the soft tissue may be curved.
[0011] In certain implementations, the means for identifying the oblique plane includes means for drawing a curve substantially parallel to the longitudinal axis of the tissue, and means for evaluating an oblique plane normal to the curve at the point where the oblique plane intersects the curve. In certain implementations a plurality of reformatted images is generated, each reformatted image being normal to the curve at the point where each reformatted image intersects the curve. The means for generating the 3-D marking may include means for generating a 3-D surface model from the matrix. In certain
implementations, the system further includes means for smoothing the 3-D marking.
[0012] In certain implementations, the 3-D model is a model of a patient's bone. The 3-D model may be a model of a patient's tibia. The 3-D model may be displayed with the image data to orient a user. In certain implementations, the means for generating the set of reformatted images includes means for computing four corners of the oblique plane, means for generating texture coordinates using the corner locations, and means for interpolating pixel intensities for the texture coordinates using the image data. In certain implementations the soft tissue is an ACL, a MCL, a PCL, or a patellar tendon of a patient.
[0013] According to one aspect, a method for locating anatomical landmarks in image data, includes receiving, at a processor, image data representative of a patient's joint, identifying an oblique plane intersecting a soft tissue represented in the image data and normal to a longitudinal axis of the soft tissue, generating, with the processor, a reformatted image, from the image data, sampled along the oblique plane, tracing a feature of the soft tissue in the reformatted image, and generating, with the processor, a 3-D marking representing the traced feature on a 3-D model. In some implementations, the image data is raw image data.
[0014] Variations and modifications will occur to those of skill in the art after reviewing this disclosure. The disclosed features may be implemented, in any combination and subcombination (including multiple dependent combinations and subcombinations), with one or more other features described herein. The various features described or illustrated above, including any components thereof, may be combined or integrated in other systems.
Moreover, certain features may be omitted or not implemented.
Brief Description of the Drawings
[0015] The foregoing and other objects and advantages will be apparent upon consideration of the following detailed description, taken in conjunction with the accompanying drawings, in which like reference characters refer to like parts throughout, and in which:
[0016] FIG. 1 shows an illustrative flowchart for identifying anatomical landmarks in image data;
[0017] FIG. 2 shows an illustrative flowchart for identifying a base oblique plane;
[0018] FIG. 3 shows an illustrative graphical user interface for identifying the base oblique plane of FIG. 2;
[0019] FIG. 4 A shows an illustrative graphical user interface for identifying a first point and a second point to identify a line parallel to a longitudinal axis of soft tissue;
[0020] FIG. 4B shows an illustrative preview window that displays a candidate base oblique plane;
[0021] FIG. 4C shows an illustrative graphical user interface for identifying a curve parallel to a longitudinal axis of soft tissue;
[0022] FIG. 5 shows an illustrative flowchart for generating reformatted images sampled along an adjusted set of orthogonal planes;
[0023] FIGS. 6A and 6B show screenshots of an illustrative graphical user interface displaying reformatted images sampled along the adjusted set of orthogonal planes;
[0024] FIG. 7 shows an illustrative graphical user interface for performing a tracing process using multiple views provided by the adjusted set of orthogonal planes; and
[0025] FIG. 8 shows illustrative 3-D markings generated from tracing data. Detailed Description
[0026] To provide an overall understanding of the systems, devices, and methods described herein, certain illustrative embodiments will be described. Although the embodiments and features described herein are specifically described for use in connection with identifying anatomical landmarks in image data for a patient's knee joint, it will be understood that the approaches disclosed are applicable to other anatomical joints as well. Moreover, the approaches outlined below may be applied to preparing patient-matched medical devices and implants to be used in surgical procedures including, but not limited to, knee, acetabular, spine arthroplasty, cranio-maxillofacial surgical procedures, shoulder arthroplasty, as well as foot, ankle, hand, and other extremity procedures.
[0027] The following disclosure provides systems, devices, and methods for identifying anatomical landmarks in image data. Various anatomical landmarks, such as ligament attachment points to the tibia, may be located in the image data using a set of reformatted images that provide a customized view of an anatomical landmark. It will be understood that the set of reformatted images may include two or more reformatted images. A ligament attachment point for which this approach is well-suited is the attachment point between the anterior cruciate ligament (ACL) and the tibia. Any other suitable soft tissue, soft tissue attachment points, or other features of a patient's joint or bone may be identified including, for example, the medial sulcus of the tibia plateau or the lateral sulcus of the tibia plateau as well as attachment points between the medial collateral ligament (MCL) and the tibia, between the posterior cruciate ligament (PCL) and the tibia, between the patellar tendon and the tibia, or any other suitable feature or landmark.
[0028] FIG. 1 shows a flowchart 100 for identifying anatomical landmarks in image data according to certain embodiments. In step 102, image data representative of a patient's joint is received. The representation of the patient's joint includes soft tissue, such as a ligament, and bone. The soft tissue (e.g. ligament) attaches to the bone at an attachment point.
Accurately identifying such an attachment point on a patient's bone provides surgeons with accurate preoperative information about the patient's anatomy and accurate models of a patient's joint for patient-matched instrumentation. The image data received in step 102 can be raw image data. For example, the image data may be received from an imaging machine without significant post-processing. The image data can be obtained using any suitable medical imaging technique including, for example, MRI, CT scan, ultrasound, X-ray, or any other suitable technique or combination thereof. The image data may be topographic image data comprising multiple image slices taken in orthogonal planes. Orthogonal planes for tomographic imaging may include medial-lateral planes (i.e. sagittal planes), anterior- posterior planes (i.e. coronal planes), and transverse planes (i.e. axial planes). These views, however, may not provide the best available view for identifying a particular ligament attachment point. To create a customized set of views that facilitate landmark identification, a user can identify a base oblique plane, oblique to one or more of the orthogonal planes in step 104. The customized set of views may include two or more views. The base oblique plane identified in step 104 is normal to the longitudinal axis of the soft tissue and intersects the soft tissue. Various methods may be used to identify a base oblique plane and are discussed in detail in relation to FIGS. 2-4.
[0029] After the base oblique plane is identified, an adjusted set of orthogonal planes that are orthogonal to the base oblique plane is generated in step 105. The adjusted set of orthogonal planes may provide views that effectively rotate the coronal, sagittal, and axial views. An adjusted set of orthogonal planes may be, for example, two or more planes, including the base plane. In step 106, reformatted images are sampled along the adjusted set of orthogonal planes generated in step 105. The reformatted images provide a customized set of views that facilitate identification of soft tissue attachment points on a patient's bone. Additionally, the reformatted images allow the landmark identification to be performed without relying on cross-sectional views of a 3-D model generated from segmentation of the image data. As a result, errors associated with 3-D model generation (e.g. segmentation errors and smoothing errors) are not introduced into the landmark identification process. Thus, accuracy may be increased by limiting the number of processing steps applied to the image data prior to landmark identification.
[0030] The generation of the reformatted images is further described in relation to FIG. 5 and an illustrative graphical user interface displaying a set of reformatted images is shown in FIGS 6B and 7. The generation of the reformatted images may be performed by a graphics processing unit (GPU) to reduce processing time, thereby allowing dynamic, "on-the-fly," generation of views more conducive to landmark identification. While a central processing unit (CPU) may also be used, a GPU may be more computationally efficient due to the highly parallel architecture of GPUs. Additionally, generating the reformatted images along the adjusted set of orthogonal planes does not necessarily require the processing capabilities needed to form 3-D models of the joint. Once the reformatted images have been generated, a user traces a feature of the soft tissue in one or more reformatted images in step 108. The process for tracing the soft tissue feature is further described in relation to FIG. 7. Because the user traces this feature from reformatted image data without relying on cross-sectional views of 3-D models generated from segmentation data, segmentation errors and smoothing errors associated with 3-D model generation may be avoided. In step 110, a computer processor uses the tracings to generate a 3-D marking representing the traced feature on a 3-D model. If desired, the 3-D marking can be smoothed as further described in relation to FIG. 8.
[0031] FIG. 2 shows a flowchart 200 for identifying a base oblique plane according to certain embodiments. For example, the steps of flowchart 200 are related to step 104 of process 100. The steps of process 200 may be performed with reference to the graphical user interface of FIGS. 3 and 4A. In step 202, the user identifies a first point and a second point along the longitudinal axis of the soft tissue. The first and second points may be identified in a first orthogonal slice of raw image data. For identifying ligament attachment points in the knee joint, the orthogonal plane used may be a medial-lateral plane. For example, as shown in FIG. 4A, a first point 410a and a second point 410b are selected in a medial-lateral plane. It may be desirable to locate the first point about where the soft tissue attaches to a patient's bone because the first point may be used later to position a base oblique plane that provides a top view of the attachment point. The second point may be located at about the center of the soft tissue. Such a placement of the first point 410a and the second point 410b is shown in FIG. 4A.
[0032] In step 204, the processor evaluates a line defined by the first and second points. For example, as shown in FIG 4 A, line 412 is evaluated based on the location of point 410a and point 410b. The line is then used in step 206 to evaluate a base oblique plane normal to the line. For any given line there is a set of many planes normal to the line and parallel to each other along the length of the line. Therefore, in step 208 one of the planes from the set of planes is identified as the base oblique plane. In some embodiments, the selected base oblique plane is the plane from the set of planes that contains the first point. For example, in FIG. 4B, the selected base oblique plane 424 contains first point 410a. Any other suitable plane from the set of planes that adequately shows the soft tissue attachment point may be selected.
[0033] FIG. 4C shows a graphical user interface for identifying a curve parallel to a longitudinal axis of soft tissue according to certain embodiments. Because the longitudinal axis of the soft tissue may be curved, a user may be able to draw a curve parallel to the longitudinal axis of the soft tissue instead of a line. When a line is drawn to define the base oblique plane, planes offset from the oblique plane will also be normal to the line. However, when a curve is used to define the base oblique plane, planes offset from that base oblique plane may in some cases not be normal to the curve at the point that the offset plane intersects the curve. Therefore, to better align the adjusted set of orthogonal planes to the curve, multiple base oblique planes may be generated, each base oblique plane being normal to the curve at the point that each plane intersects the curve.
[0034] FIG. 3 shows a graphical user interface 300 for identifying the base oblique plane of FIG. 2 according to certain embodiments. Graphical user interface 300 contains an image window 302, and buttons 306, 308, and 310. The image window displays an orthogonal slice of raw image data. In this example, a medial-lateral slice of a patient's knee joint is shown. Located near the center of the knee joint is a soft tissue 304. The soft tissue shown in the window is the anterior cruciate ligament (ACL). The ACL is a narrow soft tissue structure which attaches to the tibia and the femur. The ACL may not be adequately contained in the orthogonal planes. For better visualization of the ACL and its attachment points, an adjusted set of orthogonal views based on a base oblique plane that intersects the ACL and is normal to the longitudinal axis of the ACL is desired.
[0035] Before selecting the first and second points to evaluate the base oblique plane, the image quality may be assessed. Next case button 306 is provided so that a case can be skipped if the image quality is found to be inadequate for analysis. An image can be inadequate for a number of reasons such as low image resolution, poor visibility of the structures to be identified, or absence of segmentation data. The graphical user interface may provide the user with the ability to annotate an inadequate image to note the reason for skipping it. If, however, the image quality is found to be acceptable, the user can press "set rotation" button 308 to initiate the process of identifying an appropriate oblique plane as further described in FIG. 4A. Save button 310 is provided so that users can save their progress at any point during the landmark identification process.
[0036] FIG. 4A shows a graphical user interface for identifying a first point and a second point to identify a line parallel to a longitudinal axis of soft tissue according to certain embodiments. After the user has pressed button 308, line 412 and points 410a and 410b appear on the graphical user interface. The user can drag and drop the points 410a and 410b to align line 412 with the longitudinal axis of the soft tissue 304. It may be desirable for line 412 to be placed on the anterior side of the ACL. The first point 410a may be selected to coincide with the attachment point of the ACL, and the second point 410b may be selected to coincide with the midpoint of the ACL although any suitable point along the ACL may be selected. The first point 410a may be used later to select a plane from a set of planes that are normal to line 412. The second point 410b may be used primarily for establishing the slope of line 412. In such an implementation, the exact position of point 410a may be less important than the position of point 410b or the slope of line 412. That is, the direction of the line may be more important than its length.
[0037] FIG. 4B shows a preview window 420 that displays a candidate base oblique plane 424 according to certain embodiments. The candidate base oblique plane 424 provides a preview of the base oblique plane that would be generated from the user's current selection of points 410a and 410b. The candidate base oblique plane 424 may be shown in the preview window 420 while the user positions the points 410a and 410b. Such a preview window 420 may aid the user in the process of identifying the first and second points. The candidate base oblique plane is the oblique plane selected from the set of oblique planes normal to line 412 that contains first point 410a.
[0038] The preview window 420 also displays a reformatted image on the candidate base oblique plane 424. The user may accept or reject the candidate base oblique plane based on whether the reformatted image provides an adequate view of a ligament attachment point. The reformatted image on the candidate base oblique plane 424 may be generated using a graphics processing unit (GPU) to reduce processing time. Additionally, the preview window 410 displays a 3-D view of both the candidate base oblique plane 424 and an orthogonal slice 422 together so that the user can verify the relative orientation of the candidate base oblique plane. A 3-D model of the patient's bone 428 may also be displayed with the image data to orient the user. If the candidate base oblique plane 424 is found to be inadequate, the user may reposition points 410a and 410b to generate a new candidate base oblique plane 424. If the candidate base oblique plane 424 is found to be acceptable, the user may indicate the acceptance of the plane by clicking a button, and a separate window displaying the reformatted image may appear in the graphical user interface. Example reformatted image windows are shown in FIGS. 6B and 7.
[0039] FIG. 5 shows a flowchart 500 for generating reformatted images sampled along the adjusted set of orthogonal planes according to certain embodiments. For example, the steps of flowchart 500 are related to step 106 of process 100. In step 502, the location and orientation of the base oblique plane are received. This location and orientation information defines the base oblique plane. An adjusted set of orthogonal planes may be generated based on the base oblique plane. Each plane in the adjusted set of orthogonal planes is either parallel or orthogonal to the base oblique plane. The oblique planes in the set may be equally spaced and may span the entire volume of the raw image data. Therefore, when the reformatted images are sampled along the oblique planes, the raw image data may be effectively rotated, resulting in a new stack of image slices oriented with the base oblique plane. In the case of the oblique planes generated along a curve, there may be multiple base oblique planes. In such a case, planes orthogonal to each base oblique plane may be generated to provide a set of adjusted orthogonal views that follows the curve. It is also envisioned that the steps here can be performed without identifying a set of orthogonal planes. In such embodiments, a reformatted image may be sampled along the base oblique plane only. [0040] Four corners of the planes in the adjusted set of planes are computed in step 504. The corners of the oblique planes are determined by the intersection of the infinite plane and the bounding volume of the raw image data. For example, as shown in FIG. 4A, candidate base oblique plane 424 has corner 426 at the intersection of the oblique plane and the bounding volume of the raw image data. In step 506, the corners are used to generate texture coordinates for the planes in the adjusted set of orthogonal planes. The texture coordinates are given pixel intensity values for displaying the reformatted image on the oblique plane. Pixel intensity values may represent, for example, radiodensity measurements, attenuation coefficients, or tissue relaxation time such as Tl, T2, and proton density. In the texture mapping process of step 508, the pixel intensity values for texture coordinates that do not coincide with points in the raw image data are inferred using interpolation. Possible interpolation methods include nearest-neighbor interpolation, bilinear interpolation, trilinear interpolation, cubic interpolation, and any other suitable technique for interpolation or any combination thereof. One of the differences among the types of interpolation is the assumption of the continuity of the spatial distribution of the pixel values. Linear
interpolation, for example, assumes the intensity varies linearly so the unknown value can be determined from its four neighbors as a weighted average. Cubic interpolation, on the other hand, has a higher requirement on the continuity: not only the intensity value but also the first and second order derivatives must be continuous. Cubic interpolation thus requires more neighboring points beyond the immediate four neighbors, but produces results that are smoother. Step 508 may be performed using the graphical processing unit (GPU) or the central processing unit (CPU). Using the GPU instead of the CPU may have performance benefits due to the highly parallel architecture of GPUs. The output of the interpolation step may be a matrix of pixel intensity values, with each value in the matrix corresponding to a generated texture coordinate. The GPU or the CPU sends the output to the display
environment in step 510.
[0041] FIGS. 6 A and 6B show screenshots of a graphical user interface displaying reformatted images sampled along the adjusted set of orthogonal planes. Window 606 shows the reformatted image on the base oblique plane. The image shown in window 606 is the same reformed image that was previewed on candidate base oblique plane 424. Windows
602a and 604 show the corresponding views orthogonal to the base oblique plane. Therefore, windows 602a, 604 and 606 together provide a set of adjusted orthogonal views. Thus, by generating reformatted images along the adjusted set of orthogonal planes, the raw image data is effectively rotated to provide a new stack of image slices oriented with respect to the base oblique plane. The user may be able to advance through the reformatted images in the same way that the user would normally advance through orthogonal slices of the raw image data. Additionally, the user may be provided with a tool bar having toggle buttons that allow different views to be turned on or off.
[0042] FIG. 7 shows a graphical user interface for performing a tracing process using multiple views provided by the adjusted set of orthogonal planes according to certain embodiments. Reformatted image window 708 shows an enlarged reformatted image of a patient's knee joint and 2-D marking 710c. 2-D markings may be created by tracing the ACL attachment point in image windows 702, 704, or 708. For example, 2-D marking 710c may be created by tracing the attachment point in image window 708. The 2-D markings may be traced using a cursor controlled by computer mouse. For example, a user may click and hold a mouse button while tracing the ligament attachment point with the cursor to create the 2-D marking. Additionally, a user may advance through reformatted images during tracing. For example, the user may advance through reformatted images normal to the longitudinal axis of the soft tissue beginning just before the ligament attachment point and may stop after the end of the ligament. Advancing through images in this way may facilitate identification of the end of a ligament, and therefore the ligament attachment point, because there may be a sharp change in the grayscale of the cross-sections as the end of the ligament is passed.
Alternatively, a user may start just before the starting point of a ligament and advance through reformatted images until the end of the ligament is reached or just passed. Advancing in this way may allow user to identify attachment points on both ends of the ligament. For example, advancing from the start of the ACL to its end may facilitate identification of both the attachment point to the femur and the attachment point to the tibia.
[0043] As the user advances through reformatted images, the user may outline the cross- sections of a ligament in any of the adjusted orthogonal planes. Additionally, a user may also be provided with the ability to erase all or a part of a 2-D marking made in any view. It may be desirable to advance through 8 to 10 reformatted images during tracing to adequately mark an attachment point. As the user traces the attachment point in the cross-sections, a matrix representing the tracing data may be generated. For example, the data representing the points selected during tracing may be stored in a binary map. This tracing data may be used to generate a 3-D marking as discussed further in relation to FIG. 8.
[0044] The graphical user interface of FIG. 7 also contains surface model window 706 which displays a 3-D surface model 712 representing a patient's tibia. The 3-D surface model may help orient the user to the patient's anatomy during tracing. The 3-D surface model may be created prior to tracing based on segmentation of the raw image data. The same 2-D marking traced in the adjusted orthogonal planes is displayed as marking 710b on the 3-D surface model 712. Thus, the graphical user interface may facilitate creation of 2-D markings representing the ligament attachment point.
[0045] FIG. 8 shows 3-D markings generated from tracing data according to certain embodiments. The 3-D markings may be generated using the 2-D markings created in the tracing process described in FIG. 7. For example, 3-D marking 802, which represents the ACL attachment point to the tibia, is generated from 2-D markings 710a and 710c. In certain implementations, 3-D marking 802 may be generated from the tracing data using the marching cubes algorithm. 3-D marking 804, which represents the patellar tendon attachment point to the anterior side of the tibia, can be created by a similar process as that used to create 3-D marking 802. Alternatively, point markings 806 and 808 may each be created from the identification of a single landmark point in one or more of the adjusted orthogonal planes. The point marking 806 is selected to correspond to the medial sulcus of the tibia plateau and point marking 808 is selected to correspond to the lateral sulcus of the tibia plateau. The point markings are selected to correspond to the lowest points on the tibia condyles. Additionally, smoothing may be performed on the surface model or on the tracing data before generation of the surface model. Such smoothing may be performed for processing efficiency because smoothed 3-D markings may not require subsequent smoothing when imported into a computer aided design (CAD) environment. Thus, the tracing process described in relation to FIG. 7 may be used to generate 3-D markings representing anatomical landmarks, such as ligament attachment points.
[0046] The foregoing is merely illustrative of the principles of the disclosure, and the systems, devices, and methods can be practiced by other than the described embodiments, which are presented for purposes of illustration and not of limitation. It is to be understood that the systems, devices, and methods disclosed herein, while shown for use in knee arthroplasty systems, may be applied to systems, devices, and methods to be used in other surgical procedures including, but not limited to, acetabular, spine arthroplasty, cranio- maxillofacial surgical procedures, shoulder arthroplasty, as well as foot, ankle, hand, and extremities procedures.
[0047] Example Embodiments
Al . A method for locating anatomical landmarks in image data, comprising:
receiving, at a processor, image data representative of a patient's joint; identifying an oblique plane intersecting a soft tissue represented in the image data and normal to a longitudinal axis of the soft tissue;
identifying, with the processor, a set of planes comprising the oblique plane and at least one other plane parallel or orthogonal to the oblique plane;
generating, with the processor, a set of reformatted images, from the image data, sampled along the set of planes;
tracing a feature of the soft tissue in the set of reformatted images; and
generating, with the processor, a 3-D marking representing the traced feature on a 3-D model.
A2. The method of embodiment Al, wherein tracing the feature comprises tracing the feature in two or more reformatted images parallel to the oblique plane.
A3. The method of embodiment Al or A2, wherein tracing the feature comprises: displaying a first reformatted image on a screen; and
displaying a second reformatted image on the screen, wherein the second reformatted image is parallel to and offset from the first reformatted image.
A4. The method of any of embodiments A1-A3, wherein tracing the feature comprises outlining a cross-section of the feature in a reformatted image from the set of reformatted images.
A5. The method of any of embodiments A1-A4, wherein tracing the feature comprises storing data representing the tracing in a matrix.
A6. The method of any of embodiments A1-A5, wherein generating the set of reformatted images comprises rotating the image data.
A7. The method of any of embodiments A1-A6, wherein identifying the oblique plane comprises:
identifying a first point and a second point along the longitudinal axis of the soft tissue;
evaluating, with the processor, a line defined by the first and second points; and evaluating, with the processor, a plane normal to the line. A8. The method of any of embodiments A1-A7, wherein the first point is located about where the soft tissue attaches to a patient's bone.
A9. The method of any of embodiments A1-A8, wherein the second point is located at about the center of the soft tissue.
A10. The method of any of embodiments A1-A7, wherein the first point is located at about the start of the soft tissue and the second point is located at about the end of the soft tissue.
Al 1. The method of any of embodiments A1-A10, wherein the longitudinal axis of the soft tissue is curved.
A12. The method of embodiment Al l, wherein identifying the oblique plane comprises:
drawing a curve substantially parallel to the longitudinal axis of the tissue; and evaluating, with the processor, an oblique plane normal to the curve at the point where the oblique plane intersects the curve.
A13. The method of embodiment Al 1 or A12, wherein a plurality of reformatted images is generated, each reformatted image being normal to the curve at the point where each reformatted image intersects the curve.
A14. The method of any of embodiments A1-A13, wherein generating the 3-D marking comprises generating a 3-D surface model from the matrix.
A15. The method of any of embodiments A1-A14, further comprising smoothing the 3-D marking.
A16. The method of any of embodiments A1-A15, wherein the 3-D model is a model of a patient's bone. A17. The method of any of embodiments A1-A16, wherein the 3-D model is a model of a patient's tibia.
A18. The method of any of embodiments A1-A17 wherein the 3-D model is displayed with the image data to orient a user.
A19. The method of any of embodiments A1-A18, wherein generating the set of reformatted images comprises:
identifying, with the processor, four corners of the oblique plane; generating, with the processor, texture coordinates using the corner locations; and
interpolating, with the processor, pixel intensities for the texture coordinates using the image data. A20. The method of any of embodiments A1-A19, wherein the soft tissue is an ACL, a MCL, a PCL, or a patellar tendon of a patient.
Bl . A system for locating anatomical landmarks in image data, comprising:
means for receiving image data representative of a patient's joint;
means for identifying an oblique plane intersecting a soft tissue represented in the image data and normal to a longitudinal axis of the soft tissue;
means for identifying a set of planes parallel or orthogonal to the oblique plane;
means for generating a set of reformatted images, from the image data, sampled along the set of planes;
means for tracing a feature of the soft tissue in the set of reformatted images; and means for generating a 3-D marking representing the traced feature on a 3-D model.
B2. The system of embodiment B 1 , wherein the means for tracing the feature comprises means for tracing the feature in two or more reformatted images parallel to the oblique plane.
B3. The system of embodiment Bl or B2, wherein the means for tracing the feature comprises means for advancing through parallel reformatted images selected from the set of reformatted images. B4. The system of any of embodiments B1-B3, wherein the means for tracing the feature comprises means for outlining a cross-section of the feature in a reformatted image from the set of reformatted images.
B5. The system of any of embodiments B1-B4, wherein the means for tracing the feature comprises means for storing data representing the tracing in a matrix.
B6. The system of any of embodiments B1-B5, wherein the means for generating the set of reformatted images comprises means for rotating the image data.
B7. The system of any of embodiments B1-B6, wherein the means for identifying the oblique plane comprises:
means for identifying a first point and a second point along the longitudinal axis of the soft tissue;
means for evaluating a line defined by the first and second points; and
means for evaluating a plane normal to the line.
B8. The system of any of embodiments B1-B7, wherein the first point is located about where the soft tissue attaches to a patient's bone.
B9. The system of any of embodiments B1-B8, wherein the second point is located at about the center of the soft tissue. BIO. The system of any of embodiments B1-B7, wherein the first point is located at about the start of the soft tissue and the second point is located at about the end of the soft tissue.
Bl 1. The system of any of embodiments B1-B10, wherein the longitudinal axis of the soft tissue is curved.
B12. The system of embodiment Bl 1, wherein the means for identifying the oblique plane comprises: means for drawing a curve substantially parallel to the longitudinal axis of the tissue; and
means for evaluating an oblique plane normal to the curve at the point where the oblique plane intersects the curve.
B13. The system of embodiment Bl 1 or B12, wherein a plurality of reformatted images is generated, each reformatted image being normal to the curve at the point where each reformatted image intersects the curve. B14. The system of any of embodiments B1-B13, wherein the means for generating the 3-D marking comprises means for generating a 3-D surface model from the matrix.
B15. The system of any of embodiments B1-B14, further comprising means for smoothing the 3-D marking.
B16. The system of any of embodiments B1-B15, wherein the 3-D model is a model of a patient's bone.
B17. The system of any of embodiments B1-B16, wherein the 3-D model is a model of a patient's tibia.
B18. The system of any of embodiments B1-B17 wherein the 3-D model is displayed with the image data to orient a user. B19. The system of any of embodiments B1-B18, wherein the means for generating the set of reformatted images comprises:
means for computing four corners of the oblique plane;
means for generating texture coordinates using the corner locations; and means for interpolating pixel intensities for the texture coordinates using the image data.
B20. The system of any of embodiments B1-B19, wherein the soft tissue is an ACL, a MCL, a PCL, or a patellar tendon of a patient. CI . A method for locating anatomical landmarks in image data, comprising:
receiving, at a processor, image data representative of a patient's joint;
identifying an oblique plane intersecting a soft tissue represented in the image data and normal to a longitudinal axis of the soft tissue;
generating, with the processor, a reformatted image, from the image data, sampled along the oblique plane;
tracing a feature of the soft tissue in the reformatted image; and
generating, with the processor, a 3-D marking representing the traced feature on a 3-D model.
[0048] Variations and modifications will occur to those of skill in the art after reviewing this disclosure. The disclosed features may be implemented, in any combination and subcombination (including multiple dependent combinations and subcombinations), with one or more other features described herein. The various features described or illustrated above, including any components thereof, may be combined or integrated in other systems.
Moreover, certain features may be omitted or not implemented.
[0049] Examples of changes, substitutions, and alterations are ascertainable by one skilled in the art and could be made without departing from the scope of the information disclosed herein. All references cited herein are incorporated by reference in their entirety and made part of this application.

Claims

What is Claimed is:
1. A method for locating anatomical landmarks in image data, comprising:
receiving, at a processor, image data representative of a patient's joint;
identifying an oblique plane intersecting a soft tissue represented in the image data and normal to a longitudinal axis of the soft tissue;
identifying, with the processor, a set of planes comprising the oblique plane and at least one other plane parallel or orthogonal to the oblique plane;
generating, with the processor, a set of reformatted images, from the image data, sampled along the set of planes;
tracing a feature of the soft tissue in the set of reformatted images; and
generating, with the processor, a 3-D marking representing the traced feature on a 3-D model.
2. The method of claim 1, wherein tracing the feature comprises tracing the feature in two or more reformatted images parallel to the oblique plane.
3. The method of claim 1 or 2, wherein tracing the feature comprises:
displaying a first reformatted image on a screen; and
displaying a second reformatted image on the screen, wherein the second reformatted image is parallel to and offset from the first reformatted image.
4. The method of any of claims 1-3, wherein tracing the feature comprises outlining a cross-section of the feature in a reformatted image from the set of reformatted images.
5. The method of any of claims 1-4, wherein tracing the feature comprises storing data representing the tracing in a matrix.
6. The method of any of claims 1-5, wherein generating the set of reformatted images comprises rotating the image data.
7. The method of any of claims 1-6, wherein identifying the oblique plane comprises: identifying a first point and a second point along the longitudinal axis of the soft tissue;
evaluating, with the processor, a line defined by the first and second points; and evaluating, with the processor, a plane normal to the line.
8. The method of any of claims 1-7, wherein the longitudinal axis of the soft tissue is curved.
9. The method of claim 8, wherein identifying the oblique plane comprises:
drawing a curve substantially parallel to the longitudinal axis of the tissue; and evaluating, with the processor, an oblique plane normal to the curve at the point where the oblique plane intersects the curve.
10. The method of claim 8 or 9, wherein a plurality of reformatted images is generated, each reformatted image being normal to the curve at the point where each reformatted image intersects the curve.
11. The method of any of claims 1-10, wherein the 3-D model is a model of a patient's bone.
12. The method of any of claims 1-11, wherein the 3-D model is a model of a patient's tibia.
13. The method of any of claims 1-12 wherein the 3-D model is displayed with the image data to orient a user.
14. The method of any of claims 1-13, wherein generating the set of reformatted images comprises:
identifying, with the processor, four corners of the oblique plane; generating, with the processor, texture coordinates using the corner locations; and
interpolating, with the processor, pixel intensities for the texture coordinates using the image data.
15. The method of any of claims 1-14, wherein the soft tissue is an ACL, a MCL, a PCL, or a patellar tendon of a patient.
PCT/US2015/019615 2014-03-10 2015-03-10 Methods and systems for identifying anatomical landmarks in image data WO2015138397A1 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US15/124,988 US10235757B2 (en) 2014-03-10 2015-03-10 Methods and systems for identifying anatomical landmarks in image data
EP15761756.4A EP3116389B1 (en) 2014-03-10 2015-03-10 Method for identifying anatomical landmarks in image data
US16/152,211 US10275884B2 (en) 2014-03-10 2018-10-04 Methods and systems for identifying anatomical landmarks in image data

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201461950756P 2014-03-10 2014-03-10
US61/950,756 2014-03-10

Related Child Applications (2)

Application Number Title Priority Date Filing Date
US15/124,988 A-371-Of-International US10235757B2 (en) 2014-03-10 2015-03-10 Methods and systems for identifying anatomical landmarks in image data
US16/152,211 Continuation US10275884B2 (en) 2014-03-10 2018-10-04 Methods and systems for identifying anatomical landmarks in image data

Publications (1)

Publication Number Publication Date
WO2015138397A1 true WO2015138397A1 (en) 2015-09-17

Family

ID=54072313

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2015/019615 WO2015138397A1 (en) 2014-03-10 2015-03-10 Methods and systems for identifying anatomical landmarks in image data

Country Status (3)

Country Link
US (2) US10235757B2 (en)
EP (1) EP3116389B1 (en)
WO (1) WO2015138397A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10235757B2 (en) 2014-03-10 2019-03-19 Smith & Nephew, Inc. Methods and systems for identifying anatomical landmarks in image data
CN110913790A (en) * 2019-09-23 2020-03-24 深圳市第二人民医院 Method for positioning anatomical insertion point of medial patellar ligament

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11903691B2 (en) * 2017-07-18 2024-02-20 Siemens Healthineers Ag Combined steering engine and landmarking engine for elbow auto align
CN109102501B (en) * 2018-08-06 2021-05-28 北京市商汤科技开发有限公司 Joint image processing method and image processing equipment
WO2020086912A1 (en) * 2018-10-26 2020-04-30 Intuitive Surgical Operations, Inc. Mixed reality systems and methods for indicating an extent of a field of view of an imaging device
EP3682804A1 (en) * 2019-01-18 2020-07-22 Samsung Electronics Co., Ltd. X-ray imaging apparatus and control method thereof
CN111466954A (en) * 2020-05-25 2020-07-31 武汉中旗生物医疗电子有限公司 Three-dimensional ultrasonic stereo dissection chart generation method and device

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040122309A1 (en) * 2002-12-19 2004-06-24 Deller Timothy W. Systems and methods for creating reformatted computed tomography images
US20050148850A1 (en) * 2003-11-13 2005-07-07 Lahm Ryan P. Clinical tool for structure localization
WO2006126150A2 (en) * 2005-05-25 2006-11-30 Koninklijke Philips Electronics N.V. Automatic extraction and display of the tracheobronchial plane
US20070253610A1 (en) * 1995-12-29 2007-11-01 Pieper Steven D Anatomical visualization and measurement system
US20110255764A1 (en) * 2010-04-15 2011-10-20 Roger Lin Orientating an oblique plane in a 3d representation

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5086401A (en) 1990-05-11 1992-02-04 International Business Machines Corporation Image-directed robotic system for precise robotic surgery including redundant consistency checking
US6205411B1 (en) 1997-02-21 2001-03-20 Carnegie Mellon University Computer-assisted surgery planner and intra-operative guidance system
EP2371291A1 (en) 2010-03-31 2011-10-05 Dr. med. Roman Winkler Imaging measuring method and device
EP3116389B1 (en) 2014-03-10 2021-07-14 Smith & Nephew, Inc. Method for identifying anatomical landmarks in image data

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070253610A1 (en) * 1995-12-29 2007-11-01 Pieper Steven D Anatomical visualization and measurement system
US20040122309A1 (en) * 2002-12-19 2004-06-24 Deller Timothy W. Systems and methods for creating reformatted computed tomography images
US20050148850A1 (en) * 2003-11-13 2005-07-07 Lahm Ryan P. Clinical tool for structure localization
WO2006126150A2 (en) * 2005-05-25 2006-11-30 Koninklijke Philips Electronics N.V. Automatic extraction and display of the tracheobronchial plane
US20110255764A1 (en) * 2010-04-15 2011-10-20 Roger Lin Orientating an oblique plane in a 3d representation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
WILSON: "Estimation of Tendon Moment Arms from Three-Dimensional Magnetic Resonance Images", ANNALS OF BIOMEDICAL ENGINEERING, vol. 27, 1999, pages 247 - 256, XP008060726, DOI: 10.1114/1.180

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10235757B2 (en) 2014-03-10 2019-03-19 Smith & Nephew, Inc. Methods and systems for identifying anatomical landmarks in image data
US10275884B2 (en) 2014-03-10 2019-04-30 Smith & Nephew, Inc. Methods and systems for identifying anatomical landmarks in image data
CN110913790A (en) * 2019-09-23 2020-03-24 深圳市第二人民医院 Method for positioning anatomical insertion point of medial patellar ligament
WO2021056149A1 (en) * 2019-09-23 2021-04-01 深圳市第二人民医院 Anatomical endpoint positioning method for medial patellofemoral ligament

Also Published As

Publication number Publication date
EP3116389A1 (en) 2017-01-18
US10235757B2 (en) 2019-03-19
US20190035081A1 (en) 2019-01-31
EP3116389B1 (en) 2021-07-14
US10275884B2 (en) 2019-04-30
EP3116389A4 (en) 2017-11-15
US20170076453A1 (en) 2017-03-16

Similar Documents

Publication Publication Date Title
US10275884B2 (en) Methods and systems for identifying anatomical landmarks in image data
US10675063B2 (en) System and method for image segmentation in generating computer models of a joint to undergo arthroplasty
US11704872B2 (en) Surface and image integration for model evaluation and landmark determination
JP5684702B2 (en) A method for image segmentation in generating a computer model of a joint undergoing arthroplasty
US11478207B2 (en) Method for visualizing a bone
Xiao et al. Sparse Dictionary Learning for 3D Craniomaxillofacial Skeleton Estimation Based on 2D Face Photographs
Hancharenka et al. Preoperative planning of pelvic and lower limbs surgery by CT image processing

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 15761756

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 15124988

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

REEP Request for entry into the european phase

Ref document number: 2015761756

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 2015761756

Country of ref document: EP