汉泰.doc=>力科.trc转换程序

含光左卫

第一级代码(by Python):

1. # code by fengxh.
   
2. # transfer 汉泰.doc => Lecroy 二进制.trc格式。
   
3. # python code is created at Jun30,2023
   
4. 
   
5. from scipy.optimize import curve_fit
   
6. import numpy as np
   
7. import csv
   
8. import sys
   
9. import operator
   
10. import subprocess
    
11. import re
    
12. 
    
13. 
    
14. # 取csv文件的第N列，抛弃掉标题栏，返回数组
    
15. def load_csv_data(filename, idxBase0):
    
16.     data = []  # 存放结果的数组
    
17. 
    
18. # 打开CSV文件，并读取数据
    
19.     flagDataStart = False;
    
20.     with open(filename, 'r') as f:
    
21.         reader = csv.reader(f)
    
22.         i = 0;
    
23.         for row in reader:
    
24.             i = i+1;
    
25.             #print(i);
    
26.             if len(row) == 0:
    
27.                 if(not flagDataStart):
    
28.                     flagDataStart = True;
    
29.                     continue;
    
30.                 else:
    
31.                     break; #done!
    
32.             if(not flagDataStart): continue;
    
33.             try:
    
34.                data.append(float(row[idxBase0]))  # 取第一列数据，并转为float类型
    
35.             except:
    
36.                 continue;
    
37.     return data;
    
38. 
    
39. # 取csv文件的第N列，抛弃掉标题栏，返回数组
    
40. def get_time_scale(filename):
    
41.     flagDataStart = False;
    
42.     #string = "#CLOCK=0.004000000S"
    
43.     pattern = r"CLOCK=(\d+\.\d+)S"
    
44.     with open(filename, 'r') as f:
    
45.         reader = csv.reader(f)
    
46.         for row in reader:
    
47.             match = re.search(pattern, str(row))
    
48.             if match:
    
49.                 result = float(match.group(1));
    
50.                 return result;
    
51.     return 100e-6;
    
52. 
    
53. 
    
54. def getRmsOfFreq(arFreq, arRms, freq):
    
55.     for i in range(0, len(arFreq)):
    
56.         if(arFreq[i] == freq): return arRms[i];
    
57. 
    
58. def csv2trc(csvfile, trcfile, ratio):
    
59. # 要执行的命令
    
60. #Usage: Csv2Trc <src.csv> <tgt.trc> [<sample_rate=100e-6> <volPerAd>=1]   
    
61.     command = r"C:\Users\fengxh\Desktop\工具集\dumb\Csv2Trc_V1.0.20230531\SampleData_CSV2LecoryTrc";
    
62.     command = command + " " +csvfile + ' ' + trcfile + ' ' + str(ratio);
    
63.     # 执行命令并获取输出
    
64.     output = subprocess.check_output(command, shell=True)
    
65.     output = output.decode()  # 如果你想要将输出作为字符串处理，可以解码为返回的字节字符串
    
66.     print(output)
    
67. 
    
68. def toInteger(rms):
    
69.     vals = [];
    
70.     for i in rms:
    
71.         if(not i in vals):
    
72.             vals.append(i);
    
73. 
    
74.     sorted_vals = sorted(vals);
    
75. 
    
76.     my_dict = {sorted_vals[0]:1};
    
77.     idx = 1;
    
78.     for v in sorted_vals:
    
79.         my_dict[v] = idx;
    
80.         idx = idx +1;
    
81. 
    
82.     ret = [];
    
83.     for i in rms:
    
84.         ret.append(my_dict[i]);
    
85.             
    
86.     return ret;
    
87. 
    
88. 
    
89. def zoom_out_lecory_wave_data(filename, ratio=10):
    
90.     rms = load_csv_data(filename, 0);  #rms
    
91.     dt = get_time_scale(filename);
    
92.     print("数据的时间跨度", dt, "秒, 数据点数:", len(rms));
    
93.     rms = toInteger(rms);
    
94.     data = rms[:];
    
95.     for i in range(0, len(rms)):
    
