Hi Mathias, I managed to solve my problems by looking into the SphereLandmarkProjector class. The key was to initialize the plane geometry appropriately. In particular, it seems that the parameter plane needs to be defined by the y and z coordinates instead of x and y which I was assuming...
The 3D visualization problem does not hurt me right now since I am switching the 4th window to a 2D view when showing the CPR. For a reproduction of the 3D problem you should try your example code with the SphereLandmarkProjector instead of the PlaneLandmarkProjector. Since I am using a pretty old version of MITK it could well be that this bug has already been fixed in the current version, however. Thanks for your support again! Best, Michael. > -----Original Message----- > From: Seitel Mathias [mailto:[email protected]] > Sent: Monday, November 23, 2009 11:34 AM > To: [email protected] > Subject: Re: [mitk-users] Curved MPR > > Hi Michael, > > Apologies for the delayed reply. > > > > Unfortunately, the constructed vtkThinSplineTransform does not > > describe the curved plane that I would like to obtain. I > would like to > > define a curved intersection plane that is planar in one > direction and > > curved only in the other. > > I'm not sure if I understand you correctly: do you mean that > the resulting plane should resemble a height-field, i.e. when > looking at it exactly from above it would still appear as a > perfect planar grid? > > > > I therefore tried to use mitk::ExternAbstractTransformGeometry to > > define my own transform. I tried to understand the MITK > code in order > > to mimick mitk::ThinPlateSplineCurvedGeometry. > > What type of VTK transform are you using together with > mitk::ExternAbstractTransformGeometry? I'm afraid that the > class is currently not in active use within MITK, so I cannot > promise that it still works as expected or correctly as it is. > > The following publication describes in detail curved MPR in > MITK as it was originally implemented: > > Wolf I., Hastenteufel M., Wegner I., et al. Curved > reformations using the Medical Imaging Interaction Toolkit > (MITK). Proc SPIE Medical Imaging 2005, Volume 5744, pp831-838. > > > > To my understanding, the mapper request an intensity value at a 2D > > coordinate which is mapped by a plane geometry to a 3D coordinate > > which is mapped by the VTK transform to another 3D coordinate from > > which the image intensity is obtained. Still it is not clear to me, > > what the ranges of the respective coordinates should be and > how to set > > up the plane in order to get a correctly rendered result. > > A good place to start exploring this would be to use > mitk::ThinPlateSplineCurvedGeometry in order to automatically > initialize the plane and the internal VTK transform using the > landmark based approach. You could then determine their > parameters (via GetPlane() and > GetThinPlateSplineTransform()). To my understanding, the > plane extent should be large enough to cover the extent of > the 3D volume you are going to intersect as it is projected > onto the plane. > > > > Right now, I can only see parts of the curved plane and if I rotate > > the 3D view, things vanish and reappear. If, in the data manager, I > > click the "global reinit" button the rendering is right > again but the > > additional plane in the 3D view is straight again and at > the bottom of > > the image. Apparently, I am doing something wrong. > > This really looks like a bug in the rendering engine. Could > you post the code you used for this, so that we can try to > reproduce it? > > Thanks! > > > Mathias > > > > -------------------------------------------------------------- > ---------------- > Let Crystal Reports handle the reporting - Free Crystal > Reports 2008 30-Day > trial. Simplify your report design, integration and > deployment - and focus on > what you do best, core application coding. Discover what's new with > Crystal Reports now. http://p.sf.net/sfu/bobj-july > _______________________________________________ > mitk-users mailing list > [email protected] > https://lists.sourceforge.net/lists/listinfo/mitk-users > ------------------------------------------------------------------------------ Let Crystal Reports handle the reporting - Free Crystal Reports 2008 30-Day trial. Simplify your report design, integration and deployment - and focus on what you do best, core application coding. Discover what's new with Crystal Reports now. http://p.sf.net/sfu/bobj-july _______________________________________________ mitk-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/mitk-users
