$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Tre set di volume dati raccolti da tre diverse tecniche (microCT, cryoET e cryoSXT) sono stati utilizzati per dimostrare tre importanti caratteristiche del banco da lavoro SuRVoS: modello di formazione, Super-regione segmentazione ed etichetta spaccare. I set di dati rappresentano un gruppo eterogeneo di risultati sperimentali, per ognuno dei quali parametri di elaborazione completo sono forniti (tabella 1).
Per dimostrare il training del modello utilizzando SuRVoS Workbench, è stato scelto un set di dati relativamente ad alto contrasto con la definizione dei limiti di regione. Questo set di dati del frutto di Galium aparine, o goosegrass, sono stati raccolti usando la tomografia a raggi x fase contrasto sulla I13-2 Diamond-Manchester Imaging Beamline presso Diamond Light Source, Chilton, Oxfordshire, Regno Unito. Il campione fresco è stato montato in aria su un goniometro base sopra la fase di rotazione, ad una distanza di campione-rivelatore di 30 mm. esposizione volte erano 0.10 s con lo spettro di raggio rosa che ha un'energia media di circa 22 keV. Le proiezioni sono stati raccolti attraverso 180° con un passo di 0,1 °. Tomografiche ricostruzioni sono state eseguite utilizzando Savu30,31 con il filtro di Paganin per contrasto di fase di propagazione-basato immagini32 seguita da ricostruzione di proiezione posteriore filtrata in ASTRA toolkit33 , 34. questi dati e ' stato poi ridimensionati con 2 x 2 x 2 binning per ridurre la dimensione del file prima in ingresso in SuRVoS Workbench.
In primo luogo, i dati di input (Figura 3A) è stati filtrati e bloccati (per rimuovere i valori di intensità superiore ed inferiore nei dati) (Figura 3B). In questo modo, lo sfondo e primo piano sono state fatte più facilmente distinguibili e la consistenza della struttura interna del frutto è stata accentuata. Successivamente, supervoxels sono stati costruiti sulla cima di dataset filtrato (Figura 3). Per valutare la qualità della supervoxels, sono stati visualizzati senza i dati per verificare che le informazioni pertinenti del dataset erano ben rappresentate da supervoxels (Figura 3D). Successiva, manuale annotazioni utilizzando supervoxels sono state fornite come dati di training su tre fette del volume (Figura 3E, colori scuri). Questi dati di formazione erano sufficienti a formare il classificatore per predire (colori chiari) aree corrispondenti al fondo (verde), setola di frutta (rosso), materiale di seme (viola) e che circonda la carne (blu). Perfezionamenti morfologici sono stati usati per ripulire le segmentazioni di fori di riempimento, aumento o la riduzione come necessari (Figura 3F). Il tempo totale trascorso per determinare i parametri appropriati e segmentare questo dataset era 2 h.
Per dimostrare la segmentazione Super regione utilizzando SuRVoS Workbench, un set di dati complessi e rumoroso è stato scelto15. Questo dataset è stato raccolto utilizzando cryoET presso il centro nazionale per l'Imaging macromolecolare presso il Baylor College of Medicine, Houston, TX, Stati Uniti d'America. Brevemente, le piastrine erano tuffo congelato il bagliore scaricata e griglie TEM oro carbonio holey fiducial-trattati. Serie di inclinazione sono stati raccolti da ± 65 ° con un incremento di 2°. La serie di inclinazione allora è stata ricostruita utilizzando la proiezione posteriore ponderata in IMOD35.
Dopo il caricamento dei dati in SuRVoS (Figura 4A), è stata selezionata un'area di interesse e un set di filtro appropriato è stato applicato. In questo caso, un filtro gaussiano seguito da un filtro di variazione totale con il contrasto bloccato è stato utilizzato per accentuare i bordi e le texture dei dati (Figura 4B). Successivamente, training del modello con input dell'utente minima basata su supervoxel è stato utilizzato per segmentare la piastrina dal fondo ghiaccio e carbonio. Quindi, segmentazione semi-manuale con megavoxels e supervoxels è stata utilizzata per segmentare gli organelli. Infine, il parametro di origine supervoxel è stato cambiato in un più debole filtro denoising e supervoxel forma è stata fatta più piccolo (Vedi tabella 1) al fine di preservare meglio i confini di microtubuli per segmentazione (Figura 4). Per entrambi gli organelli e i microtubuli, sono state utilizzate annotazioni manuale rapide ogni 5-10 fette di selezionare il supervoxels che descrivono la funzionalità di interesse (Figura 4 & 4E). Il tempo totale trascorso per determinare i parametri appropriati e segmentare la regione di interesse presentato era di 6 h.
Per dimostrare l'etichetta spaccare utilizzando SuRVoS Workbench, è stato scelto un set di dati con molti, diversi organelli. Questo set di dati raccolti utilizzando cryoSXT su beamline B24 a Diamond Light Source, Chilton, Oxfordshire, UK,36. Brevemente, cellule HEK293 sono state coltivate sulle griglie di cercatore d'oro, dimensioni appropriate fiducials oro sono stati aggiunti e la griglia è stato congelato utilizzando un EM con retro-retro macchiare con trampolino. Serie di inclinazione sono stati poi raccolti su un microscopio da ± 65 ° con un incremento di 0,5 °. La serie di inclinazione allora è stata ricostruita utilizzando la proiezione posteriore ponderata in IMOD35.
Dopo il caricamento dei dati in SuRVoS (Figura 5A), è stata selezionata un'area di interesse e un filtro appropriato variazione totale è stato utilizzato per migliorare i confini degli organelli in tutto il volume (figura 5B). Successivamente, organelli erano semi-manually segmentato utilizzando megavoxels e supervoxels e poi raffinato con fori di riempimento, chiusura e dilatazione per smussare gli angoli (Figura 5). Il tempo totale per determinare i parametri appropriati e segmentare la regione di interesse presentato era di 4 h. Una volta che la segmentazione è stata finalizzata, lo Splitter di etichetta è stato utilizzato per visualizzare ogni organello come oggetto nel dataset (Figura 5) e varie caratteristiche su ogni oggetto nella trama dei dati (Figura 5E). L'interfaccia di etichetta Splitter è interattiva, aggiornando il colore associato con ogni nuova classe di etichetta sia la visualizzazione che la trama di dati. Questo permette la creazione di varie regole in base alle caratteristiche intrinseche nei dati che possono essere utilizzati per separare gli oggetti in classi utili (Figura 5F).

