ITBO20130422A1 - Metodo e apparato per generare una radiografia panoramica - Google Patents

Description

DESCRIZIONE

di una invenzione industriale dal titolo:

METODO E APPARATO PER GENERARE UNA RADIOGRAFIA PANORAMICA

L’invenzione riguarda un metodo ed un’apparecchiatura per radiografia dentale extraorale in grado di generare immagini tomografiche simili, per le loro caratteristiche, ad una immagine radiologica tomografica panoramica (nel seguito immagine panoramica). In particolare, l’invenzione riguarda un metodo per ottenere un’immagine panoramica analoga a quella di una radiografia panoramica ottenuta tramite raggi X, elaborando un volume digitale 3D in voxel ottenibile tramite tomografia computerizzata.

Ciascun voxel, nel volume 3D, presenta un rispettivo valore scalare che indica un assorbimento ai raggi X, ed à ̈ direttamente correlato alla densità dei materiali attraversati dai raggi X stessi. Usualmente, il valore di ciascun voxel in un volume 3D à ̈ rappresentato in un’immagine in livelli di grigio.

Le immagini panoramiche, a volte anche chiamate ortopantomografie, sono ben note nel campo dentale: si tratta di immagini radiografiche bidimensionali di una superficie curva predefinita approssimata alle mascelle di un paziente; tale immagine individua uno strato limitato in cui le strutture anatomiche sono a fuoco e al di fuori del quale le strutture anatomiche sono sfocate.

Le principali caratteristiche di una immagine panoramica sono:

- lo spessore dell’area sfocata intorno allo strato a fuoco;

- l’angolo di acquisizione e cioà ̈ l’angolo definito tra una direzione di acquisizione ed una direzione ortogonale allo strato a fuoco.

Queste caratteristiche possono anche variare da un pixel ad un altro dell’immagine panoramica stessa.

Un’immagine panoramica à ̈ normalmente prodotta usando uno specifico apparato per panoramiche simile a quello descritto ad esempio in EP2223653, il quale espone una testa di un paziente ai raggi X emessi da un emettitore di raggi X i quali sono ricevuti da un sensore, dotato di una matrice di ricevitori (o pixel) di poche colonne in larghezza, il quale à ̈ contraffacciato all’emettitore. Il sensore crea un’immagine radiologica parziale della porzione esposta ai raggi X della testa del paziente.

Facendo muovere trasmettitore e ricevitore lungo una traiettoria predefinita, l’immagine panoramica completa à ̈ prodotta acquisendo una pluralità di tali immagini radiologiche parziali lungo tale traiettoria ed inoltre sommando tali immagini parziali tra di loro, sovrapposte e spostate orizzontalmente di una opportuna distanza.

In alternativa, à ̈ anche possibile produrre un’immagine panoramica, che per sua stessa natura à ̈ bidimensionale, elaborando dati volumetrici 3D acquisiti con un apparato CBCT (Cone Beam Computerized Tomography).

In questo caso viene simulato un sensore di poche colonne in larghezza, tramite il quale può essere acquisita l’immagine parziale, e viene ripetuto virtualmente lo stesso metodo di acquisizione con spostamento del sensore e dell’emettitore di raggi X sulla traiettoria predefinita utilizzato nell’apparecchiatura per panoramiche descritta in precedenza. La differenza risiede nel fatto che ciascuna immagine parziale à ̈ calcolata in modo virtuale a partire dai dati in voxel del volume 3D.

Tale metodo à ̈ noto dal brevetto EP 2254475 Sirona, in cui viene simulato il metodo tomografico di sfocatura di un ortopantomografo e l’immagine panoramica à ̈ calcolata dal volume 3D precedentemente acquisito, facendo corrispondere ogni punto dell’immagine panoramica ad un voxel del volume 3D acquisito che viene selezionato come a fuoco.

Un problema delle immagini panoramiche acquisite sia tramite apparati per panoramiche sia tramite simulazione a partire da un volume 3D come in EP 2254475 SIRONA, à ̈ connesso al fatto che l’immagine panoramica che deve essere prodotta giace su una curva.

