X2D1Task.cs 14 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284
  1. using System;
  2. using System.Collections.Generic;
  3. using System.Globalization;
  4. using System.IO;
  5. using System.Linq;
  6. using System.Threading;
  7. using System.Threading.Tasks;
  8. using XdCxRhDW.Dto;
  9. namespace XdCxRhDW.TaskServer.Task
  10. {
  11. public class X2D1Task : HistoryTaskI
  12. {
  13. public override void Start(HistoryTaskProcessingDto dto)
  14. {
  15. IsRuning = true;
  16. TaskDto = dto;
  17. System.Threading.Tasks.Task.Run(async () =>
  18. {
  19. DateTime preTime = dto.StartTime;
  20. while (IsRuning || dto.EndTime.HasValue && dto.EndTime <= DateTime.Now)
  21. {
  22. string capfile = Path.Combine(dto.CapDir, $"{preTime.ToString(dto.DateDirFormat)}");
  23. IEnumerable<string> files;
  24. try
  25. {
  26. files = Directory.EnumerateFiles(capfile, "*.dat");
  27. }
  28. catch (Exception ex)
  29. {
  30. LogHelper.Info($"查询历史任务:{dto.ID}采集数据异常:{ex.Message}");
  31. Thread.Sleep(1000);
  32. continue;
  33. }
  34. var groups = files.Select(f => FileToHistoryFile(f)).GroupBy(m => m.CapTime).OrderBy(m => m.Key);
  35. foreach (var item in groups)
  36. {
  37. if (!IsRuning)
  38. {
  39. break;
  40. }
  41. var finfos = item.ToList();
  42. if (finfos.Count < 3)
  43. {
  44. LogHelper.Info($"历史任务:{dto.ID}采集时刻:{finfos.First().CapTime:yyyy-MM-dd HH:mm:ss}文件数量只有{finfos.Count}个,不满足");
  45. continue;
  46. }
  47. //超短波信号
  48. var dinfo = finfos.FirstOrDefault(m => m.Ch == 1);
  49. if (dinfo == null)
  50. {
  51. LogHelper.Info($"历史任务:{dto.ID}采集时刻:{finfos.First().CapTime:yyyy-MM-dd HH:mm:ss}未找到超短波信号ch1文件,不满足");
  52. continue;
  53. }
  54. //主星
  55. var minfo = finfos.FirstOrDefault(m => m.Ch == 2);
  56. if (minfo == null)
  57. {
  58. LogHelper.Info($"历史任务:{dto.ID}采集时刻:{finfos.First().CapTime:yyyy-MM-dd HH:mm:ss}未找到主星信号ch2文件,不满足");
  59. continue;
  60. }
  61. //邻1
  62. var ninfo = finfos.FirstOrDefault(m => m.Ch == 3);
  63. if (ninfo == null)
  64. {
  65. LogHelper.Info($"历史任务:{dto.ID}采集时刻:{finfos.First().CapTime:yyyy-MM-dd HH:mm:ss}未找到邻星信号ch3文件,不满足");
  66. continue;
  67. }
  68. try
  69. {
  70. var txInfos = await GetTxInfoAsync();
  71. var satTx = txInfos.First(t => t.TxType == EnumTxType.Rec);
  72. var cdbTx = txInfos.First(t => t.TxType == EnumTxType.Cdb);
  73. var refTx = txInfos.First(t => t.TxType == EnumTxType.Ref);
  74. //根据信号类型执行检测或参数估计
  75. if (dto.SigType == EnumSigTypeDto.Normal)
  76. {
  77. var xd = await GPUCalcAsync(minfo.FilePath, dinfo.FilePath, minfo.FsHz, 260000, 40000);
  78. var sx = await GPUCalcAsync(minfo.FilePath, ninfo.FilePath, minfo.FsHz, 0, 40000);
  79. X2D1NoXlPosDto x2D1 = new X2D1NoXlPosDto()
  80. {
  81. SigTime = minfo.CapTime,
  82. MainCode = minfo.SatId,
  83. AdjaCode = ninfo.SatId,
  84. SxDto = sx.Dt,
  85. XdDto = xd.Dt,
  86. SatTxLon = satTx.Lon,
  87. SatTxLat = satTx.Lat,
  88. CdbTxLon = cdbTx.Lon,
  89. CdbTxLat = cdbTx.Lat,
  90. RefLon = refTx.Lon,
  91. RefLat = refTx.Lat,
  92. FreqDown = minfo.FreqHz * 1e-6,
  93. FreqUp = dinfo.FreqHz * 1e-6,
  94. XdDfo = xd.Df,
  95. XdSnr = xd.Snr,
  96. SxDfo = sx.Df,
  97. SxSnr = sx.Snr,
  98. };
  99. await X2D1NoXlAsync(x2D1);
  100. }
  101. else
  102. {
  103. var mDetect = await DAMAAsync(dto.TaskCheckType, minfo.FsHz, minfo.FilePath);
  104. var dfile = await ToResampleAsync((int)dinfo.FsHz, dinfo.FilePath);
  105. var nfile = await ToResampleAsync((int)ninfo.FsHz, ninfo.FilePath);
  106. if (mDetect.Count() <= 0)
  107. {
  108. LogHelper.Info($"历史任务:{dto.ID}文件{minfo.FilePath}未检测出结果");
  109. continue;
  110. }
  111. //主星文件
  112. string mfile = mDetect.First().File1;
  113. //时隙数
  114. var smps = mDetect.Select(m => new SmpPosition() { smpStart = m.Start, smpCount = m.Length }).ToList();
  115. int dtRange = 40000;
  116. //提取时隙有效数据
  117. //var mdfile = await ExtractMergeAsync(mfile, nfile.OutFsHz, smps, 0, dtRange, minfo.Ch);
  118. //var ddfile = await ExtractMergeAsync(dfile.File, nfile.OutFsHz, smps, 260000, dtRange, dinfo.Ch);
  119. //var ndfile = await ExtractMergeAsync(nfile.File, nfile.OutFsHz, smps, 0, dtRange, ninfo.Ch);
  120. //var xds = await CPUCalcAsync(mdfile, ddfile, dfile.OutFsHz, smps, 260000, dtRange);
  121. //var sxs = await CPUCalcAsync(mdfile, ndfile, nfile.OutFsHz, smps, 0, dtRange);
  122. var xds = await CPUCalcAsync(mfile, dfile.File, dfile.OutFsHz, smps, 260000, dtRange);
  123. var sxs = await CPUCalcAsync(mfile, nfile.File, nfile.OutFsHz, smps, 0, dtRange);
  124. foreach (var sxitem in sxs)
  125. {
  126. var xd = xds.FirstOrDefault(x => x.Smpstart == sxitem.Smpstart && x.Smplen == sxitem.Smplen);
  127. var deitem = mDetect.FirstOrDefault(x => x.Start == sxitem.Smpstart && x.Length == sxitem.Smplen);
  128. if (!IsRuning)
  129. {
  130. break;
  131. }
  132. try
  133. {
  134. X2D1NoXlPosDto x2D1 = new X2D1NoXlPosDto()
  135. {
  136. TaskID = dto.ID,
  137. SigTime = minfo.CapTime.AddSeconds(sxitem.Smpstart / nfile.OutFsHz),
  138. MainCode = minfo.SatId,
  139. AdjaCode = ninfo.SatId,
  140. SxDto = sxitem.Dt,
  141. XdDto = xd == null ? 0 : xd.Dt,
  142. SatTxLon = satTx.Lon,
  143. SatTxLat = satTx.Lat,
  144. CdbTxLon = cdbTx.Lon,
  145. CdbTxLat = cdbTx.Lat,
  146. RefLon = refTx.Lon,
  147. RefLat = refTx.Lat,
  148. FreqDown = minfo.FreqHz * 1e-6,
  149. FreqUp = dinfo.FreqHz * 1e-6,
  150. XdDfo = xd == null ? 0 : xd.Df,
  151. XdSnr = xd == null ? 0 : xd.Snr,
  152. SxDfo = sxitem.Df,
  153. SxSnr = sxitem.Snr,
  154. CheckRes = new CheckResDto()
  155. {
  156. FileName = deitem.File1,
  157. ModRate = deitem.ModRate,
  158. ModType = deitem.ModType,
  159. SmpCount = deitem.Length,
  160. SmpStart = deitem.Start,
  161. UserName = deitem.UserName,
  162. PosCheckType = deitem.DmcType.GetEnumByDisplayName<EnumPosCheckTypeDto>(),
  163. }
  164. };
  165. await X2D1NoXlAsync(x2D1);
  166. }
  167. catch (Exception ex)
  168. {
  169. LogHelper.Info($"历史任务:{dto.ID}异常:{ex.Message}");
  170. continue;
  171. }
  172. }
  173. }
  174. }
  175. catch (Exception ex)
  176. {
  177. LogHelper.Info($"执行历史任务:{dto.ID}异常:{ex.Message}");
  178. continue;
  179. }
  180. }
  181. preTime = preTime.AddHours(1);
  182. }
  183. });
  184. }
  185. private async Task<bool> X2D1NoXlAsync(X2D1NoXlPosDto x2D1)
  186. {
  187. var result = await HttpHelper.PostRequestAsync<PosResDto>(baseUrl + "Pos/PosX2D1NoXlAsync", x2D1);
  188. if (result.code != 200)
  189. {
  190. LogHelper.Error($"推送历史任务定位结果:{x2D1.TaskID}异常:{result.msg}");
  191. return false;
  192. }
  193. return true;
  194. }
  195. private async Task<List<TxInfo>> GetTxInfoAsync()
  196. {
  197. List<TxInfo> infos = new List<TxInfo>();
  198. var result = await HttpHelper.GetRequestAsync<(TxInfo recTx, TxInfo cdbTx, TxInfo cxTx, TxInfo refLoc)>(baseUrl + "Tx/GetTxInfoAsync");
  199. if (result.code != 200)
  200. {
  201. infos.Add(new TxInfo() { TxType = EnumTxType.Rec });
  202. infos.Add(new TxInfo() { TxType = EnumTxType.Cdb });
  203. infos.Add(new TxInfo() { TxType = EnumTxType.Cx });
  204. infos.Add(new TxInfo() { TxType = EnumTxType.Ref });
  205. LogHelper.Error($"历史任务获取天线信息失败");
  206. return infos;
  207. }
  208. infos.Add(result.data.recTx == null ? new TxInfo() { TxType = EnumTxType.Rec } : result.data.recTx);
  209. infos.Add(result.data.cdbTx == null ? new TxInfo() { TxType = EnumTxType.Cdb } : result.data.cdbTx);
  210. infos.Add(result.data.cxTx == null ? new TxInfo() { TxType = EnumTxType.Cx } : result.data.cxTx);
  211. infos.Add(result.data.refLoc == null ? new TxInfo() { TxType = EnumTxType.Ref } : result.data.refLoc);
  212. return infos;
  213. }
  214. /// <summary>
  215. /// 获取采集文件信息
  216. /// </summary>
  217. /// <param name="filePath"></param>
  218. /// <returns></returns>
  219. private HistoryFile FileToHistoryFile(string filePath)
  220. { //读取采集文件
  221. //2024_01_31_10_01_51_000000000_ch11_-1__Nxx.xxxxxx_Exx.xxxxxx_xxxxx.xxxHz_xxx.xxxMHz_ch1_xd1.dat
  222. HistoryFile historyFile = new HistoryFile();
  223. historyFile.FilePath = filePath;
  224. var fileName = Path.GetFileNameWithoutExtension(filePath);
  225. var strs = fileName.Split(new string[] { "_", "MHz", "Hz", "ch" }, StringSplitOptions.RemoveEmptyEntries);
  226. //采集时间
  227. var datestr = string.Join("_", strs.Take(6));
  228. DateTime dateTime;
  229. bool istime = DateTime.TryParseExact(datestr, "yyyy_MM_dd_HH_mm_ss", null, DateTimeStyles.None, out dateTime);
  230. if (!istime) return historyFile;
  231. historyFile.CapTime = dateTime;
  232. //采集通道
  233. var chstr = strs.Skip(strs.Length - 2).Take(1).First();
  234. int ch;
  235. int.TryParse(chstr, out ch);
  236. //采集频点
  237. var freqstr = strs.Skip(strs.Length - 3).Take(1).First();
  238. double freqMHz;
  239. double.TryParse(freqstr, out freqMHz);
  240. //采样率
  241. var fsstr = strs.Skip(strs.Length - 4).Take(1).First();
  242. double fsHz;
  243. double.TryParse(fsstr, out fsHz);
  244. historyFile.FreqHz = freqMHz * 1e6;
  245. historyFile.Ch = ch;
  246. historyFile.FsHz = fsHz;
  247. if (ch > 1)
  248. {
  249. historyFile.SatId = GetSatId(historyFile.FreqHz);
  250. }
  251. return historyFile;
  252. }
  253. }
  254. }