Figura 1. Il Layout e le caratteristiche generali del banco da lavoro SuRVoS.
L'interfaccia GUI è sulla sinistra, mentre il riquadro di visualizzazione è sulla destra. Queste due zone sono separate da una colonna di strumenti e scorciatoie. La GUI è predisposta per guidare l'utente attraverso le fasi principali nella pre-elaborazione dei dati, scelta dei parametri di supervoxel e/o megavoxel, segmentare i dati e, se necessario utilizzando il training del modello, prima di esportare le segmentazioni. Il riquadro di visualizzazione può essere utilizzato in tre modalità: Visualizzazione base e segmentazione per visualizzare i datie qualsiasi filtri applicati e segmentare i dati, etichetta splitter per classificare gli oggetti in nuove etichette basate su aspetti inerenti ai dati e infine etichetta statistiche per misurare e visualizzare le caratteristiche degli oggetti segmentati. Per ognuna di queste modalità, il menu a discesa in alto a sinistra quali dati vengono visualizzati i controlli di angolo, e il dispositivo di scorrimento nella parte superiore controlla l'asse z. Le scorciatoie di strumento permette facile accesso controllo contrasto, trasparenza dei livelli, lo zoom, panning e tornare a "casa" nel riquadro di visualizzazione e strumenti di apertura per l'annotazione come descritto nel protocollo. Clicca qui per visualizzare una versione più grande di questa figura.

