tf_cnn.py 4.0 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293
  1. #!/usr/bin/env python
  2. # -*- coding:utf-8 -*-
  3. # @FileName :nn_bp.py
  4. # @Time :2025/2/12 10:41
  5. # @Author :David
  6. # @Company: shenyang JY
  7. from tensorflow.keras.layers import Input, Dense, LSTM, concatenate, Conv1D, Conv2D, MaxPooling1D, Reshape, Flatten
  8. from tensorflow.keras.models import Model, load_model
  9. from tensorflow.keras.callbacks import ModelCheckpoint, EarlyStopping, TensorBoard, ReduceLROnPlateau
  10. from tensorflow.keras import optimizers, regularizers
  11. from models_processing.losses.loss_cdq import rmse
  12. import numpy as np
  13. from common.database_dml import *
  14. from threading import Lock
  15. import argparse
  16. model_lock = Lock()
  17. class CNNHandler(object):
  18. def __init__(self, logger, args):
  19. self.logger = logger
  20. self.opt = argparse.Namespace(**args)
  21. self.model = None
  22. def get_model(self, args):
  23. """
  24. 单例模式+线程锁,防止在异步加载时引发线程安全
  25. """
  26. try:
  27. with model_lock:
  28. # NPHandler.model = NPHandler.get_keras_model(opt)
  29. self.model = get_h5_model_from_mongo(args, {'rmse': rmse})
  30. except Exception as e:
  31. self.logger.info("加载模型权重失败:{}".format(e.args))
  32. @staticmethod
  33. def get_keras_model(opt):
  34. # db_loss = NorthEastLoss(opt)
  35. # south_loss = SouthLoss(opt)
  36. l1_reg = regularizers.l1(opt.Model['lambda_value_1'])
  37. l2_reg = regularizers.l2(opt.Model['lambda_value_2'])
  38. nwp_input = Input(shape=(opt.Model['time_step'], opt.Model['input_size']), name='nwp')
  39. con1 = Conv1D(filters=64, kernel_size=1, strides=1, padding='valid', activation='relu', kernel_regularizer=l2_reg)(nwp_input)
  40. d1 = Dense(32, activation='relu', name='d1', kernel_regularizer=l1_reg)(con1)
  41. nwp = Dense(8, activation='relu', name='d2', kernel_regularizer=l1_reg)(d1)
  42. output = Dense(1, name='d5')(nwp)
  43. output_f = Flatten()(output)
  44. model = Model(inputs=nwp_input, outputs=output_f)
  45. adam = optimizers.Adam(learning_rate=opt.Model['learning_rate'], beta_1=0.9, beta_2=0.999, epsilon=1e-7, amsgrad=True)
  46. reduce_lr = ReduceLROnPlateau(monitor='val_loss', factor=0.01, patience=5, verbose=1)
  47. model.compile(loss=rmse, optimizer=adam)
  48. return model
  49. def train_init(self):
  50. try:
  51. if self.opt.Model['add_train']:
  52. # 进行加强训练,支持修模
  53. base_train_model = get_h5_model_from_mongo(vars(self.opt), {'rmse': rmse})
  54. base_train_model.summary()
  55. self.logger.info("已加载加强训练基础模型")
  56. else:
  57. base_train_model = self.get_keras_model(self.opt)
  58. return base_train_model
  59. except Exception as e:
  60. self.logger.info("加强训练加载模型权重失败:{}".format(e.args))
  61. def training(self, train_and_valid_data):
  62. model = self.train_init()
  63. # tf.reset_default_graph() # 清除默认图
  64. train_x, train_y, valid_x, valid_y = train_and_valid_data
  65. print("----------", np.array(train_x[0]).shape)
  66. print("++++++++++", np.array(train_x[1]).shape)
  67. model.summary()
  68. early_stop = EarlyStopping(monitor='val_loss', patience=self.opt.Model['patience'], mode='auto')
  69. history = model.fit(train_x, train_y, batch_size=self.opt.Model['batch_size'], epochs=self.opt.Model['epoch'], verbose=2, validation_data=(valid_x, valid_y), callbacks=[early_stop], shuffle=False)
  70. loss = np.round(history.history['loss'], decimals=5)
  71. val_loss = np.round(history.history['val_loss'], decimals=5)
  72. self.logger.info("-----模型训练经过{}轮迭代-----".format(len(loss)))
  73. self.logger.info("训练集损失函数为:{}".format(loss))
  74. self.logger.info("验证集损失函数为:{}".format(val_loss))
  75. return model
  76. def predict(self, test_x, batch_size=1):
  77. result = self.model.predict(test_x, batch_size=batch_size)
  78. self.logger.info("执行预测方法")
  79. return result
  80. if __name__ == "__main__":
  81. run_code = 0