]> Creatis software - clitk.git/commitdiff
Cosmetic
authortbaudier <thomas.baudier@creatis.insa-lyon.fr>
Fri, 11 Jan 2019 12:43:06 +0000 (13:43 +0100)
committertbaudier <thomas.baudier@creatis.insa-lyon.fr>
Fri, 11 Jan 2019 12:43:06 +0000 (13:43 +0100)
Remove trailing spaces and indent correctly

tools/clitkDicomRTPlan2Gate.cxx

index 3f07c8ea44ef1d760ff7549e897630bdc96eab4e..c60803405ea1fbfa24a139c0014b5d8b7c389d20 100644 (file)
@@ -38,8 +38,7 @@ int main(int argc, char * argv[])
 
   gdcm::Reader reader;
   reader.SetFileName(input_file);
-  if( !reader.Read() )
-  {
+  if(!reader.Read()) {
     std::cerr << "Failed to read: " << input_file << std::endl;
     return 1;
   }
@@ -49,81 +48,77 @@ int main(int argc, char * argv[])
   //Check file type
   gdcm::MediaStorage ms;
   ms.SetFromFile(reader.GetFile());
-  if( ms != gdcm::MediaStorage::RTIonPlanStorage )
-  {
+  if(ms != gdcm::MediaStorage::RTIonPlanStorage) {
     std::cerr << "The file <" << input_file << "> is not a RT Ion Plan." << std::endl;
     exit(0);
   }
-  
+
   std::ofstream output;
   output.open(output_file);
   output << "#TREATMENT-PLAN-DESCRIPTION" << "\n";
-  
+
   gdcm::Attribute<0x300a,0x2> label;
   label.SetFromDataSet(ds);
   gdcm::Attribute<0x300a,0x3> planName;
   planName.SetFromDataSet(ds);
-  
+
   output << "#PlanName\n"<< planName.GetValue() << "\n";
-  
+
   gdcm::Tag tFractionGroupSQ(0x300a,0x70);
   gdcm::Attribute<0x300a,0x71> fractionID;
   gdcm::Attribute<0x300a,0x78> numberOfFractions;
   gdcm::Attribute<0x300a,0x80> numberOfBeams;
-  
+
   gdcm::Tag tReferencedBeamSQ(0x300c,0x4);
   gdcm::Attribute<0x300a,0x86> meterSet;
   gdcm::Attribute<0x300c,0x6> beamID;
-  
+
   const gdcm::DataElement &FractionGroupSQ = ds.GetDataElement( tFractionGroupSQ );
   gdcm::SmartPointer<gdcm::SequenceOfItems> sqf = FractionGroupSQ.GetValueAsSQ();
   gdcm::Item &ReferencedBeamItem = sqf->GetItem(1);
   gdcm::DataSet &ReferencedBeamData = ReferencedBeamItem.GetNestedDataSet();
   numberOfFractions.SetFromDataSet(ReferencedBeamData);
-  
+
   output << "#NumberOfFractions\n"<< numberOfFractions.GetValue() << "\n";
-  
-  for( unsigned int fidx = 1; fidx <= sqf->GetNumberOfItems(); ++fidx )
-  {
+
+  for(unsigned int fidx = 1; fidx <= sqf->GetNumberOfItems(); ++fidx) {
     gdcm::Item &ReferencedBeamItem = sqf->GetItem(fidx);
     gdcm::DataSet &ReferencedBeamData = ReferencedBeamItem.GetNestedDataSet();
-    
+
     fractionID.SetFromDataSet(ReferencedBeamData);
     numberOfFractions.SetFromDataSet(ReferencedBeamData);
     numberOfBeams.SetFromDataSet(ReferencedBeamData);
     output << "##FractionID\n"<<  fractionID.GetValue() << "\n";
     output << "##NumberOfFields\n"<<  numberOfBeams.GetValue() << "\n";
-        
+
     const gdcm::DataElement &ReferencedBeamSQ = ReferencedBeamData.GetDataElement(tReferencedBeamSQ);
-    gdcm::SmartPointer<gdcm::SequenceOfItems> sqb = ReferencedBeamSQ.GetValueAsSQ();   
-    for( unsigned int bidx = 1; bidx <= sqb->GetNumberOfItems(); ++bidx )
-    {
+    gdcm::SmartPointer<gdcm::SequenceOfItems> sqb = ReferencedBeamSQ.GetValueAsSQ();
+    for(unsigned int bidx = 1; bidx <= sqb->GetNumberOfItems(); ++bidx) {
       gdcm::Item &ReferencedBeamGroupItem = sqb->GetItem(bidx);
       gdcm::DataSet &ReferencedBeamGroupData = ReferencedBeamGroupItem.GetNestedDataSet();
-      
+
       meterSet.SetFromDataSet(ReferencedBeamGroupData);
       beamID.SetFromDataSet(ReferencedBeamGroupData);
       output << "###FieldsID\n"<<  beamID.GetValue() << "\n";
     }
   }
-  
+
   gdcm::Tag tPatientSetupSQ(0x300a,0x180);
   gdcm::Attribute<0x0018,0x5100> position;
   gdcm::Attribute<0x300a,0x1b0> setupTechnique;
   gdcm::Attribute<0x300a,0x182> setupID;
-  
+
   const gdcm::DataElement &PatientSetupSQ = ds.GetDataElement( tPatientSetupSQ );
   gdcm::SmartPointer<gdcm::SequenceOfItems> sqs = PatientSetupSQ.GetValueAsSQ();
-  for( unsigned int sidx = 1; sidx <= sqs->GetNumberOfItems(); ++sidx )
-  {
+  for(unsigned int sidx = 1; sidx <= sqs->GetNumberOfItems(); ++sidx) {
     gdcm::Item &PatientSetupItem = sqs->GetItem(sidx);
     gdcm::DataSet &PatientSetupData = PatientSetupItem.GetNestedDataSet();
-    
+
     position.SetFromDataSet(PatientSetupData);
     setupTechnique.SetFromDataSet(PatientSetupData);
     setupID.SetFromDataSet(PatientSetupData);
   }
-  
+
   gdcm::Tag tIonBeamSQ(0x300a,0x3a2);
   gdcm::Attribute<0x300a,0xb2> machine;
   gdcm::Attribute<0x300a,0xc6> radiationType;
@@ -139,10 +134,10 @@ int main(int argc, char * argv[])
   gdcm::Attribute<0x300a,0x312> numberOfRangeShifters;
   gdcm::Attribute<0x300a,0x330> numberOfLateralSpreadingDevices;
   gdcm::Attribute<0x300a,0x340> numberOfRangeModulators;
-  
+
   //gdcm::Tag tSnoutSQ(0x300a,0x30c);
   //gdcm::Attribute<0x300a,0x30f> snoutID;
-  
+
   gdcm::Tag tControlPointSQ(0x300a,0x3a8);
   gdcm::Attribute<0x300a,0x15> energyUnit;
   gdcm::Attribute<0x300a,0x114> energyValue;
@@ -150,36 +145,34 @@ int main(int argc, char * argv[])
   gdcm::Attribute<0x300a,0x134> cumulativeMetersetWeight;
   gdcm::Attribute<0x300a,0x392> numberOfScannedSpots;
   gdcm::Attribute<0x300a,0x390> spotTunnedID;
-  
+
   gdcm::Tag tSpotPositionsData(0x300a,0x394);
   gdcm::Attribute<0x300a,0x394> spotPositions;
   gdcm::Tag tSpotWeightsData(0x300a,0x396);
   gdcm::Attribute<0x300a,0x396> spotWeights;
-  
+
   gdcm::Attribute<0x300a,0x11e> gantryAngle;
   gdcm::Attribute<0x300a,0x122> patientSupportAngle;
   gdcm::Attribute<0x300a,0x12c> isocenterPosition;
-  
+
   long totalCumMeterSet = 0;
-  
+
   const gdcm::DataElement &IonBeamSQ = ds.GetDataElement( tIonBeamSQ );
   gdcm::SmartPointer<gdcm::SequenceOfItems> sqi = IonBeamSQ.GetValueAsSQ();
-  
-  for( unsigned int iidx = 1; iidx <= sqi->GetNumberOfItems(); ++iidx )
-  {
+
+  for(unsigned int iidx = 1; iidx <= sqi->GetNumberOfItems(); ++iidx) {
     gdcm::Item &IonBeamItem = sqi->GetItem(iidx);
     gdcm::DataSet &IonBeamData = IonBeamItem.GetNestedDataSet();
     finalCumulativeMeterSetWeight.SetFromDataSet(IonBeamData);
     totalCumMeterSet += finalCumulativeMeterSetWeight.GetValue();
   }
-  
+
   output << "#TotalMetersetWeightOfAllFields\n" <<  totalCumMeterSet << "\n";
-  
-  for( unsigned int iidx = 1; iidx <= sqi->GetNumberOfItems(); ++iidx )
-  {
+
+  for(unsigned int iidx = 1; iidx <= sqi->GetNumberOfItems(); ++iidx) {
     gdcm::Item &IonBeamItem = sqi->GetItem(iidx);
     gdcm::DataSet &IonBeamData = IonBeamItem.GetNestedDataSet();
-    
+
     machine.SetFromDataSet(IonBeamData);
     radiationType.SetFromDataSet(IonBeamData);
     treatmentType.SetFromDataSet(IonBeamData);
@@ -194,28 +187,27 @@ int main(int argc, char * argv[])
     numberOfRangeShifters.SetFromDataSet(IonBeamData);
     numberOfLateralSpreadingDevices.SetFromDataSet(IonBeamData);
     numberOfRangeModulators.SetFromDataSet(IonBeamData);
-    
+
     /*
     const gdcm::DataElement &SnoutSQ = IonBeamData.GetDataElement(tSnoutSQ);
     gdcm::SmartPointer<gdcm::SequenceOfItems> sqsn = SnoutSQ.GetValueAsSQ();
-    for( unsigned int snidx = 1; snidx <= sqsn->GetNumberOfItems(); ++snidx )
-    {
+    for(unsigned int snidx = 1; snidx <= sqsn->GetNumberOfItems(); ++snidx) {
       gdcm::Item &SnoutItem = sqsn->GetItem(snidx);
       gdcm::DataSet &SnoutData = SnoutItem.GetNestedDataSet();
-      
+
       snoutID.SetFromDataSet(SnoutData);
     }
     */
-    
+
     const gdcm::DataElement &ControlPointSQ = IonBeamData.GetDataElement(tControlPointSQ);
     gdcm::SmartPointer<gdcm::SequenceOfItems> sqcp = ControlPointSQ.GetValueAsSQ();
     gdcm::Item &ControlPointItem = sqcp->GetItem(1);
     gdcm::DataSet &ControlPointData = ControlPointItem.GetNestedDataSet();
-    
+
     gantryAngle.SetFromDataSet(ControlPointData);
     patientSupportAngle.SetFromDataSet(ControlPointData);
     isocenterPosition.SetFromDataSet(ControlPointData);
-    
+
     output << "\n#FIELD-DESCRIPTION" << "\n";
     output << "###FieldID\n"<< ionBeamID.GetValue() << "\n";
     output << "###FinalCumulativeMeterSetWeight\n"<<  finalCumulativeMeterSetWeight.GetValue() << "\n";
@@ -226,59 +218,51 @@ int main(int argc, char * argv[])
     output << " " << isocenterPosition.GetValue(2) << "\n";
     output << "###NumberOfControlPoints\n" << numberOfControlPoints.GetValue() << "\n\n";
     output << "#SPOTS-DESCRIPTION" << "\n";
-    
-    for( unsigned int cpidx = 1; cpidx <= sqcp->GetNumberOfItems(); ++cpidx )
-    {
+
+    for(unsigned int cpidx = 1; cpidx <= sqcp->GetNumberOfItems(); ++cpidx) {
       gdcm::Item &ControlPointItem = sqcp->GetItem(cpidx);
       gdcm::DataSet &ControlPointData = ControlPointItem.GetNestedDataSet();
-      
-      if ( cpidx == 1 )
-      {
-       gantryAngle.SetFromDataSet(ControlPointData);
-       patientSupportAngle.SetFromDataSet(ControlPointData);
-       isocenterPosition.SetFromDataSet(ControlPointData);
+
+      if(cpidx == 1) {
+        gantryAngle.SetFromDataSet(ControlPointData);
+        patientSupportAngle.SetFromDataSet(ControlPointData);
+        isocenterPosition.SetFromDataSet(ControlPointData);
       }
-      
-      if ( cpidx < sqcp->GetNumberOfItems() )
-      {
-       energyUnit.SetFromDataSet(ControlPointData);
-       energyValue.SetFromDataSet(ControlPointData);
-       controlPointIndex.SetFromDataSet(ControlPointData);
-       cumulativeMetersetWeight.SetFromDataSet(ControlPointData);
-       numberOfScannedSpots.SetFromDataSet(ControlPointData);
-       spotTunnedID.SetFromDataSet(ControlPointData);
-       
-       const gdcm::DataElement &spotPositionsData = ControlPointData.GetDataElement(tSpotPositionsData);
-       spotPositions.SetFromDataElement(spotPositionsData);
-       const gdcm::DataElement &spotWeightsData = ControlPointData.GetDataElement(tSpotWeightsData);
-       spotWeights.SetFromDataElement(spotWeightsData);
+
+      if(cpidx < sqcp->GetNumberOfItems()) {
+        energyUnit.SetFromDataSet(ControlPointData);
+        energyValue.SetFromDataSet(ControlPointData);
+        controlPointIndex.SetFromDataSet(ControlPointData);
+        cumulativeMetersetWeight.SetFromDataSet(ControlPointData);
+        numberOfScannedSpots.SetFromDataSet(ControlPointData);
+        spotTunnedID.SetFromDataSet(ControlPointData);
+
+        const gdcm::DataElement &spotPositionsData = ControlPointData.GetDataElement(tSpotPositionsData);
+        spotPositions.SetFromDataElement(spotPositionsData);
+        const gdcm::DataElement &spotWeightsData = ControlPointData.GetDataElement(tSpotWeightsData);
+        spotWeights.SetFromDataElement(spotWeightsData);
       }
-      
-      if ( cpidx == sqcp->GetNumberOfItems() )
-      {
-       energyUnit.SetValue( "0" );
-       energyValue.SetValue( 0 );
-       controlPointIndex.SetFromDataSet(ControlPointData);
-       cumulativeMetersetWeight.SetFromDataSet(ControlPointData);
-       numberOfScannedSpots.SetValue( 0 );
-       spotTunnedID.SetValue( "0" );
+
+      if(cpidx == sqcp->GetNumberOfItems()) {
+        energyUnit.SetValue( "0" );
+        energyValue.SetValue( 0 );
+        controlPointIndex.SetFromDataSet(ControlPointData);
+        cumulativeMetersetWeight.SetFromDataSet(ControlPointData);
+        numberOfScannedSpots.SetValue( 0 );
+        spotTunnedID.SetValue( "0" );
       }
-      
+
       output << "####ControlPointIndex\n" << controlPointIndex.GetValue() << "\n";
       output << "####SpotTunnedID\n" << spotTunnedID.GetValue() << "\n";
       output << "####CumulativeMetersetWeight\n" << cumulativeMetersetWeight.GetValue() << "\n";
       output << "####Energy (MeV)\n"<< energyValue.GetValue() << "\n";
       output << "####NbOfScannedSpots\n" << numberOfScannedSpots.GetValue() << "\n";
       output << "####X Y Weight\n";
-      for ( unsigned int sp=0; sp < numberOfScannedSpots.GetValue(); ++sp)
-      {
-       output << spotPositions.GetValue(2*sp) << " " << spotPositions.GetValue(2*sp+1) << " " << spotWeights.GetValue(sp) << "\n";
-      }
-      
+      for( unsigned int sp=0; sp < numberOfScannedSpots.GetValue(); ++sp)
+        output << spotPositions.GetValue(2*sp) << " " << spotPositions.GetValue(2*sp+1) << " " << spotWeights.GetValue(sp) << "\n";
     }
-    
   }
-    
+
   return 0;
 }