Myu: ESP32灯油残量センサー　実験

まだ、測定値が安定しないので、新たに容量カウント補正用のセンサ（タンク入口までの短い線）を追加して、しばらくデータをとってみることにした。追加PINはGPIO15。esp32から外に出ている線がぬれるだけで、測定値に影響するので、補正できないか試す予定。

*******esp32
  #include <WiFi.h>
#include <WiFiClient.h>
#include "esp_sleep.h"
#include <Wire.h>
#include "SparkFunBME280.h"

#define TOUCH_PIN 12
#define TOUCH_PIN2 15
#define LED_PIN 2
#define SLEEP_SECONDS 21600   // 6時間

#define SERIAL_EN 0   // 1=debug, 0=off
#if SERIAL_EN
  #define S_BEGIN(...)   Serial.begin(__VA_ARGS__)
  #define S_PRINT(...)   Serial.print(__VA_ARGS__)
  #define S_PRINTLN(...) Serial.println(__VA_ARGS__)
    #define S_PRINTF(...)  Serial.printf(__VA_ARGS__)
  #define S_END(...) Serial.end()
#else
  #define S_BEGIN(...)
  #define S_PRINT(...)
  #define S_PRINTLN(...)
    #define S_PRINTF(...)
  #define S_END(...)
#endif

const char* ssid     = "***";
const char* password = "********";

const char* host = "*.*.*.*";  // Raspberry Pi IP
const int   port = 8888;
const char* path = "/py/cap";

String url;
String cprint;

BME280 bme;

void setup() {

S_BEGIN(115200);

delay(100);
   uint16_t capValue= readTouchStable(TOUCH_PIN);
   delay(300);
   uint16_t capValue2= readTouchStable(TOUCH_PIN2);
   delay(100);

Wire.begin(21, 22);  // SDA=21, SCL=22(SCKは間違い)
    bme.setI2CAddress(0x76);

if (bme.beginI2C() == false) {
       S_PRINTLN("BME280 not detected!");
    } else {
      
      S_PRINTLN("BME280 connected.");
    }

WiFi.begin(ssid, password);
  
   S_PRINTLN("Connecting to WiFi");
    while (WiFi.status() != WL_CONNECTED) {
        delay(500);
   
      S_PRINTLN(".");
    }

S_PRINTLN("\nWiFi connected");

// uint16_t capValue = touchRead(TOUCH_PIN);

S_PRINTF("Touch value: %d\n", capValue);
  S_PRINTF("Touch2 value: %d\n", capValue2);

float tempC = bme.readTempC();
    float humidity = bme.readFloatHumidity();
    float pressure_hPa = bme.readFloatPressure() / 100.0;

S_PRINTF("Temp: %.2f C, Humidity: %.2f %%, Pressure: %.2f hPa\n", tempC, humidity, pressure_hPa);
 
    url  = String(path)
         + "?cap=" + String(capValue)
         + "&t="   + String(tempC, 1)
         + "&h="   + String(humidity, 1)
         + "&p="   + String(pressure_hPa, 1)
         + "&cap2=" + String(capValue2);

cprint = String("GET ") + url + " HTTP/1.1\r\n" +
             "Host: " + host + "\r\n" +
             "Connection: close\r\n\r\n";

WiFiClient client;
    if (client.connect(host, port)) {
        
        S_PRINTLN("Connected to server");
        client.print(cprint);

while (client.connected() || client.available()) {
            if (client.available()) {
                String line = client.readStringUntil('\n');
             
              S_PRINTLN(line);
            }
        }
        client.stop();
    } else {
      
      S_PRINTLN("Connection failed url=" + url);
    }

S_PRINTLN("Going to sleep...");
    delay(100);

WiFi.disconnect(true);
    WiFi.mode(WIFI_OFF);

digitalWrite(LED_PIN, LOW);

S_END();

esp_sleep_enable_timer_wakeup(SLEEP_SECONDS * 1000000ULL);
    esp_deep_sleep_start();
}

void loop() {

}

int readTouchStable(int pin) {
   // チャネル切替直後の捨て読み
  (void)touchRead(pin);
  delay(5);

int sum = 0;
    for (int i = 0; i < 8; i++) {
        sum += touchRead(pin);
        delay(5); // ほんの少し間隔
    }
    return sum / 8;
}

