利用python進行經緯度軌跡展示

嘿！各位好久不見，距離第一次發部落格已經過去兩年多了，本人也從本科生變成了研究生，好了書歸正傳，最近在做一個關於航班滑行路徑軌跡的專案，目的是將航班的經緯度資料在地圖上顯現出來並生成一條路徑，以方便日後的滑行路徑優化與分析。本文所用的語言為python，使用的是folium包，資料在flightaware網站上可以找到，使用這個包之前還是需要先進行pip install folium

folium的基本用法

1. folium.Map([緯度，經度]，zoom start，width，height，control_scale，no_touch)   #zoom start為地圖的初始縮放尺寸，數值越大地圖放大程度就會越大，width為畫素或百分比(int或者str兩種形式)，height地圖高度，control_scale為是否需要在地圖上顯示比例尺，預設為false，no_touch為是否禁止觸碰地圖或拖拽地圖，預設為不禁止。這些功能都是可自己選擇的，一般來說只需要經度緯度和zoom_start就好了
2. folium.Marker([緯度，經度]).add_to(m)   #在地圖上設定一個標誌符號
3. folium.PolyLine(location = [[緯度1，經度1],[緯度2，經度2]...]，weight，color，opacity).add_to(m)   #在地圖上劃線將所輸入座標連線起來

python利用folium繪圖程式碼

在裝好folium包以後呢，我們就可以使用起來啦，程式碼如下所示：

from folium import plugins
import folium
import os
m = folium.Map([39.1289, 117.3539], zoom_start=10)  #中心區域的確定
location =[[39.1289, 117.3539], [39.1277262, 117.3542938], [39.1277275, 117.3543001], [39.1277262, 117.3542938], 
          [39.1277275, 117.3543001], [39.1277262, 117.3542938], [39.1277262, 117.3542938], 
          [39.1271896, 117.3541359], [39.127121, 117.354126], [39.127121, 117.354126], 
          [39.1269348, 117.3541107], [39.1268692, 117.3541061], [39.1263994, 117.3540649], 
          [39.1257591, 117.3540165], [39.125608, 117.3540192], [39.1251984, 117.3539717], 
          [39.1250038, 117.3539568], [39.1246886, 117.3539276], [39.1246033, 117.3539269], 
          [39.1244316, 117.353912], [39.1242828, 117.353912], [39.1241112, 117.3538971], 
          [39.1238623, 117.3538666], [39.1233153, 117.3538361], [39.1232643, 117.3538374], 
          [39.1230354, 117.3537478], [39.1229895, 117.353714], [39.1228638, 117.3535239], [39.122818, 117.3534493], 
          [39.1227334, 117.353241], [39.1226985, 117.3531494], [39.122652, 117.3530273], [39.122652, 117.3529968], 
          [39.1225821, 117.352829], [39.1225239, 117.3526764], [39.1224861, 117.3525835], [39.1224774, 117.3525391], 
          [39.1224657, 117.3525238], [39.1224174, 117.3523745], [39.1221886, 117.3517625], [39.1221771, 117.3517327], 
          [39.1221399, 117.3516388], [39.1221199, 117.3515834], [39.1220169, 117.3512998], [39.1219769, 117.3512115], 
          [39.1219482, 117.3511057], [39.1219188, 117.3510437], [39.121814, 117.3507996], [39.1217791, 117.3507385], 
          [39.1217558, 117.350708], [39.1215935, 117.3505982], [39.121558, 117.3505859], [39.1213417, 117.3506131], 
          [39.1211014, 117.3507475], [39.121011, 117.3507996], [39.1209528, 117.3508301], [39.120883, 117.3508759], 
          [39.1208481, 117.3509064], [39.1207352, 117.3509714], [39.1204834, 117.3511356], [39.120369, 117.3511953], 
          [39.120369, 117.3511953], [39.1202774, 117.351255], [39.1186409, 117.3522551], [39.1185321, 117.3523254], 
          [39.1182976, 117.3524641], [39.1181374, 117.3525686], [39.1173566, 117.3530426], [39.1171188, 117.3531955], 
          [39.1168213, 117.3533746], [39.1163325, 117.3536682], [39.115867, 117.3539581], [39.1156691, 117.3540802], 
          [39.1156342, 117.3540955], [39.1156342, 117.3540955], [39.1144981, 117.3547927], [39.113551, 117.3553772], 
          [39.113551, 117.3553772], [39.1125069, 117.3559869], [39.1120846, 117.3557129], [39.1120846, 117.3557129], 
          [39.1119118, 117.3553002], [39.1118546, 117.355136], [39.111782, 117.3549652], [39.1115456, 117.3542404], 
          [39.1115958, 117.3540649], [39.1115958, 117.3540649], [39.1128411, 117.3532562], [39.1131786, 117.3530426], 
          [39.113807, 117.3526459], [39.113807, 117.3526459], [39.1191902, 117.3493593], [39.1197308, 117.3490295], 
          [39.1206386, 117.3484802], [39.1208713, 117.3483276], [39.1214676, 117.347971], [39.1214676, 117.347971], 
          [39.1214676, 117.347971], [39.1214676, 117.347971]]   #輸入座標點（注意）folium包要求座標形式以緯度在前，經度在後
route = folium.PolyLine(    #polyline方法為將座標用線段形式連線起來
    location,    #將座標點連線起來
    weight=3,  #線的大小為3
    color='orange',  #線的顏色為橙色
    opacity=0.8    #線的透明度
).add_to(m)    #將這條線新增到剛才的區域m內
m.save(os.path.join(r'C:\Users\Desktop', 'Heatmap1.html'))  #將結果以HTML形式儲存到桌面上
 

 成果展示（這個應該是要在聯網的情況下進行的，有時候影象載入的有一些慢）