LDicomWaveformChannel::SetChannelSamples8

/******************************************************************/ 
/*   This is a comprehensive sample, which shows how to insert a  */ 
/*   waveform group with one ECG channel into a dataset.          */ 
/*   The main function is InsertECGCWaveform ,and the rest of the */ 
/*   functions are helping functions.                             */ 
/******************************************************************/ 
 
// Function prototypes 
BOOL SetChannelSourceAndSensitivity(LDicomWaveformChannel*  pECGChannel);  
BOOL SetChannelAnnotations(LDicomWaveformChannel*  pECGChannel);  
BOOL InsertECGCChannel( LDicomWaveformGroup* pECGWaveformGroup,  
                        L_INT16* pSamples,  
                        L_UINT32 uSampleCount);  
 
// The main function that creates the  
// waveform group and adds to the dataset 
L_INT LDicomWaveformChannel_SetChannelSamples8Example(LDicomDS &InputDS, L_INT16* pSamples, L_UINT32 uSampleCount) 
{ 
   L_INT nRet; 
   // Sanity check 
   if(!pSamples || (uSampleCount == 0) ) 
   { 
      return -1; 
   } 
 
   // Our new waveform group 
   LDicomWaveformGroup  ECGWaveformGroup;  
 
   // Reset the waveform group, we don't  
   // really need to call this!  
   ECGWaveformGroup.Reset(); 
 
   // Set the number of samples per channel. You can call  
   // LDicomWaveformGroup::GetNumberOfSamplesPerChannel to get  
   // the number of samples 
   if(!ECGWaveformGroup.SetNumberOfSamplesPerChannel(uSampleCount))  
   { 
      return -1;  
   } 
 
   // Set the sampling frequency. You can call  
   // LDicomWaveformGroup::GetSamplingFrequency to  
   // get the sampling frequency 
   ECGWaveformGroup.SetSamplingFrequency(240.00);  
   // Set sample interpretation 
   if(!ECGWaveformGroup.SetSampleInterpretation(DICOM_SAMPLE_INTERPRETATION_SS)) 
   { 
      return -1;  
   } 
   // Just in case!  
   if(ECGWaveformGroup.GetSampleInterpretation() != DICOM_SAMPLE_INTERPRETATION_SS)  
   { 
      return -1;  
   } 
 
   // No Multiplex Group Time Offset 
   // You can call LDicomWaveformGroup::GetMultiplexGroupTimeOffset  
   // to get the time offset 
   if(!ECGWaveformGroup.SetMultiplexGroupTimeOffset(FALSE, 0.0))  
   { 
      return -1;  
   } 
 
   // No Trigger Time Offset 
   // You can call LDicomWaveformGroup::GetTriggerTimeOffset  
   // to get the trigger time offset 
   if(!ECGWaveformGroup.SetTriggerTimeOffset(FALSE, 0.0))  
   { 
      return -1; 
   } 
 
   // No Trigger Sample Position 
   // You can call LDicomWaveformGroup::GetTriggerSamplePosition 
   // to get the trigger sample position 
   if(!ECGWaveformGroup.SetTriggerSamplePosition(FALSE, 0))  
   { 
      return -1;  
   } 
 
   // Waveform originality is original    
   ECGWaveformGroup.SetWaveformOriginality(DICOM_WAVEFORM_ORIGINALITY_ORIGINAL);  
 
   // Just in case 
   if(ECGWaveformGroup.GetWaveformOriginality() != DICOM_WAVEFORM_ORIGINALITY_ORIGINAL)  
   { 
      return -1; 
   } 
 
   // Set the multiplex group label 
   if(!ECGWaveformGroup.SetMultiplexGroupLabel(TEXT("SCPECG Waveform"))) 
   { 
      return -1;  
   } 
   if(ECGWaveformGroup.GetMultiplexGroupLabel() == NULL)  
   { 
      return -1;  
   } 
 
   // Set the Waveform padding value 
   // You can call LDicomWaveformGroup::GetWaveformPaddingValue 
   // to get the waveform padding value 
   if(!ECGWaveformGroup.SetWaveformPaddingValue(TRUE,32768))  
   { 
      return -1; 
   } 
 
   if (!InsertECGCChannel(&ECGWaveformGroup , pSamples, uSampleCount))  
   { 
      return -1;  
   } 
 
   // Delete any waveform groups that already exist in the dataset 
   L_UINT32 uWaveformGroupCount = InputDS.GetWaveformGroupCount(); 
   for(L_UINT32 uIndex = 0 ; uIndex < uWaveformGroupCount ; uIndex++) 
   { 
      nRet = InputDS.DeleteWaveformGroup(uIndex,0);  
      if(nRet != DICOM_SUCCESS) 
         return nRet; 
   } 
   // Insert the new waveform group into the dataset 
   nRet = InputDS.AddWaveformGroup(&ECGWaveformGroup, 0, ELEMENT_INDEX_MAX); 
   if(nRet != DICOM_SUCCESS) 
      return nRet; 
 
   return DICOM_SUCCESS; 
} 
 
