$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Le protocole d’imagerie 3D PREFUL le plus simple comprend deux séquences. Un localisateur anatomique + une mesure de la respiration libre à l’aide de l’acquisition d’un empilement d’étoiles 3D à angle d’or. La trajectoire de l’empilement d’étoiles est représentée en Figure 1. Cette trajectoire combine l’échantillonnage radial (dans le plan) et l’échantillonnage cartésien linéaire (dans la direction z). Initialement, des échantillons dans la direction de la partition kz sont acquis, suivis de la rotation de l’angle d’or pour échantillonner le plan kx-k y. Le schéma d’acquisition de l’angle d’or favorise l’auto-déclenchement avec une couverture uniforme de l’espace k, facilitant la reconstruction d’images avec des données fortement sous-échantillonnées. L’ensemble de l’examen ne dépasse pas 10 minutes. De plus amples détails sur l’acquisition de 3D PREFUL sont présentés dans la Table des matériaux.
Figure 1 : Trajectoire de l’empilement d’étoiles tourné par le nombre d’or. Notez l’échantillonnage radial appliqué le long de la dimension dans le plan (kx-ky) et l’échantillonnage cartésien appliqué le long de la direction de la tranche (kz). Ce schéma se traduit par une couverture cylindrique de l’espace k acquis. Dans cet exemple, trois partitions, dont 32 projections radiales par partition, sont représentées. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Une liste des paramètres IRM typiques utilisés pour l’imagerie 3D PREFUL sur un scanner 1,5 T est présentée dans le Tableau 1B. Pour le scanner 3T, les paramètres de séquence énumérés dans le Tableau 2 sont proposés.
Une fois l’acquisition et le transfert de données terminés, la partie reconstruction de l’image commence. Cette étape est effectuée par un code MATLAB qui appartient à la propriété intellectuelle (PI) d’une entreprise (BioVisioneers GmbH) et ne peut pas être partagé ouvertement. Le script est entièrement automatisé une fois que le chemin d’accès aux données brutes est fourni. Une vue d’ensemble schématique de la procédure de reconstruction est illustrée dans Figure 2.
Figure 2 : Vue d’ensemble schématique de la méthode IRM PREFUL (Phase Resolved Functional Lung 3D). Dans un premier temps, les données sont acquises à l’aide d’une acquisition IRM en respiration libre de 8 minutes avec une trajectoire d’empilement d’étoiles. Une fois les données transférées du scanner IRM à une unité de calcul, des images 3D à basse résolution avec une résolution temporelle d’environ 100 ms sont reconstruites. Le parenchyme pulmonaire de chaque image est segmenté et le volume pulmonaire est calculé. L’information sur le volume pulmonaire est ensuite utilisée comme signal de déclenchement. En fonction de l’amplitude et de la phase du signal de déclenchement, les projections radiales sont triées dans des bacs respiratoires qui couvrent un cycle respiratoire. La partie binning est suivie de la reconstruction dynamique d’images en pleine résolution, qui sont ensuite enregistrées au niveau inspiratoire final. Après plusieurs étapes de post-traitement, le cycle de ventilation régionale (RVent) est calculé, et la partie analyse d’image, y compris l’extraction des paramètres, est effectuée. Le cycle RVent est évalué en calculant la boucle débit-volume (FVL) pour chaque voxel. Les cartes de ventilation FVL-CM sont extraites par une évaluation de chaque FVL de voxel par rapport à une FVL de référence saine, la similarité étant évaluée à l’aide d’une métrique de corrélation croisée. Pour RVent et FVL-CM, les valeurs globales du pourcentage total de défauts de ventilation (VDP) sont quantifiées. De plus, la dynamique du cycle RVent est analysée par une analyse du temps jusqu’au pic, ce qui permet d’obtenir une carte des paramètres du temps jusqu’au pic de ventilation (VTTP). De plus, l’écart de la ventilation de pointe attendue à 50 % du cycle RVent est quantifié dans la carte VTTPDev. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Après la procédure de reconstruction d’image, toutes les images sont spatialement alignées sur un seul état fixe. L’enregistrement lui-même est effectué sur le processeur à l’aide des outils de normalisation avancés (ANTs27) ou du package d’enregistrement Forsberg28. Alors que le package ANTs reste une référence pour les tâches de recalage d’images en IRM et en imagerie CT, le recalage Forsberg facilite une procédure de recalage 10 fois plus rapide jusqu’à 9 minutes avec des résultats comparables 19. Le package d’inscription peut être choisi en fonction des priorités de l’utilisateur et des unités de calcul disponibles. Une fois l’enregistrement terminé, les images morphologiques enregistrées sont à nouveau inversées de sorte que le poumon apparaît sombre dans l’image en niveaux de gris.
Dans Figure 3, l’étape suivante du pipeline, qui implique la segmentation du parenchyme pulmonaire, est représentée. Tout d’abord, le parenchyme pulmonaire est segmenté à partir de l’image inspiratoire finale à l’aide d’un réseau neuronal convolutif avec l’architecture nnUnet architecture29. La segmentation du parenchyme pulmonaire est suivie d’une reconnaissance des vaisseaux30, qui est exclue du masque de segmentation final.
Figure 3 : Résultats exemplaires de la segmentation basée sur l’apprentissage profond pour un sujet masculin de 32 ans. La rangée du haut affiche des images morphologiques de huit coupes coronales représentatives. Dans les deuxième et troisième rangées, on peut observer le masque de parenchyme pulmonaire correspondant et le masque final à l’exclusion des vaisseaux, respectivement. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Un exemple de segmentation incorrecte du parenchyme pulmonaire est illustré dans Figure 4. Il est essentiel d’inspecter visuellement les segmentations basées sur l’apprentissage profond et, si elles sont jugées insatisfaisantes, des corrections manuelles doivent être envisagées pour améliorer la précision du masque final du parenchyme pulmonaire.
Figure 4 : Résultats exemplaires d’une segmentation incorrecte basée sur l’apprentissage profond pour un sujet masculin de 57 ans. La rangée du haut affiche des images morphologiques de huit coupes coronales représentatives. Dans les deuxième et troisième rangées, on peut observer le masque de parenchyme pulmonaire correspondant et le masque final à l’exclusion des vaisseaux, respectivement. Comme on peut le constater, plusieurs régions fibrotiques sont reconnues à tort comme des vaisseaux ou des structures non pulmonaires. Ces inexactitudes ont été corrigées manuellement, comme le montre la quatrième rangée. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Après plusieurs étapes de filtrage décrites dans le protocole, les substituts de ventilation sont calculated. 3D PREFUL MRI génère des cartes quantitatives pour la ventilation régionale statique (RVent), des métriques de corrélation dynamique débit-volume (FVL-CM) et deux paramètres basés sur l’analyse du temps de ventilation jusqu’au pic (VTTP), comme illustré dans Figure 5. Cette figure présente huit tranches coronales d’un homme en bonne santé de 32 ans. Notez la distribution homogène attendue de tous les paramètres de ventilation.
Figure 5 : IRM préliminaire 3D d’un volontaire en bonne santé (mâle de 32 ans). Cartes morphologiques représentatives (rangée du haut) et des paramètres de ventilation IRM 3D Preful (de la deuxième à la cinquième rangée) pour un volontaire sain (mâle de 32 ans). La ventilation régionale statique est représentée par la ventilation régionale (RVent), tandis que la dynamique de ventilation est évaluée à l’aide de la métrique de corrélation de boucle débit-volume (FVL-CM), le temps de ventilation jusqu’au pic (VTTP) et l’écart de VTTP (VTTPDev) représentent les paramètres de ventilation évaluant la dynamique de ventilation. Comme prévu, des valeurs de ventilation homogènes sont observées pour toutes les cartes de paramètres de ventilation. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Pour simplifier les cartes de ventilation, des cartes de défauts de ventilation (VD) sont dérivées pour RVent et FVL-CM, ce qui permet une interprétation plus rapide des résultats. Les exemples de cartes VD sont présentés en Figure 6. Pour VDRVent et VDFVL-CM, les valeurs de VDP étaient respectivement de 3,6 % et 3,0 %, ce qui se situe dans la plage normale saine. Idéalement, et en fonction de l’âge des volontaires sains, la valeur VDP ne devrait pas dépasser 10 %. Les paramètres mentionnés ci-dessus (RVent, FVL-CM et leurs cartes VD) ont été validés dans plusieurs études11,12,35,36 et sont sensibles dans la détection des maladies ainsi que dans la détection des effets induits par le traitement20,34,37,38,39.
Figure 6 : Cartes représentatives RVent et FVL-CM, y compris leurs cartes de défauts de ventilation. Les cartes RVent (rangée du haut) et FVL-CM (deuxième rangée), y compris leurs cartes des défauts de ventilation (troisième et quatrième rangées) obtenues par IRM 3D Preful pour un volontaire en bonne santé (homme de 32 ans), sont représentées. Sur les cartes VD, les régions saines sont vertes et les zones déficitaires en ventilation sont marquées en rouge. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Plusieurs études d’IRM pulmonaire ont été rapportées à des intensités de champ magnétique de 1,5 T et de 3 T. Bien qu’il y ait l’avantage théorique de 3T en raison de l’augmentation du rapport signal/bruit (SNR), cet avantage peut être compensé par des effets de susceptibilité magnétique plus prononcés à 3T. L’influence exacte de l’intensité du champ magnétique sur la qualité d’image des cartes de ventilation 3D PREFUL est actuellement inconnue. Ici, nous présentons les résultats de faisabilité (Figure 7) obtenus pour un volontaire sain (mâle de 35 ans) à l’aide d’un scanner IRM 3 T. Notons l’aspect plus hétérogène des paramètres de ventilation à 3 T par rapport à 1,5T (Figure 5).
Figure 7 : IRM préful 3D d’un volontaire en bonne santé (mâle de 35 ans). Cartographies morphologiques représentatives (rangée du haut) et des paramètres de ventilation IRM 3D Preful (de la deuxième à la cinquième rangée) pour un volontaire sain (mâle de 35 ans) à 3T. La ventilation régionale statique est représentée par la ventilation régionale (RVent), tandis que la dynamique de ventilation est évaluée à l’aide de la métrique de corrélation de boucle débit-volume (FVL-CM), le temps de ventilation jusqu’au pic (VTTP) et l’écart de VTTP (VTTPDev) représentent les paramètres de ventilation évaluant la dynamique de ventilation. Veuillez cliquer ici pour voir une version agrandie de cette figure.
La comparaison/validation de PREFUL 3D dérivé avec des mesures plus directes n’est actuellement pas publiée. Plusieurs études ont rapporté une corrélation positive entre les valeurs VDP et le chevauchement spatial des défauts de ventilation entre la technique 2D PREFUL et 129Xe11,12,14 et 19F MRI35. Une comparaison récente entre la ventilation dérivée de PREFUL 3D et les mesures de ventilation directe à l’aide de l’IRM 19F chez des patients atteints de bronchopneumopathie chronique obstructive (BPCO), d’asthme et de volontaires sains a révélé une corrélation modérée à forte au niveau mondial 40. Il s’agit d’un exemple de comparaison (Figure 8) pour un patient diagnostiqué avec une maladie pulmonaire obstructive chronique (femme de 54 ans, VEMS1 = valeur prédite de 42 %), qui a été examiné à l’aide de 3D PREFUL et d’une IRM de lavage à 19F.
Figure 8 : IRM préliminaire 3D d’une patiente de 54 ans atteinte de BPCO. Des images morphologiques (rangée du haut), des paramètres de ventilation et des cartes des défauts de ventilation respectifs d’une patiente de 54 ans atteinte de BPCO (VEMS1 = 42 % de préd., CVF = 102 % de préd.) obtenues par IRM 3D Preful (2e et 3e rangées) et IRM 19F (4e et 5e rangées) sont présentées, ainsi qu’une comparaison des cartes de défauts de ventilation des deux méthodes (dernière rangée). Dans l’ensemble, sur l’ensemble des tranches, le coefficient de Sørensen-Dice était de 25,5 % dans les zones défectueuses et de 80,6 % dans les régions saines, ce qui a donné un chevauchement spatial total de 69,2 %. Notamment, il existe une corrélation visuelle entre les zones de ventilation appariées dans les régions saines et défectueuses (représentées en vert foncé). Veuillez cliquer ici pour voir une version agrandie de cette figure.
3D PREFUL IRM peut être utilisé pour mesurer les réponses régionales à therapy20. À titre d’exemple, > montre une comparaison de trois tranches du poumon d’un sujet atteint de mucoviscidose (une femme de 43 ans) au départ (à gauche, VEMS1 = valeur prédite à 94 %) et après un traitement par modulateur CFTR (à droite, VEMS1 = valeur prédite à 112 %). Étant donné que les deux mesures sont appariées spatialement, 3D PREFUL permet l’analyse régionale de la réponse au traitement, comme illustré au bas de Figure 9. Veuillez noter l’augmentation des valeurs de la métrique de corrélation débit-volume-boucle ainsi que la réduction des défauts de ventilation après le traitement.
Figure 9 : Mesures 3D PREFUL d’une patiente de 43 ans atteinte de mucoviscidose. Cartes exemplaires des marqueurs de ventilation des mesures 3D PREFUL de base (à gauche) et post-traitement (à droite) d’une patiente de 43 ans atteinte de mucoviscidose. LaVDP FVL-CM a diminué, passant de 18,0 % (référence) à 3,6 % (post-traitement). VEMS1 % préd. de référence : 94 %, VEMS1 % préd. post-traitement : 112 %. ICV de référence : 10,48, post-traitement ICV : 9,39. Les cartes de réponse au traitement correspondantes sont illustrées en bas. Veuillez noter les zones vertes indiquant la ventilation résolue après le traitement. Veuillez cliquer ici pour voir une version agrandie de cette figure.
Tableau 1 : Liste des paramètres IRM typiques utilisés pour l’acquisition d’un radiophare d’alignement de piste et d’un PREFUL 3D sur un scanner 1.5T. (A) Alignement de piste et (B) Acquisition 3D PREFUL. Veuillez cliquer ici pour télécharger ce tableau.
Tableau 2 : Paramètres IRM proposés pour l’acquisition de PREFUL 3D sur un scanner 3T. Veuillez cliquer ici pour télécharger ce tableau.
Tableau 3 : Rapport statistique exemplaire des paramètres 3D PREFUL pour un volontaire sain (homme de 35 ans) Veuillez cliquer ici pour télécharger ce tableau.