--- /dev/null
+
+ /*=========================================================================
+
+ Program: gdcm
+ Module: $RCSfile: exPresentationState.cxx,v $
+ Language: C++
+ Date: $Date: 2009/05/20 09:45:30 $
+ Version: $Revision: 1.1 $
+
+ Copyright (c) CREATIS (Centre de Recherche et d'Applications en Traitement de
+ l'Image). All rights reserved. See Doc/License.txt or
+ http://www.creatis.insa-lyon.fr/Public/Gdcm/License.html for details.
+
+ This software is distributed WITHOUT ANY WARRANTY; without even
+ the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
+ PURPOSE. See the above copyright notices for more information.
+
+=========================================================================*/
+
+#include "gdcmFile.h"
+// ----------
+//#include "gdcmValEntry.h"
+#include "gdcmDataEntry.h"
+// ----------
+#include "gdcmSeqEntry.h"
+#include "gdcmSQItem.h"
+#include "gdcmDocEntrySet.h"
+#include "gdcmSerieHelper.h"
+#include "gdcmDirList.h"
+#include "gdcmUtil.h"
+
+#include "gdcmArgMgr.h"
+
+ GDCM_NAME_SPACE::SeqEntry *CheckIfSequenceExists(GDCM_NAME_SPACE::File *fPR, uint16_t gr, uint16_t el);
+ GDCM_NAME_SPACE::SeqEntry *CheckIfSequenceExists(GDCM_NAME_SPACE::SQItem *si, uint16_t gr, uint16_t el);
+ bool dealWithTopLevelItem(GDCM_NAME_SPACE::SQItem* currentItem);
+ bool dealWithEndLevelItem(GDCM_NAME_SPACE::SQItem* currentItem);
+ void displaySeekResult(GDCM_NAME_SPACE::SeqEntry* currentItem, uint16_t g, uint16_t e);
+ void dealWithSequence(GDCM_NAME_SPACE::SeqEntry* se);
+// -------------------------------------------------------------------------------------------------------
+
+int main(int argc, char *argv[])
+{
+ START_USAGE(usage)
+ "\n exPresentationState :\n ",
+ "Extracts and displays the Graphic annotation / Text Objet Sequences of gdcm-parsable Dicom file",
+ " ",
+ "usage: exPresentationState {filein=inputFileName|dirin=inputDirectoryName}",
+ " [ { [noshadowseq] | [noshadow][noseq] } ] [debug] ",
+ " filein : Name of the image file ",
+ " PRFile : Name of the PresentationState file : ",
+ " noshadowseq: user doesn't want to load Private Sequences ",
+ " noshadow : user doesn't want to load Private groups (odd number) ",
+ " noseq : user doesn't want to load Sequences ",
+ " verbose : developper wants to run the program in 'verbose mode' ",
+ " debug : developper wants to run the program in 'debug mode' ",
+ FINISH_USAGE
+
+ // ----- Initialize Arguments Manager ------
+
+ GDCM_NAME_SPACE::ArgMgr *am = new GDCM_NAME_SPACE::ArgMgr(argc, argv);
+
+ if (am->ArgMgrDefined("usage") || argc == 1)
+ {
+ am->ArgMgrUsage(usage); // Display 'usage'
+ delete am;
+ return 0;
+ }
+
+ if (am->ArgMgrDefined("debug"))
+ GDCM_NAME_SPACE::Debug::DebugOn();
+
+ bool verbose = am->ArgMgrDefined("verbose");
+
+ const char *fileName = am->ArgMgrGetString("filein");
+ const char *PRName = am->ArgMgrGetString("PRfile");
+
+ /* if unused Param we give up */
+ if ( am->ArgMgrPrintUnusedLabels() )
+ {
+ am->ArgMgrUsage(usage);
+ delete am;
+ return 0;
+ }
+ delete am; // ------ we don't need Arguments Manager any longer ------
+
+// ============================================================
+// Read the input image.
+// ============================================================
+
+ GDCM_NAME_SPACE::File *f = GDCM_NAME_SPACE::File::New( );
+
+/*
+ //f->SetLoadMode(GDCM_NAME_SPACE::LD_NOSEQ | GDCM_NAME_SPACE::LD_NOSHADOW);
+ f->SetFileName( fileName );
+ f->SetMaxSizeLoadEntry(0xffff);
+ bool res = f->Load();
+
+ if( GDCM_NAME_SPACE::Debug::GetDebugFlag())
+ {
+ std::cout << "---------------------------------------------" << std::endl;
+ f->Print();
+ std::cout << "---------------------------------------------" << std::endl;
+ }
+ if (!res) {
+ std::cerr << "Sorry, " << fileName << " not a gdcm-readable "
+ << "DICOM / ACR File"
+ << std::endl;
+ f->Delete();
+ return 1;
+ }
+ std::cout << " ... is readable " << std::endl;
+*/
+
+// =================================================================================
+
+ GDCM_NAME_SPACE::File *fPR = GDCM_NAME_SPACE::File::New( );
+ fPR->SetFileName( PRName );
+ fPR->SetMaxSizeLoadEntry(0xffff);
+ bool res2 = fPR->Load();
+
+ if (!res2) {
+ std::cout << "Sorry, " << PRName << " not a gdcm-readable "
+ << "DICOM / ACR File"
+ << std::endl;
+ f->Delete();
+ return 1;
+ }
+ std::cout << " ... is readable " << std::endl;
+
+ GDCM_NAME_SPACE::SeqEntry *se;
+
+
+ se = CheckIfSequenceExists( fPR, 0x0070, 0x0001);
+ //displaySeekResult(se, 0x0070, 0x0001);
+ if (!se)
+ {
+ std::cout << "Hopeless : " << std::hex << 0x0070 << "|" << 0x0001 << std::dec << " doesn't exist..." <<std::endl;
+ exit (0);
+ }
+
+ std::cout << "[" << PRName << "] is a gdcm-readable PresentationState file, that holds one or more 'ROI'\n" <<std::endl;
+
+ GDCM_NAME_SPACE::SQItem* currentItem = se->GetFirstSQItem(); // Get the first 'ROI'
+ int i =0;
+ while (currentItem != NULL)
+ {
+ std::cout << "======== Deal With 'ROI' n." << i << std::endl;
+
+ // do what you want with the current 'ROI'
+ dealWithTopLevelItem(currentItem);
+ //...
+
+ currentItem = se->GetNextSQItem(); // Get the next 'ROI'
+ i++;
+ }
+
+ /*
+ bool res3 = TestPresentationState(f, fPR);
+
+ if (res3)
+ {
+ std::cout << "[" << PRName << "] is a gdcm-readable PresentationState file" <<std::endl;
+ }
+ else
+ {
+ std::cout << "Sorry, [" << PRName << "] is not a gdcm-readable PresentationState file" <<std::endl;
+ }
+*/
+
+
+ f->Delete();
+ fPR->Delete();
+return 0;
+}
+
+//----------------------------------------------------------------------------------------------------
+
+bool dealWithTopLevelItem(GDCM_NAME_SPACE::SQItem* currentItem)
+{
+
+ int tabElement[] = { 0x0008, 0x0009, 0x005a, 0x0060, 0x0086,
+ 0x0308, 0x0309, 0x030A, 0x030d, 0x0311, 0x0314,
+ 0x0318, 0x031c, 0x031e, 0x0402, 0x0404, 0x0000 };
+
+ bool res = false;
+ GDCM_NAME_SPACE::SeqEntry *se;
+
+ se = CheckIfSequenceExists(currentItem, 0x0008, 0x1140);
+ displaySeekResult(se, 0x0008, 0x1140);
+ if (se != 0)
+ {
+ res = true;
+ dealWithSequence(se);
+ }
+
+ for(int i=0; tabElement[i]!=0x0000; i++)
+ {
+ se = CheckIfSequenceExists(currentItem, 0x0070, tabElement[i]);
+ displaySeekResult(se, 0x0070, tabElement[i]);
+ if (se != 0)
+ {
+ res = true;
+ dealWithSequence(se);
+ }
+ }
+ return (res);
+}
+
+
+//----------------------------------------------------------------------------------------------------
+
+bool dealWithEndLevelItem(GDCM_NAME_SPACE::SQItem* currentItem)
+{
+
+ int tabElement[] = { 0x0008, 0x0009, 0x005a, 0x0060, 0x0086,
+ 0x0308, 0x0309, 0x030A, 0x030d, 0x0311, 0x0314,
+ 0x0318, 0x031c, 0x031e, 0x0402, 0x0404, 0x0000 };
+
+ bool res = false;
+ GDCM_NAME_SPACE::SeqEntry *se;
+
+ se = CheckIfSequenceExists(currentItem, 0x0008, 0x1140);
+ displaySeekResult(se, 0x0008, 0x1140);
+ if (se != 0)
+ {
+ res = true;
+ dealWithSequence(se);
+ }
+
+ for(int i=0; tabElement[i]!=0x0000; i++)
+ {
+ se = CheckIfSequenceExists(currentItem, 0x0070, tabElement[i]);
+ displaySeekResult(se, 0x0070, tabElement[i]);
+ if (se != 0)
+ {
+ res = true;
+ dealWithSequence(se);
+ }
+ }
+ return (res);
+}
+//----------------------------------------------------------------------------------------------------
+
+void dealWithSequence(GDCM_NAME_SPACE::SeqEntry* se)
+{
+ uint16_t g = se->GetGroup();
+ uint16_t e = se->GetElement();
+ std::cout << std::hex << "\n------------------------ deal with " << g <<"|" << e << std::dec
+ << " " << se->GetName() << std::endl;
+
+ GDCM_NAME_SPACE::SQItem *si = se->GetFirstSQItem();
+ if (!si) {
+ std::cout << "Sequence " << std::hex << g <<"|" << e << std::dec << "has no Item ?!?" << std::endl;
+ return;
+ }
+
+ if (g == 0x0008) {
+
+ si->Print(std::cout);
+ } else if (g == 0x0070) {
+
+ si->Print(std::cout);
+ } else {
+ std::cout << "Unexpected Group " << std::hex << g << std::hex << std::endl;
+ }
+
+ si = se->GetNextSQItem();
+ if (si)
+ {
+ std::cout << "Sequence " << std::hex << g <<"|" << e << std::dec << "should have only ONE Item ?!?" << std::endl;
+ si->Print(std::cout);
+ return;
+ }
+}
+
+//----------------------------------------------------------------------------------------------------
+
+void displaySeekResult(GDCM_NAME_SPACE::SeqEntry* se, uint16_t g, uint16_t e)
+{
+ if (se)
+ {
+ // std::cout << std::hex << g << "|" << e << std::dec << " [" << se->GetName() << "] exists" <<std::endl;
+ }
+ else
+ {
+ std::cout << " No " << std::hex << g << "|" << e << std::dec << " found" <<std::endl;
+ }
+}
+
+//----------------------------------------------------------------------------------------------------
+
+bool TestPresentationState(GDCM_NAME_SPACE::File *f, GDCM_NAME_SPACE::File *fPR)
+{
+
+/*------------------------------------------------------
+
+0070 0001 SQ 1 Graphic Annotation Sequence
+0070 0008 SQ 1 Text Object Sequence
+0070 0009 SQ 1 Graphic Object Sequence
+0070 005a SQ 1 Displayed Area Selection Sequence
+0070 0060 SQ 1 Graphic Layer Sequence
+0070 0086 SQ 1 Content Creator's Identification Code Sequence
+0070 0308 SQ 1 Registration Sequence
+0070 0309 SQ 1 Matrix Registration Sequence
+0070 030a SQ 1 Matrix Sequence
+0070 030d SQ 1 Registration Type Code Sequence
+0070 0311 SQ 1 Fiducial Identifier Code Sequence
+0070 0314 SQ 1 Used Fiducials Sequence
+0070 0318 SQ 1 Graphic Coordinates Data Sequence
+0070 031c SQ 1 Fiducial Set Sequence
+0070 031e SQ 1 Fiducial Sequence
+0070 0402 SQ 1 Blending Sequence
+0070 0404 SQ 1 Referenced Spatial Registration Sequence
+------------------------------------------------------- */
+
+/*------------------------------------------------------
+Relevant part of Dicom V3 Dict
+
+0070 0000 UL 1 Group Length
+0070 0001 SQ 1 Graphic Annotation Sequence
+0070 0002 CS 1 Graphic Layer
+0070 0003 CS 1 Bounding Box Annotation Units
+0070 0004 CS 1 Anchor Point Annotation Units
+0070 0005 CS 1 Graphic Annotation Units
+0070 0006 ST 1 Unformatted Text Value
+0070 0008 SQ 1 Text Object Sequence
+0070 0009 SQ 1 Graphic Object Sequence
+0070 0010 FL 2 Bounding Box Top Left Hand Corner
+0070 0011 FL 2 Bounding Box Bottom Right Hand Corner
+0070 0012 CS 1 Bounding Box Text Horizontal Justification
+0070 0014 FL 2 Anchor Point
+0070 0015 CS 1 Anchor Point Visibility
+0070 0020 US 1 Graphic Dimensions
+0070 0021 US 1 Number of Graphic Points
+0070 0022 FL 2-2n Graphic Data
+0070 0023 CS 1 Graphic Type
+0070 0024 CS 1 Graphic Filled
+0070 0041 CS 1 Image Horizontal Flip
+0070 0042 US 1 Image Rotation
+0070 0052 SL 2 Displayed Area Top Left Hand Corner
+0070 0053 SL 2 Displayed Area Bottom Right Hand Corner
+0070 005a SQ 1 Displayed Area Selection Sequence
+0070 0060 SQ 1 Graphic Layer Sequence
+0070 0062 IS 1 Graphic Layer Order
+0070 0066 US 1 Graphic Layer Recommended Display Grayscale Value
+0070 0067 US 3 Graphic Layer Recommended Display RGB Value (RET)
+0070 0068 LO 1 Graphic Layer Description
+0070 0080 CS 1 Content Label
+0070 0081 LO 1 Content Description
+0070 0082 DA 1 Presentation Creation Date
+0070 0083 TM 1 Presentation Creation Time
+0070 0084 PN 1 Content Creator's Name
+0070 0086 SQ 1 Content Creator's Identification Code Sequence
+0070 0100 CS 1 Presentation Size Mode
+0070 0101 DS 2 Presentation Pixel Spacing
+0070 0102 IS 2 Presentation Pixel Aspect Ratio
+0070 0103 FL 1 Presentation Pixel Magnification Ratio
+0070 0306 CS 1 Shape Type
+0070 0308 SQ 1 Registration Sequence
+0070 0309 SQ 1 Matrix Registration Sequence
+0070 030a SQ 1 Matrix Sequence
+0070 030c CS 1 Frame of Reference Transformation Matrix Type
+0070 030d SQ 1 Registration Type Code Sequence
+0070 030f ST 1 Fiducial Description
+0070 0310 SH 1 Fiducial Identifier
+0070 0311 SQ 1 Fiducial Identifier Code Sequence
+0070 0312 FD 1 Contour Uncertainty Radius
+0070 0314 SQ 1 Used Fiducials Sequence
+0070 0318 SQ 1 Graphic Coordinates Data Sequence
+0070 031a UI 1 Fiducial UID
+0070 031c SQ 1 Fiducial Set Sequence
+0070 031e SQ 1 Fiducial Sequence
+0070 0401 US 3 Graphic Layer Recommended Display CIELab Value
+0070 0402 SQ 1 Blending Sequence
+0070 0403 FL 1 Relative Opacity
+0070 0404 SQ 1 Referenced Spatial Registration Sequence
+0070 0405 CS 1 Blending Position
+------------------------------------------------------- */
+
+}
+
+
+// ----------------------------------------------------------------------------------
+
+
+
+// ----------------------------------------------------------------------------------
+
+
+GDCM_NAME_SPACE::SeqEntry *CheckIfSequenceExists( GDCM_NAME_SPACE::File *fPR, uint16_t gr, uint16_t el)
+{
+ GDCM_NAME_SPACE::SeqEntry *se= fPR->GetSeqEntry(gr, el);
+ return se;
+}
+
+// ----------------------------------------------------------------------------------
+
+GDCM_NAME_SPACE::SeqEntry *CheckIfSequenceExists( GDCM_NAME_SPACE::SQItem *si, uint16_t gr, uint16_t el)
+{
+ GDCM_NAME_SPACE::SeqEntry *se= si->GetSeqEntry(gr, el);
+ return se;
+}
projects / gdcm.git / commitdiff