Hi Luis,<br><br>I was getting a segmentation fault because I had forgotten to change the Dimension from 2 to 3 in DicomReadPrintTag.cxx.<br><br>For the input image, I have:<br>(0020|0037) Image Orientation (Patient) = 0.0\0.0\-1.0\0.0\1.0\0.0 <br>
<br>For the output image, I have:<br>(0020|0037) Image Orientation (Patient) = 0.000000\0.000000\1.000000\0.000000\1.000000\0.000000<br><br>There is a difference. Since the program is changing the orientation when processed through the pipeline, is the best approach to manually change the orientation tag of the dictionary and then copy the revised dictionary to the output file?<br>
<br>john<br><br>Below are the full tags:<br><br>[jdrozd@trumpet DicomImageReadPrintTags]$ ./DicomImageReadPrintTags correctedsubject5.dcm <br>(0002|0000) Group Length = 214 <br>
(0002|0001) File Meta Information Version = AAE= <br>(0002|0002) Media Storage SOP Class UID = 1.2.840.10008.5.1.4.1.1.2 <br>(0002|0003) Media Storage SOP Instance UID = 1.2.826.0.1.3680043.2.1125.1.10607669833050788267094246636093811<br>
(0002|0010) Transfer Syntax UID = 1.2.840.10008.1.2.1 <br>(0002|0012) Implementation Class UID = 147.144.143.155 <br>
(0002|0013) Implementation Version Name = ITK/GDCM 1.2.4 <br>(0002|0016) Source Application Entity Title = NOTSPECIFIED <br>
(0008|0008) Image Type = DERIVED\PRIMARY <br>(0008|0012) Instance Creation Date = 20091029 <br>
(0008|0013) Instance Creation Time = 135224 <br>(0008|0016) SOP Class UID = 1.2.840.10008.5.1.4.1.1.2 <br>
(0008|0018) SOP Instance UID = 1.2.826.0.1.3680043.2.1125.1.10607669833050788267094246636093811 <br>(0008|0020) Study Date = 20081030 <br>
(0008|0030) Study Time = 164348.940 <br>(0008|0050) Accession Number = <br>
(0008|0060) Modality = CT <br>(0008|0064) Conversion Type = WSD <br>
(0008|0070) Manufacturer = Manifacturer <br>(0008|0080) Institution Name = GDCM Hospital <br>
(0008|0090) Referring Physician's Name = Refering Phisician <br>(0010|0010) Patient's Name = Patient 188858520 <br>
(0010|0020) Patient ID = 1747233212 <br>(0010|0030) Patient's Birth Date = 19500101 <br>
(0010|0040) Patient's Sex = M <br>(0018|0088) Spacing Between Slices = 1.207500 <br>
(0018|1164) Imager Pixel Spacing = 0.945750\0.945750 <br>(0020|000d) Study Instance UID = 1.2.826.0.1.3680043.2.1125.1.53653479342656887425609263146205769 <br>
(0020|000e) Series Instance UID = 1.2.826.0.1.3680043.2.1125.1.47155241092896696198844159625679986 <br>(0020|0010) Study ID = 1533117581 <br>
(0020|0011) Series Number = 2135500125 <br>(0020|0013) Instance Number = 165 <br>
(0020|0020) Patient Orientation = L\P <br>(0020|0032) Image Position (Patient) = 199.237496852874\0.0\0.0 <br>
(0020|0037) Image Orientation (Patient) = 0.0\0.0\-1.0\0.0\1.0\0.0 <br>(0020|4000) Image Comments = NOT FOR CLINICAL USE <br>
(0028|0002) Samples per Pixel = 1 <br>(0028|0004) Photometric Interpretation = MONOCHROME2 <br>
(0028|0008) Number of Frames = 166 <br>(0028|0010) Rows = 256 <br>
(0028|0011) Columns = 256 <br>(0028|0030) Pixel Spacing = 0.945750\0.945750 <br>
(0028|0034) Pixel Aspect Ratio = 1\1 <br>(0028|0100) Bits Allocated = 16 <br>
(0028|0101) Bits Stored = 16 <br>(0028|0102) High Bit = 15 <br>
(0028|0103) Pixel Representation = 0 <br>(0028|1052) Rescale Intercept = 0.0 <br>
(0028|1053) Rescale Slope = 1.0 <br>(0028|1054) Rescale Type = US <br>
(7fe0|0000) Group Length = 21757964 <br>Patient's Name (0010|0010) is: Patient 188858520 <br>
Performing Physician's Name (0008|1050): (No Value Found in File) <br>PixelType: scalar <br>
Component Type: unsigned_short <br> <br>[jdrozd@trumpet DicomImageReadPrintTags]$ ./DicomImageReadPrintTags outsubject5.dcm <br>
(0002|0000) Group Length = 194 <br>(0002|0001) File Meta Information Version = AAE= <br>(0002|0002) Media Storage SOP Class UID = 1.2.840.10008.5.1.4.1.1.7 <br>
(0002|0003) Media Storage SOP Instance UID = 1.2.826.0.1.3680043.2.1143.8995209879494441053058872310646575724<br>(0002|0010) Transfer Syntax UID = 1.2.840.10008.1.2.1 <br>
(0002|0012) Implementation Class UID = 147.144.143.155 <br>(0002|0013) Implementation Version Name = ITK/GDCM 1.2.4 <br>
(0008|0012) Instance Creation Date = 20091126 <br>(0008|0013) Instance Creation Time = 150853 <br>
(0008|0016) SOP Class UID = 1.2.840.10008.5.1.4.1.1.7 <br>(0008|0018) SOP Instance UID = 1.2.826.0.1.3680043.2.1143.8995209879494441053058872310646575724 <br>
(0008|0020) Study Date = 20091126 <br>(0008|0030) Study Time = 150853 <br>
(0008|0050) Accession Number = <br>(0008|0060) Modality = OT <br>
(0008|0064) Conversion Type = SYN <br>(0008|0070) Manufacturer = GDCM Factory<br>(0008|0080) Institution Name = GDCM Hospital<br>(0008|0090) Referring Physician's Name =<br>
(0010|0010) Patient's Name = GDCM^Patient<br>(0010|0020) Patient ID = GDCM ID<br>(0010|0030) Patient's Birth Date =<br>(0010|0040) Patient's Sex =<br>(0018|0088) Spacing Between Slices = 1.207500<br>(0018|1164) Imager Pixel Spacing = 0.945750\0.945750<br>
(0020|000d) Study Instance UID = 1.2.826.0.1.3680043.2.1125.1.99663966693826559960117773595375861<br>(0020|000e) Series Instance UID = 1.2.826.0.1.3680043.2.1125.1.80656207154489419332286356071443484<br>(0020|0010) Study ID =<br>
(0020|0011) Series Number =<br>(0020|0013) Instance Number =<br>(0020|0020) Patient Orientation = H\P<br>(0020|0032) Image Position (Patient) = 199.237503\0.000000\0.000000<br>(0020|0037) Image Orientation (Patient) = 0.000000\0.000000\1.000000\0.000000\1.000000\0.000000<br>
(0028|0002) Samples per Pixel = 1<br>(0028|0004) Photometric Interpretation = MONOCHROME2<br>(0028|0008) Number of Frames = 166<br>(0028|0010) Rows = 256<br>(0028|0011) Columns = 256<br>(0028|0030) Pixel Spacing = 0.945750\0.945750<br>
(0028|0100) Bits Allocated = 16<br>(0028|0101) Bits Stored = 16<br>(0028|0102) High Bit = 15<br>(0028|0103) Pixel Representation = 1<br>(7fe0|0000) Group Length = 21757964<br>Patient's Name (0010|0010) is: GDCM^Patient<br>
Performing Physician's Name (0008|1050): (No Value Found in File)<br>PixelType: scalar<br>Component Type: short<br><br><div class="gmail_quote">On Thu, Nov 26, 2009 at 2:44 PM, John Drozd <span dir="ltr"><<a href="mailto:john.drozd@gmail.com">john.drozd@gmail.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi Luis:<br><br>Sorry for taking some time to get back to you. I took a break from my computer to refresh my mind so I did an hour of ice skating at our university hockey arena. I only fell a few times.<br>
<br>Here is the output that you requested:<br>
<br>From the added std::cout statements to my code, I have:<br><br>output from reader->GetOutput()->GetDirection()<br>0 0 1<br>0 1 0<br>-1 0 0<br><br>output from castFilter->GetOutput()->GetDirection()<br>0 0 1<br>
0 1 0<br>-1 0 0<br><br>output from connectedThreshold->GetOutput()->GetDirection()<br>0 0 1<br>0 1 0<br>-1 0 0<br><br>output from caster->GetOutput()->GetDirection()<br>0 0 1<br>0 1 0<br>-1 0 0<br><br>Also here are the tag statements from running that you requested:<br>
<i><b>(I get a segmentation fault when I run this program on the output image<br>and I'm trying to figure out why it seg faults on the output image)</b></i><br><br> <br>[jdrozd@trumpet DicomImageReadPrintTags]$ ./DicomImageReadPrintTags correctedsubject5.dcm<br>
(0002|0000) Group Length = 214 <br>(0002|0001) File Meta Information Version = AAE= <br>(0002|0002) Media Storage SOP Class UID = 1.2.840.10008.5.1.4.1.1.2 <br>
(0002|0003) Media Storage SOP Instance UID = 1.2.826.0.1.3680043.2.1125.1.10607669833050788267094246636093811<br>(0002|0010) Transfer Syntax UID = 1.2.840.10008.1.2.1 <br>
(0002|0012) Implementation Class UID = 147.144.143.155 <br>(0002|0013) Implementation Version Name = ITK/GDCM 1.2.4 <br>
(0002|0016) Source Application Entity Title = NOTSPECIFIED <br>(0008|0008) Image Type = DERIVED\PRIMARY <br>
(0008|0012) Instance Creation Date = 20091029 <br>(0008|0013) Instance Creation Time = 135224 <br>
(0008|0016) SOP Class UID = 1.2.840.10008.5.1.4.1.1.2 <br>(0008|0018) SOP Instance UID = 1.2.826.0.1.3680043.2.1125.1.10607669833050788267094246636093811 <br>
(0008|0020) Study Date = 20081030 <br>(0008|0030) Study Time = 164348.940 <br>
(0008|0050) Accession Number = <br>(0008|0060) Modality = CT <br>
(0008|0064) Conversion Type = WSD <br>(0008|0070) Manufacturer = Manifacturer <br>
(0008|0080) Institution Name = GDCM Hospital <br>(0008|0090) Referring Physician's Name = Refering Phisician <br>
(0010|0010) Patient's Name = Patient 188858520 <br>(0010|0020) Patient ID = 1747233212 <br>
(0010|0030) Patient's Birth Date = 19500101 <br>(0010|0040) Patient's Sex = M <br>
(0018|0088) Spacing Between Slices = 1.207500<br>(0018|1164) Imager Pixel Spacing = 0.945750\0.945750<br>(0020|000d) Study Instance UID = 1.2.826.0.1.3680043.2.1125.1.53653479342656887425609263146205769<br>(0020|000e) Series Instance UID = 1.2.826.0.1.3680043.2.1125.1.47155241092896696198844159625679986<br>
(0020|0010) Study ID = 1533117581<br>(0020|0011) Series Number = 2135500125<br>(0020|0013) Instance Number = 165<br>(0020|0020) Patient Orientation = L\P<br>(0020|0032) Image Position (Patient) = 199.237496852874\0.0\0.0<br>
(0020|0037) Image Orientation (Patient) = 0.0\0.0\-1.0\0.0\1.0\0.0<br>(0020|4000) Image Comments = NOT FOR CLINICAL USE<br>(0028|0002) Samples per Pixel = 1<br>(0028|0004) Photometric Interpretation = MONOCHROME2<br>(0028|0008) Number of Frames = 166<br>
(0028|0010) Rows = 256<br>(0028|0011) Columns = 256<br>(0028|0030) Pixel Spacing = 0.945750\0.945750<br>(0028|0034) Pixel Aspect Ratio = 1\1<br>(0028|0100) Bits Allocated = 16<br>(0028|0101) Bits Stored = 16<br>(0028|0102) High Bit = 15<br>
(0028|0103) Pixel Representation = 0<br>(0028|1052) Rescale Intercept = 0.0<br>(0028|1053) Rescale Slope = 1.0<br>(0028|1054) Rescale Type = US<br>(7fe0|0000) Group Length = 21757964<br>Patient's Name (0010|0010) is: Patient 188858520<br>
Performing Physician's Name (0008|1050): (No Value Found in File)<br>PixelType: scalar<br>Component Type: unsigned_short<br><br>[jdrozd@trumpet DicomImageReadPrintTags]$ ./DicomImageReadPrintTags outsubject5.dcm<br>Segmentation fault<br>
[jdrozd@trumpet DicomImageReadPrintTags]$<br><br>Thanks,<br><font color="#888888">john</font><div><div></div><div class="h5"><br><br><br><div class="gmail_quote">On Thu, Nov 26, 2009 at 1:27 PM, Luis Ibanez <span dir="ltr"><<a href="mailto:luis.ibanez@kitware.com" target="_blank">luis.ibanez@kitware.com</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">Hi John,<br>
<br>
Please use the program:<br>
<br>
Insight/Examples/IO/DicomImageReadPrintTags.cxx<br>
<br>
To print the DICOM tags of both the input and output images.<br>
<br>
Look for the string:<br>
<br>
(0020|0037) Image Orientation (Patient)<br>
<br>
and post the values (for both images) to the mailing list.<br>
<br>
They should have a look similar to:<br>
<br>
(0020|0037) Image Orientation (Patient) =<br>
0.9983222\0.02102422\0.05395257\-0.02105488\0.9997783\1.490116e-08<br>
<br>
----<br>
<br>
<br>
Also add the following print out statements to your code:<br>
<br>
std::cout << reader->GetOutput()->GetDirection() << std::endl;<br>
std::cout << casterFilter->GetOutput()->GetDirection() << std::endl;<br>
std::cout << connectedThreshold->GetOutput()->GetDirection() << std::endl;<br>
std::cout << caster->GetOutput()->GetDirection() << std::endl;<br>
<br>
after the call to<br>
<br>
writer->Update();<br>
<br>
<br>
In this way we could sort out at what level of the process<br>
the Direction information is being lost.<br>
<br>
<br>
Thanks,<br>
<br>
<br>
Luis<br>
<br>
<br>
----------------------------------------------------------------------------<br>
<div><div></div><div>On Thu, Nov 26, 2009 at 12:57 PM, John Drozd <<a href="mailto:john.drozd@gmail.com" target="_blank">john.drozd@gmail.com</a>> wrote:<br>
> Hello,<br>
><br>
> Can anyone please tell me why my connect threshold segmentation code, where<br>
> I input a 3d dicom brain volume and output a 3d dicom ventricle<br>
> segmentation, gives me my segmentation "upside down" (relative to the<br>
> original inputted image) when I view it in 3D Slicer? When I output the<br>
> code in any other format eg. .mhd, the segmentation is right side up.<br>
><br>
> Thank you.<br>
><br>
> john<br>
><br>
> Below is my uncommented code:<br>
><br>
> /*<br>
> to run type:<br>
> ./ConnectedThresholdImageFilter correctedsubject5.dcm outsubject5.dcm 103<br>
> 142 95 17100 17300<br>
> */<br>
><br>
> #if defined(_MSC_VER)<br>
> #pragma warning ( disable : 4786 )<br>
> #endif<br>
><br>
> #ifdef __BORLANDC__<br>
> #define ITK_LEAN_AND_MEAN<br>
> #endif<br>
><br>
><br>
> #include "itkConnectedThresholdImageFilter.h"<br>
><br>
> #include "itkImage.h"<br>
> #include "itkCastImageFilter.h"<br>
><br>
> #include "itkCurvatureFlowImageFilter.h"<br>
><br>
> #include "itkImageFileReader.h"<br>
> #include "itkImageFileWriter.h"<br>
><br>
> #include "itkGDCMImageIO.h"<br>
><br>
> #include "itkVersion.h"<br>
><br>
> #include "itkOrientedImage.h"<br>
> #include "itkMinimumMaximumImageFilter.h"<br>
><br>
> #include "itkGDCMImageIO.h"<br>
> #include "itkGDCMSeriesFileNames.h"<br>
> #include "itkNumericSeriesFileNames.h"<br>
><br>
> #include "itkImageSeriesReader.h"<br>
> #include "itkImageSeriesWriter.h"<br>
><br>
> #include "itkResampleImageFilter.h"<br>
> #include "itkShiftScaleImageFilter.h"<br>
><br>
> #include "itkIdentityTransform.h"<br>
> #include "itkLinearInterpolateImageFunction.h"<br>
><br>
> #include <itksys/SystemTools.hxx><br>
><br>
> #include "gdcm/src/gdcmFile.h"<br>
> #include "gdcm/src/gdcmUtil.h"<br>
><br>
> #include <string><br>
><br>
> int main( int argc, char *argv[])<br>
> {<br>
> if( argc < 7 )<br>
> {<br>
> std::cerr << "Missing Parameters " << std::endl;<br>
> std::cerr << "Usage: " << argv[0];<br>
> std::cerr << " inputImage outputImage seedX seedY seedZ lowerThreshold<br>
> upperThreshold" << std::endl;<br>
><br>
> return 1;<br>
> }<br>
><br>
> typedef float InternalPixelType;<br>
><br>
> const unsigned int Dimension = 3;<br>
><br>
> typedef itk::Image< InternalPixelType, Dimension > InternalImageType;<br>
><br>
> typedef signed short OutputPixelType;<br>
><br>
> typedef itk::Image< OutputPixelType, Dimension > OutputImageType;<br>
> typedef itk::Image< float, Dimension > OutputImageType2;<br>
> typedef itk::CastImageFilter< InternalImageType, OutputImageType ><br>
> CastingFilterType;<br>
> CastingFilterType::Pointer caster = CastingFilterType::New();<br>
><br>
> const unsigned int ImageDimension = 3;<br>
> typedef signed short PixelType;<br>
><br>
> typedef itk::Image< PixelType, ImageDimension > FixedImageType;<br>
> typedef itk::Image< float, ImageDimension > FloatImageType;<br>
><br>
> typedef itk::ImageFileReader< FixedImageType > ReaderType;<br>
> typedef itk::ImageFileWriter< OutputImageType > WriterType;<br>
> typedef itk::ImageFileWriter< FloatImageType > WriterType2;<br>
><br>
> ReaderType::Pointer reader = ReaderType::New();<br>
> WriterType::Pointer writer = WriterType::New();<br>
> WriterType2::Pointer writer2 = WriterType2::New();<br>
><br>
> typedef itk::GDCMImageIO ImageIOTypefixed;<br>
> ImageIOTypefixed::Pointer gdcmImageIOfixed = ImageIOTypefixed::New();<br>
> reader->SetImageIO( gdcmImageIOfixed );<br>
><br>
> typedef itk::GDCMImageIO ImageIOTypefixed2;<br>
> ImageIOTypefixed2::Pointer gdcmImageIOfixed2 = ImageIOTypefixed2::New();<br>
><br>
> reader->SetFileName( argv[1] );<br>
><br>
> reader->Update();<br>
><br>
> typedef itk::CurvatureFlowImageFilter< InternalImageType,<br>
> InternalImageType ><br>
> CurvatureFlowImageFilterType;<br>
><br>
> CurvatureFlowImageFilterType::Pointer smoothing =<br>
> CurvatureFlowImageFilterType::New();<br>
><br>
> typedef itk::ConnectedThresholdImageFilter< InternalImageType,<br>
> InternalImageType > ConnectedFilterType;<br>
><br>
> ConnectedFilterType::Pointer connectedThreshold =<br>
> ConnectedFilterType::New();<br>
><br>
> typedef signed short InputAPixelType;<br>
> typedef float OutputBPixelType;<br>
><br>
> typedef itk::Image< InputAPixelType, 3 > InputAImageType;<br>
> typedef itk::Image< OutputBPixelType, 3 > OutputBImageType;<br>
><br>
> typedef itk::CastImageFilter< InputAImageType, OutputBImageType ><br>
> CastFilterType;<br>
><br>
> CastFilterType::Pointer castFilter = CastFilterType::New();<br>
><br>
><br>
> castFilter->SetInput( reader->GetOutput() );<br>
><br>
><br>
> connectedThreshold->SetInput( castFilter->GetOutput() );<br>
><br>
> caster->SetInput( connectedThreshold->GetOutput() );<br>
><br>
><br>
> smoothing->SetNumberOfIterations( 20 ); //was 5<br>
> smoothing->SetTimeStep( 0.125 );<br>
><br>
> const InternalPixelType lowerThreshold = atof( argv[6] );<br>
> const InternalPixelType upperThreshold = atof( argv[7] );<br>
><br>
> connectedThreshold->SetLower( lowerThreshold );<br>
> connectedThreshold->SetUpper( upperThreshold );<br>
><br>
> connectedThreshold->SetReplaceValue( 255 );<br>
><br>
> InternalImageType::IndexType index;<br>
><br>
> index[0] = atoi( argv[3] );<br>
> index[1] = atoi( argv[4] );<br>
><br>
> //added<br>
> index[2] = atoi( argv[5] );<br>
><br>
> std::cout << index << std::endl;<br>
><br>
> // Software Guide : BeginCodeSnippet<br>
> connectedThreshold->SetSeed( index );<br>
><br>
> //obtain a 5 x 5 bounding region of seeds<br>
> int ii, jj, kk;<br>
><br>
> ii = index[0];<br>
> jj = index[1];<br>
> kk = index[2];<br>
><br>
> for (int i = ii; i < ii + 5; i++)<br>
> for (int j = jj; j < jj + 5; j++)<br>
> for (int k = kk; k < kk + 5; k++)<br>
> {<br>
><br>
> index[0] = i;<br>
> index[1] = j;<br>
> index[2] = k;<br>
> connectedThreshold->AddSeed( index );<br>
> }<br>
><br>
> for (int i = ii; i > ii - 5; i--)<br>
> for (int j = jj; j > jj - 5; j--)<br>
> for (int k = kk; k > kk - 5; k--)<br>
> {<br>
><br>
> index[0] = i;<br>
> index[1] = j;<br>
> index[2] = k;<br>
> connectedThreshold->AddSeed( index );<br>
> }<br>
><br>
> connectedThreshold->Print(std::cout,17100);<br>
><br>
><br>
><br>
> typedef itk::MetaDataDictionary DictionaryType;<br>
><br>
> DictionaryType inputdict = reader->GetMetaDataDictionary();<br>
><br>
> writer->SetMetaDataDictionary( inputdict );<br>
><br>
> writer->SetFileName( argv[2] );<br>
><br>
> writer->SetInput( caster->GetOutput() );<br>
><br>
> try<br>
> {<br>
> writer->Update();<br>
> }<br>
> catch( itk::ExceptionObject & excep )<br>
> {<br>
> std::cerr << "Exception caught !" << std::endl;<br>
> std::cerr << excep << std::endl;<br>
> }<br>
><br>
> return 0;<br>
> }<br>
><br>
><br>
</div></div><div><div></div><div>> _____________________________________<br>
> Powered by <a href="http://www.kitware.com" target="_blank">www.kitware.com</a><br>
><br>
> Visit other Kitware open-source projects at<br>
> <a href="http://www.kitware.com/opensource/opensource.html" target="_blank">http://www.kitware.com/opensource/opensource.html</a><br>
><br>
> Kitware offers ITK Training Courses, for more information visit:<br>
> <a href="http://www.kitware.com/products/protraining.html" target="_blank">http://www.kitware.com/products/protraining.html</a><br>
><br>
> Please keep messages on-topic and check the ITK FAQ at:<br>
> <a href="http://www.itk.org/Wiki/ITK_FAQ" target="_blank">http://www.itk.org/Wiki/ITK_FAQ</a><br>
><br>
> Follow this link to subscribe/unsubscribe:<br>
> <a href="http://www.itk.org/mailman/listinfo/insight-users" target="_blank">http://www.itk.org/mailman/listinfo/insight-users</a><br>
><br>
><br>
</div></div></blockquote></div><br>
</div></div></blockquote></div><br>