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; 
}