data_new = np.array([]) data = file_data['面积(㎡)'].values for i in data: data_new = np.append(data_new,np.array(i[:-2])) data = data_new.astype(np.float64) file_data.loc[:'面积(㎡)']=data file_data报错修改

时间: 2023-11-28 21:06:22 浏览: 70
根据您提供的代码,我看到您可能是想将 `file_data` 中的 `面积(㎡)` 字段的数值后面的单位(即 '㎡')去掉,并将剩余的数据转换为浮点数,然后更新回 `file_data` 中相应的字段。但是您的代码有一些问题,可能导致报错。 首先,您初始化 `data_new` 的时候使用了一个空的 numpy 数组,这样的话在循环过程中每次都需要将 `data_new` 扩展一次,这样的效率非常低下。建议您在初始化 `data_new` 的时候指定一个合适的大小。 其次,您在循环过程中使用了 `np.append()`,这个函数的效率也比较低下,因为它每次都会创建一个新的数组,并将原数组的数据复制到新数组中。建议您使用 Python 列表进行存储,然后在循环结束后一次性转换为 numpy 数组。 最后,您在更新 `file_data` 的时候使用了 `loc`,但是 `loc` 的用法是从行标签或列标签中选择行或列,而不是对整个 DataFrame 进行赋值。建议您使用 `iloc` 或直接对 `file_data['面积(㎡)']` 进行赋值。 下面是修改后的代码示例: ``` data = file_data['面积(㎡)'].values data_new = [float(i[:-2]) for i in data] file_data['面积(㎡)'] = data_new ``` 这个代码将 `data` 中的每个元素去掉最后两个字符(即 '㎡'),然后转换为浮点数,存储到 `data_new` 列表中。最后将 `data_new` 赋值给 `file_data` 中的 '面积(㎡)' 列。

相关推荐

import pyntcloud from scipy.spatial import cKDTree import numpy as np def pass_through(cloud, limit_min=-10, limit_max=10, filter_value_name="z"): """ 直通滤波 :param cloud:输入点云 :param limit_min: 滤波条件的最小值 :param limit_max: 滤波条件的最大值 :param filter_value_name: 滤波字段(x or y or z) :return: 位于[limit_min,limit_max]范围的点云 """ points = np.asarray(cloud.points) if filter_value_name == "x": ind = np.where((points[:, 0] >= limit_min) & (points[:, 0] <= limit_max))[0] x_cloud = pcd.select_by_index(ind) return x_cloud elif filter_value_name == "y": ind = np.where((points[:, 1] >= limit_min) & (points[:, 1] <= limit_max))[0] y_cloud = cloud.select_by_index(ind) return y_cloud elif filter_value_name == "z": ind = np.where((points[:, 2] >= limit_min) & (points[:, 2] <= limit_max))[0] z_cloud = pcd.select_by_index(ind) return z_cloud # -------------------读取点云数据并可视化------------------------ # 读取原始点云数据 cloud_before=pyntcloud.PyntCloud.from_file("./data/pcd/000000.pcd") # 进行点云下采样/滤波操作 # 假设得到了处理后的点云(下采样或滤波后) pcd = o3d.io.read_point_cloud("./data/pcd/000000.pcd") filtered_cloud = pass_through(pcd, limit_min=-10, limit_max=10, filter_value_name="x") # 获得原始点云和处理后的点云的坐标值 points_before = cloud_before.points.values points_after = filtered_cloud.points.values # 使用KD-Tree将两组点云数据匹配对应,求解最近邻距离 kdtree_before = cKDTree(points_before) distances, _ = kdtree_before.query(points_after) # 计算平均距离误差 ade = np.mean(distances) print("滤波前后的点云平均距离误差为:", ade) o3d.visualization.draw_geometries([filtered_cloud], window_name="直通滤波", width=1024, height=768, left=50, top=50, mesh_show_back_face=False) # 创建一个窗口,设置窗口大小为800x600 vis = o3d.visualization.Visualizer() vis.create_window(width=800, height=600) # 设置视角点 ctr = vis.get_view_control() ctr.set_lookat([0, 0, 0]) ctr.set_up([0, 0, 1]) ctr.set_front([1, 0, 0])这段程序有什么问题吗