Nella figura 2. Super-regione gerarchia riduce la complessità di segmentazione di immagini.
Un'immagine del set di dati per la segmentazione della Berkeley (BSDS50037) è stata utilizzata per dimostrare le proprietà e gli effetti delle Super-regioni. L'immagine originale (a sinistra) è composto da migliaia di voxel, che sono poi riuniti in raggruppamenti adiacente, simili per creare qualche centinaio supervoxels (centro). Supervoxels possono essere raccolte anche in raggruppamenti adiacente, simili per creare alcune decine di megavoxels (a destra). Con ogni raggruppamento, la complessità del compito di segmentazione è diminuita, per risorse computazionali e di manuale. Cosa importante, ecco un esempio 2D, tuttavia, sia supervoxels che megavoxels sono 3D. Clicca qui per visualizzare una versione più grande di questa figura.

Nella figura 3. Elaborazione di un oggetto Dataset utilizzando una strategia di segmentazione di formazione modello microCT.
R. una sola fetta 2D dei dati grezzi. B. applicazione di un filtro bloccato variazione totale ai dati grezzi migliorato i confini tra i vari aspetti del corpo fruttifero. C. sono stati scelti i parametri appropriati supervoxel. D. una regione di interesse (scatola rossa in C) è indicata per dimostrare che i confini dei dati sono presenti nel supervoxels essi stessi. E. tre fette del volume con le annotazioni manuali delle varie aree del set di dati visualizzati in colori scuri (verde, rossi, blu e viola) e previsioni dopo aver eseguito il training del modello visualizzati negli stessi colori luce. F. la stesse tre fette con la segmentazione finale dopo le previsioni di formazione del modello sono state accettate. Barre della scala sono 1 mm. per favore clicca qui per visualizzare una versione più grande di questa figura.

Nella figura 4. Elaborazione di un oggetto Dataset utilizzando una strategia di segmentazione della Super-regione di cryoET.
R. una sola fetta 2D dei dati grezzi. B. una regione di interesse (riquadro rosso in A) con un insieme di strati filtro applicato per accentuare i confini degli organelli. C. esempio di annotazione di un organello utilizzando Super-regioni. Un singolo organello è mostrato con supervoxels sovrapposto con l'annotazione manuale utente visualizzato in nero (a sinistra) e il supervoxels scelto con tale annotazione in blu (a destra). D. esempio di annotazione di un microtubulo utilizzando Super-regioni. Una singola regione di microtubuli è mostrata con l'annotazione manuale utente visualizzato in nero (a sinistra) e il supervoxels scelto con tale annotazione in verde (a destra). E. la segmentazione finale con la piastrina segmentato da sfondo utilizzando il modello di formazione (Vedi tabella 1 per i dettagli) e vari organelli e microtubuli segmentati utilizzando una strategia di segmentazione della Super-regione. Colori non indicano tipi specifici organelli che sono mostrati qui prima classificazione. Scala bar in A, B ed E sono 1 μm e in C e D sono 0,5 μm. Clicca qui per visualizzare una versione più grande di questa figura.