Ne segue che in ciascuna immagine parziale acquisita dal sensore di ricezione (reale o virtuale) solo una colonna di tale sensore, tipicamente la colonna centrale, risulterà inquadrare correttamente il punto voluto su tale curva, le altre colonne avranno invece un errore di posizionamento che, quando viene effettuata la somma delle immagini parziali, porta ad una leggera sfocatura anche nella curva stessa che si desidera essere a fuoco.

Un altro problema à ̈ che le principali caratteristiche di una immagine panoramica ottenuta sono collegate, per ciascun pixel, a quelle dei pixel adiacenti in quanto le immagini parziali sono sommate tra loro sovrapposte.

Un metodo alternativo per elaborare dati volumetrici 3D ed ottenere un’immagine 2D bidimensionale “simile†, come di seguito vedremo, ad una immagine panoramica à ̈ noto da EP1973075 CEFLA oppure da US8325874 Imaging Sciences International.

EP1973075 CEFLA illustra un metodo di elaborazione che prevede di estrarre uno strato di volume 3D di spessore definito giacente tra due diverse curve, una più interna alla testa del paziente ed una più esterna, e di calcolare un assorbimento in tale strato di volume 3D, elaborando valori di assorbimento in una direzione ortogonale ad una curva centrale tra la curva esterna e la curva interna. Grazie a questo metodo di elaborazione non c’à ̈ sfocatura, e quindi l’immagine 2D ottenibile à ̈ “simile†ad un immagine panoramica nel senso che à ̈ più “pulita†rispetto ad una immagine panoramica tradizionale in quanto alcune strutture anatomiche non risultano visibili.

Questo costituisce un problema, dal momento che gli operatori medici non sono familiari con un tipo di immagine panoramica 2D senza tali strutture anatomiche e preferiscono un’immagine panoramica tradizionale.

Scopo della presente invenzione à ̈ fornire un metodo di elaborazione ed un’apparecchiatura per ottenere immagini panoramiche 2D che abbiano una accuratezza migliorata partendo da acquisizioni 3D.

Questo scopo à ̈ ottenuto con un metodo ed un’apparecchiatura che hanno le caratteristiche delle rivendicazioni indipendenti. Forme realizzative vantaggiose e raffinamenti sono specificati nelle rivendicazioni dipendenti da queste.

Ulteriori vantaggi e proprietà della presente invenzione sono descritte nella descrizione seguente, in cui forme realizzative esemplari della presente invenzione sono spiegate in dettaglio sulla base dei disegni:

Figura 1 mostra un volume 3D acquisito tramite una tomografia computerizzata a fascio conico (CBCT) ed una superficie curva selezionata in tale volume 3D, che sarà la superficie a fuoco dell’immagine panoramica, in cui à ̈ mostrato un punto P identificato in tale superficie curva ed un insieme di rette passanti per tale punto e presentanti ciascuna un diverso angolo di incidenza, contenute entro un angolo solido avente il punto come vertice; Figura 2 mostra una illustrazione schematica di una sezione assiale dell’immagine di Figura 1, immersa in una griglia rappresentante una pluralità di voxel ed in cui à ̈ mostrata in pianta la superficie curva selezionata, una pluralità di primi punti identificati in tale superficie curva ed anche un ingrandimento di una porzione di tale superficie curva con tali primi punti identificati. Con riferimento a Figura 1, con 1 viene indicato un volume digitale 3D espresso in voxel il quale à ̈ acquisito tramite CBCT (Cone Beam Computerized Tomography). Ciascun voxel presenta un rispettivo valore di assorbimento ai raggi X espresso da un numero scalare il quale à ̈ visualizzato in Figura 1 come livelli di grigio. Si precisa che un volume 3D 1 analogo potrebbe essere anche acquisito tramite tecniche diverse rispetto alla CBCT, ad esempio tramite tomografia assiale computerizzata o TAC con fascio a ventaglio (fan beam) (in inglese CAT da computed axial tomography).

