@ohuoy/easymap 1.0.29 → 1.1.0

package/data/sys_log.sql

@@ -0,0 +1,41 @@
+/*
+ Navicat Premium Data Transfer
+
+ Source Server         : 80
+ Source Server Type    : PostgreSQL
+ Source Server Version : 120000 (120000)
+ Source Host           : 192.168.11.80:5432
+ Source Catalog        : sdgydb
+ Source Schema         : public
+
+ Target Server Type    : PostgreSQL
+ Target Server Version : 120000 (120000)
+ File Encoding         : 65001
+
+ Date: 16/09/2025 13:57:56
+*/
+
+
+-- ----------------------------
+-- Table structure for sys_log
+-- ----------------------------
+DROP TABLE IF EXISTS "public"."sys_log";
+CREATE TABLE "public"."sys_log" (
+  "sl_id" int4 NOT NULL DEFAULT nextval('sys_log_sl_id_seq'::regclass),
+  "sl_time" timestamptz(6) NOT NULL,
+  "sl_machine_name" varchar(200) COLLATE "pg_catalog"."default" NOT NULL,
+  "sl_process_name" varchar(200) COLLATE "pg_catalog"."default" NOT NULL,
+  "sl_level" varchar(20) COLLATE "pg_catalog"."default" NOT NULL,
+  "sl_logger" varchar(8000) COLLATE "pg_catalog"."default",
+  "sl_message" varchar(8000) COLLATE "pg_catalog"."default",
+  "sl_call_site" varchar(8000) COLLATE "pg_catalog"."default",
+  "sl_thread" varchar(50) COLLATE "pg_catalog"."default",
+  "sl_exception" varchar(8000) COLLATE "pg_catalog"."default",
+  "sl_stacktrace" text COLLATE "pg_catalog"."default"
+)
+;
+
+-- ----------------------------
+-- Primary Key structure for table sys_log
+-- ----------------------------
+ALTER TABLE "public"."sys_log" ADD CONSTRAINT "sys_log_pkey" PRIMARY KEY ("sl_id");

package/dist/example/create_report.py