Nella figura 5. Analisi di un Dataset utilizzando lo strumento di Splitter etichetta cryoSXT.
R. una sola fetta 2D dei dati grezzi. B. una regione di interesse (scatola rossa nella A) con un filtro di variazione totale applicato per accentuare gli organelli. C. la segmentazione finale con supervoxels sovrapposti. D. La parte di visualizzazione della barra di divisione etichetta con gli organelli classificati utilizzando le regole visualizzate in F. E. la parte di trama dello splitter etichetta visualizzare l'intensità media all'interno di ogni oggetto, con le regole visualizzate in F applicata. Ogni linea verticale lungo l'asse x rappresenta un unico oggetto ed è con codifica a colori per abbinare la classe a che è stato assegnato. F. esempi di regole di classificazione per separare i vari oggetti in base alle loro proprietà intrinseche. Barre della scala sono 1 μm. Clicca qui per visualizzare una versione più grande di questa figura.
| Nome / Dataset | Fonte | P1 | P2 | P3 | P4 |
| Filtro gaussiano | | Sigma | | | |
| Gamma / predefinito | - | [0,5, 10] / 1 | | | |
| CryoET (G1) | Dati grezzi | 1 | | | |
| CryoET (G2) | Dati grezzi | 2 | | | |
| Variazione totale | | Lambda | Spaziatura | # Iter | Morsetto |
| Gamma / predefinito | - | [0.1, 30] / 10 | [0.1, 10] / 1 | [50, 500] / 100 | - |
| MicroCT (TV1) | Dati grezzi | 10 | 1 | 100 | (1,-) |
| CryoET (TV2) | G1 | 7 | 1 | 200 | - |
| CryoET (TV3) | G2 | 10 | 1 | 100 | - |
| CryoSXR (TV4) | Dati grezzi | 7 | 1 | 100 | - |
| Sogliatura | | Vmin | Vmax | | |
| Gamma / predefinito | | | |
CryoET (TH1)
TV3
0
-
Centraggio gaussiana
Sigma
Gamma / predefinito
-
[0,5, 10] / 2
MicroCT (GC1)
TV1
2
Gaussiana normalizzazione
Sigma
Gamma / predefinito
-
[0,5, 10] / 2
MicroCT (GN1)
TV1
2
Laplaciano di gaussiana
Sigma
Trebbiare
Risposta
Gamma / predefinito
-
[0,5, 10] / 2
[Sì/No] / No
[Luminoso/scuro] / Bright
MicroCT (LG1)
TV1
2
No
Luminoso
Differenza di gaussiane
Sigma Init
Rapporto di Sigma
Gamma / predefinito
-
[0,5, 10] / 2
[1,1, 3] / 1.6
MicroCT (DG1)
TV1
2
1.6
CryoET (DG2)
TV3
2
1.6
Tensore di struttura part.
Sigma1
Zona di Sigma
Gamma / predefinito
-
[0,5, 10] / 2
[0,5, 10] / 2
CryoET (ST1)
TV3
2
2
Supervoxels
Forma
Spaziatura
Compattezza
Gamma / predefinito
-
[1, 10] / 10
[0.1, 5] / 1
[1, 200] / 20
MicroCT (SV1)
TV1
(10, 10, 10)
(1, 1, 1)
30
CryoET (SV2)
TV3
(10, 10, 10)
(1, 1, 1)
50
CryoET (SV3)
TV2
(3, 5, 5)
(1, 1, 1)
50
CryoSXT (SV4)
TV4
(10, 10, 10)
(1, 1, 1)
30
Megavoxels
Lambda
# Bidoni
Gamma
Gamma / predefinito
-
[0.01, 1] 0,1
[10, 200] / 20
Nessuno, auto o [0, 1] / nessuno
CryoET (MV1)
SV2
0.1
50
Automatico
TV1
CryoSXT (MV2)
SV4
0.4
50
Nessuno
TV4
Training del modello
Regione
Classificatore
Raffinatezza
Disponibile / Default
[voxel / supervoxel]
[Ensembles, SVM, modelli lineari Online] /
[Nessuno, Potts, aspetto] / aspetto
Ensemble - RF
microCT
TV3
SV1
Foresta di casuale:
Aspetto
TH1
GC1
-Albero #:
-Lambda:
GN1
[10, 100] / 100
[1, 500] / 10
LG1
100
50
DG1
cryoET
TV2 DG2 ST1
SV2
Foresta di casuale:
Aspetto
-Albero #:
-Lambda:
[10, 100] / 100
[1, 500] / 10
100
50
Raffinatezza di annotazione
Raggio
Gamma / predefinito
[1, 20] / 1
microCT
Apertura
5
Fori di riempimento
1
Dilatazione
2
cryoET
Apertura
3
Fori di riempimento
1
Dilatazione
2
cryoSXT
Apertura
3
Fori di riempimento
1
Dilatazione
2
Tabella 1. Parametri ottimizzati utilizzato per elaborare tutte le tre serie di dati (microCT, cryoET e cryoSXT).>
Per ogni parametro predefinito e una gamma di uso generale vengono fornite. In molti casi, i dati filtrati viene utilizzati come fonte per l'elaborazione a valle. In questi casi, l'abbreviazione è usata per denotare il nuovo dataset di origine. Ad esempio, G1 (i dati grezzi cryoET filtrata gaussiana) utilizzati come input durante un filtro di variazione totale per creare TV2. Informazioni viene presentate solo per aspetti del banco di lavoro che sono stati utilizzati per ogni set di dati di processo. Ad esempio, Training del modello non è stato utilizzato durante l'elaborazione di dataset cryoSXT presentato qui, quindi parametri non sono dato per questo.