第一级代码(by Python):- # code by fengxh.
- # transfer 汉泰.doc => Lecroy 二进制.trc格式。
- # python code is created at Jun30,2023
- from scipy.optimize import curve_fit
- import numpy as np
- import csv
- import sys
- import operator
- import subprocess
- import re
- # 取csv文件的第N列,抛弃掉标题栏,返回数组
- def load_csv_data(filename, idxBase0):
- data = [] # 存放结果的数组
- # 打开CSV文件,并读取数据
- flagDataStart = False;
- with open(filename, 'r') as f:
- reader = csv.reader(f)
- i = 0;
- for row in reader:
- i = i+1;
- #print(i);
- if len(row) == 0:
- if(not flagDataStart):
- flagDataStart = True;
- continue;
- else:
- break; #done!
- if(not flagDataStart): continue;
- try:
- data.append(float(row[idxBase0])) # 取第一列数据,并转为float类型
- except:
- continue;
- return data;
- # 取csv文件的第N列,抛弃掉标题栏,返回数组
- def get_time_scale(filename):
- flagDataStart = False;
- #string = "#CLOCK=0.004000000S"
- pattern = r"CLOCK=(\d+\.\d+)S"
- with open(filename, 'r') as f:
- reader = csv.reader(f)
- for row in reader:
- match = re.search(pattern, str(row))
- if match:
- result = float(match.group(1));
- return result;
- return 100e-6;
- def getRmsOfFreq(arFreq, arRms, freq):
- for i in range(0, len(arFreq)):
- if(arFreq[i] == freq): return arRms[i];
- def csv2trc(csvfile, trcfile, ratio):
- # 要执行的命令
- #Usage: Csv2Trc <src.csv> <tgt.trc> [<sample_rate=100e-6> <volPerAd>=1]
- command = r"C:\Users\fengxh\Desktop\工具集\dumb\Csv2Trc_V1.0.20230531\SampleData_CSV2LecoryTrc";
- command = command + " " +csvfile + ' ' + trcfile + ' ' + str(ratio);
- # 执行命令并获取输出
- output = subprocess.check_output(command, shell=True)
- output = output.decode() # 如果你想要将输出作为字符串处理,可以解码为返回的字节字符串
- print(output)
- def toInteger(rms):
- vals = [];
- for i in rms:
- if(not i in vals):
- vals.append(i);
- sorted_vals = sorted(vals);
- my_dict = {sorted_vals[0]:1};
- idx = 1;
- for v in sorted_vals:
- my_dict[v] = idx;
- idx = idx +1;
- ret = [];
- for i in rms:
- ret.append(my_dict[i]);
-
- return ret;
- def zoom_out_lecory_wave_data(filename, ratio=10):
- rms = load_csv_data(filename, 0); #rms
- dt = get_time_scale(filename);
- print("数据的时间跨度", dt, "秒, 数据点数:", len(rms));
- rms = toInteger(rms);
- data = rms[:];
- for i in range(0, len(rms)):
- data[i//round(ratio) + round(i%(round(ratio))*len(rms)/ratio)] = rms[i];
- data_2d = [[element] for element in data]
- filename = filename + '.zoom_out' + str(ratio) + ".csv";
- with open(filename, 'w', newline='') as csvfile:
- writer = csv.writer(csvfile)
- writer.writerows(data_2d)
- csv2trc(filename, filename+".trc", dt*ratio/len(rms));
- print("done!");
- #调用示例
- filename = "35.25_40K.doc";
- zoom_out_lecory_wave_data(filename,1);
Copy the Code 第二级代码,SampleData_CSV2LecoryTrc,writened in C#, build by VS2022
- /*
- * lecroy_trc.cs
- *
- * A csv to lecory .trc bin file format helper modules.
- *
- * thanks to Lecroy, for its powerful offline wavefrom analyse tools.
- * and Surrea1@CSDN for provide the .trc file format.
- *
- * ----
- * V1.0.20230530 by twicave
- * created. the first verison, coding in vs2022(c#)
- *
- */
- using System;
- using System.Collections.Generic;
- using System.Data.SqlTypes;
- using System.Diagnostics.SymbolStore;
- using System.Linq;
- using System.Reflection.PortableExecutable;
- using System.Security.Cryptography;
- using System.Text;
- using System.Threading;
- using System.Threading.Tasks;
- using System.Runtime.InteropServices;
- using static System.Net.Mime.MediaTypeNames;
- using System.ComponentModel;
- namespace SampleData_CSV2LecoryTrc
- {
- internal class lecroy_trc
- {
- static byte[] lecroy_trc_header = new byte[]{
- 0x23,0x39,0x30,0x30,0x32,0x30,0x30,0x30,0x31,0x34,0x34,0x57,0x41,0x56,0x45,0x44,
- 0x45,0x53,0x43,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x4C,0x45,0x43,0x52,0x4F,
- 0x59,0x5F,0x32,0x5F,0x33,0x00,0x00,0x00,0x00,0x00,0x00,0x01,0x00,0x01,0x00,0x5A,
- 0x01,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
- 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xB6,0x83,0x1E,0x00,0x00,0x00,0x00,0x00,0x00,
- 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x4C,0x45,0x43,0x52,0x4F,0x59,0x57,0x53,0x34,
- 0x32,0x34,0x00,0x00,0x00,0x00,0x00,0x4B,0x45,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
- 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xDB,0x41,0x0F,0x00,0xDB,
- 0x41,0x0F,0x00,0x40,0x42,0x0F,0x00,0x00,0x00,0x00,0x00,0xDA,0x41,0x0F,0x00,0x00,
- 0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x01,0x00,0x00,0x00,0x01,
- 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x0D,0x00,0x80,0x39,0x00,0x00,0x00,0xBF,0x00,
- 0x00,0xFE,0x46,0x00,0x00,0x00,0xC7,0x0B,0x00,0x01,0x00,0xBD,0x37,0x06,0x35,0x1E,
- 0xF4,0xC8,0xB6,0x68,0x78,0xD1,0xBF,0x46,0xB6,0xF3,0xFD,0xD4,0x78,0xD1,0xBF,0x56,
- 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
- 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
- 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x53,
- 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
- 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,
- 0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0xFF,
- 0xEB,0xAF,0x2C,0x99,0xBC,0x0B,0xD8,0x31,0x27,0x39,0x40,0x02,0x0F,0x1A,0x04,0xE7,
- 0x07,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x00,0x01,0x00,0x20,
- 0x00,0x02,0x00,0x00,0x00,0x80,0x3F,0x13,0x00,0x00,0x00,0x00,0x00,0x80,0x3F,0x00,
- 0x00,0x00,0xBF,0x00,0x00 };
- public class trc_field_define
- {
- public string standardName;
- public string field_name;
- public int offetIdxbase0;
- public int len;
- public string type;
- public string memo;
- public string example;
- public trc_field_define(string standardName, string field_name, int offsetIdxbase0, int len, string type, string memo, string example)
- {
- this.standardName = standardName;
- this.field_name = field_name;
- this.offetIdxbase0 = offsetIdxbase0;
- this.len = len;
- this.type = type;
- this.memo = memo;
- this.example = example;
- }
- };
- static List<trc_field_define> trc_header_fields = new List<trc_field_define>{
- new trc_field_define("totalLength","文件总长字段", 0, 11, "string", "%文件头,主要用于指示文件大小,数值为9e10+文件字节长度(不包括此文件头的11个字节)", "#9020000350"),
- new trc_field_define("no_use","大小端", 11+34, 2, "u16", "数据大小端", "=1小端"),
- new trc_field_define("no_use", "文件头长", 11+36, 4, "u32", "%文件头WAVEDSEC块的长度,4B, long, 36-39", "可认为定长:346,+11字节头部=0x165 == sizeof(lecroy_trc_header)"),
- new trc_field_define("allDataLength","波形数据总长", 11+60, 4, "u32", "WAVE_ARRAY_1,%4B,long,60-63,数据列1长度", "20000004,一般是采样点数总长=时间+数据"),
- new trc_field_define("timeDataLength","时间数据总长", 11+112, 4, "u32", "", "10000002, 为波形数据总长/2"),
- new trc_field_define("sampleDataLength","采样数据总长", 11+116, 4, "u32", "", "10000002,为波形数据总长/2"),
- new trc_field_define("pntsPerScreen","PNTS_PER_SCREEN", 11+120, 4, "u32", "区间数", "为波形数据总长/2-2"),
- new trc_field_define("lastvalidPnts","LAST_VALID_PNT", 11+128, 4, "u32", "区间最后一个数据点", "为波形数据总长/2-1"),
- new trc_field_define("vertical_gain","VERTICAL_GAIN", 11+156, 4, "float", "最终的数值=采样值*Gain+offset", "1.373120030621067e-04"),
- new trc_field_define("vertical_offset","VERTICAL_OFFSET", 11+160, 4, "float", "垂直零点", "0"),
- new trc_field_define("max","MAX_VALUE", 11+164, 4, "u32", "最大值", "28875"),
- new trc_field_define("min","MIN_VALUE", 11+168, 4, "u32", "最小值", "-29131"),
- new trc_field_define("horizInterval","HORIZ_INTERVAL", 11+176, 4, "float", "采样率", "1e-10"),
- new trc_field_define("horizOffset","HORIZ_OFFSET", 11+180, 8, "double", "水平偏移", "-5e-4"),
- new trc_field_define("pixelOffset","PIXEL_OFFSET", 11+188, 8, "u32", "屏幕当前偏移?", "-5e-4"),
- new trc_field_define("horiz_uncertainty","HORIZ_UNCERTAINTY", 11+292, 4, "float", "区间数", "1e-12,采样率的1%"),
- new trc_field_define("timeBase","TIMEBASE", 11+324, 2, "u16", "水平栅格时长?这是enum", "24"),
- new trc_field_define("fixed_Vert_Gain","FIXED_VERT_GAIN", 11+332, 2, "u16", "垂直固定增益?enum", "18")
- };
- public struct trcHeaderFields
- {
- public string totalLength; //0x9e10+0x0b+lenOfAd*2 "#...."
- public int allDataLength; //2*lenOfAd(bytes)
- public int timeDataLength; //lenOfAd;
- public int sampleDataLength; //lenOfAd;
- public int pntsPerScreen; //lenfOfAd-2;
- public int lastvalidPnts; //lenOfAd-1
- public float vertical_gain; //垂直增益
- public UInt32 vertical_offset; //垂直偏移
- public UInt32 max;
- public UInt32 min;
- public float horizInterval; //采样率
- public double horizOffset; //中心点
- public double pixelOffset; //又一个水平偏移
- public float horiz_uncertainty; //samplerate/4
- public UInt16 timeBase; //水平档位
- public UInt16 fixed_Vert_Gain; //垂直档位
- };
- public static Dictionary<string, trc_field_define> dicOfFields = new Dictionary<string, trc_field_define>() { };
- static public bool WriteTrcFromCsv(float sampleRate, float volPerAd, uint[] sampleData, FileInfo fi)
- {
- bool ret = false;
- /*
- * %Header
- %用户文本块
- TEXT % 视上面定义的长度而定
- % TRIGTIME
- TRIGGER_T % 8B f 原名结尾是TIME,但与前面定义重复,故改为T
- TRIGGER_OFFSET % 8B f
- % RISTIME
- RIS_OFFSET % 8B f
- % DATA_ARRAY_1
- DATA_ARRAY_1
- % DATA_ARRAT_2
- DATA_ARRAY_2
- % SIMPLE
- SIMPLE int16模式
- % DUAL
- 版权声明:本文为CSDN博主「Surrea1」的原创文章,遵循CC 4.0 BY - SA版权协议,转载请附上原文出处链接及本声明。
- 原文链接:https://blog.csdn.net/Surrea1/article/details/126807087
- */
- try
- {
- int totallength = lecroy_trc_header.Length + 2 * sampleData.Length - 11; //no 11 more
- StringBuilder sb = new StringBuilder();
- sb.AppendFormat("#9{0:000000000}", totallength);
- trcHeaderFields headerFields;
- headerFields.totalLength = sb.ToString();
- headerFields.timeDataLength = sampleData.Length;
- headerFields.sampleDataLength = sampleData.Length;
- headerFields.allDataLength = sampleData.Length * 2; //in bytes.
- headerFields.horizInterval = sampleRate;
- headerFields.horizOffset = 0;
- headerFields.horiz_uncertainty = sampleRate / 4;
- headerFields.pixelOffset = 0;
- headerFields.pntsPerScreen = sampleData.Length;
- headerFields.lastvalidPnts = sampleData.Length-1;
- headerFields.vertical_gain = volPerAd; //垂直增益
- headerFields.vertical_offset = 0; //垂直偏移
- headerFields.timeBase = 15;
- headerFields.fixed_Vert_Gain = 1; //垂直增益
- byte[] header = new byte[lecroy_trc_header.Length];
- for (int i = 0; i < header.Length; ++i)
- {
- header[i] = 0;
- }
- lecroy_trc_header.CopyTo(header, 0);
- if (dicOfFields.Count == 0)
- {
- foreach (trc_field_define i in trc_header_fields)
- {
- if (dicOfFields.ContainsKey(i.standardName)) continue;
- dicOfFields.Add(i.standardName, i);
- }
- }
- Encoding.ASCII.GetBytes(headerFields.totalLength).CopyTo(header, dicOfFields["totalLength"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.timeDataLength).CopyTo(header, dicOfFields["timeDataLength"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.sampleDataLength).CopyTo(header, dicOfFields["sampleDataLength"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.allDataLength).CopyTo(header, dicOfFields["allDataLength"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.horizInterval).CopyTo(header, dicOfFields["horizInterval"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.horizOffset).CopyTo(header, dicOfFields["horizOffset"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.horiz_uncertainty).CopyTo(header, dicOfFields["horiz_uncertainty"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.pixelOffset).CopyTo(header, dicOfFields["pixelOffset"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.pntsPerScreen).CopyTo(header, dicOfFields["pntsPerScreen"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.lastvalidPnts).CopyTo(header, dicOfFields["lastvalidPnts"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.vertical_gain).CopyTo(header, dicOfFields["vertical_gain"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.vertical_offset).CopyTo(header, dicOfFields["vertical_offset"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.timeBase).CopyTo(header, dicOfFields["timeBase"].offetIdxbase0);
- BitConverter.GetBytes(headerFields.fixed_Vert_Gain).CopyTo(header, dicOfFields["fixed_Vert_Gain"].offetIdxbase0);
- FileStream fs = new FileStream(fi.FullName, FileMode.OpenOrCreate, FileAccess.Write);
- BinaryWriter bw = new BinaryWriter(fs);
- /*(byte[] zero24 = new byte[24];
- for (int i = 0; i < zero24.Length; ++i)
- {
- zero24[i] = 0;
- }
- bw.Write(zero24);*/
-
- UInt16 min = UInt16.MaxValue;
- UInt16 max = UInt16.MinValue;
- foreach (uint data in sampleData)
- {
- UInt16 d16 = (UInt16)data;
- if(min>d16) min = d16;
- if(max<d16) max = d16;
- }
- BitConverter.GetBytes((float)min).CopyTo(header, dicOfFields["min"].offetIdxbase0);
- BitConverter.GetBytes((float)max).CopyTo(header, dicOfFields["max"].offetIdxbase0);
- bw.Write(header);
- foreach (uint data in sampleData)
- {
- UInt16 d16 = (UInt16)data;
- bw.Write(BitConverter.GetBytes(d16));
- }
- bw.Flush();
- fs.Close();
- ret = true;
- }
- catch (Exception ex)
- {
- Console.WriteLine(ex.Message);
- return false;
- }
- return ret;
- }
- }
- }
Copy the Code
python代码中,波形抽样,扩容的代码是第一版的代码,后续改进的部分欠奉。
C#代码中.csv的读取代码欠奉。
最终在力科分析软件中运行的效果如下:
时基的部分做了调整,但是ad值部分,没有对照档位做适配。有兴趣的同志可以在这些代码的基础上继续工作。
|