重庆梁平网站制作公司,p2p的网站建设,怎样用ps设计网站模板,网络推广怎么干目前金属矿开采#xff0c;爆破还是主要的破岩方式#xff0c;为了保证巷道采场的安全#xff0c;需要对爆破震动进行监测#xff0c;获取的监测数据如附件#xff0c;第1列数据为震动的序号#xff0c;第2、3、4列为x,y,z三个方向的震动速度#xff0c;往往由于各种因素…目前金属矿开采爆破还是主要的破岩方式为了保证巷道采场的安全需要对爆破震动进行监测获取的监测数据如附件第1列数据为震动的序号第2、3、4列为x,y,z三个方向的震动速度往往由于各种因素如环境、采集设备故障等导致采集的数据存在缺失、异常、冗余等现象这会降低系统分析和预测能力因此如何将获取“脏数据”进行科学合理的清洗是亟待解决的问题。一般说来脏数据主要存在异常、缺失和冗余这三类如某一行的数据不是4列或者x,y,z缺失或者格式为字符不是数字等。对数据进行清洗后请计算出爆破振动的合速度建立新的文件保存在第5列在文件开头保存程序运行时的系统时间和程序开发者姓名并绘制出合速度的波形图,绘制波形使用turtle库或matplotlib库。 X(S),X(cm/s),Y(cm/s),Z(cm/s) -2023.2048,-0.040775947,-0.026564745,-0.01956255 -0.2046,-0.040775947,-0.026564745,-0.01956255 -0.2044,-0.043995101,-0.025502155,-0.023680981 -0.2042,-0.045068152,-0.027627335,-0.023680981 -0.204,-0.046141203,-0.032940284,-0.020592158 -0.204,-0.046141203,-0.032940284,-0.020592158,ttytt 110.204,-0.046141203,-0.032940284,-0.020592158 -0.204,,-0.032940284,-0.020592158 -0.204,dog,-0.032940284,-0.020592158 -0.2038,-0.048287306,-0.035065463,-0.024710589 -0.2036,-0.045068152,-0.031877694,-0.023680981 -0.2034,-0.046141203,-0.024439565,-0.018532942 -0.2032,-0.035410691,-0.023376976,-0.022651373 -0.203,-0.035410691,-0.022314386,-0.01956255 -0.2028,-0.035410691,-0.021251796,-0.020592158 -0.2026,-0.032191537,-0.021251796,-0.01956255 -0.2024,-0.033264589,-0.018064027,-0.021621766 · · · 实现代码如下
import csv
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
import datetimewith open(test.txt, w) as f:f.write(ID:20220000,Name:XXX,time:%s\n % datetime.datetime.now().strftime(%Y-%m-%d %H:%M:%S))out open(data.csv, w, newline, encodingutf-8)
csv_writer csv.writer(out, dialectexcel)f open(data.txt, r)
flag 0
for line in f.readlines():lst list(line.strip(\n).split(,)) # 将字符串转为列表从而可以按单元格写入csvfor item in lst:if item X(S) or item X(cm/s) or item Y(cm/s) or item Z(cm/s):flag 4breakelif item.islower():flag 0breakelif item.find(.) 3 and item.find(.) 0:flag flag 1continueif flag 4:csv_writer.writerow(lst)flag 0# 读取数据文件
reader pd.read_csv(data.csv, header[0])
X reader.iloc[:, 1]
Y reader.iloc[:, 2]
Z reader.iloc[:, 3]
T reader.iloc[:, 0]
# 计算合速度
V np.sqrt(X ** 2 Y ** 2 Z ** 2) * (X / abs(X))
# 将V存入第五列
reader[V] V
reader.to_csv(data.csv, indexFalse)
# 生成结果表单
out open(result.csv, w, newline, encodingutf-8)
csv_writer2 csv.writer(out, dialectexcel)
f1 open(test.txt)
for line in f1.readlines():lst list(line.strip(\n).split(,))csv_writer2.writerow(lst)
f2 open(data.csv)
for line in f2.readlines():lst list(line.strip(\n).split(,))csv_writer2.writerow(lst)
# 绘制波形图
plt.figure(num1, figsize(8, 4))
plt.plot(T, V, colorgreen)
plt.xlabel(time)
plt.ylabel(velocity)
plt.title(Blasting vibration composite velocity waveform graph)
plt.grid()
plt.show() 效果 生成的数据文件会自动和代码文件储存在一起。
本代码在pycharm中编写在其他编译器中如有问题概不负责。