Una apparecchiatura (non illustrata) viene fornita per l’elaborazione del volume digitale 3D 1 in voxel, in modo da ottenere un’immagine panoramica analoga a quella di una radiografia panoramica ottenuta direttamente da un paziente tramite emissione di raggi X.

L’apparecchiatura comprende un’unità di elaborazione ed opzionalmente mezzi di acquisizione del volume 3D 1 stesso, che nel caso sopra esaminato possono essere di tipo CBCT oppure CAT.

In altre parole, l’apparecchiatura può comprendere sia i mezzi di acquisizione di un volume 3D 1 da un paziente, sia l’unità di elaborazione per elaborare il volume 3D 1 direttamente acquisito. In alternativa, l’apparecchiatura può essere solo dotata dell’unità di elaborazione, in quanto il volume 3D 1 da elaborare à ̈ già stato precedentemente acquisito. In quest’ultimo caso, l’unità di elaborazione potrebbe essere collocata anche a grande distanza dai mezzi di acquisizione del volume 3D 1 stesso.

In uso, l’unità di elaborazione dell’apparecchiatura esegue un metodo di elaborazione secondo la presente invenzione che consente di ottenere un’immagine panoramica analoga a quella di una radiografia panoramica ottenuta direttamente da un paziente tramite emissione di raggi X.

Il metodo della presente invenzione comprende la fase a) di selezionare una superficie curva 2 nel volume 3D 1.

Per selezionare tale superficie curva 2, può essere consentito ad un operatore umano di operare la selezione tramite un’opportuna interfaccia grafica, e/o prevedere una selezione automatica di tale superficie curva 2 mediante un algoritmo di ricerca di strutture anatomiche caratteristiche, e/o prevedere una forma predefinita della superficie curva 2.

Tutte le sopra indicate possibilità di selezione possono anche essere combinate in modo opportuno, ad esempio ad un operatore può essere consentito di modificare, tramite interfaccia operatore, una forma predefinita della superficie curva 2 e/o definita tramite algoritmo di ricerca automatica.

Inoltre, prima di tale fase a) di selezione della superficie curva 2, può essere opzionalmente previsto di rimuovere eventuali strutture indesiderate dal volume 3D acquisito 1, in modo da semplificare il volume 3D 1 stesso ai fini della sua successiva elaborazione.

Dopo aver selezionato la superficie curva 2, il metodo comprende la fase b) di identificare una pluralità di primi punti 3 sulla superficie curva 2 (Figura 2) e di collegare ciascun punto P di tali primi punti 3 identificati con un corrispondente pixel dell’immagine panoramica da ottenere, al fine di mappare i primi punti 3 della superficie curva 2 su corrispondenti pixel dell’immagine panoramica da ottenere.

Solo a scopo illustrativo, la superficie curva 2 di Figura 1 ha un’altezza approssimativa di 1.500 punti ed una larghezza di 3.000 punti. Con il presente metodo devono essere mappati 1.500 x 3.000 = 4.500.000 punti della superficie curva 2 su corrispondenti pixel dell’immagine panoramica da ottenere.

Il metodo comprende inoltre la fase c) di associare, per ciascun punto P di tale pluralità di primi punti 3 un rispettivo insieme di rette 4 passanti per tale punto P identificato, ciascuna retta presentando un diverso angolo di incidenza, tali rette essendo contenute entro un rispettivo angolo solido 5 avente il punto identificato P come vertice. Per ciascuna retta, il rispettivo angolo di incidenza à ̈ definito come angolo tra la retta stessa ed una retta perpendicolare X a tale superficie curva 2 passante per il punto identificato P.

Per variare le caratteristiche di qualità dell’immagine panoramica punto per punto oppure in porzioni specifiche, à ̈ possibile selezionare per ciascun punto identificato P oppure per zone determinate di tale superficie curva 2 il numero delle rette di tale insieme 4 ed un’angolazione specifica per ciascuna retta.

Ad esempio, una distribuzione delle rette 4 entro l’angolo solido 5 determina una profondità di campo dell’immagine panoramica la quale à ̈ inversamente proporzionale all’ampiezza di tale distribuzione.