*******bottle

@app.route('/py/cap') 
@local_only
def cap_handler():

cap_value = request.query.cap  # ?cap=123 または ?cap=list10   list10だと最新の１０個データ表示
    if cap_value:
        if cap_value.startswith("list"):
            # ?cap=listN の場合
            try:
                n = int(cap_value[4:])
                lines = read_file(LOG_FILE_CAP)
                start = max(0, len(lines) - n)
                latest_lines = lines[start:]
                response.content_type = 'application/json'
                # JSON形式で返す
                result = []
                for line in latest_lines:
                    dt_str,value_str,value_t,value_p,value_h,value_str2 = line.strip().split(',', 5)
                    result.append({
                     "datetime": dt_str,
                     "cap": int(value_str),
                     "temp": float(value_t),
                     "hpa": float(value_p),
                     "hum": float(value_h),
                     "cap2":int(value_str2)
                     })
                return json.dumps(result)
            except ValueError:
                response.status = 400
                return "Invalid list parameter"
        else:
            # 数値の場合
            try:
                t_value= request.query.t
                p_value= request.query.p
                h_value= request.query.h
                cap2_value= request.query.cap2

latest_cap = int(cap_value)
                latest_t=float(t_value)
                latest_p=float(p_value)
                latest_h=float(h_value)
                #latest_cap2=int(cap2_value)
                latest_cap2 = int(cap2_value) if cap2_value else 0
                now_str = datetime.now().strftime("%Y-%m-%d %H:%M:%S")
                line = f"{now_str},{latest_cap},{latest_t},{latest_p},{latest_h},{latest_cap2}\n"
                
                # ログを読み込んで追記
                lines = read_file(LOG_FILE_CAP)
                lines.append(line)
                # 最大行数を維持
                if len(lines) > MAX_LINES:
                    lines = lines[-MAX_LINES:]
                write_log(lines)

return f"Cap value {latest_cap} saved.\n"
            except ValueError:
                response.status = 400
                return "Invalid value"
    else:
        # パラメータなし → 最新値返す
        if latest_cap is not None:
            response.content_type = 'application/json'
            return json.dumps({"latest_cap": latest_cap})
        else:
            response.status = 404
            return "No data yet"

@app.route('/py/captable')
@local_only
def captable():
    return """
<!DOCTYPE html>
<html lang="ja">
<head>
<meta charset="UTF-8">
<meta name="viewport" content="width=device-width, initial-scale=1"> 
<title>灯油タンク残量概算</title>
<style>
body {
  font-family: sans-serif;
  margin: 1em;
}

/* 入力エリア */
.controls {
  display: flex;
  flex-wrap: wrap;
  gap: 0.5em;
  align-items: center;
}
.controls label, .controls button {
  font-size: 1em;
}

/* テーブル */
table {
  border-collapse: collapse;
  width: 100%;
  margin-top: 1em;
  font-size: 0.95em;
}
th, td {
  border: 1px solid #999;
  padding: 0.5em;
  text-align: center;
}

/* スマホ向け：フォント少し小さく、テーブル横スクロール可 */
@media (max-width: 600px) {
  body { font-size: 0.9em; }
  .controls { flex-direction: column; align-items: flex-start; }
  table { display: block; overflow-x: auto; white-space: nowrap; }
  th, td { padding: 0.4em 0.6em; }
}
</style>
</head>
<body>
<h3>残量データ一覧</h3>

let prevFiltered = null;
let alpha = 0.07;
let threshold = 20;

let aval=1.4;
let eval=2.05;
let hval=51;
let jval=0.6;
let lval="45";

async function loadParams(){
    const response = await fetch("/api/get_params");
    const data = await response.json();

return {
        aval: Number(data.aval)  || 1.4,
        eval: Number(data.eval)  || 2.05,
        hval: Number(data.hval)  || 51,
        jval: Number(data.jval)  || 0.6,
        alpha: Number(data.alpha) || 0.07,
        threshold: Number(data.threshold) || 20,
        lval: data.lval || "45"
    };
}

function capToLevel(cap, cap2, ondo, hum) {
    if (cap === 0) return 0.0;

//let rawValue = 1.4 * (1/cap - 2.05/cap2) * 100000 + 51 - 0.6 * hum;
    let rawValue = aval * (1/cap - eval/cap2) * 100000 + hval - jval * hum;

// 初回はそのまま
    if (prevFiltered === null) {
        prevFiltered = rawValue;
        return rawValue;
    }

//let alpha = 0.07;

// 給油などの急増は即追従
    if (rawValue - prevFiltered > threshold) {
        prevFiltered = rawValue;
    } else {
        prevFiltered = alpha * rawValue + (1 - alpha) * prevFiltered;
    }

return prevFiltered;
}

async function loadData() {
  prevFiltered = null;

const params = await loadParams();

const lval2 = params.lval;

const res = await fetch('/py/cap?cap=list'+lval2);

const data = await res.json();
  const tbody = document.querySelector('#capTable tbody');
  tbody.innerHTML = '';

// 最新を上に
  data.slice().forEach(item => {
    const tr = document.createElement('tr');
    const conv = capToLevel(item.cap,item.cap2,item.temp,item.hum);
    tr.innerHTML = `
      <td>${item.datetime}</td>
      <td>${item.cap}</td>
      <td>${conv.toFixed(2)}</td>
      <td>${item.temp}</td>
      <td>${item.hpa}</td>
      <td>${item.hum}</td>
      <td>${item.cap2}</td>
    `;
    tbody.appendChild(tr);
  });
}

loadData();
</script>
</body>
</html>
"""

