$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Le immagini μCT longitudinali per questo studio sono state raccolte alle settimane 0 e 6 da difetti radiali critici di 3 mm in ratti adulti Charles River SASCO-SD trattati con uno scaffold a base di policaprolattone (PCL). Tutto l'uso degli animali è stato eseguito in conformità con i protocolli approvati dal Comitato per le risorse animali (UCAR) dell'Università di Rochester. La raccolta delle immagini μCT è stata eseguita utilizzando Scanco Medical VivaCT 40.
NOTA: Le fasi principali di questo protocollo sono suddivise in segmentazione dell'immagine μCT, allineamento, selezione e ritaglio del ROI, analisi e visualizzazione (Figura 1). Il protocollo per la segmentazione delle immagini μCT è adattato da Kenney et al. (2022)11.

Figura 1: Diagramma riassuntivo del flusso di lavoro. Le fasi del protocollo sono principalmente suddivise in segmentazione delle immagini, allineamento del modello, selezione del ROI e analisi del volume. Clicca qui per visualizzare una versione più grande di questa figura.
1. Segmentazione dell'immagine μCT
- Apertura di immagini all'interno del software
- Avvia il software Amira (di seguito denominato software di visualizzazione) e seleziona Open Data. All'interno della finestra pop-up, vai alla directory della cartella desiderata per selezionare il punto temporale iniziale .dcm file. Premere CTRL+A per selezionare tutte le immagini e fare clic su Apri immagini.
- Una finestra pop-up mostrerà le informazioni sul set di dati dell'immagine aperto; fare clic su OK per procedere.
- Apparirà un avviso sul ridimensionamento. Fare clic su Converti in float (l'opzione predefinita) per procedere.
- L'icona del set di dati apparirà all'interno della sezione Vista progetto . Per modificare il nome del set di dati, fare clic una volta sull'icona del set di dati e premere F2; una volta inserito il nome desiderato, fare clic su OK. Stabilire un modello di denominazione per identificare la scansione (ad esempio, 'Materiale dell'impalcatura ossea' - 'WKX').
- Filtraggio e soglia dell'immagine μCT
- Con le immagini caricate nella sezione Vista progetto , fare clic sul visualizzatore 2D predefinito Ortofetta per visualizzare gli orientamenti e le sezioni 2D. Una visualizzazione 2D apparirà nella finestra di visualizzazione a destra.
- Per visualizzare le immagini in 3D, fare clic con il pulsante destro del mouse sul set di dati, cercare Rendering volume e selezionare. All'interno della sezione Proprietà per il rendering del volume, regolare la soglia inferiore della mappa dei colori su 2500 e premere Invio.
NOTA: 2D e 3D possono essere attivati e disattivati facendo clic sul quadrato blu accanto a Ortogonale o Rendering volume, rispettivamente. Per la finestra di visualizzazione a destra, sono disponibili più opzioni di visualizzazione lungo la parte superiore della barra. Questi includono diversi strumenti per spostare, ruotare o ingrandire/rimpicciolire l'immagine; strumenti per modificare l'orientamento dell'immagine in viste XY, XZ e YZ e strumenti per misurazioni semplici.
- Se necessario, alcuni set di dati potrebbero richiedere la rimozione di artefatti o dati esterni all'area di interesse utilizzando una modifica del volume. Per aggiungere una modifica del volume, fare clic con il pulsante destro del mouse sul set di dati, cercare Modifica volume e selezionare.
- Fare clic su Modifica volume nella sezione Visualizzazione progetto . Nella sezione Proprietà , seleziona TabBox dalla prima finestra a discesa a destra di Strumento.
- Quando si selezionano le regioni da mantenere, regolare la TabBox facendo clic e trascinando gli angoli verdi utilizzando il cursore Interagisci (l'icona del cursore appuntita lungo la barra superiore della finestra di visualizzazione). La regolazione di TabBox potrebbe richiedere l'utilizzo di visualizzazioni alternative (XY, XZ o YZ) per identificare le aree da rimuovere.
- Una volta che la TabBox circonda l'area da mantenere, selezionare Taglia all'esterno per rimuovere l'artefatto o i dati all'esterno della TabBox. Nella vista del progetto verrà generato un nuovo set di dati modificato.
- Fare clic e trascinare la stringa di rendering del volume dal set di dati originale al set di dati modificato. In questo modo l'immagine visualizzata nella finestra di visualizzazione verrà sostituita con il set di dati modificato. All'interno della sezione Proprietà per il rendering del volume, regolare la soglia inferiore della mappa dei colori su 2500 e premere Invio.
- Fare clic con il pulsante destro del mouse sui dati modificati, cercare Filtro mediano e selezionare. Nella sezione Proprietà per il filtro mediano, selezionare 3D per l'interpretazione e fare clic su Applica. In questo modo verrà creato un nuovo set di dati filtrato (.filtered) nella sezione Visualizzazione progetto .
- Fare clic con il pulsante destro del mouse sui dati filtrati, cercare Soglia interattiva e selezionare. Nella sezione Proprietà per la soglia interattiva, selezionare 3D per Tipo di anteprima, regolare il valore inferiore per l'intervallo di intensità su 2500 e fare clic su Applica. In questo modo verrà creato un nuovo set di dati con soglia (.thresholded) nella sezione Visualizzazione progetto .
- Fare clic e trascinare la stringa di rendering del volume dal set di dati modificato al set di dati con soglia. Disattiva la soglia interattiva e attiva il lacerazione del volume (se disattivata).
- Segmentazione dell'immagine μCT
- Fare clic sul set di dati filtrato, quindi passare alla scheda Segmentazione facendo clic su Segmentazione sotto la barra dei menu. All'interno della scheda Segmentazione , i controlli sono simili a quelli della scheda Progetto .
NOTA: Ci sono alcuni strumenti aggiuntivi rilevanti per la finestra di visualizzazione all'interno della scheda Segmentazione . C'è un pulsante di zoom +/- nell'angolo superiore della finestra di visualizzazione che può essere utilizzato in modo che l'immagine riempia la finestra di visualizzazione. L'immagine può essere centrata nella finestra di visualizzazione utilizzando la barra di scorrimento sinistra/destra lungo la parte inferiore della finestra di visualizzazione. Sopra questa barra di scorrimento, è presente una seconda barra di scorrimento per spostarsi tra le fette di dati μCT.
- Per garantire che la scheda Segmentazione sia impostata correttamente, nella sezione Editor di segmentazione , l'immagine è impostata sul set di dati filtrato e il campo Etichetta è impostato sul set di dati con soglia. Nella sezione Materiali , fai doppio clic sul materiale denominato Materiale e rinominalo in Osso.
- Fai clic su Nuovo sotto l'Editor di segmentazione per creare un nuovo campo etichetta. Nella finestra pop-up, rinominarlo utilizzando una convenzione di denominazione simile (ad esempio, 'Materiale dell'impalcatura ossea' - 'WKX' - Marker) come prima e come Marker e fare clic su OK.
- Regolare la finestra di visualizzazione 2D utilizzando viste alternative (XY, XZ o YZ) per visualizzare una vista del piano sagittale o coronale dell'immagine.
- Usa pennello, lazo, bacchetta e soglia all'interno della sezione Selezione per identificare le aree di interesse. Per questo protocollo, selezionare lo strumento pennello . Regola le dimensioni dello strumento pennello ; In genere, un puntatore di medie dimensioni funziona bene.
- Regola la pila di immagini fino alla prima sezione in cui l'osso del radio inizia a essere visibile. Ciò potrebbe richiedere di spostarsi più volte tra le sezioni per determinare quali volumi appartengono al radio e all'ulna e quando iniziano.
- Una volta identificata la fetta, usa il cursore del pennello per disegnare facendo clic e trascinando il cursore attorno a questo segmento di osso. Apparirà una linea rossa nel punto in cui il cursore è stato disegnato. Disegna all'interno della sezione dell'osso, non all'esterno dell'osso; Questo aiuterà il software di visualizzazione a differenziare tra le due ossa.
- Continuate a spostarvi tra le fette e a disegnare attorno al raggio. Generalmente, attraverso la sezione centrale della pila di immagini, i disegni possono essere realizzati ogni 20-30 sezioni di immagine. Lungo l'estremità della pila di immagini o intorno al difetto di dimensioni critiche, crea disegni ogni 5-10 sezioni di immagine per aiutare il software di visualizzazione mentre l'osso svanisce in entrata o in uscita.
- Una volta completata l'intera lunghezza del raggio, fare clic su Selezione, Riempimento, quindi su Tutte le sezioni nel menu a discesa della barra dei menu. Questo riempirà tutti i disegni realizzati con lo strumento pennello. Questi appariranno ora ombreggiati con il rosso anziché con i contorni.
- Fare clic su Selezione e Interpolazione nel menu a discesa della barra dei menu. In questo modo verranno riempite le regioni ombreggiate del raggio osseo in tutte le sezioni in base ai disegni realizzati.
NOTA: A volte mancheranno alcune sezioni di osso. Prima di procedere scorrendo tutte le fette per assicurarsi di aver generalmente identificato l'osso. Se una sezione di osso è stata mancante, utilizzare lo strumento pennello per girare intorno a questa regione su alcune fette che circondano la sezione mancante e ripetere i due passaggi precedenti (1.3.9 e 1.3.10).
- Una volta che il raggio è stato completamente identificato, nella sezione Materiali , fai doppio clic sul materiale denominato Interno e rinominalo in Raggio.
- Fare clic sul simbolo Aggiungi (+) situato sotto la sezione Selezione per aggiungerlo come materiale. Il raggio è ora delineato nella finestra di visualizzazione. Fare clic sull'icona Blocca per il materiale Raggio nella sezione Materiali prima di procedere; Ciò garantirà che non vengano apportate modifiche a questo materiale.
- Completa lo stesso processo per l'osso dell'ulna. Nella sezione Materiali , fare clic su Aggiungi per creare un nuovo materiale. Fai doppio clic sul nuovo materiale e rinominalo Ulna.
- Ripetere i passaggi da 1.3.6 a 1.3.12 per l'osso dell'ulna.
- Una volta completati i passaggi per le ossa del raggio e dell'ulna, sblocca il raggio facendo clic sull'icona del lucchetto per il materiale del raggio nella sezione Materiali . Passa alla scheda Progetto facendo clic su Progetto sotto la barra dei menu. A questo punto è presente un set di dati denominato Marcatori applicato al set di dati filtrato nella sezione Visualizzazione progetto.
- Fare clic con il pulsante destro del mouse sul set di dati filtrato e cercare Marker Based Watershed Inside Mark (Image Segmentation) e selezionare. All'interno della sezione Proprietà , facendo clic sui menu a discesa, impostare i dati Dati come set di dati filtrati, i marcatori come set di dati marcatori e la maschera binaria come set di dati con soglia e fare clic su Applica.
- Verrà generato un file .grown all'interno della sezione Project View . Fare clic e trascinare la stringa di rendering del volume dal set di dati con soglia al set di dati .grown per visualizzare la segmentazione del radio e dell'osso dell'ulna all'interno della finestra di visualizzazione.
- Per convertire il file .grown in 8-bit, fai clic sul set di dati .grown e cerca Converti tipo immagine e seleziona. Il file .grown verrà convertito da 16 bit a 8 bit.
2. Allineamento
- Allineamento degli assi ed estrazione del radio e dell'ulna
NOTA: Se il modello attuale è il punto temporale iniziale, tutti i passaggi devono essere completati. In caso contrario, iniziare questa sezione dal passaggio 2.1.4. Per consentire una sezione trasversale del raggio, il modello del punto temporale iniziale deve essere allineato perpendicolarmente al piano della sezione orto; Questo modello servirà come base per l'allineamento delle settimane successive.
- Nella sezione Vista progetto , attivare la sezione Ortogonale (se disattivata) e quindi fare clic sul set di dati .grown. Nella sezione Proprietà , fare clic sull'icona Editor di trasformazioni . Utilizzando il cursore Interagisci, regola l'angolo del raggio e dell'ulna facendo clic sui punti di allineamento dell'asse verde in modo che la sezione ortogonale crei una sezione trasversale attraverso il raggio.
- Una volta che questo è stato allineato, è necessario creare un nuovo set di dati per salvare questa trasformazione. Fare clic con il pulsante destro del mouse sul set di dati .grown, cercare Ricampiona immagine trasformata e selezionare.
- Fare clic su Ricampiona immagine trasformata e, all'interno della sezione Proprietà, impostare il valore Dati sul set di dati .grown, Interpolazione sul vicino più vicino, Modalità su esteso, Mantieni su Dimensioni Voxel e Padding su 0. Fare clic su Applica per creare un nuovo set di dati .transformed.
- Per estrarre il raggio e l'ulna dal modello di segmentazione combinato, fare clic con il pulsante destro del mouse sul set di dati .transformed, cercare Estrai etichetta e selezionare.
- Fare clic su Estrai etichetta, all'interno della sezione Proprietà , impostare le etichette sul set di dati .transformed , l'ID etichetta su 1 e l'opzione Esporta in binario selezionata. Fare clic su Applica per creare un set di dati Risultato.
- Fare clic e trascinare la stringa di rendering del volume dal set di dati .grown al set di dati Result . Fare clic sul set di dati dei risultati e premere F2 per rinominare il file (ad esempio, 'Materiale dell'impalcatura ossea' - 'WKX' - Raggio).
- Ripetere i passaggi da 2.1.4 a 2.1.6 per l'ulna impostando l'ID etichetta su 2 e rinominando il file dei risultati per Ulna.
- Salvataggio file e apertura immagini longitudinali
NOTA: Se il modello attuale è il punto temporale iniziale, tutti i passaggi devono essere completati. In caso contrario, iniziare questa sezione dal passaggio 2.2.2.
- Per salvare il file del punto temporale iniziale, fare clic su File, Salva progetto con nome e Imposta posizione directory file nel menu a discesa. Impostare il tipo di salvataggio come file Amira Project e file di dati (pack &go) (*.hx) e denominare il file come versione master (ad esempio, 'Bone Scaffold Material'_ 'WKX' MASTER).
- Apri il file principale del punto temporale iniziale (se non è già aperto) e salva una copia facendo clic su File, Salva progetto con nome e Imposta posizione directory file nel menu a discesa. Impostare il tipo di salvataggio come file Amira Project e file di dati (pack & go) (*.hx) e denominare il file (ad esempio, 'Bone Scaffold Material'_ 'WK0' e 'WKX'). Questo file verrà utilizzato per confrontare il punto temporale iniziale con le settimane successive senza sovrascrivere il file master.
- All'interno del nuovo file di confronto, aprire le immagini per il punto temporale di confronto facendo clic su Apri dati nella sezione Vista progetto . All'interno della finestra pop-up, vai alla directory della cartella desiderata per selezionare il punto temporale iniziale .dcm file. Premere CTRL+A per selezionare tutte le immagini e fare clic su Apri immagini.
- Una finestra pop-up fornirà informazioni relative al set di dati dell'immagine aperto; fare clic su OK per procedere.
- Apparirà un avviso sul ridimensionamento. Fare clic su Converti in float (l'opzione predefinita) per procedere.
- L'icona del set di dati apparirà all'interno della sezione Vista progetto . Per modificare il nome del set di dati, fare clic una volta sull'icona del set di dati e premere F2; una volta inserito il nome desiderato, fare clic su OK. È possibile stabilire un modello di denominazione per identificare la scansione (ad esempio, "Materiale dell'impalcatura ossea" - "WKX").
- Allineamento del modello
NOTA: Questo processo discuterà l'allineamento delle segmentazioni ossee del radio. Lo stesso processo può essere applicato all'ulna, se necessario. Non allineare entrambe le ossa contemporaneamente; Ciò causerà problemi con il corretto allineamento.
- Disattiva qualsiasi punto temporale iniziale Rendering volume. Non viene visualizzato nulla all'interno della finestra di visualizzazione.
- Fare clic con il pulsante destro del mouse sul set di dati di confronto (settimana successiva), cercare Ortho Slice e selezionare.
- Ripetere la sezione 1 (a partire dal passaggio 1.2), la sezione 2 (passaggio 2.1) seguendo eventuali note aggiuntive per i numeri di passaggio forniti. Al termine di queste sezioni, riprendere dal passaggio 2.3.4. Il completamento di queste sezioni produrrà la segmentazione del set di dati della settimana successiva ed estrarrà le ossa del radio e dell'ulna.
- Fare clic con il pulsante destro del mouse sul set di dati di confronto (settimana successiva) per l'osso del radio estratto, cercare la procedura guidata di registrazione dell'immagine e selezionare. Nella sezione Proprietà , impostare Dati come set di dati della settimana di confronto per l'osso del radio estratto, Riferimento come set di dati del punto temporale iniziale per l'osso del radio estratto.
- Per la sezione Azioni della procedura guidata di registrazione delle immagini , fare clic su Ignora per il passaggio 1 di 4. Per i passaggi 2 e 3 di 4, utilizzando il cursore Interagisci, regolare la TabBox sull'area comune tra il punto temporale iniziale e i set di dati della settimana di confronto facendo clic su Applica in Azione dopo ogni passaggio. Per il passaggio 4 di 4, impostare Metrica come Correlazione, Trasformazione come Rigida e Pre-allineamento come Allinea assi principali e fare clic su Applica in Azione.
- Una volta allineati i set di dati, fare clic con il pulsante destro del mouse sul set di dati di confronto (settimana successiva) per l'osso del radio estratto, cercare Ricampiona immagine trasformata e selezionare. Nella sezione Proprietà, impostare i dati Dati come set di dati della settimana di confronto per l'osso del radio estratto, l'Interpolazione al vicino più vicino, la Modalità su Estesa, Mantieni su Dimensione Voxel e il Valore di spaziatura interna su 0 e fare clic su Applica. Verrà creato un nuovo set di dati .transformed per il set di dati di confronto.
3. Selezione e ritaglio del ROI
NOTA: Completare prima il ritaglio ROI per determinare il numero di fette attorno alla frattura di dimensioni critiche. Una volta determinati questi numeri di fetta, questi possono essere utilizzati nello stesso momento per l'ulna, se necessario. Per questo processo, il ritaglio non può essere annullato una volta applicato al modello.
- Regione di interesse coltura
- Fare clic su e attivare la sezione ortogonale per il punto temporale iniziale e impostare i dati sul set di dati del punto temporale iniziale per il raggio estratto. Impostate l'orientamento in modo che il piano produca un taglio trasversale attraverso l'osso del radio.
- Utilizzando la barra di scorrimento Numero sezione nella sezione Proprietà , modificare il numero della sezione per determinare le posizioni delle sezioni prossimali e distali che circondano il difetto di dimensioni critiche. Per entrambe le fette, determinare la fetta più distale e prossimale dove la frattura incontra la diafisi del radio osseo. Documentare il numero della sezione; Questo numero varia per ogni modello.
- Fare clic su e attivare la sezione ortogonale per la settimana di confronto e impostare i dati sul set di dati del punto temporale iniziale per il raggio estratto. Impostare l'orientamento in modo che il piano produca un taglio trasversale attraverso l'osso del radio.
- Utilizzando la barra di scorrimento Numero sezione nella sezione Proprietà e con il set di dati del punto temporale iniziale che mostra la sezione ortogonale distale, modificare il numero della sezione della settimana di confronto in modo che si allinei con la sezione distale del punto temporale iniziale. Il software di visualizzazione consente un controllo visivo dell'allineamento; Quando le sezioni sono sovrapposte, appariranno come una singola sezione nella finestra di visualizzazione.
- Annotare il numero della sezione distale del set di dati della settimana di confronto. Ripetere il passaggio 3.1.4 per la fetta prossimale.
NOTA: Come ulteriore controllo, determinare la differenza tra la sezione prossimale e distale sia per il set di dati iniziale che per quello della settimana di confronto. Questi valori devono essere equivalenti, illustrando un ROI di dimensioni coerenti.
- Fare clic sul punto temporale iniziale per il raggio estratto. Nella sezione Proprietà , fare clic sullo strumento Editor di ritaglio .
- All'interno del pop-up dell'editor di ritaglio , inserisci i valori minimo e massimo nel campo X, Y o Z appropriato (la finestra di visualizzazione modificherà il ROI man mano che vengono inseriti). Una volta inseriti i valori, fare clic su OK. Il set di dati verrà ritagliato in base al ROI di input.
- Ripetere il passaggio 3.1.7 per il set di dati della settimana di confronto.
4. Analisi e visualizzazione
- Analisi dei volumi
- Per determinare il volume del set di dati del punto temporale iniziale, fare clic con il pulsante destro del mouse sul set di dati del punto temporale iniziale trasformato per il raggio estratto, cercare Statistiche materiale e selezionare. Nella sezione Proprietà , impostare Dati come punto temporale iniziale trasformato per il raggio estratto, impostare Seleziona come materiali e fare clic su Applica. Un nuovo . Il set di dati MaterialStatistics verrà creato per il set di dati del punto temporale iniziale.
- Fare clic sul pulsante . Set di dati MaterialStatistics per il punto temporale iniziale per il raggio estratto e all'interno della finestra Proprietà , fare clic su Mostra foglio di calcolo. Fare clic sulla scheda Tabelle sopra la finestra Proprietà . All'interno di questa scheda, viene visualizzato il volume per il set di dati del punto temporale iniziale ritagliato per il raggio estratto.
- Ripetere i passaggi 4.1.1 e 4.1.2 per il set di dati della settimana di confronto per il raggio estratto. All'interno della scheda Tabella , il punto temporale iniziale e la settimana di confronto avranno schede separate in cui è possibile fare clic.
- Visualizzazione di set di dati
- Per visualizzare la variazione del volume osseo, fare clic con il pulsante destro del mouse sul set di dati trasformato della settimana di confronto per il raggio estratto, cercare Aritmetica e selezionare. Fare clic su Aritmetica e all'interno della finestra Proprietà , impostare l'input A come set di dati trasformato nella settimana di confronto per il raggio estratto, impostare l'input B come set di dati del punto temporale iniziale per il raggio estratto, impostare l'input C come NESSUNA ORIGINE, impostare il tipo di risultato come input A, lasciare l'opzione deselezionata, impostare i canali dei risultati come l'input A, e impostare Expression su A-B.
- Fare clic sul set di dati dei risultati e premere F2 per rinominare il file (ad esempio, 'Bone Scaffold Material' - 'WKX' - Bone Change). Fare clic con il pulsante destro del mouse su questo set di dati dei risultati, cercare Genera superficie e selezionare. Fare clic su Genera superficie e all'interno della finestra Proprietà , fare clic su Applica e, nella finestra pop-up, fare clic su Continua. Verrà creato un nuovo set di dati .surf.
- Fare clic con il pulsante destro del mouse sul set di dati .surf, cercare Vista superficie e selezionare. Una vista di superficie del risultato aritmetico verrà mostrata nella finestra di visualizzazione.
- Per modificare il colore della vista Superficie, fare clic sulla vista Superficie nella finestra Vista progetto , quindi all'interno della finestra Proprietà , fare clic sul menu a discesa per Colori e selezionare Costante. Fare clic sulla mappa dei colori e assegnare il colore preferito.
- Per visualizzare la variazione del volume osseo nel set di dati della settimana iniziale, fare clic con il pulsante destro del mouse sul set di dati .transformed, cercare Estrai etichetta e selezionare. Fai clic su Estrai etichetta e, nella sezione Proprietà , imposta le etichette sul set di dati .transformed , l'ID etichetta su 1 e l'opzione Esporta in binario selezionata. Fare clic su Applica per creare un set di dati Risultato . Fare clic sul set di dati dei risultati e premere F2 per rinominare il file (ad esempio, 'Materiale dell'impalcatura ossea' - 'WKX' - Completo - Raggio).
- Fare clic con il pulsante destro del mouse su questo set di dati dei risultati, cercare Genera superficie e selezionare. Fare clic su Genera superficie e, all'interno della finestra Proprietà , fare clic su Applica e nella finestra pop-up, fare clic su Continua. Verrà creato un nuovo set di dati .surf.
- Fare clic con il pulsante destro del mouse sul set di dati .surf, cercare Vista superficie e selezionare. Una vista di superficie del risultato aritmetico verrà mostrata nella finestra di visualizzazione.
- Per modificare il colore della vista Superficie, fare clic sulla vista Superficie nella finestra Vista progetto , quindi all'interno della finestra Proprietà, fare clic sul menu a discesa per Colori e selezionare Costante. Fare clic sulla mappa dei colori e assegnare il colore preferito.
- Per continuare con altre settimane, tornare al passaggio 2.2.2.