Ciò consente di migliorare la qualità complessiva dell’immagine panoramica ottenuta in quanto à ̈ possibile una maggiore flessibilità nella elaborazione della stessa. La possibilità di variare punto per punto sulla superficie curva il numero di rette dell’insieme 4 ed il loro angolo d’incidenza, nonché opzionalmente la loro distribuzione, consente ad esempio di elaborare la superficie curva 2 in corrispondenza della zona di transizione tra la zona frontale e la zona dei canini in modo diverso rispetto alle rimanenti zone laterali dei molari. Tale zona di transizione à ̈ critica in quanto l’arcata dentale da acquisire à ̈ qui fortemente curva e quindi richiede un’attenzione diversa rispetto alla zona dei molari-premolari o degli incisivi.

Le porzioni di immagine panoramica corrispondenti a tali zone della superficie curva presenteranno quindi caratteristiche di qualità diverse.

Il metodo comprende inoltre la fase d) di calcolare un assorbimento lungo ciascuna di dette rette nel volume 3D, e cioà ̈ internamente al volume 3D, e la fase e) di attribuire al punto P identificato dei primi punti 3 un valore di assorbimento, calcolato come funzione degli assorbimenti calcolati per ciascuna delle rette dell’insieme 4 nell’angolo solido 5 associato al punto P identificato.

Il metodo comprende inoltre la fase f) di attribuire al pixel dell’immagine panoramica collegato al punto P identificato dei primi punti 3 un valore espresso in livelli di grigio, il quale à ̈ calcolato come funzione del valore di assorbimento del punto P identificato. Ad esempio, così come il valore di assorbimento à ̈ un numero scalare rappresentato graficamente in un immagine di tomografia assiale con livelli di grigio, analogamente il valore del pixel dell’immagine panoramica può assumere lo stesso valore scalare del valore di assorbimento calcolato, oppure può assumere un valore diverso, ad esempio proporzionale al valore di assorbimento ma non identico.

Abbiamo detto che nella fase d) del metodo viene calcolato per ciascuna retta dell’insieme 4 di rette un assorbimento nel volume 3D.

Per il calcolo di tale assorbimento viene innanzitutto selezionato per ciascuna retta un insieme di secondi punti (non rappresentati) nel volume 3D, in una fase d1). Ciascuna retta infatti passa attraverso il punto identificato P e si estende anche all’interno del volume 3D, definito da voxel.

Considerata quindi una retta specifica, per ciascun punto selezionato dei secondi punti nel volume 3D sono individuati i voxel circostanti il punto selezionato, in una fase d2). Viene inoltre calcolata in una fase d3) una media pesata del valore di assorbimento di tali voxel circostanti e al punto selezionato dei secondi punti della retta viene attribuito un valore di assorbimento pari a tale media pesata, in una fase d4).

I valori di assorbimento calcolati per ciascun punto selezionato dei secondi punti della retta considerata sono sommati, dando loro un peso dipendente dalla distanza dal punto identificato P che identifica uno strato da mettere a fuoco, e tale somma viene attribuita a tale retta come assorbimento nel volume 3D.

La media pesata calcolata nella fase d3) à ̈ ottenuta valutando la distanza del punto selezionato dei secondi punti rispetto al centro dei voxel circostanti il punto stesso. Nel dettaglio, per ciascuno di tali voxel circostanti il quale presenta un rispettivo centro, viene considerato il rispettivo valore di assorbimento in corrispondenza di tale centro e viene valutata una distanza tra tale centro ed il punto selezionato di tali secondi punti. La media pesata à ̈ ottenuta sommando ciascun valore di assorbimento di ciascun voxel moltiplicato per un coefficiente che à ̈ funzione della distanza del centro del voxel al punto selezionato di tali secondi punti.

In tale modo, ciascun punto dell’insieme dei secondi punti della retta ha un valore di assorbimento che tiene conto della posizione del punto stesso all’interno del volume 3D, e cioà ̈ la disposizione del punto rispetto ai voxel del volume 3D che lo circondano. Una volta calcolato l’assorbimento lungo ciascuna retta nel volume 3D, abbiamo detto che il metodo prevede di calcolare nella fase e) l’assorbimento nel punto P del primo insieme di punti 3 della superficie curva 2 come funzione di tutti gli assorbimenti calcolati per ciascuna delle rette.