// Add an ECG channel to the group 
BOOL InsertECGCChannel (LDicomWaveformGroup* pECGWaveformGroup, L_INT16* pSamples, L_UINT32 uSampleCount)  
{ 
   LDicomWaveformChannel*  pECGChannel= NULL;  
   L_UINT32 uNumberOfSamples  =0;  
   L_INT32  nChannelNumber =  -1;  
   L_INT32  nMinValue      =  0;  
   L_INT32  nMaxValue      =  0;  
 
   if(pECGWaveformGroup == NULL)  
   { 
      return FALSE; 
   } 
 
   // Add a channel to the group 
   pECGChannel = pECGWaveformGroup->AddChannel(/*ELEMENT_INDEX_MAX*/); 
   if(pECGChannel == NULL)  
   { 
      return FALSE;  
   } 
   // If we want to update a channel, we can first call  
   // LDicomWaveformGroup::DeleteChannel and then call  
   // LDicomWaveformGroup::AddChannel with the same  
   // index we passed to DeleteChannel 
 
 
   // Make sure that the channel really got added  
   // This is for the purposes of this sample only, because  
   // the check we did in the previous statement is enough 
   if(   (pECGWaveformGroup->GetNumberOfChannels() !=1)||  
         (pECGWaveformGroup->GetChannel(0) ==NULL) )  
   { 
      return FALSE;    
   }    
   // Sanity checks!  
   if((pECGChannel->GetWaveformGroup() != pECGWaveformGroup)||  
      (pECGChannel->GetIndex() != 0))  
   { 
      return FALSE;    
   }    
   /* 
         Set the channel samples 
         The data we are setting in here is 16 bit data.  
         We would call SetChannelSamples8 if the data were 8 bit.  
         We can also call SetChannelSamples32, which will set the data  
         as either 8-bit or 16-bit depending on the sample interpretation    
   */ 
   if(pECGChannel->SetChannelSamples16(pSamples, uSampleCount)!=uSampleCount)  
   { 
      return FALSE;  
   } 
   // Just in case 
   if((!pECGChannel->GetChannelSamples(&uNumberOfSamples))||  
      (uNumberOfSamples == 0))  
   { 
      return FALSE;    
   } 
   // Set the channel source and sensitivity 
   if(!SetChannelSourceAndSensitivity(pECGChannel))  
   { 
      return FALSE;  
   } 
   // Set the channel status 
   pECGChannel->SetChannelStatus(DICOM_CHANNEL_STATUS_OK);  
   if(pECGChannel->GetChannelStatus()!=DICOM_CHANNEL_STATUS_OK)  
  { 
      return FALSE;    
   } 
   // Set the channel time skew. You can call  
   // LDicomWaveformChannel::GetChannelTimeSkew to get the time skew 
   // You can also call LDicomWaveformChannel::SetChannelSampleSkew 
   // And LDicomWaveformChannel::GetChannelSampleSkew to set and get  
   // the channel samples skew 
   pECGChannel->SetChannelTimeSkew(0.000000);  
 
 
   // Set the waveform channel number 
   if(!pECGChannel->SetWaveformChannelNumber(TRUE, 0))  
   { 
      return FALSE;    
   } 
   if((!pECGChannel->GetWaveformChannelNumber(&nChannelNumber))||  
      (nChannelNumber !=0))  
   { 
      return FALSE;    
   } 
   // Set the channel label    
   if(!pECGChannel->SetChannelLabel(TEXT("First Channel"))) 
   { 
      return FALSE;    
   } 
   if(!pECGChannel->GetChannelLabel()) 
   { 
      return FALSE;    
   } 
   // No channel offset 
   // You can also call LDicomWaveformChannel::GetChannelOffset 
   // To get the channel offset 
   if(!pECGChannel->SetChannelOffset(FALSE, 0.0))  
   { 
      return FALSE;       
   } 
   // Set filter low frequency 
   // You can also call LDicomWaveformChannel::GetFilterLowFrequency 
   // To get the filter low frequency    
   if(!pECGChannel->SetFilterLowFrequency(TRUE, 0.050))  
   { 
      return FALSE;       
   }    
   // Set filter high frequency 
   // You can also call LDicomWaveformChannel::GetFilterHighFrequency 
   // To get the filter high frequency    
   if(!pECGChannel->SetFilterHighFrequency(TRUE, 100.00))  
   { 
      return FALSE;       
   }    
   // Set the channel minimum value 
   if(!pECGChannel->SetChannelMinimumValue(TRUE, -386))  
   { 
      return FALSE;       
   }    
   // Set the channel maximum value 
   if(!pECGChannel->SetChannelMaximumValue(TRUE, 1264))  
   { 
      return FALSE;       
   } 
   // Just in case 
   if((!pECGChannel->GetChannelMinimumValue(&nMinValue))||  
      (nMinValue !=-386))  
   { 
      return FALSE;       
   } 
   // Just in case 
   if((!pECGChannel->GetChannelMaximumValue(&nMaxValue))||  
      (nMaxValue !=1264))  
   { 
      return FALSE;       
   }           
   /* 
      You can also call the following functions to set and get 
 the notch  
      filter frequency and bandwidth 
      LDicomWaveformChannel::SetNotchFilterFrequency 
      LDicomWaveformChannel::GetNotchFilterFrequency 
      LDicomWaveformChannel::SetNotchFilterBandwidth 
      LDicomWaveformChannel::GetNotchFilterBandwidth    
   */  
 
   // Last, but not least, set the channel annotations!  
   return SetChannelAnnotations(pECGChannel);  
} 
 
