Click here to Skip to main content
Click here to Skip to main content
Add your own
alternative version

Medical image visualization using WPF

, 27 Sep 2012
The article demonstrates the visualization of medical images (DICOM) using WPF.
using System;
using System.Windows;
using System.Windows.Media;
using System.Windows.Media.Animation;
using System.Windows.Media.Imaging;
using System.Windows.Media.Media3D;
using System.Xml.Linq;

using DICOMViewer.Helper;

// Implementation of the Image Slice for the Image Flow View using Windows Presentation Foundation (WPF).
// Major parts of the code have been taken from the 'WPF Cover Flow Tutorial':
// http://d3dal3.blogspot.com/2008/10/wpf-cover-flow-tutorial-part-1.html

namespace DICOMViewer.ImageFlow
{
    class ImageSlice : ModelVisual3D
    {
        private XDocument mXDocument = null;
        private string mFileName = null;
        private string mZValue = null;
        private ModelVisual3D m3DModel = null;
        private Model3DGroup mModelGroup = null;
        private AxisAngleRotation3D mRotation = null;
        private TranslateTransform3D mTranslation = null;
        private WriteableBitmap mBitmap = null;
        private ImageBrush mImageBrush = null;

        // Create new ImageSlice
        public ImageSlice(XDocument theXDocument, string theFileName, string theZValue, ModelVisual3D the3DModel)
        {
            mXDocument = theXDocument;
            mFileName = theFileName;
            mZValue = theZValue;
            m3DModel = the3DModel;
        }

        public void Animate(double aRotationAngle, double aTranslationX, double aTranslationY, double aTranslationZ, TimeSpan aAnimationDuration)
        {
            var rotateAnimation = new DoubleAnimation(aRotationAngle, aAnimationDuration);
            var xAnimation = new DoubleAnimation(aTranslationX, aAnimationDuration);
            var yAnimation = new DoubleAnimation(aTranslationY, aAnimationDuration);
            var zAnimation = new DoubleAnimation(aTranslationZ, aAnimationDuration);

            if(mRotation != null)
                mRotation.BeginAnimation(AxisAngleRotation3D.AngleProperty, rotateAnimation);

            if (mTranslation != null)
            {
                mTranslation.BeginAnimation(TranslateTransform3D.OffsetXProperty, xAnimation);
                mTranslation.BeginAnimation(TranslateTransform3D.OffsetYProperty, yAnimation);
                mTranslation.BeginAnimation(TranslateTransform3D.OffsetZProperty, zAnimation);
            }
        }

        public void SetBitmap()
        {
            if (mBitmap == null)
            {
                CTSliceInfo CTSliceInfo = new CTSliceInfo(mXDocument, mFileName);

                mBitmap = CTSliceInfo.GetPixelBufferAsBitmap();
                mBitmap.Freeze();
                
                mImageBrush = new ImageBrush();
                mImageBrush.ImageSource = mBitmap;
                
                mModelGroup = new Model3DGroup();
                mModelGroup.Children.Add(new GeometryModel3D(Tessellate(), new DiffuseMaterial(mImageBrush)));
                this.Content = mModelGroup;
                
                m3DModel.Children.Add(this);
            }
        }

        public void ResetBitmap()
        {
            if (mImageBrush != null)
                mImageBrush.ImageSource = null;

            mImageBrush = null;
            mBitmap = null;
            m3DModel.Children.Remove(this);
        }

        new public void Transform(double aRotationAngle, double aTranslationX, double aTranslationY, double aTranslationZ)
        {
            mRotation = new AxisAngleRotation3D(new Vector3D(0, 1, 0), aRotationAngle);
            mTranslation = new TranslateTransform3D(aTranslationX, aTranslationY, aTranslationZ);

            var transformGroup = new Transform3DGroup();
            transformGroup.Children.Add(new RotateTransform3D(mRotation));
            transformGroup.Children.Add(mTranslation);

            if(mModelGroup != null)
                mModelGroup.Transform = transformGroup;
        }

        public string FileName { get { return mFileName; } }
        public string ZValue { get { return mZValue; } }

        private static Geometry3D Tessellate()
        {
            var p0 = new Point3D(-1, -1, 0);
            var p1 = new Point3D(1, -1, 0);
            var p2 = new Point3D(1, 1, 0);
            var p3 = new Point3D(-1, 1, 0);

            var q0 = new Point(0, 0);
            var q1 = new Point(1, 0);
            var q2 = new Point(1, 1);
            var q3 = new Point(0, 1);

            return BuildMesh(p0, p1, p2, p3, q0, q1, q2, q3);
        }

        private static Vector3D CalculateNormal(Point3D p0, Point3D p1, Point3D p2)
        {
            var v0 = new Vector3D(p1.X - p0.X, p1.Y - p0.Y, p1.Z - p0.Z);
            var v1 = new Vector3D(p2.X - p1.X, p2.Y - p1.Y, p2.Z - p1.Z);
            return Vector3D.CrossProduct(v0, v1);
        }

        private static MeshGeometry3D BuildMesh(Point3D p0, Point3D p1, Point3D p2, Point3D p3,
                                                Point q0, Point q1, Point q2, Point q3)
        {
            var mesh = new MeshGeometry3D();
            mesh.Positions.Add(p0);
            mesh.Positions.Add(p1);
            mesh.Positions.Add(p2);
            mesh.Positions.Add(p3);

            var normal = CalculateNormal(p0, p1, p2);
            mesh.TriangleIndices.Add(0);
            mesh.TriangleIndices.Add(1);
            mesh.TriangleIndices.Add(2);
            mesh.Normals.Add(normal);
            mesh.Normals.Add(normal);
            mesh.Normals.Add(normal);
            mesh.TextureCoordinates.Add(q3);
            mesh.TextureCoordinates.Add(q2);
            mesh.TextureCoordinates.Add(q1);

            normal = CalculateNormal(p2, p3, p0);
            mesh.TriangleIndices.Add(2);
            mesh.TriangleIndices.Add(3);
            mesh.TriangleIndices.Add(0);
            mesh.Normals.Add(normal);
            mesh.Normals.Add(normal);
            mesh.Normals.Add(normal);
            mesh.TextureCoordinates.Add(q0);
            mesh.TextureCoordinates.Add(q1);
            mesh.TextureCoordinates.Add(q2);

            mesh.Freeze();
            return mesh;
        }
    }
}

By viewing downloads associated with this article you agree to the Terms of Service and the article's licence.

If a file you wish to view isn't highlighted, and is a text file (not binary), please let us know and we'll add colourisation support for it.

License

This article, along with any associated source code and files, is licensed under The Code Project Open License (CPOL)

Share

About the Author

Bernhard Kügle
Software Developer (Senior) Siemens Healthcare
Germany Germany
Currently working as a Requirement Engineer for Siemens Healthcare in Germany.

| Advertise | Privacy | Mobile
Web01 | 2.8.140821.2 | Last Updated 27 Sep 2012
Article Copyright 2012 by Bernhard Kügle
Everything else Copyright © CodeProject, 1999-2014
Terms of Service
Layout: fixed | fluid