@app.route('/py/dashboard2')
@local_only
def dashboard2():
    return """
<!DOCTYPE html>
<html lang="ja">
<head>
<meta charset="UTF-8">
<title>Dashboard</title>

.framebox {
  display: flex;
    flex-direction: column;  /* ← iframe を縦方向に並べたい場合 */
    align-items: center;     /* ← 横方向の中央に配置 */
}
/* 画面幅に応じて伸縮し、横スクロールが出ないようにする */
.framebox iframe {
    width: 90%;
    aspect-ratio: 16 / 9;   /* ← 高さ400pxの代わりにコレ！ */
    border: 1px solid #ccc;

}
</style>
</head>

<body>
<h2>Dashboard</h2>

<form id="paramForm">
  aval:<input type="text" id="aval" name="aval">
  eval:<input type="text" id="eval" name="eval">
  hval:<input type="text" id="hval" name="hval">
  jval:<input type="text" id="jval" name="jval">
  alpha:<input type="text" id="alpha" name="alpha">
  threshold:<input type="text" id="threshold" name="threshold">
  lval:<input type="text" id="lval" name="lval">
  <br>

Conversion formula:
  aval*(1/c-eval/c2)*10^5+hval-jval*hum　
  RawVal*alpha+(1-alpha)*PrevVal　lval=listCnt
  <br>
  <button type="submit">パラメータ更新</button>
</form>

document.addEventListener("DOMContentLoaded", function() {

fetch("/api/get_params")
    .then(response => response.json())
    .then(data => {

document.getElementById("aval").value  = data.aval  || "1.4";
        document.getElementById("eval").value  = data.eval  || "2.05";
        document.getElementById("hval").value  = data.hval  || "51";
        document.getElementById("jval").value  = data.jval  || "0.6";
        document.getElementById("alpha").value = data.alpha || "0.07";
        document.getElementById("threshold").value = data.threshold || "20";
        document.getElementById("lval").value = data.lval || "45";
      })
    .catch(error => {
        console.error("読み込みエラー:", error);
    });

});

document.getElementById("paramForm").addEventListener("submit", function(e) {
    e.preventDefault();  // フォームの通常送信を止める

const data = {
        aval: document.getElementById("aval").value,
        eval: document.getElementById("eval").value,
        hval: document.getElementById("hval").value,
        jval: document.getElementById("jval").value,
        alpha: document.getElementById("alpha").value,
        threshold: document.getElementById("threshold").value,
        lval: document.getElementById("lval").value
    };

fetch("/api/update_params", {
        method: "POST",
        headers: {
            "Content-Type": "application/json"
        },
        body: JSON.stringify(data)
    })
    .then(response => response.json())
    .then(result => {
        document.getElementById("result").innerHTML =
            "保存結果: " + result.message;
    })
    .catch(error => {
        document.getElementById("result").innerHTML =
            "エラーが発生しました";
        console.error(error);
    });
});
</script>

</body>
</html>
"""

