網站首頁 編程語言 正文
1.運行環境及數據
Python3.7、PyCharm Community Edition 2021.1.1,win10系統。
使用的庫:matplotlib、numpy、sklearn、pandas等
數據:CSV文件,包含時間,經緯度,高程等數據
2.基于時間序列的分析2D
讀取時間列和高程做一下分析:
代碼如下:
from PIL import Image import matplotlib.pyplot as plt import numpy as np from sklearn.cluster import KMeans, MiniBatchKMeans import pandas as pd ? if __name__ == "__main__": ? ? data = pd.read_csv(r"H:\CSDN_Test_Data\UseYourTestData.csv") ? ? x, y = data['Time (sec)'], data['Height (m HAE)'] ? ? n = len(x) ? ? x = np.array(x) ? ? x = x.reshape(n, 1)#reshape 為一列 ? ? y = np.array(y) ? ? y = y.reshape(n, 1)#reshape 為一列 ? ? data = np.hstack((x, y)) #水平合并為兩列 ? ? k = 8 ?# 設置顏色聚類的類別個數(我們分別設置8,16,32,64,128進行對比) ? ? cluster = KMeans(n_clusters=k) ?# 構造聚類器 ? ? C = cluster.fit_predict(data) ? ? # C_Image = cluster.fit_predict(data) ? ? print("訓練總耗時為:%s(s)" % (Trainingtime).seconds) ? ? plt.figure() ? ? plt.scatter(data[:, 0], data[:, 1], marker='o', s=2, c=C) ? ? plt.show()
結果展示:
2.1 2000行數據結果展示
2.2 6950行數據結果展示
2.3 300M,約105萬行數據結果展示
CPU立馬90%以上了。大約1-2分鐘,也比較快了。
markersize
有些大了, 將markersize
改小一些顯示,設置為0.1,點太多還是不明顯。
?3.經緯度高程三維坐標分類顯示3D-空間點聚類
修改代碼,讀取相應的列修改為X,Y,Z坐標:如下:
from PIL import Image import matplotlib.pyplot as plt import numpy as np from sklearn.cluster import KMeans, MiniBatchKMeans import pandas as pd from mpl_toolkits.mplot3d import Axes3D ? if __name__ == "__main__": ? ? data = pd.read_csv(r"H:\CSDN_Test_Data\UseYourTestData.csv") ? ? x, y,z = data['Longitude (deg)'],data['Latitude (deg)'], ?data['Height (m HAE)'] ? ? n = len(x) ? ? x = np.array(x) ? ? x = x.reshape(n, 1)#reshape 為一列 ? ? y = np.array(y) ? ? y = y.reshape(n, 1)#reshape 為一列 ? ? z = np.array(z) ? ? z = z.reshape(n, 1) ?# reshape 為一列 ? ? data = np.hstack((x, y, z)) #水平合并為兩列 ? ? k = 8 ?# 設置顏色聚類的類別個數(我們分別設置8,16,32,64,128進行對比) ? ? cluster = KMeans(n_clusters=k) ?# 構造聚類器 ? ? C = cluster.fit_predict(data) ? ? ? # C_Image = cluster.fit_predict(data) ? ? print("訓練總耗時為:%s(s)" % (Trainingtime).seconds) ? ? fig = plt.figure() ? ? ax = Axes3D(fig) ? ? ? ax.scatter(data[:, 0], data[:, 1],data[:, 2], s=1, c=C) ? ? # 繪制圖例 ? ? ax.legend(loc='best') ? ? # 添加坐標軸 ? ? ax.set_zlabel('Z Label', fontdict={'size': 15, 'color': 'red'}) ? ? ax.set_ylabel('Y Label', fontdict={'size': 15, 'color': 'red'}) ? ? ax.set_xlabel('X Label', fontdict={'size': 15, 'color': 'red'}) ? ? plt.show()
3.1 2000行數據結果顯示
由于經度在緯度方向上在17m范圍類,所以立體效果較差,可以換其他數據測試。
3.2 300M的CSV數據計算顯示效果
105萬行數據顯示結果:
原文鏈接:https://blog.csdn.net/soderayer/article/details/123727669
相關推薦
- 2022-09-04 ffmpeg網頁視頻流m3u8?ts實現視頻下載_相關技巧
- 2023-07-26 TypeScript中的聯合類型、類型別名、接口、類型斷言
- 2022-10-01 python3中_from...import...與import?...之間的區別詳解(包/模塊)_
- 2022-04-12 Cannot read property ‘forEach‘ of undefined
- 2022-07-10 詳解BlockingQueue阻塞隊列的使用
- 2022-04-20 Python?設計模式中的創建型建造者模式_python
- 2022-07-28 Python常用Web框架Django、Flask與Tornado介紹_python
- 2023-12-02 vscode無法連接寶塔ftp排雷
- 最近更新
-
- window11 系統安裝 yarn
- 超詳細win安裝深度學習環境2025年最新版(
- Linux 中運行的top命令 怎么退出?
- MySQL 中decimal 的用法? 存儲小
- get 、set 、toString 方法的使
- @Resource和 @Autowired注解
- Java基礎操作-- 運算符,流程控制 Flo
- 1. Int 和Integer 的區別,Jav
- spring @retryable不生效的一種
- Spring Security之認證信息的處理
- Spring Security之認證過濾器
- Spring Security概述快速入門
- Spring Security之配置體系
- 【SpringBoot】SpringCache
- Spring Security之基于方法配置權
- redisson分布式鎖中waittime的設
- maven:解決release錯誤:Artif
- restTemplate使用總結
- Spring Security之安全異常處理
- MybatisPlus優雅實現加密?
- Spring ioc容器與Bean的生命周期。
- 【探索SpringCloud】服務發現-Nac
- Spring Security之基于HttpR
- Redis 底層數據結構-簡單動態字符串(SD
- arthas操作spring被代理目標對象命令
- Spring中的單例模式應用詳解
- 聊聊消息隊列,發送消息的4種方式
- bootspring第三方資源配置管理
- GIT同步修改后的遠程分支