@@ -0,0 +1,326 @@
+import numpy as np
+from sklearn.cluster import MeanShift
+import datetime
+import os
+import math
+import json
+
+def read_SN(pols):
+    if pols == '':  # 文件无法读取时，返回空
+        raise FileExistsError('该文件不存在或无法读取')
+    SN=pols[0]['lidarSN']
+    thetime=datetime.datetime.strptime(pols[0]['date'], '%Y-%m-%d %H:%M:%S')
+    DATE=thetime.strftime("%Y%m%d")
+    return SN,DATE
+
+
+#从文件中读取污染源极坐标和时间
+def read_xyz(pols,lng0,lat0):
+    if pols == '':      # 文件无法读取时，返回空
+        raise FileExistsError('该文件不存在或无法读取')
+    #MINDISTANCE=1000
+    allx = []
+    alltime = []
+    ally = []
+    for k in pols:
+        n = 0
+        if 'lng' in k:
+            lng=k['lng']
+        if 'lat' in k:
+            lat=k['lat']
+        if 'Lng' in k:
+            lng=k['Lng']
+        if 'Lat' in k:
+            lat=k['Lat']
+        if 'date' in k:
+            date=k['date']
+        if 'Date' in k:
+            date=k['Date']
+        x,y=lnglat2xy(lng,lat,lng0,lat0)
+        allx.append(x)
+        ally.append(y)
+        atime=datetime.datetime.strptime(date, '%Y-%m-%d %H:%M:%S')
+        alltime.append(atime)
+    #print('国耀量子 版权所有')
+    return allx, ally, alltime
+
+#把污染源空间聚类
+def pts_cluster(x, y, polar=False, bandwidth=500):
+    # 转为np数组
+    x = np.array(x)
+    y = np.array(y)
+    if polar:   # 输入为极坐标数据, x 为r, y为theta, 正北为0, 顺时针.
+        x0 = x * np.sin(np.radians(y))
+        y0 = x * np.cos(np.radians(y))
+        d = np.c_[x0, y0]
+    else:       # 笛卡尔坐标
+        d = np.c_[x, y]
+    #   调用sklearn中的均值偏移函数
+    ms = MeanShift(bandwidth=bandwidth, bin_seeding=True)
+    # 训练数据
+    ms.fit(d)
+    # 数据集中每个点被各自bandwidth内中心点的访问次数来进行分类
+    output_label = ms.labels_
+    # 中心点坐标集合,该集合为数据形式[x,y]
+    cluster_centers = ms.cluster_centers_
+    if polar:   # 输入为极坐标, 输出时将center转为极坐标输出
+        r = np.sqrt(cluster_centers[:, 0]**2 + cluster_centers[:, 1]**2)
+        theta = np.degrees(np.arctan2(cluster_centers[:, 0], cluster_centers[:, 1]))
+        theta[theta < -90] = theta[theta < -90] + 360   # [-180, 180] -> [-90, 270]
+        cluster_centers = np.c_[r, theta].round(3)
+    return output_label, cluster_centers
+
+#把聚类后的污染区域的时间排出来
+def label_time(label,time):
+    combinenumber=max(label)
+    polcombine=[[]for i in range(combinenumber+1)]
+
+    i=0
+    for polpoint in label:
+        polcombine[polpoint].append(time[i])
+        i=i+1
+    i=0
+    for polpoint in polcombine:
+        polcombine[i]=sorted(polpoint)
+        i=i+1
+    polcombine=sorted(polcombine,key=len,reverse=True)
+
+
+
+    return polcombine
+
+#生成每个污染区域的污染时间段
+def pol_time(polcombine,maxdiff=3600,mintime=1800):
+    pn=len(polcombine)
+    poltimequantum = [[] for i in range(pn)]
+    polintensity = [[] for i in range(pn)]
+    polintensity_fr = [[] for i in range(pn)]
+    poltime = [[] for i in range(pn)]
+    poltimequantum_fr = [[] for i in range(pn)]
+    poltime_fr = [[] for i in range(pn)]
+    j=0;
+    for point in polcombine:
+        i=0
+        for atime in point:
+            #poltime_fr[j].append(atime.strftime('%Y%m%d%H%M'))
+            if i==0:
+                k = 1  # 初始化排放量
+                start=atime
+            if i>0:
+                d=atime-lasttime
+                k=k+1
+                if i == len(point) - 1:
+                    if d.seconds<=maxdiff:
+                        lasttime = atime
+                        k=k+1
+
+                    if (lasttime-start).seconds>mintime or k>3:
+
+                        polstart = start-datetime.timedelta(minutes=(start.minute)%5+30)
+                        polstop = lasttime+datetime.timedelta(minutes=5-(lasttime.minute)%5)
+                        poltimequantum[j].append([polstart,polstop])
+                        polintensity[j].append(k-1)
+                    start = atime
+                    k = 1
+
+                elif d.seconds>=maxdiff:
+                    if (lasttime-start).seconds>mintime or k>3:
+                        polstart = start-datetime.timedelta(minutes=(start.minute)%5+30)
+                        polstop = lasttime+datetime.timedelta(minutes=5-(lasttime.minute)%5)
+                        poltimequantum[j].append([polstart,polstop])
+                        polintensity[j].append(k-1)
+
+
+                    start=atime
+                    k = 1
+            i=i+1
+            lasttime = atime
+        poltimequantum[j],poltimequantum_fr[j],polintensity_fr[j]=time_combine(poltimequantum[j],polintensity[j])
+        [poltime[j].append(x) for x in polcombine[j] if time_just(poltimequantum[j],x)]
+        [poltime_fr[j].append(x.strftime('%Y%m%d%H%M%S')) for x in polcombine[j] if time_just(poltimequantum[j],x)]
+        j=j+1
+    return poltimequantum,polintensity_fr,poltime,poltimequantum_fr,poltime_fr
+
+#生成每个污染区域的污染时间段
+def pol_time2(polcombine):
+    pn = len(polcombine)#聚合数量
+    poltimequantum = [[] for i in range(pn)]
+    polintensity = [[] for i in range(pn)]
+    polintensity_fr = [[] for i in range(pn)]
+    poltime = [[] for i in range(pn)]
+    poltimequantum_fr = [[] for i in range(pn)]
+    poltime_fr = [[] for i in range(pn)]
+    j=0
+    for point in polcombine:
+
+        n2=len(point)#聚合内报警次数
+        start=point[0]
+        stop=point[-1]
+        polstart = start-datetime.timedelta(minutes=(start.minute)%5+30)
+        polstop = stop+datetime.timedelta(minutes=5-(stop.minute)%5)
+        poltimequantum[j].append([polstart,polstop])
+        polintensity[j].append(n2)
+        poltimequantum[j],poltimequantum_fr[j],polintensity_fr[j]=time_combine(poltimequantum[j],polintensity[j])
+        [poltime[j].append(x) for x in polcombine[j] if time_just(poltimequantum[j],x)]
+        [poltime_fr[j].append(x.strftime('%Y%m%d%H%M%S')) for x in polcombine[j] if time_just(poltimequantum[j],x)]
+        j=j+1
+    return poltimequantum,polintensity_fr,poltime,poltimequantum_fr,poltime_fr
+
+
+def time_just(timelist,atime):
+    just=0
+    for timeq in timelist:
+        if atime>timeq[0] and atime<timeq[1]:
+            just=1
+    return just
+
+#多个时间段若密集则合并
+def time_combine(timelist,polintensity):
+    n=len(timelist)
+    if n==0:
+        return [],[],[]
+    combine=[timelist[0]]
+    combine_fr=[[timelist[0][0].strftime('%Y%m%d%H%M%S'),timelist[0][1].strftime('%Y%m%d%H%M%S')]]
+    polintensity_fr=[polintensity[0]]
+    if n<2:
+        return combine,combine_fr,polintensity
+    j=0
+    for i in range(n-1):
+        just=time_combine_just(timelist[i],timelist[i+1])
+        if just:
+            combine[j][1]=timelist[i+1][1]
+            combine_fr[j][1]=timelist[i+1][1].strftime('%Y%m%d%H%M%S')
+            polintensity_fr[j]=polintensity_fr[j]+polintensity[i+1]
+        else:
+            combine.append(timelist[i+1])
+            combine_fr.append([timelist[i + 1][0].strftime('%Y%m%d%H%M%S'),timelist[i + 1][1].strftime('%Y%m%d%H%M%S')])
+            polintensity_fr.append(polintensity[i+1])
+            j=j+1
+    return combine,combine_fr,polintensity_fr
+
+
+
+
+#判断是否要合并区间
+def time_combine_just(time1,time2):
+    start1 = time1[0]
+    stop1 = time1[1]
+    start2 = time2[0]
+    stop2 = time2[1]
+    if (stop1-start1+stop2-start2)>0.8*(stop2-start1):
+        return True
+    return False
+
+
+#角度转弧度
+def radian_list(degree):
+    radian=[]
+    for i in degree:
+        radian.append(math.radians(i))
+    return radian
+
+#生成报告
+def final_report(center,poltimequantum,polintensity,poltime,lng0,lat0):
+    i=0
+    finalreport=[]
+
+    for point in center:
+        if poltimequantum[i]:
+            polone={}
+            xy=point.tolist()
+            x=xy[0]
+            y=xy[1]
+            lng,lat=xy2lnglat(x,y,lng0,lat0)
+            polone['x'] = x
+            polone['y'] = y
+            polone['lng']=lng
+            polone['lat']=lat
+            polone['intensity']=polintensity[i]
+            polone['timequantum']=poltimequantum[i]
+            polone['time']=poltime[i]
+            #print(json.dumps(polreport))
+            finalreport.append(polone)
+
+        i=i+1
+
+
+    #with open("write_json.json", "w", encoding='utf-8') as f:
+    #    json.dump(finalreport, f, ensure_ascii=False)
+    #json.dumps(finalreport)
+    return finalreport
+
+
+
+
+
+
+
+# def create_path(ppath):
+#     if os.path.exists(ppath):
+#         pass
+#     else:
+#         os.mkdir(ppath)
+#     return 1
+
+def lnglat2xy(lng,lat,lng0,lat0):
+    #lng0=117.148978#雷达经度
+    #lat0=36.689944#雷达纬度
+    ra = 6378140  # ⾚道半径
+    rb = 6356755  # 极半径
+    #flatten = (ra - rb) / ra  # Partial rate of the earth
+    # change angle to radians
+    radlat = math.radians(lat)
+    radlng = math.radians(lng)
+    radlat0 = math.radians(lat0)
+    radlng0 = math.radians(lng0)
+    x=ra*math.cos(radlat0)*(radlng-radlng0)
+    y=ra*(radlat-radlat0)
+    return x,y
+
+def xy2lnglat(x,y,lng0,lat0):
+    #lng0=117.148978#雷达经度
+    #lat0=36.689944#雷达纬度
+    ra = 6378140  # ⾚道半径
+    radlat0 = math.radians(lat0)
+    lat=lat0+math.degrees(y/ra)
+    lng=lng0+math.degrees(x/ra/math.cos(radlat0))
+    lng6=round(lng,6)
+    lat6=round(lat,6)
+    return lng6,lat6
+
+def save_data(finalreport,filename='write_json.json'):
+
+    with open(filename, "w", encoding='utf-8') as f:
+        json.dump(finalreport, f, ensure_ascii=False)
+    return 1
+
+
+def main_loop(lng0, lat0, windfilter, jsonstr=""):
+    raw = json.loads(jsonstr)
+    allx, ally, time = read_xyz(raw,lng0,lat0)
+    #SN,DATE=read_SN(raw)
+    #print(SN+DATE)
+    alllabel, allcenter = pts_cluster(allx, ally, bandwidth=500, polar=False)
+    polcombine = label_time(alllabel, time)
+    # ax = plt.subplot(111, projection='polar') \
+    # rootpath = 'E:\\version3.7test\POL_report2\\'
+    # lidarpath = rootpath + SN + '\\'
+    # datepath = lidarpath + DATE + '\\'
+    # create_path(rootpath)
+    # create_path(lidarpath)
+    # create_path(datepath)
+
+    poltimequantum, polintensity_fr,poltime, poltimequantum_fr,poltime_fr= pol_time(polcombine)
+    #finalreport = final_report(allcenter, poltimequantum, poltime)
+    finalreport_fr = final_report(allcenter, poltimequantum_fr, polintensity_fr,poltime_fr,lng0,lat0)
+    #printfinalreport=json.dumps(finalreport_fr)
+    #print(printfinalreport)
+    #save_data(finalreport_fr)
+    printfinalreport = json.dumps(finalreport_fr)
+    if len(finalreport_fr) == 0 and not windfilter:
+        poltimequantum, polintensity_fr, poltime, poltimequantum_fr, poltime_fr = pol_time2(polcombine)
+        # finalreport = final_report(allcenter, poltimequantum, poltime)
+        finalreport_fr = final_report(allcenter, poltimequantum_fr, polintensity_fr, poltime_fr,lng0,lat0)
+        return json.dumps(finalreport_fr)        
+    else:
+        return printfinalreport

package/index.js

@@ -124,6 +124,11 @@ export class EasyMap {
         })       
     }
 
+
+    setStyle(url){
+      this.map.setStyle(url)
+      this.map.theme.value = this.map.style.stylesheet.theme
+    }
     addScale(postion="bottom-right"){
       if(this.ScaleControl){
         this.map.addControl(this.ScaleControl,postion);

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@ohuoy/easymap",
-  "version": "1.0.29",
+  "version": "1.1.0",
   "description": "self map easy use",
   "main": "main.js",
   "scripts": {