import numpy as np import matplotlib.pyplot as plt import pickle as pkl import pandas as pd import tensorflow.keras from tensorflow.keras.models import Sequential, Model, load_model from tensorflow.keras.layers import LSTM, GRU, Dense, RepeatVector, TimeDistributed, Input, BatchNormalization, \ multiply, concatenate, Flatten, Activation, dot from sklearn.metrics import mean_squared_error,mean_absolute_error from tensorflow.keras.optimizers import Adam from tensorflow.python.keras.utils.vis_utils import plot_model from tensorflow.keras.callbacks import EarlyStopping from keras.callbacks import ReduceLROnPlateau df = pd.read_csv('lorenz.csv') signal = df['signal'].values.reshape(-1, 1) x_train_max = 128 signal_normalize = np.divide(signal, x_train_max) def truncate(x, train_len=100): in_, out_, lbl = [], [], [] for i in range(len(x) - train_len): in_.append(x[i:(i + train_len)].tolist()) out_.append(x[i + train_len]) lbl.append(i) return np.array(in_), np.array(out_), np.array(lbl) X_in, X_out, lbl = truncate(signal_normalize, train_len=50) X_input_train = X_in[np.where(lbl <= 9500)] X_output_train = X_out[np.where(lbl <= 9500)] X_input_test = X_in[np.where(lbl > 9500)] X_output_test = X_out[np.where(lbl > 9500)] # Load model model = load_model("model_forecasting_seq2seq_lstm_lorenz.h5") opt = Adam(lr=1e-5, clipnorm=1) model.compile(loss='mean_squared_error', optimizer=opt, metrics=['mae']) #plot_model(model, to_file='model_plot.png', show_shapes=True, show_layer_names=True) # Train model early_stop = EarlyStopping(monitor='val_loss', patience=20, verbose=1, mode='min', restore_best_weights=True) #reduce_lr = ReduceLROnPlateau(monitor='val_loss', factor=0.2, patience=9, verbose=1, mode='min', min_lr=1e-5) #history = model.fit(X_train, y_train, epochs=500, batch_size=128, validation_data=(X_test, y_test),callbacks=[early_stop]) #model.save("lstm_model_lorenz.h5") # 对测试集进行预测 train_pred = model.predict(X_input_train[:, :, :]) * x_train_max test_pred = model.predict(X_input_test[:, :, :]) * x_train_max train_true = X_output_train[:, :] * x_train_max test_true = X_output_test[:, :] * x_train_max # 计算预测指标 ith_timestep = 10 # Specify the number of recursive prediction steps # List to store the predicted steps pred_len =2 predicted_steps = [] for i in range(X_output_test.shape[0]-pred_len+1): YPred =[],temdata = X_input_test[i,:] for j in range(pred_len): Ypred.append (model.predict(temdata)) temdata = [X_input_test[i,j+1:-1],YPred] # Convert the predicted steps into numpy array predicted_steps = np.array(predicted_steps) # Plot the predicted steps #plt.plot(X_output_test[0:ith_timestep], label='True') plt.plot(predicted_steps, label='Predicted') plt.legend() plt.show()

import pandas as pd import os from scipy import integrate, signal import numpy as np import matplotlib import matplotlib.pyplot as plt matplotlib.rcParams['font.sans-serif'] = ['SimHei'] # 显示中文 matplotlib.rcParams['axes.unicode_minus'] = False # 用来正常显示正负号 #y = pd.read_excel(r'C:\Users\ppddcsm\Desktop\第一批数据拆分\第一批1号1振\A1-1-600.xlsx', usecols=[1],index_col=False, header=None ,skiprows=[0]) folder_path = r'C:\Users\ppddcsm\Desktop\第二批数据拆分\第二批1号1振' file_names = os.listdir(folder_path) for file_name in file_names: file_path = os.path.join(folder_path, file_name) y = pd.read_excel(file_path) N = len(y) fs = 1280 dt = 1/fs #t_axis = [i * dt for i in range(len(y))] # 时间轴 t_axis = [i * dt for i in range(len(y))] y1 = y.swapaxes(0, 1) # 矩阵转置 data = y1.fillna(-1).values #获取数据,将缺失值标记设置为-1,并转换为NumPy数组对象 t = data.flatten() # 展平数组 a = np.array(t) # 梯形法 cumtrapz累计计算积分,cumtrapz(y, x=None, dx=1.0, axis=-1, initial=None)。y: 需要被积分的数值序列;x: y中元素的间距,积分变量,若为空,则y元素的间距默认为dx; # 续:dx: 如果x为空,y中元素的间距由dx给出;axis: 确定积分轴;initial: 如果提供,则用该值作为返回值的第一个数值。 #y_int = integrate.cumtrapz(np.array(a), x=None, dx=0.00078125, initial=0)*1000 # m到mm转换要乘1000 #y_int = np.multiply(integrate.cumtrapz(np.array(a), x=None, dx=0.00078125, initial=0), 1000) Y = integrate.cumtrapz(np.array(a), x=None, dx=0.00078125, initial=0) y_int = np.multiply(Y, 1000)

最新推荐

recommend-type

rockyou.txt

rockyou
recommend-type

ASP+ACCESS网上人才信息管理系统(源代码+论文)【ASP】.zip