//Set the channel source and sensitivity 
BOOL SetChannelSourceAndSensitivity(LDicomWaveformChannel*  pECGChannel)  
{ 
       
   // Let's use the DICOM context group tables!  
 
   ///////////////////////Channel Source/////////////////////////    
   // Load the ECG Leads table  
   LDicomContextGroup::Load(CID_3001);    
   pDICOMCONTEXTGROUP pECGLeadsGroup = LDicomContextGroup::Find(CID_3001);  
   if(pECGLeadsGroup)  
   { 
      // 5.6.3-9-1 is Lead I (Einthoven)  
      pDICOMCODEDCONCEPT pCodedConcept = LDicomContextGroup::FindCodedConcept(pECGLeadsGroup,  
                                                                              TEXT("SCPECG"), 
                                                                              TEXT("5.6.3-9-1")); 
      if(pCodedConcept)  
      { 
         DICOMCODESEQUENCEITEM   DicomSourceSequenceItem;  
         L_UINT16                uRet;  
 
         memset(&DicomSourceSequenceItem , 0 , sizeof(DICOMCODESEQUENCEITEM));  
         DicomSourceSequenceItem.uStructSize = sizeof(DICOMCODESEQUENCEITEM);          
         
         DicomSourceSequenceItem.pszCodeValue              = pCodedConcept->pszCodeValue                 ;  
         DicomSourceSequenceItem.pszCodingSchemeDesignator = pCodedConcept->pszCodingSchemeDesignator    ;  
         DicomSourceSequenceItem.pszCodingSchemeVersion    = pCodedConcept->pszCodingSchemeVersion       ;  
         DicomSourceSequenceItem.pszCodeMeaning            = pCodedConcept->pszCodeMeaning               ;  
 
         // Set the channel source          
         uRet = pECGChannel->SetChannelSource(&DicomSourceSequenceItem);  
         if(uRet!= DICOM_SUCCESS)  
         { 
            return FALSE;  
         } 
         // You can call LDicomWaveformChannel::GetChannelSource () 
         // to get the channel source 
      } 
      else 
      { 
        return FALSE;  
      } 
   } 
   else 
   { 
      return FALSE;  
   } 
 
   ///////////////////////Channel Sensitivity//////////////////////// 
    
   // Cardiology Units of Measurement 
   LDicomContextGroup::Load(CID_3082);    
   pDICOMCONTEXTGROUP pCardiologyUnitsGroup = LDicomContextGroup::Find(CID_3082);  
   if(pCardiologyUnitsGroup)  
   { 
      pDICOMCODEDCONCEPT pCodedConcept = LDicomContextGroup::FindCodedConcept(pCardiologyUnitsGroup,  
                                                                              TEXT("UCUM"), 
                                                                              TEXT("mV")); 
      if(pCodedConcept)  
      { 
         DICOMCODESEQUENCEITEM   DicomSourceSequenceItem;  
         L_UINT16                uRet;  
 
         memset(&DicomSourceSequenceItem , 0 , sizeof(DICOMCODESEQUENCEITEM));  
         DicomSourceSequenceItem.uStructSize = sizeof(DICOMCODESEQUENCEITEM);          
         
         DicomSourceSequenceItem.pszCodeValue              = pCodedConcept->pszCodeValue                 ;  
         DicomSourceSequenceItem.pszCodingSchemeDesignator = pCodedConcept->pszCodingSchemeDesignator    ;  
         DicomSourceSequenceItem.pszCodingSchemeVersion    = pCodedConcept->pszCodingSchemeVersion       ;  
         DicomSourceSequenceItem.pszCodeMeaning            = pCodedConcept->pszCodeMeaning               ;  
 
         // Set the channel sensitivity 
         uRet = pECGChannel->SetChannelSensitivity(TRUE,  
                                                   0.001220,  
                                                   &DicomSourceSequenceItem);  
         if(uRet  != DICOM_SUCCESS)  
         { 
            return FALSE;  
         }          
         // You can call LDicomWaveformChannel::GetChannelSensitivity() 
         // to get the channel sensitivity 
      } 
      else 
      { 
         return FALSE;  
      } 
   } 
   else 
   { 
      return FALSE;  
   } 
 
   return TRUE;  
} 
 
