OTSEDSSim.cpp 7.3 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286
  1. #include "stdafx.h"
  2. #include "OTSEDSSim.h"
  3. namespace OTSController {
  4. COTSEDSSim::COTSEDSSim()
  5. {
  6. m_bDoQuantification = FALSE;
  7. }
  8. COTSEDSSim::~COTSEDSSim()
  9. {
  10. }
  11. void COTSEDSSim::FinishedInstance()
  12. {
  13. return;
  14. }
  15. // initialization
  16. BOOL COTSEDSSim::Init()
  17. {
  18. // simulation, never be wrong
  19. return TRUE;
  20. }
  21. // collect spectrum at the given position (to controller buffer)
  22. BOOL COTSEDSSim::CollectSpectrum(DWORD a_nMilliseconds, const CPoint& a_oPoint)
  23. {
  24. // read simulation spectrum into buffer
  25. memset(m_nRayData, 0, sizeof(DWORD) * (long)EDSConst::XANA_CHANNELS);
  26. if (!CollectASpectrumFromTxtFile(m_nRayData, (DWORD)EDSConst::XANA_CHANNELS))
  27. {
  28. // failed to call CollectASpectrumFromTxtFile method
  29. LogTrace(__FILE__, __LINE__, _T("COTSEDSSim::CollectSpectrum: failed to call CollectASpectrumFromTxtFile method."));
  30. }
  31. // always return TRUE
  32. return TRUE;
  33. }
  34. // collects spectrum (to controller buffer)
  35. BOOL COTSEDSSim::CollectSpectrum(DWORD a_nMilliseconds)
  36. {
  37. // read simulation spectrum into buffer
  38. memset(m_nRayData, 0, sizeof(DWORD) * (long)EDSConst::XANA_CHANNELS);
  39. if (!CollectASpectrumFromTxtFile(m_nRayData, (DWORD)EDSConst::XANA_CHANNELS))
  40. {
  41. // failed to call CollectASpectrumFromTxtFile method
  42. LogTrace(__FILE__, __LINE__, _T("COTSEDSSim::CollectSpectrum: failed to call CollectASpectrumFromTxtFile method."));
  43. }
  44. // always return TRUE
  45. return TRUE;
  46. }
  47. // collects spectrum (to given buffer)
  48. BOOL COTSEDSSim::CollectSpectrum(DWORD a_nMilliseconds, long* a_pCounts, DWORD a_nBufferSize)
  49. {
  50. // input check
  51. ASSERT(a_pCounts);
  52. if (!a_pCounts)
  53. {
  54. // invalid input data buffer
  55. LogTrace(__FILE__, __LINE__, _T("COTSEDSSim::CollectSpectrum: invalid input data buffer."));
  56. return FALSE;
  57. }
  58. // read simulation spectrum into buffer
  59. memset(a_pCounts, 0, sizeof(long) * (long)a_nBufferSize);
  60. if (!CollectASpectrumFromTxtFile((DWORD*)a_pCounts, a_nBufferSize))
  61. {
  62. // failed to call CollectASpectrumFromTxtFile method
  63. LogTrace(__FILE__, __LINE__, _T("COTSEDSSim::CollectSpectrum: failed to call CollectASpectrumFromTxtFile method."));
  64. }
  65. // always return TRUE
  66. return TRUE;
  67. }
  68. // get live time
  69. float COTSEDSSim::GetLiveTime()
  70. {
  71. // is simulation?
  72. return (float)0.1;
  73. }
  74. BOOL COTSEDSSim::IsSupportQuantification()
  75. {
  76. return TRUE;
  77. }
  78. BOOL COTSEDSSim::GetQuantificationMethods(std::vector<CString>& a_vMethods)
  79. {
  80. return TRUE;
  81. }
  82. BOOL COTSEDSSim::QuantifyXrayPoint(CPosXray* a_pXRayPoint, LPCTSTR a_sMethodName)
  83. {
  84. return TRUE;
  85. }
  86. BOOL COTSEDSSim::QuantifySpectrumFile(LPCTSTR a_sFilePathName, LPCTSTR a_sMethodName, CElementChemistriesList& a_listElementChemistry)
  87. {
  88. return TRUE;
  89. }
  90. BOOL COTSEDSSim::QuantifySpectrumOut(DWORD a_nMilliseconds, long* a_pCounts, DWORD a_nBufferSize, CElementChemistriesList& a_listElementChemistry)
  91. {
  92. return TRUE;
  93. }
  94. BOOL COTSEDSSim::GetXRayByPoints(std::vector<CPosXrayPtr>& a_vXRayPoints, const DWORD a_nXRayAQTime)
  95. {
  96. // read simulation spectrum into buffer
  97. memset(m_nRayData, 0, sizeof(DWORD) * (long)EDSConst::XANA_CHANNELS);
  98. if (!CollectASpectrumFromTxtFile(m_nRayData, (DWORD)EDSConst::XANA_CHANNELS))
  99. {
  100. // failed to call CollectASpectrumFromTxtFile method
  101. LogTrace(__FILE__, __LINE__, _T("COTSEDSSim::GetXRayByPoints: failed to call CollectASpectrumFromTxtFile method."));
  102. }
  103. for (auto pXrayPoi : a_vXRayPoints)
  104. {
  105. pXrayPoi->SetXrayData(m_nRayData);
  106. if (m_bDoQuantification)
  107. {
  108. CElementChemistriesList listElementQuantifyData;
  109. //CString sResult = _T("Quant=O,K-Serie,28.5\nQuant=Al,L-Serie,26.4\nQuant=S,L-Serie,1.6\nQuant=Ca,L-Serie,15.8\nQuant=Fe,L-Serie,9.5");
  110. CString sResult = _T("Quant=O,K-Serie,28.5\nQuant=S,L-Serie,10\nQuant=Mn,L-Serie,20\nQuant=Fe,L-Serie,9.5");
  111. CElementPtr pElement = CElementPtr(new CElement());
  112. listElementQuantifyData = pElement->ExtractElementChemistrys(sResult);
  113. pXrayPoi->SetElementQuantifyData(listElementQuantifyData);
  114. }
  115. }
  116. Sleep(500);
  117. // always return TRUE
  118. return TRUE;
  119. }
  120. BOOL COTSEDSSim::GetXRayByFeatures(std::vector<CPosXrayPtr>& a_listXRayPoints,
  121. std::vector<std::vector<BrukerSegment>>& a_listFeatures,
  122. const DWORD a_nXRayAQTime)
  123. {
  124. // read simulation spectrum into buffer
  125. memset(m_nRayData, 0, sizeof(DWORD) * (long)EDSConst::XANA_CHANNELS);
  126. if (!CollectASpectrumFromTxtFile(m_nRayData, (DWORD)EDSConst::XANA_CHANNELS))
  127. {
  128. // failed to call CollectASpectrumFromTxtFile method
  129. LogTrace(__FILE__, __LINE__, _T("COTSEDSSim::GetXRayByFeatures: failed to call CollectASpectrumFromTxtFile method."));
  130. }
  131. for (auto pXrayPoi : a_listXRayPoints)
  132. {
  133. pXrayPoi->SetXrayData(m_nRayData);
  134. if (m_bDoQuantification)
  135. {
  136. CElementChemistriesList listElementQuantifyData;
  137. //CString sResult = _T("Quant=O,K-Serie,28.5\nQuant=Al,L-Serie,26.4\nQuant=S,L-Serie,1.6\nQuant=Ca,L-Serie,15.8\nQuant=Fe,L-Serie,9.5");
  138. CString sResult = _T("Quant=O,K-Serie,28.5\nQuant=Fe,L-Serie,49.5");
  139. CElementPtr pElement = CElementPtr(new CElement());
  140. listElementQuantifyData = pElement->ExtractElementChemistrys(sResult);
  141. pXrayPoi->SetElementQuantifyData(listElementQuantifyData);
  142. }
  143. }
  144. // always return TRUE
  145. return TRUE;
  146. }
  147. //BOOL COTSEDSSim::GetXRayByFeaturesFromMultiPoint(std::vector<CPosXrayPtr>& a_listXRayPoints,
  148. // std::vector<std::vector<BrukerSegment>>& a_listFeatures,
  149. // const DWORD a_nXRayAQTime)
  150. //{
  151. // int nPointNum = 0;
  152. // for (auto pFeature: a_listFeatures)
  153. // {
  154. // std::vector<CPosXrayPtr> listXRayPoints;
  155. // std::vector<BrukerSegment> listSegment = pFeature;
  156. // int nSize = listSegment.size();
  157. // for (int i = 0; i < nSize -2; i= i+2)
  158. // {
  159. // BrukerSegment pSegment = listSegment[i];
  160. // int nLength = pSegment.XCount;
  161. // int nStart = pSegment.XStart;
  162. // CPoint pt;
  163. // pt.x = nStart;
  164. // pt.y = nStart + nLength / 2 ;
  165. // CPosXrayPtr pXray = CPosXrayPtr(new CPosXray());
  166. // pXray->SetPosition(pt);
  167. // listXRayPoints.push_back(pXray);
  168. // }
  169. // if (!GetXRayByPoints(listXRayPoints, a_nXRayAQTime))
  170. // {
  171. // LogErrorTrace(__FILE__, __LINE__, _T("GetXRayByFeaturesFromMultiPoint: Can't get xray from points."));
  172. // return FALSE;
  173. // }
  174. // int nNum = 0;
  175. // DWORD XRayData[GENERALXRAYCHANNELS];
  176. // memset(XRayData, 0, sizeof(DWORD) *GENERALXRAYCHANNELS);
  177. // for (auto pXRayPoint : listXRayPoints)
  178. // {
  179. //
  180. // DWORD* pXRayData;
  181. // pXRayData = pXRayPoint->GetXrayData();
  182. //
  183. // for (int m = 0; m < GENERALXRAYCHANNELS; m++)
  184. // {
  185. // XRayData[m] += pXRayData[m];
  186. // }
  187. // nNum++;
  188. // }
  189. // if (nNum == 0)
  190. // {
  191. // LogErrorTrace(__FILE__, __LINE__, _T("GetXRayByFeaturesFromMultiPoint:there is points."));
  192. // return FALSE;
  193. // }
  194. // for (int k = 0; k < GENERALXRAYCHANNELS; k++)
  195. // {
  196. // XRayData[k] = (int)((double)XRayData[k]/(double)nNum + 0.5);
  197. // }
  198. // CPosXrayPtr pXrayPoint = CPosXrayPtr(new CPosXray());
  199. // a_listXRayPoints[nPointNum]->SetXrayData(XRayData);
  200. // nPointNum++;
  201. // }
  202. //
  203. // // always return TRUE
  204. // return TRUE;
  205. //}
  206. void COTSEDSSim::SetQuantification(BOOL a_bQuantification)
  207. {
  208. m_bDoQuantification = a_bQuantification;
  209. return;
  210. }
  211. BOOL COTSEDSSim::GetQuantification()
  212. {
  213. return TRUE;
  214. }
  215. // Get number of channels
  216. DWORD COTSEDSSim::GetNumberOfChannels(void)
  217. {
  218. return (DWORD)2000;
  219. }
  220. // Get the x-Ray data
  221. DWORD* COTSEDSSim::GetXRayData()
  222. {
  223. return m_nRayData;
  224. }
  225. }