[vtkusers] Volume Rendering of 3D CT Scan similar to 2D X Ray image

Arnaud Gelas arnaudgelas at gmail.com
Tue Jan 28 08:12:54 EST 2014


Hi Bill,

I have been playing with existing ITK examples, but results are not good enough for my case…

Would you mind sharing the code, please?

Thanks!!

Cheers,
Arnaud

On Jan 24, 2014, at 4:33 PM, Bill Lorensen <bill.lorensen at gmail.com> wrote:

> I have some code we wrote years ago. It used
> typedef itk::ExpNegativeImageFilter<InputImageType,InternalImageType>
>      ExpNegativeType;
>    ExpNegativeType::Pointer expNegative = ExpNegativeType::New();
>    expNegative->SetInput( thresh->GetOutput());
>    expNegative->SetFactor( -.0011 );
> 
> followed by
>   // grab a row for the 'DRR'
>    typedef itk::AccumulateImageFilter<InternalImageType,InternalImageType>
> AccumulateType;
>    AccumulateType::Pointer APAccumulate = AccumulateType::New();
>    AccumulateType::Pointer lateralAccumulate = AccumulateType::New();
> 
>    APAccumulate->SetAverage(false);
>    lateralAccumulate->SetAverage(false);
>    APAccumulate->SetAccumulateDimension(1);
>    lateralAccumulate->SetAccumulateDimension(0);
> 
>    APAccumulate->SetInput ( expNegative->GetOutput() );
>    lateralAccumulate->SetInput ( expNegative->GetOutput() );
>    APAccumulate->Update();
>    lateralAccumulate->Update();
> 
> I'll be happy to share the code with the caveat that it has lots of
> application specific code in it.
> 
> There is also the ITK Example:
> Examples/Filtering/DigitallyReconstructedRadiograph1.cxx
> 
> Bill
> 
> On Fri, Jan 24, 2014 at 10:23 AM, Arnaud Gelas <arnaudgelas at gmail.com> wrote:
>> Hi all,
>> 
>> I would like to render as realistic as possible 1 3D CT Scan, as if it was 2D X Ray Image.
>> 
>> I have tried Average Intensity Projection, but results are much smoother than what it should look like.
>> 
>> Has anyone any experience on the question? Any advice appreciated?
>> 
>> Do I need to generate a specific transfer function for this CT? Could I apply the same transfer function to another one?
>> Has anyone got such a transfer function that could be shared?
>> 
>> Thanks in advance,
>> 
>> Best regards
>> Arnaud Gelas
>> _______________________________________________
>> Powered by www.kitware.com
>> 
>> Visit other Kitware open-source projects at http://www.kitware.com/opensource/opensource.html
>> 
>> Please keep messages on-topic and check the VTK FAQ at: http://www.vtk.org/Wiki/VTK_FAQ
>> 
>> Follow this link to subscribe/unsubscribe:
>> http://www.vtk.org/mailman/listinfo/vtkusers
> 
> 
> 
> -- 
> Unpaid intern in BillsBasement at noware dot com



More information about the vtkusers mailing list