// Add annotations for the channel 
BOOL SetChannelAnnotations(LDicomWaveformChannel*  pECGChannel)  
{ 
 
   // Delete any existing channel annotations 
   for(L_UINT32 uIndex = 0 ; uIndex < 
pECGChannel->GetAnnotationCount();uIndex++) 
   { 
      if(pECGChannel->GetAnnotation(uIndex))  
      { 
         pECGChannel->DeleteAnnotation(uIndex);  
      } 
   } 
 
   DICOMWAVEFORMANNOTATION ECGWaveformAnnotation;  
 
   memset(&ECGWaveformAnnotation,0, sizeof(DICOMWAVEFORMANNOTATION));  
   ECGWaveformAnnotation.uStructSize = sizeof(DICOMWAVEFORMANNOTATION);  
 
   // Over here we are populating the DICOMCODESEQUENCEITEM structure 
   // directly. What we can also do is fill it with values from 
   // the LDicomContextGroup table which lists all the context 
   // group tables defined by the DICOM standard.  You can look at 
   // the SetChannelSourceAndSensitivity above for an Example on 
   // how to do that 
 
   // 
   DICOMCODESEQUENCEITEM CodedName;  
   memset(&CodedName,0, sizeof(DICOMCODESEQUENCEITEM));  
   CodedName.uStructSize               = sizeof(DICOMCODESEQUENCEITEM);  
   CodedName.pszCodeValue              = TEXT("8867-4"); 
   CodedName.pszCodeMeaning            = TEXT("Heart rate"); 
   CodedName.pszCodingSchemeDesignator = TEXT("LN"); 
   CodedName.pszCodingSchemeVersion    = TEXT("19971101"); 
 
   ECGWaveformAnnotation.pCodedName = &CodedName;  
 
   DICOMCODESEQUENCEITEM MeasurementUnits;  
   memset(&MeasurementUnits,0, sizeof(DICOMCODESEQUENCEITEM));  
   MeasurementUnits.uStructSize               = 
sizeof(DICOMCODESEQUENCEITEM);  
   MeasurementUnits.pszCodeValue              = TEXT("{H.B.}/min"); 
   MeasurementUnits.pszCodeMeaning            = TEXT("Heart beat per minute"); 
   MeasurementUnits.pszCodingSchemeDesignator = TEXT("UCUM"); 
   MeasurementUnits.pszCodingSchemeVersion    = TEXT("1.4"); 
 
   ECGWaveformAnnotation.pMeasurementUnits = &MeasurementUnits;  
 
   L_DOUBLE dNumericValue = 69.00;  
   ECGWaveformAnnotation.pNumericValue = &dNumericValue ;  
   ECGWaveformAnnotation.uNumericValueCount = 1;  
 
   if(pECGChannel->AddAnnotation(&ECGWaveformAnnotation) != DICOM_SUCCESS)  
   { 
      return FALSE;  
   } 
 
   return TRUE;  
}