pv-current/deviation/covar.hy

" 实现标准输入格式到算子输入的转换
保证不论输入的是一个计算单元，还是多个计算单元，不影响算子实现

大写字段名用于输入、输出，小写字段名用于计算过程中的存储字段"

(import [pandas :as pd] 
        logging
        json
        [numpy :as np]
        [typing [Tuple]])

(setv START-TIME "13:00"
      END-TIME "17:00"
      TF 30
      I0 1
      IM 12
      STATION-ID "STATION_ID"      ; 对应输入数据表中的场站标识
      DEVICE-ID "DEVICE_ID"        ; 对应输入数据表中的设备标识
      TIME-HEADER "MONITOR_TIME"
      REPORT-TIME "ALERT_TIME"
      V-HEADER "VOLTAGE_VALUE"
      CURRENT-HEADER "CURRENT_VALUE"
      GRP-ID "grp_id"              ; 分组标识符，每组包含属于一个计算单元的多个组串
      STR-HEADER "str_no"
      DEVIATION "DEVIATION"        ; 离散度计算结果列名
      STR-STATUS "STR_STATUS"      ; 组串状态

      COM-FAIL -1                  ; 通讯故障标志
      UNIT-INVALID -2              ; 计算单元无有效值标志
      STR-INVALID -3              ; 组串无有效值标志
      
      COM-FAIL-DF (.DataFrame pd
                              {DEVIATION [COM-FAIL] STR_STATUS ["-"]}
                              :index (.Index pd [END_TIME])) 
      UNIT-INVALID-DF (.DataFrame pd
                              {DEVIATION [UNIT_INVALID] STR_STATUS ["-"]}
                              :index (.Index pd [END_TIME])))

(defn get-valid-zone [indf thr] "获取有效计算区间 [t1, t2]"
  (setv above-i0-vals (. indf loc [(.any (>= (. indf [CURRENT-HEADER]) thr) :axis 1)]))
  (assert not (.empty above-i0-flag)) ; 通过 step 2 的验证，本步中至少应该有一个大于 I0 的电流值
  (. above-i0-vals index [0]))

(defn filter-window [indf start-index] "按照文档第3步要求对过滤和校正组串电流值"
  (setv df2 (. indf [start-index:]))
  (setv df3 ))