ASP+ACCESS网上人才信息管理系统(源代码+论文)【ASP】
recommend-type

河北金融学院经济大数据课设2024年 软科学校爬虫课设

河北金融学院经济大数据课设2024年 软科学校爬虫课设
recommend-type

widgetsnbextension-4.0.0b0-py3-none-any.whl

Python库是一组预先编写的代码模块,旨在帮助开发者实现特定的编程任务,无需从零开始编写代码。这些库可以包括各种功能,如数学运算、文件操作、数据分析和网络编程等。Python社区提供了大量的第三方库,如NumPy、Pandas和Requests,极大地丰富了Python的应用领域,从数据科学到Web开发。Python库的丰富性是Python成为最受欢迎的编程语言之一的关键原因之一。这些库不仅为初学者提供了快速入门的途径,而且为经验丰富的开发者提供了强大的工具,以高效率、高质量地完成复杂任务。例如,Matplotlib和Seaborn库在数据可视化领域内非常受欢迎,它们提供了广泛的工具和技术,可以创建高度定制化的图表和图形,帮助数据科学家和分析师在数据探索和结果展示中更有效地传达信息。
recommend-type

zigbee-cluster-library-specification

最新的zigbee-cluster-library-specification说明文档。
recommend-type

管理建模和仿真的文件

管理Boualem Benatallah引用此版本:布阿利姆·贝纳塔拉。管理建模和仿真。约瑟夫-傅立叶大学-格勒诺布尔第一大学,1996年。法语。NNT:电话:00345357HAL ID:电话:00345357https://theses.hal.science/tel-003453572008年12月9日提交HAL是一个多学科的开放存取档案馆,用于存放和传播科学研究论文,无论它们是否被公开。论文可以来自法国或国外的教学和研究机构,也可以来自公共或私人研究中心。L’archive ouverte pluridisciplinaire
recommend-type

MATLAB结构体与对象编程:构建面向对象的应用程序,提升代码可维护性和可扩展性

![MATLAB结构体与对象编程:构建面向对象的应用程序,提升代码可维护性和可扩展性](https://picx.zhimg.com/80/v2-8132d9acfebe1c248865e24dc5445720_1440w.webp?source=1def8aca) # 1. MATLAB结构体基础** MATLAB结构体是一种数据结构,用于存储和组织相关数据。它由一系列域组成,每个域都有一个名称和一个值。结构体提供了对数据的灵活访问和管理,使其成为组织和处理复杂数据集的理想选择。 MATLAB中创建结构体非常简单,使用struct函数即可。例如: ```matlab myStruct
recommend-type

详细描述一下STM32F103C8T6怎么与DHT11连接

STM32F103C8T6可以通过单总线协议与DHT11连接。连接步骤如下: 1. 将DHT11的VCC引脚连接到STM32F103C8T6的5V电源引脚; 2. 将DHT11的GND引脚连接到STM32F103C8T6的GND引脚; 3. 将DHT11的DATA引脚连接到STM32F103C8T6的GPIO引脚,可以选择任一GPIO引脚,需要在程序中配置; 4. 在程序中初始化GPIO引脚,将其设为输出模式,并输出高电平,持续至少18ms,以激活DHT11; 5. 将GPIO引脚设为输入模式,等待DHT11响应,DHT11会先输出一个80us的低电平,然后输出一个80us的高电平,
recommend-type

JSBSim Reference Manual

JSBSim参考手册,其中包含JSBSim简介,JSBSim配置文件xml的编写语法,编程手册以及一些应用实例等。其中有部分内容还没有写完,估计有生之年很难看到完整版了,但是内容还是很有参考价值的。
recommend-type

"互动学习:行动中的多样性与论文攻读经历"

多样性她- 事实上SCI NCES你的时间表ECOLEDO C Tora SC和NCESPOUR l’Ingén学习互动,互动学习以行动为中心的强化学习学会互动,互动学习,以行动为中心的强化学习计算机科学博士论文于2021年9月28日在Villeneuve d'Asq公开支持马修·瑟林评审团主席法布里斯·勒菲弗尔阿维尼翁大学教授论文指导奥利维尔·皮耶昆谷歌研究教授:智囊团论文联合主任菲利普·普雷教授,大学。里尔/CRISTAL/因里亚报告员奥利维耶·西格德索邦大学报告员卢多维奇·德诺耶教授,Facebook /索邦大学审查员越南圣迈IMT Atlantic高级讲师邀请弗洛里安·斯特鲁布博士,Deepmind对于那些及时看到自己错误的人...3谢谢你首先,我要感谢我的两位博士生导师Olivier和Philippe。奥利维尔,"站在巨人的肩膀上"这句话对你来说完全有意义了。从科学上讲,你知道在这篇论文的(许多)错误中,你是我可以依