Dear Ahmed,
really thanks a lot for getting back to the list and posting the solution,
this is highly appreciated.
I'm happy you found the solution and can proceed with your OsiriX integration.
Actually, it would be really interesting to see it running, feel free to keep
me posted.
Best,
Luca
On Feb 23, 2013, at 6:43 AM, Ahmed Harouni wrote:
> Dear Luca
>
> I managed to solve my problem. I don't want to wast your time responding. I
> converted the point co-ordinates from RAS to IJK. Also I wasn't passing the
> point correctly. Now it is working great.
>
> Thanks
>
> int VMTKseed1[3],VMTKseed2[3]; // fill it with points in IJK coordinates
>
> vtkIdList *my_VTKseedIdList1=vtkIdList::New();
> vtkIdList *my_VTKseedIdList2=vtkIdList::New();
> vtkIdType pttyp1,pttyp2;
> pttyp1=myVTKspeedImage->ComputePointId(VMTKseed1);
> pttyp2=myVTKspeedImage->ComputePointId(VMTKseed2);
> /**************/
> my_VTKseedIdList1->InsertNextId(pttyp1);
> my_VTKseedIdList2->InsertNextId(pttyp2);
> vtkvmtkCollidingFrontsImageFilter *myVMTKcollidingFronts =
> vtkvmtkCollidingFrontsImageFilter::New();
> myVMTKcollidingFronts->SetInput(myVTKspeedImage);
> myVMTKcollidingFronts->SetSeeds1(my_VTKseedIdList1);
> myVMTKcollidingFronts->SetSeeds2(my_VTKseedIdList2);
> myVMTKcollidingFronts->ApplyConnectivityOn();
> myVMTKcollidingFronts->StopOnTargetsOn();
> myVMTKcollidingFronts->Update();
>
> Ahmed A. Harouni
>
>
> On Thu, Feb 21, 2013 at 12:30 PM, Ahmed Harouni <ahmed.haro...@gmail.com>
> wrote:
> Dear Luca
>
> Thanks for your previous help. I have managed to use
> vtkvmtkVesselnessMeasureImageFilter class and am displaying the results
> successfully to Osirix output. Now I trying to use
> vtkvmtkCollidingFrontsImageFilter. I am following the pythons scripts and
> getting the expected results.
> I can't determine ow to pass the two seed points to the
> vtkvmtkCollidingFrontsImageFilter class using the SetSeeds1. I think I am
> using the wrong co-ordinate system. I capture the ROI from Osirix which is
> given to me in RAS system in mm. My understanding that VMTK uses IJK system
> but it already has the origin and spacing values in the vtkImageData.
> Should I convert my point coordinate manually using matrix transformation or
> am I missing something? (Bellow is my code for
> vtkvmtkCollidingFrontsImageFilter)
>
> Thanks for your help.
>
>
> vtkIdList *my_VTKseedIdList1=vtkIdList::New();
> vtkIdList *my_VTKseedIdList2=vtkIdList::New();
> // vtkSmartPointer<vtkPoints> point11 =vtkSmartPointer<vtkPoints>::New();
> vtkPoints *point11 =vtkPoints::New();
> vtkPoints *points22=vtkPoints::New();
> vtkIdType pttyp1,pttyp2;
>
>
> pttyp1=point11->InsertNextPoint(5,5,5);//(15.3,110.8,946.2);//{0,90,922};//{283,242,251};//
> {18.5,-148.0,115};
> pttyp2=points22->InsertNextPoint(15,15,15);//(-10.8,95.8,931);//
> {80,146,951};//{ 338,276, 218};//{44.0, -155.0, 110.5};
> my_VTKseedIdList1->InsertNextId(pttyp1);
> my_VTKseedIdList2->InsertNextId(pttyp2);
>
> vtkvmtkCollidingFrontsImageFilter *myVMTKcollidingFronts =
> vtkvmtkCollidingFrontsImageFilter::New();
> myVMTKcollidingFronts->SetInput(myVTKspeedImage); //myVTKspeedImage
> is: vtkImageData * myVTKspeedImage and I call it's update function and I
> see the expected results when displaying it.
> myVMTKcollidingFronts->SetSeeds1(my_VTKseedIdList1);
> myVMTKcollidingFronts->SetSeeds2(my_VTKseedIdList2);
> myVMTKcollidingFronts->ApplyConnectivityOn();
> myVMTKcollidingFronts->StopOnTargetsOn();
> myVMTKcollidingFronts->Update(); // I get an empty image !
>
>
> Ahmed A. Harouni
>
> ------------------------------------------------------------------------------
> Everyone hates slow websites. So do we.
> Make your web apps faster with AppDynamics
> Download AppDynamics Lite for free today:
> http://p.sf.net/sfu/appdyn_d2d_feb_______________________________________________
> vmtk-users mailing list
> vmtk-users@lists.sourceforge.net
> https://lists.sourceforge.net/lists/listinfo/vmtk-users
------------------------------------------------------------------------------
Everyone hates slow websites. So do we.
Make your web apps faster with AppDynamics
Download AppDynamics Lite for free today:
http://p.sf.net/sfu/appdyn_d2d_feb
_______________________________________________
vmtk-users mailing list
vmtk-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/vmtk-users
