Details:
EPID information:
Image Dimension: 611 576 1
Image Spacing: 0.25mm 0.25mm 1.00mm
contourDRR Information:
Image Dimension: 270 256 1
Image Spacing: 0.5859mm 0.5859mm 1.0000mm
With these specifactions, at firstly I made the same type these volumes (EPID and contourDRR) to Double. Then I applied geodesicActiveContourLevelSetImageFilter. I get the following error message:
Exception thrown in SimpleITK GeodesicActiveContourLevelSetImageFilter_Execute: /work/Preview/Slicer-0-build/SimpleITK/Code/BasicFilters/src/sitkImageFilter.cxx:63:
sitk::ERROR: Input “featureImage” for “GeodesicActiveContourLevelSetImageFilter” has size of [ 270, 256, 1 ] which does not match the primary input’s size of [ 611, 576, 1 ]!
So I use the module of General Registration (ANTs) to unify the specifications of these volumes (EPID and contourDRR).
Inputs:
Fixed Image: EPID
Moving Image: contourDRR
Stages (Presets): Rigid
Setting:
Computation Precision: double
Outputs:
Transformed Volume: ContourDRRTRansformed
I get the following result:
EPID information:
Image Dimension: 611 576 1
Image Spacing: 0.25mm 0.25mm 1.00mm
ContourDRRTRansformed Information:
Image Dimension: 611 576 1
Image Spacing: 0.25mm 0.25mm 1.0000mm
At now, I run geodesicActiveContourLevelSetImageFilter on these volumes: EPID and ContourDRRTRansformed nodes.I get the following error message:
Exception thrown in SimpleITK GeodesicActiveContourLevelSetImageFilter_Execute: /work/Preview/Slicer-0-build/ITK/Modules/Filtering/ImageFilterBase/include/itkRecursiveSeparableImageFilter.hxx:227:
ITK ERROR: RecursiveGaussianImageFilter(0x7f8d78567a40): The number of pixels along direction 2 is less than 4. This filter requires a minimum of four pixels along the dimension to be processed.
I do not know to add three slices (exactly one or copy) in the direction of IS in 3DSlicer. For doing it, I save these volumes ( EPID and ContourDRRTRansformed ) in nifti format and use medcon for example with following commands:
$ cp ContourDRRTRansformed.nii ContourDRRTRansformed1.nii
$ cp ContourDRRTRansformed.nii ContourDRRTRansformed2.nii
$ cp ContourDRRTRansformed.nii ContourDRRTRansformed3.nii
$ cp ContourDRRTRansformed.nii ContourDRRTRansformed3.nii
$ medcon -f ContourDRRTRansformed1.nii ContourDRRTRansformed2.nii ContourDRRTRansformed3.nii ContourDRRTRansformed4.nii -stack3d -c nii
...
$ cp EPID.nii EPID1.nii
$ cp EPID.nii EPID2.nii
$ cp EPID.nii EPID3.nii
$ cp EPID.nii EPID3.nii
$ medcon -f EPID1.nii EPID2.nii EPID3.nii EPID4.nii -stack3d -c nifti
Result:
EPID information:
Image Dimension: 611 576 4
Image Spacing: 0.25mm 0.25mm 1.00mm
contourDRR Information:
Image Dimension: 611 576 4
Image Spacing: 0.25mm 0.25mm 1.00mm
After doing it, I again I run geodesicActiveContourLevelSetImageFilter on these volumes: EPID and ContourDRRTRansformed. As mentioned in my previous message, unfortunately, after using this filter, I took an output image that had all the pixels equal to one.
Please guide me to solving this problem.
Thanks.
Shahrokh