Segmentação de discos intervertebrais lombares para modelação e simulação computacional

Tese de Mestrado Integrado, Engenharia Biomédica e Biofísica (Engenharia Clínica e Instrumentação Médica), 2022, Universidade de Lisboa, Faculdade de CiênciasA lombalgia é a principal causa de incapacidade a nível mundial. A degeneração do disco intervertebral é uma das causas da lombalgia, podendo em casos avançados necessitar da remoção do disco intervertebral e substituição deste por um implante. Este implante pode consistir num dispositivo contendo enxerto ósseo (fusão espinhal) ou num disco intervertebral artificial (artroplastia discal). Ambos os métodos apresentam vantagens e desvantagens, pelo que é importante estudar, através de modelação e simulação em elementos finitos, a forma como implantes específicos afetam a biomecânica da coluna lombar antes de os inserir. Esta modelação personalizada requer a capacidade de segmentar as estruturas anatómicas relevantes a partir de imagens médicas. O presente trabalho teve como principal objetivo a implementação/desenvolvimento de um método para localizar e segmentar automaticamente discos intervertebrais lombares em 3D a partir de imagens de ressonância magnética em ponderação T2, com o intuito de auxiliar a construção de modelos de elementos finitos da coluna lombar a partir de casos reais, fornecendo informação precisa e personalizada sobre a forma dos discos intervertebrais do paciente. O desenvolvimento do método para permitir adicionalmente segmentar separadamente as duas principais estruturas do disco – anel fibroso e núcleo pulposo – e detetar automaticamente casos em que a degeneração não permite fazer esta distinção foi posteriormente seguido como objetivo secundário. O método de segmentação foi desenvolvido a partir de um método pré-existente na literatura para realização de segmentações 2D no perfil sagital, tendo este sido parcialmente implementado, modificado e adaptado para uso em 3D. O método permitiu realizar segmentações com uma exatidão média de 87.0 ± 3.7% medida pelo coeficiente de Dice em relação a segmentações manuais de referência. Esta eficácia é comparável com outros métodos de segmentação 3D na literatura. Este método apresenta a vantagem de ser significativamente mais rápido que a maioria dos métodos existentes, demorando apenas alguns segundos para completar uma segmentação dos discos lombares. O método para detetar degeneração discal identificou corretamente o estado de 96% dos discos (saudáveis e degenerados) com que foi testado.Back pain, especially in the lumbar spine, is the main cause of disability in the world. Intervertebral disc (IVD) degeneration is one of the causes of back pain. In some cases this requires the removal of the disc and its replacement with an implant. This implant may consist of either a cage containing bone graft (spinal fusion) or an artificial IVD (disc arthroplasty). Both of these treatments have advantages and disadvantages, which is why it is important to study, through computer modeling and finite element simulation, the ways in which specific implants affect the biomechanics of the lumbar spine before inserting them. This customized modeling requires the ability to segment the relevant anatomical structures from medical images. The present work had as its main objective the implementation/development of a method for localizing and automatically segmenting lumbar IVDs in 3D from T2 weighted magnetic resonance imaging, with the goal of supporting and complementing the generation of finite element models from real lumbar spines, by providing accurate and personalized information on the shape of the patient’s IVDs. The development of the method to also allow performing separate segmentations of the IVD’s two main structures – annulus fibrosus and nucleus pulposus – as well as automatically detecting degenerated IVDs where this distinction is no longer possible was later pursued as a secondary objective. The segmentation method was developed from a pre-existing method in the literature aimed at performing 2D segmentations in the sagittal profile, which was partially implemented, modified and adapted to 3D use. The method performed segmentations with a mean accuracy of 87.0 ± 3.7% as measured by the Dice coefficient in relation to manually segmented reference volumes, or ground truths. This method has the advantage of being significantly faster than most existing 3D segmentation methods, requiring only a few seconds to perform a complete segmentation of the lumbar discs. The method for detecting IVD degeneration correctly identified the status of 96% of the discs (healthy and degenerated) on which it was tested