使用Python+OpenCV统计每日用电量(附代码演练)

AI算法与图像处理

共 10288字,需浏览 21分钟

 ·

2021-03-21 20:25

点击下面卡片关注AI算法与图像处理”,选择加"星标"或“置顶”

重磅干货,第一时间送达

Hexing 预付费电表归PLN所有。这个工具可以用来代替必须由操作人员手动记录的模拟仪表。现在用预付的电表,我可以从KWH,安培,电压和我的信用卡中看到更准确的测量。
本文将是一个使用OpenCV库读取预付表测量值的示例。目前,我使用raspberry pi 4来运行我的代码,使用尺寸为4 x 3.1 x 2.5cm的相机。如果你想要扩展USB线,我建议设置USB线的最大长度为3米。超过3米会失去USB信号,导致摄像头无法连接。
将相机面对预付费电表脉冲LED灯:
设置好相机后,如上图所示,就可以开始测量了。

让我们编码!

kwh_pln_sensor.py

我用这个代码实时监控脉冲LED灯
首先,导入所需的库
from sqlalchemy import create_engine
import cv2
import numpy as np
import pandas as pd
import datetime
import time
  • sqlalchemy是连接数据库(对象关系映射器)的助手,在本例中为PostgreSQL
  • cv2(开源的计算机视觉库,或称为OpenCV):这个python库可以帮助从相机中过滤出像素颜色,用于过滤红色
  • pandas使创建矩阵列和行以将数据保存到数据库变得更容易
  • datetimetime,管理日期和时间格式的库
如果发现ModuleNotFoundError错误,则应该首先安装库
pip install opencv-python
pip install sqlalchemy
pip install pandas
创建连接数据库的引擎
conn = 
create_engine("postgresql+psycopg2://root:12345678@localhost:5432/plnstats").connect()
我用root用户名和密码12345678在本地主机服务器端口5432上创建了名为plnstats的数据库
初始化相机并捕获它。
cap = cv2.VideoCapture(0)
(ret, frame) = cap.read()
cv2.VideoCapture(0),你可以根据计算机视频设备索引将索引从0更改为1、2,…等。
cap.read()从相机读取正确的状态帧和图像帧
while True:
(ret, frame) = cap.read()
size = frame.size
image = cv2.cvtColor(frame, cv2.COLOR_BGR2HSV)

将原始帧转换为RGB颜色,image= cv2.cvtColor(frame, cv2.COLOR_BGR2HSV)

为了使用遮罩滤镜,设置RGB的上下限,我尝试确定恒定的RGB并应用代码
lower = np.array([168250])
upper = np.array([179255255])
mask = cv2.inRange(image, lower, upper)
计算红色与非红色的百分比
no_red = cv2.countNonZero(mask)
frac_red = np.divide(float(no_red), size)
percent_red = np.multiply((float(frac_red)), 100)
设置阈值,这取决于相机,我使用阈值10.0
if (percent_red >= 10.0):
保存到pd.DataFrame
data_capture = {
'color_percentage': percent_red,
'created_on': datetime.datetime.now()
}
df = pd.DataFrame(columns=['color_percentage','created_on'],data=data_capture,index=[0])
并保存到数据库中
df.to_sql('home_pln_kwh_sensor', schema='public', con=conn, if_exists='append',index=False)
这是kwh_pln_sensor.py的完整代码
from sqlalchemy import create_engine
import cv2
import numpy as np
import pandas as pd
import datetime
import time

conn = create_engine("postgresql+psycopg2://root:12345678@localhost:5432/plnstats").connect()

cap = cv2.VideoCapture(0)
(ret, frame) = cap.read()

capture_status = True



while capture_status:
    (ret, frame) = cap.read()

    size = frame.size

    image = cv2.cvtColor(frame, cv2.COLOR_BGR2HSV)
    lower = np.array([168250])
    upper = np.array([179255255])
    mask = cv2.inRange(image, lower, upper)

    no_red = cv2.countNonZero(mask)
    frac_red = np.divide(float(no_red), size)
    percent_red = np.multiply((float(frac_red)), 100)

    print(percent_red)    
    if (percent_red >= 10.0): #estimated from the red colour percentage impulse LED lamp
        data_capture = {
            'color_percentage': percent_red,
            'created_on': datetime.datetime.now()
        }

        df = pd.DataFrame(columns=['color_percentage','created_on'],data=data_capture,index=[0])
        df.to_sql('home_pln_kwh_sensor', schema='public', con=conn, if_exists='append',index=False)

cap.release()
cv2.destroyAllWindows()
可以删除第31行中的代码。该语句将帮助你检查相机是否成功捕捉到LED闪烁

那么如何总结传感器捕获?
我想提取关于小时使用情况和信用卡使用的信息。之后,我就可以估计自己在家的使用情况,并在我的信用卡被清空之前被注意到。

kwh_pln_calculate.py

