XcorrUtils.cs 11 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246
  1. using System;
  2. using System.Diagnostics;
  3. using System.IO;
  4. using System.Linq;
  5. using System.Reflection;
  6. using System.Threading.Tasks;
  7. using System.Collections;
  8. using System.Collections.Generic;
  9. using XdCxRhDW.Dto;
  10. using System.Security.Policy;
  11. namespace XdCxRhDW.Api
  12. {
  13. public class XcorrStruct
  14. {
  15. public String file1 { get; set; }
  16. public String file2 { get; set; }
  17. public long smpStart { get; set; } //开始样点
  18. public double smpCount { get; set; } //样点数
  19. public double samplingRate { get; set; } //采样率
  20. public double dtCenter { get; set; } //时差中心
  21. public double dtRange { get; set; } //时差范围
  22. public double dfRange { get; set; } //频差范围
  23. public double snrThreshold { get; set; } //信噪比门限
  24. public XcorrStruct Copy()
  25. {
  26. XcorrStruct xItem = new XcorrStruct();
  27. xItem.file1 = file1;
  28. xItem.file2 = file2;
  29. xItem.smpCount = smpCount;
  30. xItem.samplingRate = samplingRate;
  31. xItem.dtCenter = dtCenter;
  32. xItem.dtRange = dtRange;
  33. xItem.dfRange = dfRange;
  34. xItem.smpStart = smpStart;
  35. xItem.snrThreshold = snrThreshold;
  36. return xItem;
  37. }
  38. }
  39. public class XcorrUtils
  40. {
  41. Process cafp;
  42. Process dmp;
  43. // D:/data/test/r1.dat D:/data/test/r4.dat 1562500 -250000 20 16384 100000 1048576 11
  44. public async Task<CafResult> Calc(XcorrStruct xs)
  45. {
  46. CafResult res = new CafResult();
  47. cafp = new Process();
  48. await Task.Run(() =>
  49. {
  50. int flag = 0;
  51. if (xs.smpStart == 0)
  52. {
  53. flag = 2;//文件参估,避免滑动时溢出
  54. cafp.StartInfo.Arguments = $"{flag} \"{xs.file1}\" \"{xs.file2}\" {(Int64)(xs.samplingRate)} {xs.dtCenter} {xs.dtRange} {xs.dfRange} {xs.smpCount} {xs.snrThreshold}";
  55. }
  56. else
  57. {
  58. var minStart = xs.dtRange / 2 / 1e6 * xs.samplingRate;
  59. if (xs.smpStart < minStart)//滑动计算起时值必须要大于等于此数
  60. {
  61. xs.smpStart = (long)minStart;
  62. }
  63. cafp.StartInfo.Arguments = $"{flag} \"{xs.file1}\" \"{xs.file2}\" {(Int64)(xs.samplingRate)} {xs.dtCenter} {xs.dtRange} {xs.dfRange} {xs.smpStart} {xs.smpCount} {xs.snrThreshold}";
  64. }
  65. var exePath= Path.Combine(cafp.StartInfo.WorkingDirectory, "AddIns\\检测Cpu参估");
  66. cafp.StartInfo.WorkingDirectory = exePath;
  67. cafp.StartInfo.FileName = Path.Combine(exePath,"XcorrCpu.exe");
  68. cafp.StartInfo.CreateNoWindow = true;
  69. cafp.StartInfo.RedirectStandardError = true;
  70. cafp.StartInfo.RedirectStandardOutput = true;
  71. cafp.StartInfo.UseShellExecute = false;
  72. cafp.Start();
  73. Stopwatch stopWatch = new Stopwatch();
  74. stopWatch.Start();
  75. cafp.WaitForExit();
  76. stopWatch.Stop();
  77. TimeSpan ts = stopWatch.Elapsed;
  78. var str = cafp.StandardOutput.ReadToEnd();
  79. res.FromLine(str);
  80. res.file1 = xs.file1;
  81. res.file2 = xs.file2;
  82. res.tm = (int)stopWatch.Elapsed.TotalMilliseconds;
  83. res.smpstart = xs.smpStart;
  84. res.smplen = (long)xs.smpCount;
  85. });
  86. return res;
  87. }
  88. /// <summary>
  89. /// 参估画图
  90. /// </summary>
  91. /// <returns></returns>
  92. public async Task<List<ImageResult>> DrawImage(XcorrStruct xs)
  93. {
  94. cafp = new Process();
  95. List<ImageResult> list = new List<ImageResult>();
  96. await Task.Run(() =>
  97. {
  98. FileInfo file1 = new FileInfo(xs.file1);
  99. FileInfo file2 = new FileInfo(xs.file2);
  100. long totalsamp = file1.Length < file2.Length ? file1.Length / 4 : file2.Length / 4;
  101. //样点数为0时计算所有样本
  102. if (xs.smpCount == 0)
  103. {
  104. xs.smpCount = (int)totalsamp - xs.smpStart;
  105. }
  106. else if (xs.smpCount > 0 && xs.smpCount < 1)
  107. {
  108. xs.smpCount = (int)(totalsamp * xs.smpCount);
  109. }
  110. string outFile = Path.Combine(AppDomain.CurrentDomain.BaseDirectory, "CgImag.txt");
  111. int flag = 1;
  112. if (xs.smpStart == 0)
  113. {
  114. flag = 3;//文件画图,避免滑动时溢出
  115. cafp.StartInfo.Arguments = $"{flag} \"{xs.file1}\" \"{xs.file2}\" {(long)xs.samplingRate} {xs.dtCenter} {xs.dtRange} {xs.dfRange} {xs.smpCount} {xs.snrThreshold} {outFile}";
  116. }
  117. else
  118. {
  119. var minStart = xs.dtRange / 2 / 1e6 * xs.samplingRate;
  120. if (xs.smpStart < minStart)//滑动计算起时值必须要大于等于此数
  121. {
  122. xs.smpStart = (long)minStart;
  123. }
  124. cafp.StartInfo.Arguments = $"{flag} \"{xs.file1}\" \"{xs.file2}\" {(long)(xs.samplingRate)} {xs.dtCenter} {xs.dtRange} {xs.dfRange} {xs.smpStart} {xs.smpCount} {xs.snrThreshold} {outFile}";
  125. }
  126. cafp.StartInfo.FileName = Path.Combine(cafp.StartInfo.WorkingDirectory, "AddIns\\检测Cpu参估\\XcorrCpu.exe");
  127. //cafp.StartInfo.Arguments = $"1 \"{xs.file1}\" \"{xs.file2}\" {(long)(xs.samplingRate)} {xs.dtCenter} {xs.dtRange} {xs.dfRange} {xs.smpStart} {xs.smpCount} {xs.snrThreshold} {outFile}";
  128. cafp.StartInfo.CreateNoWindow = true;
  129. cafp.StartInfo.RedirectStandardError = true;
  130. cafp.StartInfo.RedirectStandardOutput = true;
  131. cafp.StartInfo.UseShellExecute = false;
  132. cafp.Start();
  133. Stopwatch stopWatch = new Stopwatch();
  134. stopWatch.Start();
  135. cafp.WaitForExit();
  136. stopWatch.Stop();
  137. TimeSpan ts = stopWatch.Elapsed;
  138. var str = cafp.StandardOutput.ReadToEnd();
  139. if (str.StartsWith("1:") && File.Exists(outFile))
  140. {
  141. var lines = File.ReadAllLines(outFile);
  142. if (lines.Length > 1)
  143. {
  144. var lineFirstArr = lines[0].Split(",".ToCharArray(), StringSplitOptions.RemoveEmptyEntries);
  145. int fsHz = Convert.ToInt32(lineFirstArr[0]);
  146. int xFlag = Convert.ToInt32(lineFirstArr[1]);
  147. int xMax = Convert.ToInt32(lineFirstArr[2]);
  148. float yMax = Convert.ToSingle(lineFirstArr[3]);
  149. for (int i = 1; i < lines.Length; i++)
  150. {
  151. var lineArr = lines[i].Split(",".ToCharArray(), StringSplitOptions.RemoveEmptyEntries);
  152. ImageResult ir = new ImageResult()
  153. {
  154. FsHz = fsHz,
  155. XFlag = xFlag,
  156. XMax = xMax,
  157. YMax = yMax,
  158. XValue = Convert.ToInt32(lineArr[0]),
  159. YValue = Convert.ToSingle(lineArr[1]),
  160. Snr = Convert.ToSingle(lineArr[2]),
  161. };
  162. list.Add(ir);
  163. }
  164. }
  165. File.Delete(outFile);
  166. }
  167. });
  168. return list;
  169. }
  170. public async Task<IEnumerable<DmcResult>> DmcCheckAsync(string fileName, double fsHz, DmcType dmcType, double? bandKHz = null)
  171. {
  172. if (bandKHz == null || bandKHz.Value <= 0) bandKHz = 25;
  173. string dmcCmd = "all";
  174. dmp = new Process();
  175. string pFileName = Path.Combine(dmp.StartInfo.WorkingDirectory, "AddIns\\检测Cpu参估\\dmc.exe");
  176. string pArguments = string.Empty;
  177. switch (dmcType)
  178. {
  179. case DmcType.DAMA:
  180. dmcCmd = "dm";
  181. if (fsHz != 96000)
  182. throw new Exception("DAMA只支持96K采样率");
  183. pArguments = $"{dmcCmd} \"{fileName}\" -c true";//-c包含ccow
  184. break;
  185. case DmcType.IBS:
  186. dmcCmd = "nd";
  187. pArguments = $"{dmcCmd} \"{fileName}\" -w {bandKHz}";// -f {fs}" -c --dmfirst";
  188. break;
  189. case DmcType.Ky5758:
  190. pFileName = Path.Combine(dmp.StartInfo.WorkingDirectory, "AddIns\\检测Cpu参估\\enc.exe");
  191. //enc.exe enc-test.dat 0.096 5 0全部文件
  192. pArguments = $"{fileName} {fsHz * 1e-6} {5} {0}";
  193. break;
  194. default:
  195. break;
  196. }
  197. return await Task.Run(() =>
  198. {
  199. dmp.StartInfo.FileName = pFileName;
  200. dmp.StartInfo.Arguments = pArguments;
  201. dmp.StartInfo.CreateNoWindow = true;
  202. dmp.StartInfo.RedirectStandardError = true;
  203. dmp.StartInfo.RedirectStandardOutput = true;
  204. dmp.StartInfo.UseShellExecute = false;
  205. dmp.Start();
  206. Stopwatch stopWatch = new Stopwatch();
  207. stopWatch.Start();
  208. dmp.WaitForExit();
  209. stopWatch.Stop();
  210. TimeSpan ts = stopWatch.Elapsed;
  211. var str = dmp.StandardOutput.ReadToEnd();
  212. return ConvertDmcResult(dmcType, str, ts.TotalMilliseconds);
  213. });
  214. }
  215. private IEnumerable<DmcResult> ConvertDmcResult(DmcType type, string res, double tm)
  216. {
  217. var lines = res.Split(Environment.NewLine.ToArray(), StringSplitOptions.RemoveEmptyEntries);
  218. foreach (var line in lines)
  219. {
  220. var items = line.Split('\t');
  221. if (items.Length < 2) continue;
  222. int start = int.Parse(items[0]);
  223. int length = int.Parse(items[1]);
  224. string userName = "";
  225. if (items.Length >= 3)
  226. userName = items[2];
  227. var item = new DmcResult(start, length, userName, (int)tm);
  228. if (type == DmcType.DAMA)
  229. item.ModType = "BPSK";
  230. else
  231. item.ModType = string.Empty;
  232. yield return item;
  233. }
  234. }
  235. }
  236. }