XcorrUtils.cs 11 KB

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