与前面一样,导入所需的库
import pandas as pd
from datetime import datetime,timedelta
from sqlalchemy import create_engine
我用timedelta来得到昨天的日期(你也可以用来计算不同的时间)
创建常量。IMPULSE_KWH = 1000表示每1000个LED闪烁表示1千瓦时,RUPIAH_PER_KWH = 1444.70表示R-1型/TR 1.301-2.200 VA为1444.70卢比
IMPULSE_KWH = 1000 
RUPIAH_PER_KWH = 1444.70
这一天的值和昨天的值,格式为年-月-日
now_date = datetime.now().strftime('%Y-%m-%d')
yesterday_date = datetime.now() - timedelta(days=1)
yesterday_date = yesterday_date.strftime('%Y-%m-%d')
连接PostgreSQL数据库
conn = create_engine("postgresql+psycopg2://root:12345678@localhost:5432/plnstats").connect()
我想总结一下昨天kWh的使用情况
df = pd.read_sql(sql='SELECT color_percentage, created_on from public.home_pln_kwh_sensor where created_on >= \'{yesterday} 00:00:00\' and created_on < \'{currdate} 00:00:00\''.format(yesterday=yesterday_date,currdate=now_date),con=conn)
将created_on列转换为pandas datetime
df['created_on'] = pd.to_datetime(df.created_on)
df['created_on'] = df['created_on'].dt.strftime('%Y-%m-%d %H:%M:%S')
规格化数据,其中一秒只有一个LED闪烁
df = df.drop_duplicates(subset=['created_on'])
每小时写一栏,以获取更多说明
df['created_on'] = pd.to_datetime(df.created_on)
df['hour'] = df['created_on'].dt.strftime('%H')
df['day'] =  df['created_on'].dt.strftime('%Y-%m-%d')
创建一天的清单以检查每天的信用卡使用情况
list_day = list(df['day'].unique())
创建一个新的数据帧进行总结
df_summary = pd.DataFrame(columns=['day','hour','sum_impulse','charges_amount'])
计算
for dday in list_day:
list_hour = list(df[df['day'] == str(dday)]['hour'].unique())
for y in list_hour:
wSumKWH = df[(df['day'] == str(dday)) & (df['hour'] == str(y))]
df_summary = df_summary.append({'hour':y,'day':dday,'sum_impulse':wSumKWH.groupby('hour').count()['day'][0],'charges_amount':(wSumKWH.groupby('hour').count()['day'][0])*RUPIAH_PER_KWH/IMPULSE_KWH},ignore_index=True)
df_summary.to_sql('home_pln_kwh_hour', schema='public', con=conn, if_exists='append', index=False)
credit_summary = df_summary[df_summary['day'] == yesterday_date].groupby('day').sum()
credit_summary.to_sql('home_pln_kwh_summary', schema='public', con=conn, if_exists='append', index=True)
你可以在数据库中查询摘要,也可以像这样从数据帧中提取摘要
图1为传感器每小时的分组,图2为千瓦时*千瓦时价格计算后的总结
然后,你可以看到每小时的最高价格。
这是kwh_pln_calculate.py的完整代码
import pandas as pd
from datetime import datetime,timedelta
from sqlalchemy import create_engine

IMPULSE_KWH = 1000 #depends on prepaid meter
RUPIAH_PER_KWH = 1444.70  # R-1/TR 1.301 – 2.200 VA

def run_cron():
    now_date = datetime.now().strftime('%Y-%m-%d')

    yesterday_date = datetime.now() - timedelta(days=1)
    yesterday_date = yesterday_date.strftime('%Y-%m-%d')

    conn = create_engine("postgresql+psycopg2://root:12345678@localhost:5432/plnstats").connect()

    df = pd.read_sql(sql='SELECT color_percentage, created_on from public.home_pln_kwh_sensor where created_on >= \'{yesterday} 00:00:00\' and created_on < \'{currdate} 00:00:00\''.format(yesterday=yesterday_date,currdate=now_date),con=conn)

    df['created_on'] = pd.to_datetime(df.created_on)

    df['created_on'] = df['created_on'].dt.strftime('%Y-%m-%d %H:%M:%S')
    df = df.drop_duplicates(subset=['created_on'])
    df['created_on'] = pd.to_datetime(df.created_on)
    df['hour'] = df['created_on'].dt.strftime('%H')
    df['day'] =  df['created_on'].dt.strftime('%Y-%m-%d')

    list_day = list(df['day'].unique())

    df_summary = pd.DataFrame(columns=['day','hour','sum_impulse','charges_amount'])

    #DETAIL PER HOUR
    for dday in list_day:
        list_hour = list(df[df['day'] == str(dday)]['hour'].unique())
        for y in list_hour:
            wSumKWH = df[(df['day'] == str(dday)) & (df['hour'] == str(y))]
            df_summary = df_summary.append({'hour':y,'day':dday,'sum_impulse':wSumKWH.groupby('hour').count()['day'][0],'charges_amount':(wSumKWH.groupby('hour').count()['day'][0])*RUPIAH_PER_KWH/IMPULSE_KWH},ignore_index=True)

    df_summary.to_sql('home_pln_kwh_hour', schema='public', con=conn, if_exists='append', index=False)
    credit_summary = df_summary[df_summary['day'] == yesterday_date].groupby('day').sum()

    credit_summary.to_sql('home_pln_kwh_summary', schema='public', con=conn, if_exists='append', index=True)


if __name__ == '__main__':
    run_cron()
谢谢你的阅读,希望你可以用它来衡量每个月的用电量


个人微信(如果没有备注不拉群!
请注明:地区+学校/企业+研究方向+昵称



下载1:何恺明顶会分享


AI算法与图像处理」公众号后台回复:何恺明,即可下载。总共有6份PDF,涉及 ResNet、Mask RCNN等经典工作的总结分析


下载2:终身受益的编程指南:Google编程风格指南


AI算法与图像处理」公众号后台回复:c++,即可下载。历经十年考验,最权威的编程规范!


下载3 CVPR2021

AI算法与图像处公众号后台回复:CVPR即可下载1467篇CVPR 2020论文 和 CVPR 2021 最新论文

点亮 ,告诉大家你也在看



浏览 53
点赞
评论
收藏
分享

手机扫一扫分享

分享
举报
评论
图片
表情
推荐
点赞
评论
收藏
分享

手机扫一扫分享

分享
举报