基于Python的地圖繪制教程
本文將介紹通過Python繪制地形圖的方法，所需第三方Python相關模塊包括 rasterio、geopandas、cartopy 等，可通過 pip 等方式安裝。
1 示例代碼
1.1 導入相關模塊
| import rasterio | |
| import geopandas as gpd | |
| import numpy as np | |
| import cartopy.crs as ccrs | |
| import matplotlib.pyplot as plt | |
| from matplotlib.colors import LinearSegmentedColormap |
1.2 設置地圖字體及字號
| plt.rcParams["font.family"] = "Times New Roman" | |
| plt.rcParams["font.size"] = 14 |
1.3 繪制地圖
| # 設置投影：墨卡托投影 | |
| # cartopy 投影說明：https://www.cnblogs.com/youxiaogang/p/14247184.html | |
| projection = ccrs.Mercator() | |
| # 繪制地圖 | |
| fig, ax = plt.subplots(figsize=(20, 10), subplot_kw={'projection': projection}) | |
| # 設置地圖范圍（數值參數為對應投影下的范圍坐標） | |
| ax.set_extent([13530000, 14630000, 4960000, 5850000], crs=projection) | |
| # 讀取矢量文件 | |
| shp = gpd.read_file("Data/Jilin_Mercator.shp") | |
| shp.plot(ax=ax, transform=projection, edgecolor="black", linewidth=1, facecolor="none") | |
| # 創建自定義顏色映射 | |
| colors = ["#369121", "#95C769", "#FFFFBF", "#E6865A", "#D14E30", "#BA1414"] | |
| n_bins = 100 # 定義色帶的顏色數量 | |
| cmap_name = "green_brown" | |
| cm = LinearSegmentedColormap.from_list(cmap_name, colors, N=n_bins) | |
| # 讀取地形柵格數據 | |
| dataset = rasterio.open("Data/DEM_Jilin_Mercator.tif") | |
| data = dataset.read(1) # 讀取第一個波段的數據 | |
| nodata_value = dataset.nodata # 獲取NoData值 | |
| # 創建掩膜，去除NoData區域 | |
| data = np.ma.masked_where(data == nodata_value, data) | |
| # 將地形柵格添加到地圖中 | |
| extent = [dataset.bounds.left, dataset.bounds.right, | |
| dataset.bounds.bottom, dataset.bounds.top] | |
| im = ax.imshow(data, origin="upper", extent=extent, | |
| transform=projection, cmap=cm) | |
| # 繪制網格線并添加標簽 | |
| gl = ax.gridlines(draw_labels=True, linestyle="--", color="#4F4F4F") | |
| gl.xlocator = plt.FixedLocator(range(120, 135, 3)) | |
| gl.ylocator = plt.FixedLocator(range(40, 50, 2)) | |
| # 添加色帶，設置色帶的縮放比例為 0.4，主圖和色帶之間的間距為 0.1 | |
| cbar = plt.colorbar(im, ax=ax, orientation="horizontal", shrink=0.4, pad=0.1) | |
| cbar.set_label("Elevation (m)", labelpad=10) # 設置色帶標簽與色帶的距離為 10 點 | |
| cbar.ax.xaxis.set_label_position("top") # 設置色帶標簽位置 | |
| plt.savefig("Pic.jpg", dpi=600) | |
| plt.show() |
2 結果圖展示
審核編輯：黃飛
聲明：本文內容及配圖由入駐作者撰寫或者入駐合作網站授權轉載。文章觀點僅代表作者本人，不代表電子發燒友網立場。文章及其配圖僅供工程師學習之用，如有內容侵權或者其他違規問題，請聯系本站處理。
舉報投訴
-
python
+關注
關注
56文章
4807瀏覽量
85037 -
地圖繪制
+關注
關注
0文章
4瀏覽量
2374
原文標題：基于Python的地圖繪制教程
文章出處：【微信號：magedu-Linux，微信公眾號：馬哥Linux運維】歡迎添加關注！文章轉載請注明出處。
發布評論請先 登錄
相關推薦
豐田汽車地圖繪制服務比谷歌地圖服務還要好
普及，谷歌地圖也正一步一步走下神壇，而這或許正是 豐田 公司想要得到的結果。 豐田將推出一款全新基于眾包模式的地圖繪制服務，通過眾包模式豐田車行地圖能夠保持不間斷的更新。在豐田看來，谷
路徑規劃算法之柵格地圖繪制
a. 本欄路徑規劃利用矩陣（二維數組）來表示柵格地圖（因為對于矩陣，無論MATLAB、C++還是Python，矩陣更適合數組的表達，更便于編程。）
用Python來繪制世界地圖
Python之所以這么流行，是因為它不僅能夠應用于科技領域，還能用來做許多其他學科的研究工具，繪制地圖便是其功能之一。 今天我們用matplot工具包之一的 mpl_toolkits 來繪制世界
Diagrams ：基于Python繪制云系統架構的模塊
Diagrams 是一個基于Python繪制云系統架構的模塊，它能夠通過非常簡單的描述就能可視化架構，并支持以下6個云產品的圖標： AWS、Azure、GCP、K8s、阿里云 和 Oracle 云
Diagrams：基于Python繪制云系統架構的模塊
Diagrams 是一個基于Python繪制云系統架構的模塊，它能夠通過非常簡單的描述就能可視化架構，并支持以下6個云產品的圖標： AWS、Azure、GCP、K8s、阿里云 和 Oracle 云
Diagrams：基于Python繪制云系統架構的模塊
Diagrams 是一個基于Python繪制云系統架構的模塊，它能夠通過非常簡單的描述就能可視化架構，并支持以下6個云產品的圖標： AWS、Azure、GCP、K8s、阿里云 和 Oracle 云
繪制同切圓python代碼怎么運行
繪制同切圓是一個很有趣的數學問題，可以使用Python語言進行實現。在這篇文章中，我們將探討同切圓的概念、繪制同切圓的算法和Python代碼的實現。 同切圓的概念 同切圓是指具有相同圓
工商網監
評論