96.         data[i//round(ratio) + round(i%(round(ratio))*len(rms)/ratio)] = rms[i];
    
97. 
    
98.     data_2d = [[element] for element in data]
    
99.     filename = filename + '.zoom_out' + str(ratio) + ".csv";
    
100.     with open(filename, 'w', newline='') as csvfile:
     
101.         writer = csv.writer(csvfile)
     
102.         writer.writerows(data_2d)
     
103.     csv2trc(filename, filename+".trc", dt*ratio/len(rms));
     
104.     print("done!");
     
105. 
     
106. #调用示例
     
107. filename = "35.25_40K.doc";
     
108. zoom_out_lecory_wave_data(filename,1);
     

Copy the Code

第二级代码，SampleData_CSV2LecoryTrc，writened in C#, build by VS2022

1. /*
   
2. *  lecroy_trc.cs
   
3. *  
   
4. *  A csv to lecory .trc bin file format helper modules.
   
5. *  
   
6. *  thanks to Lecroy, for its powerful offline wavefrom analyse tools.
   
7. *  and Surrea1@CSDN for provide the .trc file format.
   
8. *  
   
9. *  ----
   
10. *  V1.0.20230530 by twicave
    
11. *          created. the first verison, coding in vs2022(c#)
    
12. *  
    
13. */
    
14. 
    
15. using System;
    
16. using System.Collections.Generic;
    
17. using System.Data.SqlTypes;
    
18. using System.Diagnostics.SymbolStore;
    
19. using System.Linq;
    
20. using System.Reflection.PortableExecutable;
    
21. using System.Security.Cryptography;
    
22. using System.Text;
    
23. using System.Threading;
    
24. using System.Threading.Tasks;
    
25. using System.Runtime.InteropServices;
    
26. using static System.Net.Mime.MediaTypeNames;
    
27. using System.ComponentModel;
    
28. 
    
29. namespace SampleData_CSV2LecoryTrc
    
30. {
    
31. 
    
32.     internal class lecroy_trc
    
33.     {
    
34.         static byte[] lecroy_trc_header = new byte[]{
    
35.         0x23,0x39,0x30,0x30,0x32,0x30,0x30,0x30,0x31,0x34,0x34,0x57,0x41,0x56,0x45,0x44,
    
36.         0x45,0x53,0x43,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x4C,0x45,0x43,0x52,0x4F,
    
37.         0x59,0x5F,0x32,0x5F,0x33,0x00,0x00,0x00,0x00,0x00,0x00,0x01,0x00,0x01,0x00,0x5A,
    
38.         0x01,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
    
39.         0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xB6,0x83,0x1E,0x00,0x00,0x00,0x00,0x00,0x00,
    
40.         0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x4C,0x45,0x43,0x52,0x4F,0x59,0x57,0x53,0x34,
    
41.         0x32,0x34,0x00,0x00,0x00,0x00,0x00,0x4B,0x45,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
    
42.         0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xDB,0x41,0x0F,0x00,0xDB,
    
43.         0x41,0x0F,0x00,0x40,0x42,0x0F,0x00,0x00,0x00,0x00,0x00,0xDA,0x41,0x0F,0x00,0x00,
    
44.         0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x01,
    
45.         0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x0D,0x00,0x80,0x39,0x00,0x00,0x00,0xBF,0x00,
    
46.         0x00,0xFE,0x46,0x00,0x00,0x00,0xC7,0x0B,0x00,0x01,0x00,0xBD,0x37,0x06,0x35,0x1E,
    
47.         0xF4,0xC8,0xB6,0x68,0x78,0xD1,0xBF,0x46,0xB6,0xF3,0xFD,0xD4,0x78,0xD1,0xBF,0x56,
    
48.         0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
    
49.         0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
    
50.         0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x53,
    
51.         0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
    
52.         0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
    
53.         0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xFF,
    
54.         0xEB,0xAF,0x2C,0x99,0xBC,0x0B,0xD8,0x31,0x27,0x39,0x40,0x02,0x0F,0x1A,0x04,0xE7,
    
55.         0x07,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x01,0x00,0x20,
    
56.         0x00,0x02,0x00,0x00,0x00,0x80,0x3F,0x13,0x00,0x00,0x00,0x00,0x00,0x80,0x3F,0x00,
    
57.         0x00,0x00,0xBF,0x00,0x00 };
    
58. 
    
59.         public class trc_field_define
    
60.         {
    
61.             public string standardName;
    
62.             public string field_name;
    
63.             public int offetIdxbase0;
    
64.             public int len;
    
65.             public string type;
    
66.             public string memo;
    
67.             public string example;
    
68.             public trc_field_define(string standardName, string field_name, int offsetIdxbase0, int len, string type, string memo, string example)
    
69.             {
    
70.                 this.standardName = standardName;
    
71.                 this.field_name = field_name;
    
72.                 this.offetIdxbase0 = offsetIdxbase0;
    
73.                 this.len = len;
    
74.                 this.type = type;
    
75.                 this.memo = memo;
    
76.                 this.example = example;
    
77.             }
    
78.         };
    
79. 
    
80.         static List<trc_field_define> trc_header_fields = new List<trc_field_define>{
    
81.             new trc_field_define("totalLength","文件总长字段", 0, 11, "string", "%文件头，主要用于指示文件大小，数值为9e10+文件字节长度（不包括此文件头的11个字节）", "#9020000350"),
    
82.             new trc_field_define("no_use","大小端", 11+34, 2, "u16", "数据大小端", "=1小端"),
    
83.             new trc_field_define("no_use", "文件头长", 11+36, 4, "u32", "%文件头WAVEDSEC块的长度，4B, long, 36-39", "可认为定长：346，+11字节头部=0x165 == sizeof(lecroy_trc_header)"),
    
84.             new trc_field_define("allDataLength","波形数据总长", 11+60, 4, "u32", "WAVE_ARRAY_1,%4B,long,60-63，数据列1长度", "20000004，一般是采样点数总长=时间+数据"),
    
85.             new trc_field_define("timeDataLength","时间数据总长", 11+112, 4, "u32", "", "10000002, 为波形数据总长/2"),
    
86.             new trc_field_define("sampleDataLength","采样数据总长", 11+116, 4, "u32", "", "10000002，为波形数据总长/2"),
    
87.             new trc_field_define("pntsPerScreen","PNTS_PER_SCREEN", 11+120, 4, "u32", "区间数", "为波形数据总长/2-2"),
    
88.             new trc_field_define("lastvalidPnts","LAST_VALID_PNT", 11+128, 4, "u32", "区间最后一个数据点", "为波形数据总长/2-1"),
    
89.             new trc_field_define("vertical_gain","VERTICAL_GAIN", 11+156, 4, "float", "最终的数值=采样值*Gain+offset", "1.373120030621067e-04"),
    
90.             new trc_field_define("vertical_offset","VERTICAL_OFFSET", 11+160, 4, "float", "垂直零点", "0"),
    
91.             new trc_field_define("max","MAX_VALUE", 11+164, 4, "u32", "最大值", "28875"),
    
92.             new trc_field_define("min","MIN_VALUE", 11+168, 4, "u32", "最小值", "-29131"),
    
93.             new trc_field_define("horizInterval","HORIZ_INTERVAL", 11+176, 4, "float", "采样率", "1e-10"),
    
94.             new trc_field_define("horizOffset","HORIZ_OFFSET", 11+180, 8, "double", "水平偏移", "-5e-4"),
    
95.             new trc_field_define("pixelOffset","PIXEL_OFFSET", 11+188, 8, "u32", "屏幕当前偏移？", "-5e-4"),
    
96.             new trc_field_define("horiz_uncertainty","HORIZ_UNCERTAINTY", 11+292, 4, "float", "区间数", "1e-12，采样率的1%"),
    
97.             new trc_field_define("timeBase","TIMEBASE", 11+324, 2, "u16", "水平栅格时长？这是enum", "24"),
    
98.             new trc_field_define("fixed_Vert_Gain","FIXED_VERT_GAIN", 11+332, 2, "u16", "垂直固定增益？enum", "18")
    
99.         };
    
100. 
     
101.         public struct trcHeaderFields
     
102.         {
     
103.             public string totalLength; //0x9e10+0x0b+lenOfAd*2 "#...."
     
104.             public int allDataLength; //2*lenOfAd(bytes)
     
105.             public int timeDataLength; //lenOfAd;
     
106.             public int sampleDataLength; //lenOfAd;
     
107.             public int pntsPerScreen; //lenfOfAd-2;
     
108.             public int lastvalidPnts; //lenOfAd-1
     
109.             public float vertical_gain; //垂直增益
     
110.             public UInt32 vertical_offset; //垂直偏移
     
111.             public UInt32 max;
     
112.             public UInt32 min;
     
113.             public float horizInterval; //采样率
     
114.             public double horizOffset; //中心点
     
115.             public double pixelOffset; //又一个水平偏移
     
116.             public float horiz_uncertainty; //samplerate/4
     
117.             public UInt16 timeBase; //水平档位
     
118.             public UInt16 fixed_Vert_Gain; //垂直档位
     
119.         };
     
120. 
     
121.         public static Dictionary<string, trc_field_define> dicOfFields = new Dictionary<string, trc_field_define>() { };
     
122. 
     
123. 
     
124.         static public bool WriteTrcFromCsv(float sampleRate, float volPerAd, uint[] sampleData, FileInfo fi)
     
125.         {
     
126.             bool ret  = false;
     
127. 
     
128.             /*
     
129.              * %Header
     
130.                %用户文本块
     
131.                     TEXT % 视上面定义的长度而定
     
132.                % TRIGTIME
     
133.                     TRIGGER_T % 8B f 原名结尾是TIME，但与前面定义重复，故改为T
     
134.                     TRIGGER_OFFSET % 8B f
     
135.                % RISTIME
     
136.                     RIS_OFFSET % 8B f
     
137.                % DATA_ARRAY_1
     
138.                     DATA_ARRAY_1
     
139.                % DATA_ARRAT_2
     
140.                     DATA_ARRAY_2
     
141.                % SIMPLE
     
142.                     SIMPLE int16模式
     
143.                % DUAL
     
144.                     版权声明：本文为CSDN博主「Surrea1」的原创文章，遵循CC 4.0 BY - SA版权协议，转载请附上原文出处链接及本声明。
     
145.                     原文链接：https://blog.csdn.net/Surrea1/article/details/126807087
     
146.             */
     
147.             try
     
148.             {
     
149.                 int totallength = lecroy_trc_header.Length + 2 * sampleData.Length - 11; //no 11 more
     
150.                 StringBuilder sb = new StringBuilder();
     
151.                 sb.AppendFormat("#9{0:000000000}", totallength);
     
152. 
     
153.                 trcHeaderFields headerFields;
     
154.                 headerFields.totalLength = sb.ToString();
     
155.                 headerFields.timeDataLength = sampleData.Length;
     
156.                 headerFields.sampleDataLength = sampleData.Length;
     
157.                 headerFields.allDataLength = sampleData.Length * 2; //in bytes.
     
158.                 headerFields.horizInterval = sampleRate;
     
159.                 headerFields.horizOffset = 0;
     
160.                 headerFields.horiz_uncertainty = sampleRate / 4;
     
161.                 headerFields.pixelOffset = 0;
     
162.                 headerFields.pntsPerScreen = sampleData.Length;
     
163.                 headerFields.lastvalidPnts = sampleData.Length-1;
     
164.                 headerFields.vertical_gain = volPerAd; //垂直增益
     
165.                 headerFields.vertical_offset = 0; //垂直偏移
     
166.                 headerFields.timeBase = 15;
     
167.                 headerFields.fixed_Vert_Gain = 1; //垂直增益
     
168. 
     
169.                 byte[] header = new byte[lecroy_trc_header.Length];
     
170.                 for (int i = 0; i < header.Length; ++i)
     
171.                 {
     
172.                     header[i] = 0;
     
173.                 }
     
174.                 lecroy_trc_header.CopyTo(header, 0);
     
175. 
     
176. 
     
177.                 if (dicOfFields.Count == 0)
     
178.                 {
     
179.                     foreach (trc_field_define i in trc_header_fields)
     
180.                     {
     
181.                         if (dicOfFields.ContainsKey(i.standardName)) continue;
     
182.                         dicOfFields.Add(i.standardName, i);
     
183.                     }
     
184.                 }
     
185.                 Encoding.ASCII.GetBytes(headerFields.totalLength).CopyTo(header, dicOfFields["totalLength"].offetIdxbase0);
     
186.                 BitConverter.GetBytes(headerFields.timeDataLength).CopyTo(header, dicOfFields["timeDataLength"].offetIdxbase0);
     
187.                 BitConverter.GetBytes(headerFields.sampleDataLength).CopyTo(header, dicOfFields["sampleDataLength"].offetIdxbase0);
     
188.                 BitConverter.GetBytes(headerFields.allDataLength).CopyTo(header, dicOfFields["allDataLength"].offetIdxbase0);
     
189.                 BitConverter.GetBytes(headerFields.horizInterval).CopyTo(header, dicOfFields["horizInterval"].offetIdxbase0);
     
190.                 BitConverter.GetBytes(headerFields.horizOffset).CopyTo(header, dicOfFields["horizOffset"].offetIdxbase0);
     
191.                 BitConverter.GetBytes(headerFields.horiz_uncertainty).CopyTo(header, dicOfFields["horiz_uncertainty"].offetIdxbase0);
     
192.                 BitConverter.GetBytes(headerFields.pixelOffset).CopyTo(header, dicOfFields["pixelOffset"].offetIdxbase0);
     
193.                 BitConverter.GetBytes(headerFields.pntsPerScreen).CopyTo(header, dicOfFields["pntsPerScreen"].offetIdxbase0);
     
194.                 BitConverter.GetBytes(headerFields.lastvalidPnts).CopyTo(header, dicOfFields["lastvalidPnts"].offetIdxbase0);
     
195.                 BitConverter.GetBytes(headerFields.vertical_gain).CopyTo(header, dicOfFields["vertical_gain"].offetIdxbase0);
     
196.                 BitConverter.GetBytes(headerFields.vertical_offset).CopyTo(header, dicOfFields["vertical_offset"].offetIdxbase0);
     
197.                 BitConverter.GetBytes(headerFields.timeBase).CopyTo(header, dicOfFields["timeBase"].offetIdxbase0);
     
198.                 BitConverter.GetBytes(headerFields.fixed_Vert_Gain).CopyTo(header, dicOfFields["fixed_Vert_Gain"].offetIdxbase0);
     
199. 
     
200.                 FileStream fs = new FileStream(fi.FullName, FileMode.OpenOrCreate, FileAccess.Write);
     
201.                 BinaryWriter bw = new BinaryWriter(fs);
     
202. 
     
203.                 /*(byte[] zero24 = new byte[24];
     
204.                 for (int i = 0; i < zero24.Length; ++i)
     
205.                 {
     
206.                     zero24[i] = 0;
     
207.                 }
     
208.                 bw.Write(zero24);*/
     
209.                
     
210.                 UInt16 min = UInt16.MaxValue;
     
211.                 UInt16 max = UInt16.MinValue;
     
212.                 foreach (uint data in sampleData)
     
213.                 {
     
214.                     UInt16 d16 = (UInt16)data;
     
215.                     if(min>d16) min = d16;
     
216.                     if(max<d16) max = d16;
     
217.                 }
     
218.                 BitConverter.GetBytes((float)min).CopyTo(header, dicOfFields["min"].offetIdxbase0);
     
219.                 BitConverter.GetBytes((float)max).CopyTo(header, dicOfFields["max"].offetIdxbase0);
     
220.                 bw.Write(header);
     
221.                 foreach (uint data in sampleData)
     
222.                 {
     
223.                     UInt16 d16 = (UInt16)data;
     
224.                     bw.Write(BitConverter.GetBytes(d16));
     
225.                 }
     
226. 
     
227.                 bw.Flush();
     
228.                 fs.Close();
     
229.                 ret = true;
     
230.             }
     
231.             catch (Exception ex)
     
232.             {
     
233.                 Console.WriteLine(ex.Message);
     
234.                 return false;
     
235.             }
     
236. 
     
237.             return ret;
     
238. 
     
239.         }
     
240.     }
     
241. }
     

Copy the Code

python代码中，波形抽样，扩容的代码是第一版的代码，后续改进的部分欠奉。
C#代码中.csv的读取代码欠奉。

最终在力科分析软件中运行的效果如下：

时基的部分做了调整，但是ad值部分，没有对照档位做适配。有兴趣的同志可以在这些代码的基础上继续工作。