In particolare tale fase e) comprende calcolare una media pesata dei valori di assorbimento delle rette dell’insieme di rette 4.

Come detto in precedenza, ciascuna retta à ̈ contraddistinta da un angolo di incidenza, definito in precedenza rispetto alla retta perpendicolare X. Per calcolare la media pesata, sono sommati i valori di assorbimento di ciascuna retta, ciascun valore di assorbimento moltiplicato per un coefficiente funzione dell’angolo di incidenza della retta stessa. Questo ad esempio consente di attribuire coefficienti maggiori a rette con piccoli angoli di incidenza rispetto a rette più angolate, per privilegiare gli angoli di incidenza sostanzialmente nulli e cioà ̈ acquisizioni ortogonali alla superficie curva 2 nel volume 3D.

Come abbiamo detto in precedenza, la fase b) comprende identificare una pluralità di primi punti 3 sulla superficie curva 2 ed il metodo di calcolo della presente invenzione à ̈ basato su elaborazioni fatte per ciascun punto P di tali primi punti 3.

Per identificare tali primi punti 3 nella superficie curva 2, sono eseguite sezioni di tale superficie curva 2 con superfici predefinite in una fase b1) e sono ottenute linee di intersezione 6 tra la superficie curva 2 e tali superfici predefinite in una fase b2). Tali primi punti 3 sono identificati nelle linee di intersezione 6 così definite.

La forma di queste superfici predefinite può essere piana o curva e sarà curva nel caso si vogliano applicare eventuali distorsioni all’immagine o in funzione delle caratteristiche di qualità dell’immagine panoramica che si intende ottenere.

Considerando un asse longitudinale A nel volume 3D, Ã ̈ possibile eseguire le sezioni in tale volume come nella fase b1) utilizzando piani paralleli ortogonali a tale asse longitudinale A, per ottenere quindi sezioni assiali di detto volume 3D.

I punti identificati tra la pluralità di primi punti 3 nelle linee di intersezione 6 possono essere scelti equidistanti tra loro oppure a distanza variabile, quest’ultima selezionata opportunamente come funzione delle caratteristiche di qualità dell’immagine panoramica che si intende ottenere.

In Figura 2 à ̈ mostrata una linea di intersezione 6 ottenuta tramite una sezione assiale e sono mostrati a puro titolo esemplificativo alcuni dei punti identificati come pluralità di primi punti 3. Se consideriamo due punti consecutivi identificati come Z e Y, estremità di un primo tratto 7 della linea di intersezione 6 ed ulteriori due punti consecutivi K e J estremità di un secondo tratto 8 di tale linea 6, si può notare come il primo tratto 7 abbia lunghezza L1 minore della lunghezza L2 del secondo tratto 8. La distanza tra i punti Z e Y à ̈ quindi minore della distanza tra i punti K e J e ciò dipende da un angolo di osservazione medio con il quale si vuole inquadrare la superficie curva 2, in corrispondenza del primo tratto 7 oppure del secondo tratto 8.

Secondo una forma realizzativa del metodo della presente invenzione, la superficie curva 2 da analizzare nel volume 3D comprende linee generatrici tra loro parallele, le quali sono inoltre parallele all’asse longitudinale del volume 3D. In altre parole, data una qualsiasi linea curva in una sezione assiale del volume 3D, la superficie curva 2 tramite la quale mappare l’immagine panoramica à ̈ quella che viene proiettata da tale linea curva secondo la direzione normale al piano di sezione stesso.

In questo caso specifico, la creazione dell’immagine panoramica può procedere considerando tutti i punti giacenti su uno di tali segmenti proiezione e cioà ̈ su una colonna. Tutti i punti identificati di una medesima colonna della superficie curva 2 saranno pertanto associati a corrispondenti pixel dell’immagine panoramica, e per ogni punto saranno considerate rette di diversa angolazione ma contenute in un angolo piano di ampiezza determinata.