@app.post('/api/update_params')
def update_params():
    data = request.json
    # 保存（例：jsonファイルに保存）
    with open("params.json", "w") as f:
        json.dump(data, f)

return {"message": "パラメータを保存しました"}

@app.get('/api/get_params')
def get_params():
    if not os.path.exists("params.json"):
        return {
            "aval": "",
            "eval": "",
            "hval": "",
            "jval": "",
            "alpha": "",
            "threshold":"",
            "lval": ""

}

with open("params.json", "r") as f:
        return json.load(f)

@app.route('/py/capgraph')
@local_only
def capgraph():
    return """
<!DOCTYPE html>
<html lang="ja">
<head>
<meta charset="UTF-8">
<title>残量グラフ</title>

<script src="https://cdn.jsdelivr.net/npm/chart.js@3.9.1/dist/chart.min.js"></script>

let prevFiltered = null;
let alpha = 0.07;
let threshold = 20;   // 給油判定（例）

// aval*(1/c-eval/c2)*10^5+hval-jval*hum
let aval=1.4;
let eval=2.05;
let hval=51;
let jval=0.6;
let lval="45";

async function loadParams(){
    const response = await fetch("/api/get_params");
    const data = await response.json();

function capToLevel(cap, cap2, ondo, hum) {
    if (cap === 0) return 0.0;

let rawValue = aval * (1/cap - eval/cap2) * 100000 + hval - jval * hum;

// 初回はそのまま
    if (prevFiltered === null) {
        prevFiltered = rawValue;
        return rawValue;
    }

//let alpha = 0.03;

return prevFiltered;
}

function capToLevel2(cap) {
    if (cap !== 0) {
        return 56450.0 / cap;
    } else {
        return 0.0;
    }
}

let capChart = null;

async function loadGraphData() {

prevFiltered = null;

const params = await loadParams();

const lval2= params.lval

const res = await fetch('/py/cap?cap=list'+lval2);
    const data = await res.json();

// x軸ラベル
    const labels = data.map(d => d.datetime.slice(5, 16));  // "12-06 15:22"

loadGraphData();
</script>
</body>
</html>
"""

</script>
</body>
</html>
"""

※その後、測定用のCPSをc1、補正用のCPSをc2とすると、容量がカウントの逆数なので、(1/c1-2.05/c2)*10^5*係数a+切片bで近似すると、比較的うまくいくようだった。残量との差をエクセルで計算して、分散が最小になるように、aとbを調整してみた。この値をグラフにしてみると、ほぼ湿度と逆相関のようだったので、最終的に（上記の値）+51-0.6*湿度とするとさらに近似できそうだった。（このときも分散を見ながら手作業で調整。２次の項も入れてみたがかえってよくなかった。）

　これでも、変動が多いけれど、おそらく空気中の湿度というより、表面の水滴などが影響しているのでないかとも予想？されます。（湿度と相関あっても、完全な相関でないので）そこで、平均すればあるていど変動を抑えられるかと思い、移動平均的な計算を入れることにした。ChatGPTによると、 Value=RawValue*0.07+PrevValue*0.93のように前回データとその時の測定値に重みを付けながら足していくといいらしい。これで、しばらく様子をみることにした。

　今のところ、この移動平均を取り入れたら、以下のグラフ（赤い線）のようになり、わりと実測値に近くなった気がした。

　ChatGPTによると指数移動平均というらしい。低い周波数（残量の変化はかなりゆっくり進むので、あうらしい）を通すフィルターと同じ効果あるのだそう。ある意味LPFのようなものだとか。

※ただ、その後、最初のほうのデータにより、前半の状況が変わってくることもあるようなので、まだ十分とはいえないよう。

※変換式のパラメータをWebで変更できるようにして、実測値と計算値の分散が最小なるようにときどき調整している。残量が少なくなってきたときのほうが、データが多く蓄積してくるので、平均化されて誤差は減っている気もする。大事なのは、灯油がなくなりつつある頃の状況を知ることなので、そういう意味では、なんとか使える？かも。

Myu

2025年12月17日水曜日

ESP32灯油残量センサー　実験

0 件のコメント:

コメントを投稿

ブログアーカイブ

2025年12月17日水曜日

ESP32灯油残量センサー 実験

0 件のコメント:

コメントを投稿

ESP32灯油残量センサー　実験