Segmentación de un hueso de la rodilla en 3D
Para cuantificar y medir las propiedades de un componente en un volumen, la segmentación es un primer paso necesario. Para segmentar el tejido óseo en un volumen de MRT, un algoritmo de agrupamiento se utiliza para lograr una segmentación áspera y aplicar un algoritmo grow-cut para obtener el resultado final.
Out[1]= | |
El volumen es procesado usando MedianFilter para regularizar el ruido. Una segmentación inicial se realiza agrupando intensidades de vóxel utilizando ClusteringComponents. Esto divide los datos en tres regiones: espacio vacío, tejido muscular y el resto, que incluye hueso, grasa y piel.
Out[2]= | |
Extraiga el segmento con mayor intensidad media, que representa hueso, grasa y tejido cutáneo.
Out[3]= | |
Out[4]= | |
GrowCutComponents pueden proporcionar una agradable segmentación final. Con el fin de crear dos marcadores para áreas de hueso y sin hueso, uno puede usar operaciones morfológicas.
El marcador de hueso puede ser calculado al erosionar el segmento usando un radio 4, el cual elimina toda la piel delgada y las capas de grasa y también erosiona parte del hueso.
Out[5]= | |
Para los marcadores de áreas sin hueso, dilatan el centro del hueso más de lo que ha sido erosionado y extraen la superficie alrededor del volumen dilatado.
Out[6]= | |
Ejecute el algoritmo grow-cut en 3D para refinar la segmentación.
Out[7]= | |
Visualice la segmentación.
Out[8]= | |
Las medidas de tamaño y densidad pueden ser calculadas para cada segmento.
Out[9]= | |
Out[37]//TraditionalForm= |
| |