Secondo questa forma di attuazione, il metodo di elaborazione risulta computazionalmente meno gravoso.

Unitamente al metodo, viene inoltre fornito un programma comprendente codice per implementare un metodo per ottenere un’immagine panoramica analoga a quella di una radiografia panoramica ottenuta tramite raggi X, in cui il metodo elabora un volume digitale 3D secondo quanto detto in precedenza ed in cui il programma viene eseguito in un’apparecchiatura descritta come in precedenza.

Grazie all’invenzione, à ̈ possibile ottenere un’immagine panoramica migliorata in quanto per ciascun punto P della superficie curva 2 viene calcolato un assorbimento all’interno del volume 3D separatamente ed indipendentemente dai rimanenti punti. Possono inoltre essere variati, per ciascun punto P identificato o per zone predefinite della superficie curva, i parametri di elaborazione del volume 3D, ciò consentendo elevata flessibilità di acquisizione delle strutture anatomiche stesse.

Si aggiunge inoltre che essendo l’immagine panoramica elaborata a partire dai voxel del volume 3D e non solo dai voxel di uno strato di esso, l’immagine panoramica ottenuta à ̈ familiare anche agli operatori medici a cui à ̈ destinata.

Il calcolo dell’assorbimento lungo l’insieme di rette prescinde infatti dai vincoli imposti da un emettitore ed un ricevitore di raggi X di dimensioni stabilite mobili lungo una traiettoria predefinita. L’acquisizione risulta quindi non più vincolata a tali dimensioni e l’immagine panoramica risulta migliorata.

Claims (4)

1. RIVENDICAZIONI 1. Metodo per ottenere un’immagine panoramica analoga a quella di una radiografia panoramica ottenuta da un paziente tramite emissione di raggi X, elaborando un volume digitale 3D acquisito (1) in voxel ciascun voxel presentando valori di assorbimento di raggi X, comprendente le seguenti fasi: a) selezionare una superficie curva (2) nel volume 3D (1); b) identificare una pluralità di primi punti (3) su detta superficie curva (2) e collegare ciascun punto di detti primi punti (3) identificati con un corrispondente pixel dell’immagine panoramica da ottenere, per mappare detti primi punti (3) di detta superficie curva (2) su corrispondenti pixel di detta immagine panoramica; caratterizzato dal fatto di c) associare per ciascun punto identificato (P) di detti primi punti (3) un rispettivo insieme di rette (4) passanti per detto punto identificato (P) e presentanti ciascuna un diverso angolo di incidenza, dette rette essendo contenute entro un angolo solido (5) avente detto punto identificato (P) come vertice; d) calcolare un assorbimento nel volume 3D (1) lungo ciascuna di dette rette; e) attribuire a detto punto identificato (P) di detti primi punti (3) un valore di assorbimento calcolato come funzione di tutti gli assorbimenti calcolati per ciascuna di dette rette di detto insieme (4); f) attribuire a detto pixel collegato a detto punto identificato (P) di detti primi punti (3) un valore espresso in livelli di grigio calcolato come funzione del valore di assorbimento di detto punto identificato (P).
2. 2. Metodo secondo la rivendicazione 1, in cui detto associare un rispettivo insieme di rette (4) in detta fase c) comprende selezionare in modo opportuno per ciascun punto identificato (P), e/o per zone determinate di tale superficie curva, un numero predeterminato di dette rette di detto insieme (4) ed uno specifico angolo di incidenza per ciascuna retta per variare caratteristiche di qualità di detta immagine panoramica.
3. 3. Metodo secondo una delle rivendicazioni da 1 a 2, in cui detto calcolare per ciascuna retta un assorbimento nel volume 3D lungo detta stessa retta in detta fase d) comprende: d1) selezionare un insieme di secondi punti di detta retta nel volume 3D (1); d2) per ciascun punto selezionato di detti secondi punti di detta retta individuare i voxel circostanti detto punto; d3) calcolare una media pesata del valore di assorbimento di detti voxel circostanti; d4) attribuire a detto punto selezionato di detti secondi punti un valore di assorbimento pari a detta media pesata; d5) sommare tutti i valori di assorbimento calcolati per ciascuno di detti secondi punti dando loro un peso dipendente dalla distanza da detto punto predefinito (P) ed attribuire detta somma a detta retta come assorbimento nel volume 3D. 4. Metodo secondo la rivendicazione 3, in cui detto calcolare una media pesata del valore di assorbimento di detti voxel circostanti in detta fase d3) comprende: - per ciascuno di detti voxel circostanti il quale presenta un rispettivo centro, considerare il rispettivo valore di assorbimento in corrispondenza di detto centro; - valutare una rispettiva distanza tra detto centro e detto punto selezionato di detti secondi punti; - sommare ciascun valore di assorbimento di ciascun voxel moltiplicato per un coefficiente funzione di detta distanza del centro di detto voxel da detto punto selezionato di detti secondi punti. 5. Metodo secondo una delle rivendicazioni precedenti, in cui detto attribuire a detto punto identificato (P) di detti primi punti (3) un valore di assorbimento in detta fase e) comprende calcolare una media pesata dei valori di assorbimento di dette rette di detto insieme di rette (4). 6. Metodo secondo la rivendicazione 5, in cui detto calcolare detta media pesata comprende sommare ciascun valore di assorbimento di ciascuna retta di detto insieme di rette (4) moltiplicato per un coefficiente funzione di detto angolo di incidenza, definito come angolo tra detta retta ed una retta perpendicolare (X) a detta superficie curva (2) e passante per detto punto identificato (P) di detti primi punti (3). 7. Metodo secondo una delle rivendicazioni precedenti, in cui detto selezionare una superficie curva (2) nel volume 3D (1) in detta fase a) comprende: - consentire ad un operatore umano di operare detta selezione tramite un’opportuna interfaccia grafica; - e/o prevedere una selezione automatica di tale superficie curva (2) mediante un algoritmo di ricerca di strutture anatomiche caratteristiche; - e/o prevedere una forma predefinita di detta superficie curva (2). 8. Metodo secondo una delle rivendicazioni precedenti, in cui detto identificare una pluralità di primi punti (3) su detta superficie curva in detta fase b) comprende: b1) eseguire sezioni di tale superficie curva con superfici predefinite; b2) ottenere linee di intersezione (6) tra detta superficie curva (2) e dette superfici predefinite ed identificare detta pluralità di primi punti (3) in dette linee di intersezione (6). 9. Metodo secondo la rivendicazione 8, in cui detto volume 3D (1) presenta un asse longitudinale (A) ed in cui detto eseguire dette sezioni in detta fase b1) comprende utilizzare piani paralleli ortogonali a detto asse longitudinale per ottenere sezioni assiali di detto volume 3D. 10. Metodo secondo una delle rivendicazioni da 8 o 9, in cui detto identificare detta pluralità di primi punti (3) in dette linee di intersezione (6) in detta fase b2) comprende selezionare punti tra loro equidistanti oppure selezionare punti a distanza variabile, detta distanza variabile essendo funzione delle caratteristiche della immagine panoramica che si intende ottenere. 11. Metodo secondo una delle rivendicazioni precedenti, e comprendente rimuovere eventuali strutture anatomiche indesiderate da detto volume 3D (1) prima di detto selezionare una superficie curva (2) in detto volume 3D (1) in detta fase a). 12. Apparecchiatura per l’elaborazione di un volume digitale 3D acquisito in voxel, ciascun voxel presentando un rispettivo assorbimento di raggi X, atta in uso ad applicare il metodo secondo una delle rivendicazioni 1-11. 13. Apparecchiatura secondo la rivendicazione 12, e comprendente mezzi di acquisizione di detto volume 3D ed un’unità di elaborazione di detto volume 3D.
4. 4. Programma comprendente codice per implementare un metodo per ottenere un’immagine panoramica analoga a quella di una radiografia panoramica ottenuta tramite raggi X, elaborando un volume digitale 3D acquisito secondo una delle rivendicazioni da 1 a 11 quando detto programma viene eseguito in una apparecchiatura secondo le rivendicazione 12 e 13.