国产无码综合区,色欲AV无码国产永久播放,无码天堂亚洲国产AV,国产日韩欧美女同一区二区

<dl id="l555p"></dl>

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

2年前作者：洛的地理研學(xué)分類：Toy博客閱讀(17)違法舉報

這篇具有很好參考價值的文章主要介紹了基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）。希望對大家有所幫助。如果存在錯誤或未考慮完全的地方，請大家不吝賜教，您也可以點擊"舉報違法"按鈕提交疑問。

1 gma 簡介

gma 是一個基于 Python 的地理、氣象數(shù)據(jù)快速處理和數(shù)據(jù)分析函數(shù)包（Geographic and Meteorological Analysis，gma）。gma 網(wǎng)站：地理與氣象分析庫。

gma 的主要功能有哪些？

氣候氣象（例如 SPEI、SPI、ET0 等）。
遙感指數(shù)（例如 NDVI、EVI、TVDI 等）。
數(shù)學(xué)運算（例如數(shù)據(jù)平滑、評估、濾波、拉伸、增強變換等）。
系統(tǒng)交互（例如獲取路徑、重命名、壓縮等操作）。
空間雜項（例如計算空間距離、面積計算，坐標(biāo)轉(zhuǎn)換、空間插值等操作）。
柵格處理（例如柵格鑲嵌、裁剪、重采樣、重投影、格式轉(zhuǎn)換、數(shù)據(jù)融合等）。
柵格分析（例如 DEM 坡度、坡向、陰影、等值線等計算）。
矢量處理（例如矢量裁剪、擦除、交集、融合、重投影等）。
地圖工具（例如柵格、矢量數(shù)據(jù)繪圖，指北針、比例尺等生成，坐標(biāo)系定義等）

本文所用的 gma 版本及下載

目前，本文所用的 gma 版本為 gma-1.1.5a15。

系統(tǒng) (X64)： Window 10+
Python 版本： 3.10.9

下載地址：

鏈接：https://pan.baidu.com/s/12wyfOGo0IMsIQIca7iOrxg?pwd=9scj
提取碼：9scj

其他版本可使用 pip install gma 安裝

2 發(fā)展歷程

gma 繪圖功能以 ArcGIS 繪圖功能為參考，目標(biāo)是構(gòu)建一個以 Python 為基礎(chǔ)的快速地理繪圖工具。
其大約與2022 年 11 月疫情居家期間開始構(gòu)建，目前開放（指的是可以直接調(diào)用，且有使用幫助的）功能包括：

plot：主要功能是創(chuàng)建地圖框。
rcs: 為配合繪圖功能，構(gòu)建的擴展空間參考模塊。目的在于初始化已有坐標(biāo)系，并創(chuàng)建自定義坐標(biāo)系。
inres：內(nèi)置資源模塊，主要提供 gma 內(nèi)置的一些矢量數(shù)據(jù)、柵格數(shù)據(jù)，方便底圖添加。

從 1.1.2 版本開始，gma 繪圖功能的主要更新說明如下：

1.1.2：新增 map 地圖工具包，開始支持空間繪圖（基于 gma 的矢量圖層（Layer）和柵格數(shù)據(jù)集類（DataSet））。
1.1.3：添加?xùn)鸥駭?shù)據(jù)集繪制功能（離散型柵格）等。
1.1.4：添加?xùn)鸥駭?shù)據(jù)集繪制功能（重分類或唯一值柵格）等，開始內(nèi)置一些底圖數(shù)據(jù)（世界國家、世界陸地和海洋DEM）。
1.1.5（將在2023年4月的某個周末發(fā)布）：添加矢量要素（Feature）繪制功能，添加更多的地圖數(shù)據(jù)（湖泊、陸地、河流、海洋以及自然地球影像）。

其中也伴隨著眾多 BUG 修復(fù)和功能優(yōu)化，以及一些配套函數(shù)的構(gòu)建，例如（以下函數(shù)均位于gma.io 模塊下）：

CreateFeatureFromPoints：根據(jù)坐標(biāo)點創(chuàng)建 Feature（然后可按照矢量要素的方法添加）
CreateLayerFromFeature：將多個Feature 合并生成 Layer（然后可按照矢量圖層方法添加）
ReadArrayAsDataSet：將數(shù)組讀取為 DataSet（然后可按照柵格數(shù)據(jù)集的方式添加）
TranslateFeatureToDataSet：將矢量要素轉(zhuǎn)換為柵格數(shù)據(jù)集（類似于矢量轉(zhuǎn)柵格）
TranslateLayerToDataSet：將矢量圖層轉(zhuǎn)換為柵格數(shù)據(jù)集（類似于矢量轉(zhuǎn)柵格）

3 主要繪圖思路與步驟

第一步：添加一個地圖框。與 ArcGIS 地圖框類似，后續(xù)所有數(shù)據(jù)均在地圖框內(nèi)添加和繪制。
第二步：添加數(shù)據(jù)（包括矢量圖層、矢量要素、柵格數(shù)據(jù)集等）、并配置參數(shù)（顏色、大小、寬度等）。與ArcGIS 類似，這些數(shù)據(jù)在添加到地圖框后均會被重投影到地圖框配置的底圖坐標(biāo)系上；可配置的參數(shù)類似ArcGIS的符號系統(tǒng)，例如矢量數(shù)據(jù)的背景色、邊界色等，柵格數(shù)據(jù)的色帶（顏色條）、拉伸方式（例如直方圖均衡化）和數(shù)據(jù)變換（例如 gamma 變換）等。
第三步：添加地圖整飾要素（例如：指北針、比例尺和圖例等）。與 ArcGIS 的地圖整飾要素類似，可通過參數(shù)控制顏色、大小、樣式等。

4 主要繪圖功能示例（gma 1.1.5 +）

from gma.map import rcs, plot, inres

4.1 繪制包含國界線的自然地球圖

# 1.初始化一個地圖框，底圖為 WGS84 坐標(biāo)系
MapF = plot.MapFrame(Axes = None, BaseMapProj = 'WGS84', Extent = None)

# 2.將內(nèi)置的世界矢量圖層添加到地圖框
# 2.1 添加內(nèi)置國界（矢量數(shù)據(jù)），并設(shè)置填充色透明，線寬 0.2，線黑色
MapL1 = MapF.AddLayer(inres.WorldLayer.Country, FaceColor = 'none', LineWidth = 0.2, EdgeColor = 'black', Zorder = 1)
# 2.2 添加自然地球底圖（RGB 柵格，不需要配置色帶）
MapD1 = MapF.AddDataSetDiscrete(inres.WorldDataSet.NaturalEarth)

# 3.獲取經(jīng)緯網(wǎng)，但不顯示經(jīng)緯網(wǎng)，為地圖框刻度使用
GridLines = MapF.AddGridLines(LineColor = 'none')

# 4. 設(shè)置地圖框邊框
Frame = MapF.SetFrame()

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

4.2 繪制陸地、海洋、河流和湖泊數(shù)據(jù)

# 1.初始化一個地圖框，底圖為 WGS84 坐標(biāo)系
MapF = plot.MapFrame(Axes = None, BaseMapProj = 'WGS84', Extent = None)

# 2.將內(nèi)置的世界矢量圖層添加到地圖框
MapL1 = MapF.AddLayer(inres.WorldLayer.Land, FaceColor = 'white', LineWidth = 0.2, EdgeColor = 'black')
MapL2 = MapF.AddLayer(inres.WorldLayer.Lake, FaceColor = 'lightblue', LineWidth = 0.2, EdgeColor = '#BEE8FF', Zorder = 2)
MapL3 = MapF.AddLayer(inres.WorldLayer.River, LineColor = 'lightblue', LineWidth = 0.2, Zorder = 3)
MapL4 = MapF.AddLayer(inres.WorldLayer.Ocean, FaceColor = '#BEE8FF', EdgeColor = 'none')

# 3.獲取經(jīng)緯網(wǎng)，但不顯示經(jīng)緯網(wǎng)，為地圖框刻度使用
GridLines = MapF.AddGridLines(LineColor = 'none')

# 4. 設(shè)置地圖框邊框
Frame = MapF.SetFrame()

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

4.3 最好把中國放在地圖中央

# 0 使用 gma 定義一個中心經(jīng)線為 150°E 的簡易圓柱投影（即：改變了中央經(jīng)線的地理坐標(biāo)系）
Proj = rcs.CustomGCS(CentralLongitude = 150)

# 1.初始化一個地圖框，底圖 0 定義的坐標(biāo)系
MapF = plot.MapFrame(Axes = None, BaseMapProj = Proj, Extent = None)

# 2.將內(nèi)置的世界矢量圖層添加到地圖框
MapL1 = MapF.AddLayer(inres.WorldLayer.Land, FaceColor = 'white', LineWidth = 0.2, EdgeColor = 'black')
MapL2 = MapF.AddLayer(inres.WorldLayer.Lake, FaceColor = 'lightblue', LineWidth = 0.2, EdgeColor = '#BEE8FF', Zorder = 2)
MapL3 = MapF.AddLayer(inres.WorldLayer.River, LineColor = 'lightblue', LineWidth = 0.2, Zorder = 3)
MapL4 = MapF.AddLayer(inres.WorldLayer.Ocean, FaceColor = '#BEE8FF', EdgeColor = 'none')

# 3.獲取經(jīng)緯網(wǎng)，但不顯示經(jīng)緯網(wǎng)，為地圖框刻度使用
GridLines = MapF.AddGridLines(LineColor = 'none')

# 4. 設(shè)置地圖框邊框
Frame = MapF.SetFrame()

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

4.4 為添加的柵格數(shù)據(jù)集進行拉伸和變換（數(shù)據(jù)增強）

4.4.1 拉伸數(shù)據(jù)集（以百分比截斷拉伸為例）

# 0 使用 gma 定義一個中心經(jīng)線為 150°E 的簡易圓柱投影（即：改變了中央經(jīng)線的地理坐標(biāo)系）
Proj = rcs.CustomGCS(CentralLongitude = 150)

# 1.初始化一個地圖框，底圖 0 定義的坐標(biāo)系
MapF = plot.MapFrame(Axes = None, BaseMapProj = Proj, Extent = None)

# 2.將內(nèi)置的世界矢量圖層添加到地圖框
# 2.1 添加內(nèi)置國界（矢量數(shù)據(jù)），并設(shè)置填充色透明，線寬 0.2，線黑色
MapL1 = MapF.AddLayer(inres.WorldLayer.Country, FaceColor = 'none', LineWidth = 0.2, EdgeColor = 'black', Zorder = 1)
# 2.2 添加自然地球底圖，并對數(shù)據(jù)進行百分比截斷拉伸（相關(guān)參數(shù)詳見：gma.math.Stretch）
MapD1 = MapF.AddDataSetDiscrete(inres.WorldDataSet.NaturalEarth, Stretch = 'Percentage')

# 3.獲取經(jīng)緯網(wǎng)，但不顯示經(jīng)緯網(wǎng)，為地圖框刻度使用
GridLines = MapF.AddGridLines(LineColor = 'none')

# 4. 添加地圖框邊框
Frame = MapF.SetFrame()

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

4.4.2 變換數(shù)據(jù)集（以 Gamma 變換為例）

# 0 使用 gma 定義一個中心經(jīng)線為 150°E 的簡易圓柱投影（即：改變了中央經(jīng)線的地理坐標(biāo)系）
Proj = rcs.CustomGCS(CentralLongitude = 150)

# 1.初始化一個地圖框，底圖 0 定義的坐標(biāo)系
MapF = plot.MapFrame(Axes = None, BaseMapProj = Proj, Extent = None)

# 2.將內(nèi)置的世界矢量圖層添加到地圖框
# 2.1 添加內(nèi)置國界（矢量數(shù)據(jù)），并設(shè)置填充色透明，線寬 0.2，線黑色
MapL1 = MapF.AddLayer(inres.WorldLayer.Country, FaceColor = 'none', LineWidth = 0.2, EdgeColor = 'black', Zorder = 1)
# 2.2 添加自然地球底圖，并對數(shù)據(jù)進行 gamma 變換，并設(shè)置變換值（相關(guān)參數(shù)詳見：gma.math.Correction）
MapD1 = MapF.AddDataSetDiscrete(inres.WorldDataSet.NaturalEarth, Correction = 'Gamma', CorrectionPROP = dict(GammaV = 2))

# 3.獲取經(jīng)緯網(wǎng)，但不顯示經(jīng)緯網(wǎng)，為地圖框刻度使用
GridLines = MapF.AddGridLines(LineColor = 'none')

# 4. 添加地圖框邊框
Frame = MapF.SetFrame()

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

4.5 繪制單波段數(shù)據(jù)集，并為其配置色帶

# 0 使用 gma 定義一個中心經(jīng)線為 150°E 的簡易圓柱投影（即：改變了中央經(jīng)線的地理坐標(biāo)系）
Proj = rcs.CustomGCS(CentralLongitude = 150)

# 1.初始化一個地圖框，底圖 0 定義的坐標(biāo)系
MapF = plot.MapFrame(Axes = None, BaseMapProj = Proj, Extent = None)

# 2.將內(nèi)置的世界矢量圖層添加到地圖框
# 2.1 添加內(nèi)置國界（矢量數(shù)據(jù)），并設(shè)置填充色透明，線寬 0.2，線黑色
MapL1 = MapF.AddLayer(inres.WorldLayer.Country, FaceColor = 'none', LineWidth = 0.2, EdgeColor = 'black', Zorder = 1)
# 2.2 添加世界陸地和海洋 DEM ，并配置 jet 色帶
MapD1 = MapF.AddDataSetDiscrete(inres.WorldDataSet.DEM, CMap = 'jet')

# 3.獲取經(jīng)緯網(wǎng)，但不顯示經(jīng)緯網(wǎng)，為地圖框刻度使用
GridLines = MapF.AddGridLines(LineColor = 'none')

# 4. 添加地圖框邊框
Frame = MapF.SetFrame()

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

4.6 在地圖上繪制點、線、面

Points = [[112.5, 34.4], [30, 60], [2.3, 48.8], [47.5, -18.9]]
# 創(chuàng)建點、線、面 Feature
FPoints = plot.CreatePlotFeature(Points, Type = 'MultiPoint')
FLine = plot.CreatePlotFeature(Points, Type = 'Line')
FPolygon = plot.CreatePlotFeature(Points, Type = 'Polygon')

# 1.初始化一個地圖框
MapF = plot.MapFrame(Axes = None, BaseMapProj = Proj, Extent = None)

# 2.將內(nèi)置的世界矢量圖層添加到地圖框
MapD1 = MapF.AddDataSetDiscrete(inres.WorldDataSet.NaturalEarth, Zorder = 0)
# 2.2 添加面、線、點
MapL1 = MapF.AddFeature(FPolygon, FaceColor = (0.8, 0.5, 0.5, 0.5), Zorder = 1)
MapL2 = MapF.AddFeature(FLine, LineColor = 'green', LineWidth = 0.5, Zorder = 2)
MapL3 = MapF.AddFeature(FPoints, PointColor = 'red', PointSize = 2, Zorder = 3)

MapF.AddGridLines(LineColor = 'none')
Frame = MapF.SetFrame()

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

4.7 換個視角看一下

Points = [[112.5, 34.4], [30, 60], [2.3, 48.8], [47.5, -18.9]]
# 創(chuàng)建點、線、面 Feature
FPoints = plot.CreatePlotFeature(Points, Type = 'MultiPoint')
FLine = plot.CreatePlotFeature(Points, Type = 'Line')
FPolygon = plot.CreatePlotFeature(Points, Type = 'Polygon')

# 0 使用 gma 定義一個中心經(jīng)線為 112°E 的 Albers 等面積投影
Proj = rcs.AlbersEqualArea(CentralLongitude = 112)

# 1.初始化一個地圖框，底圖 0 定義的坐標(biāo)系
MapF = plot.MapFrame(Axes = None, BaseMapProj = Proj, Extent = None)

# 2.將內(nèi)置的世界矢量圖層添加到地圖框
MapD1 = MapF.AddDataSetDiscrete(inres.WorldDataSet.NaturalEarth, Zorder = 0)
# 2.2 添加面、線、點
MapL1 = MapF.AddFeature(FPolygon, FaceColor = (0.8, 0.5, 0.5, 0.5), Zorder = 1)
MapL2 = MapF.AddFeature(FLine, LineColor = 'green', LineWidth = 0.5, Zorder = 2)
MapL3 = MapF.AddFeature(FPoints, PointColor = 'red', PointSize = 2, Zorder = 3)

MapF.AddGridLines(LineColor = 'none')
Frame = MapF.SetFrame()

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

5 從空間插值到地圖繪制全過程

5.1 從反距離權(quán)重插值到繪圖

示例數(shù)據(jù)下載：

鏈接：https://pan.baidu.com/s/1-lNWdo-V5bsa1r5JS_N4wg?pwd=agzp
提取碼：agzp

5.1.1 讀取數(shù)據(jù)并插值

# import gma
import pandas as pd
from gma import io
from gma.smc import Interpolate
from gma.map import rcs, plot, inres

Data = pd.read_excel("Interpolate.xlsx")
Points = Data.loc[:, ['經(jīng)度','緯度']].values
Values = Data.loc[:, ['值']].values

# 步驟1：反距離權(quán)重插值
IDWD = Interpolate.IDW(Points, Values, Resolution = 0.01, InProjection = 'WGS84')

# 步驟2：將插值結(jié)果轉(zhuǎn)換為 DataSet 數(shù)據(jù)集
IDWDataSet = io.ReadArrayAsDataSet(IDWD.Data, Projection = 'WGS84', Transform = IDWD.Transform)

5.1.2 繪圖（離散型）

# 步驟3：繪圖
# 0 使用 gma 定義一個中心經(jīng)線為 125°E ，兩條標(biāo)準(zhǔn)緯線為 40°N, 52°N 的 Albers 等面積投影
Proj = rcs.AlbersEqualArea(CentralLongitude = 125, StandardParallels = (40, 52))

# 1.初始化一個地圖框，并配置視圖范圍 
MapF = plot.MapFrame(Axes = None, BaseMapProj = Proj, Extent = [112, 38, 138, 54])

# 2.將內(nèi)置的世界矢量圖層添加到地圖框
MapL1 = MapF.AddLayer(inres.WorldLayer.Country, FaceColor = 'none', LineWidth = 0.2, EdgeColor = 'black', Zorder = 2)
MapL2 = MapF.AddLayer(inres.WorldLayer.Ocean, FaceColor = '#BEE8FF', EdgeColor = 'none')

MapD1 = MapF.AddDataSetDiscrete(IDWDataSet, CMap = 'jet', Zorder = 1)

# 3.獲取經(jīng)緯網(wǎng)
GridLines = MapF.AddGridLines(LONRange = (100, 150, 5), LATRange = (30, 60, 5))

# 4.設(shè)置地圖框邊框
Frame = MapF.SetFrame()

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

5.1.3 地圖整飾要素（指北針、比例尺）修飾

# 0 使用 gma 定義一個中心經(jīng)線為 125°E ，兩條標(biāo)準(zhǔn)緯線為 40°N, 52°N 的 Albers 等面積投影
Proj = rcs.AlbersEqualArea(CentralLongitude = 125, StandardParallels = (40, 52))

# 1.初始化一個地圖框，并配置視圖范圍 
MapF = plot.MapFrame(Axes = None, BaseMapProj = Proj, Extent = [112, 38, 138, 54])

# 2.將內(nèi)置的世界矢量圖層添加到地圖框
MapL1 = MapF.AddLayer(inres.WorldLayer.Country, FaceColor = 'none', LineWidth = 0.2, EdgeColor = 'black', Zorder = 2)
MapL2 = MapF.AddLayer(inres.WorldLayer.Ocean, FaceColor = '#BEE8FF', EdgeColor = 'none')

MapD1 = MapF.AddDataSetDiscrete(IDWDataSet, CMap = 'jet', Zorder = 1)

# 3.獲取經(jīng)緯網(wǎng)
GridLines = MapF.AddGridLines(LONRange = (100, 150, 5), LATRange = (30, 60, 5))

# 4.添加地圖整飾要素
AddCompass = MapF.AddCompass(LOC = (0.1, 0.8), Color = 'black')
ScaleBar = MapF.AddScaleBar(LOC = (0.1, 0.02), Width = 0.3, Color = 'black', FontSize = 7)

# 5.設(shè)置地圖框邊框
Frame = MapF.SetFrame()

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

5.1.4 繪制重分類圖并添加圖例

# 0 使用 gma 定義一個中心經(jīng)線為 125°E ，兩條標(biāo)準(zhǔn)緯線為 40°N, 52°N 的 Albers 等面積投影
Proj = rcs.AlbersEqualArea(CentralLongitude = 125, StandardParallels = (40, 52))

# 1.初始化一個地圖框，并配置視圖范圍 
MapF = plot.MapFrame(Axes = None, BaseMapProj = Proj, Extent = [112, 38, 138, 54])

# 2.將內(nèi)置的世界矢量圖層添加到地圖框
MapL1 = MapF.AddLayer(inres.WorldLayer.Country, FaceColor = 'none', LineWidth = 0.2, EdgeColor = 'black', Zorder = 2)
MapL2 = MapF.AddLayer(inres.WorldLayer.Ocean, FaceColor = '#BEE8FF', EdgeColor = 'none')

MapD1 = MapF.AddDataSetClassify(IDWDataSet,
                                CMap = 'jet',
                                Remap = [[-25, 0], [-20, 1], [-15, 2], [-10, 3],[-5, 4], [100, 5]],
                                Method = 'Range', 
                                Labels = ['<= -25', '-25 ~ -20', '-20 ~ -15','-15 ~ -10','-10 ~ -5',' > -5'],
                                Zorder = 1     
                               )

# 3.獲取經(jīng)緯網(wǎng)
GridLines = MapF.AddGridLines(LONRange = (100, 150, 5), LATRange = (30, 60, 5))

# 4.添加地圖整飾要素
AddCompass = MapF.AddCompass(LOC = (0.1, 0.8), Color = 'black')
ScaleBar = MapF.AddScaleBar(LOC = (0.1, 0.02), Width = 0.3, Color = 'black', FontSize = 7)
Legend = MapF.AddLegend(LegendName = '氣溫（℃）', TitleAlignment = 'left', PlotID = [2], LOC = (0.75, 0.5), TitleFont = 'SimSun')

# 5.設(shè)置地圖框邊框
Frame = MapF.SetFrame()

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

5.2 其他插值方法（自然鄰域法、克里金法、趨勢面法、B-樣條函數(shù)法）對比

import matplotlib.pyplot as plt

TD = Interpolate.Trend(Points, Values, Resolution = 0.1)
KD = Interpolate.Kriging(Points, Values, Resolution = 0.1)
NND = Interpolate.NaturalNeighbor(Points, Values, Resolution = 0.1)
BSD = Interpolate.BSpline(Points, Values, Resolution = 0.1)

TDDataSet = io.ReadArrayAsDataSet(TD.Data, Projection = 'WGS84', Transform = TD.Transform)
KDDataSet = io.ReadArrayAsDataSet(KD.Data, Projection = 'WGS84', Transform = KD.Transform)
NNDDataSet = io.ReadArrayAsDataSet(NND.Data, Projection = 'WGS84', Transform = NND.Transform)
BSDDataSet = io.ReadArrayAsDataSet(BSD.Data, Projection = 'WGS84', Transform = BSD.Transform)

# 0 使用 gma 定義一個中心經(jīng)線為 125°E ，兩條標(biāo)準(zhǔn)緯線為 40°N, 52°N 的 Albers 等面積投影
Proj = rcs.AlbersEqualArea(CentralLongitude = 125, StandardParallels = (40, 52))

PlotDS = [TDDataSet, KDDataSet, NNDDataSet, BSDDataSet]

Titles = ['Trend', 'Kriging', 'NaturalNeighbor', 'BSpline']

fig = plt.figure(figsize = (8, 6), dpi = 300)

for i in range(4):
    Axes = plt.subplot(2, 2, i + 1)

    # 1.初始化一個地圖框，并配置視圖范圍 
    MapF = plot.MapFrame(Axes = Axes, BaseMapProj = Proj, Extent = [110, 38, 140, 54])

    # 2.將內(nèi)置的世界矢量圖層添加到地圖框
    MapL1 = MapF.AddLayer(inres.WorldLayer.Country, FaceColor = 'none', LineWidth = 0.2, EdgeColor = 'black', Zorder = 2)
    MapL2 = MapF.AddLayer(inres.WorldLayer.Ocean, FaceColor = '#BEE8FF', EdgeColor = 'none')

    MapD1 = MapF.AddDataSetClassify(PlotDS[i],
                                    CMap = 'jet',
                                    Remap = [[-25, 0], [-20, 1], [-15, 2], [-10, 3],[-5, 4], [100, 5]],
                                    Method = 'Range', 
                                    Labels = ['<= -25', '-25 ~ -20', '-20 ~ -15','-15 ~ -10','-10 ~ -5',' > -5'],
                                    Zorder = 1     
                                   )

    # 3.獲取經(jīng)緯網(wǎng)
    GridLines = MapF.AddGridLines(LONRange = (100, 150, 5), LATRange = (30, 60, 5))

    # 4.添加地圖整飾要素
    AddCompass = MapF.AddCompass(LOC = (0.1, 0.75), Color = 'black')
    ScaleBar = MapF.AddScaleBar(LOC = (0.1, 0.02), Width = 0.3, Color = 'black', FontSize = 7)
    if i == 3:
        Legend = MapF.AddLegend(LegendName = '氣溫（℃）', 
                                TitleAlignment = 'left', 
                                LabelFontSize = 5, 
                                PlotID = [2], 
                                LOC = (0.74, 0.35), 
                                HandleLength = 2.3,
                                HandleHeight = 1.5,
                                TitleFontSize = 7,
                                TitleFont = 'SimSun')

    Frame = MapF.SetFrame()
     
    Axes.set_title(Titles[i], fontsize = 8, y = 0.85, fontdict = {'family':'Times New Roman'})
    
fig.subplots_adjust(wspace = 0.02, hspace = 0.01)

基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

6 其他相關(guān)拓展

6.1 查看當(dāng)前系統(tǒng)所有色帶

DefinedCMaps = plot.GetPreDefinedCMaps(PlotCMap = True)

（部分圖）
基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）

6.2 查看當(dāng)前系統(tǒng)所有字體

SystemFonts = plot.GetSystemFonts(ShowFonts = True)

（部分圖）
基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）文章來源地址http://www.zghlxwxcb.cn/news/detail-407742.html

到了這里，關(guān)于基于 Python 的地理空間繪圖指南（以 gma 空間繪圖功能為例）的文章就介紹完了。如果您還想了解更多內(nèi)容，請在右上角搜索TOY模板網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關(guān)文章，希望大家以后多多支持TOY模板網(wǎng)！

本文來自互聯(lián)網(wǎng)用戶投稿，該文觀點僅代表作者本人，不代表本站立場。本站僅提供信息存儲空間服務(wù)，不擁有所有權(quán)，不承擔(dān)相關(guān)法律責(zé)任。如若轉(zhuǎn)載，請注明出處：如若內(nèi)容造成侵權(quán)/違法違規(guī)/事實不符，請點擊違法舉報進行投訴反饋，一經(jīng)查實，立即刪除！

分享到：

領(lǐng)支付寶紅包贊助服務(wù)器費用

地理空間分析12——地理位置數(shù)據(jù)隱私與安全
在數(shù)字化時代，地理位置數(shù)據(jù)成為了眾多應(yīng)用程序和服務(wù)不可或缺的一部分。從導(dǎo)航應(yīng)用到社交媒體，從廣告定位到城市規(guī)劃，地理位置數(shù)據(jù)的應(yīng)用范圍廣泛。然而，這些數(shù)據(jù)的收集和使用也引發(fā)了廣泛的隱私和安全擔(dān)憂。本文將探討地理位置數(shù)據(jù)隱私的挑戰(zhàn)和重要性，并介
2024年03月19日
瀏覽(28)
matlab做經(jīng)濟地理、地理距離、經(jīng)濟距離空間權(quán)重矩陣
首先講下地理加權(quán)空間權(quán)重矩陣：該矩陣的經(jīng)濟含義是通過不同點的坐標(biāo)系之間的距離遠(yuǎn)近來衡量兩地之間的關(guān)系重要程度，當(dāng)兩點之間距離較遠(yuǎn)，所占的權(quán)重越低，而距離越近，權(quán)重越高。故操作如下：首先需要導(dǎo)入坐標(biāo)數(shù)據(jù)： A=csvread(\\\'JWD.csv\\\',1,0); % JWD.csv是文件名,csvrea
2023年04月12日
瀏覽(24)
Elasticsearch 地理空間計算
? 地理位1置搜索：使用查詢API進行搜索 1、Geo Distance Query Geo Distance Query是在給定的距離范圍內(nèi)搜索數(shù)據(jù)的查詢。例如，可以搜索給定地理位置附近指定距離內(nèi)的所有文檔。 2、?Geo Bounding Box Query Geo Bounding Box Query是一種以四邊形形式搜索地理坐標(biāo)數(shù)據(jù)的查詢。它可以定義頂點
2024年02月12日
瀏覽(28)
Python地理數(shù)據(jù)處理 22：基于arcpy批量操作（四）
代碼描述：遍歷a文件夾下的所有tif影像，并使用每個a文件夾中的tif影像對b文件夾下的所有tif影像進行裁剪。裁剪后的柵格將以兩個tif文件進行組合命名，并保存到另一個文件夾中。獲取柵格數(shù)據(jù)的平均值，并輸出程序運行進度：程序運行進度：某文件夾中包含多個子文件
2024年02月04日
瀏覽(21)
Elasticsearch 地理空間搜索 - 遠(yuǎn)超 OpenSearch
作者：來自 Elastic?Nathan_Reese 2021 年，OpenSearch 和 OpenSearch Dashboards 開始作為 Elasticsearch 和 Kibana 的分支。盡管 OpenSearch 和 OpenSearch Dashboards 具有相似的血統(tǒng)，但它們不提供相同的功能。在分叉時，只能克隆開源許可的功能。這給 OpenSearch 留下了一部分功能。自分叉以來，這種
2024年02月01日
瀏覽(26)
Python地理數(shù)據(jù)處理二十一：基于arcpy批量操作（三）
實現(xiàn)將給定的 .shp 文件中的所有省份作為裁剪范圍，對給定的 .tif 文件進行裁剪，輸出所有省份的單獨 .tif 文件：實現(xiàn)對文件名前14個字符一樣的tif圖像進行柵格運算求和：如：XXXX_XXX_2003.M01_Mean、XXXX_XXX_2003.M02_Mean、XXXX_XXX_2003.M03_Mean；XXXX_XXX_2004.M01_Mean、XXXX_XXX_2004.M02_Mean、
2024年02月01日
瀏覽(22)
kaggle學(xué)習(xí)筆記-情感和地理空間分析
秘魯食品評論中的情緒和地理空間分析自然語言處理 (NLP) 是人工智能的一個分支，致力于讓計算機能夠像人類一樣理解文本和口語單詞。另一方面，地理空間分析是對圖像、GPS、衛(wèi)星攝影和歷史數(shù)據(jù)的收集、顯示和操作，這些數(shù)據(jù)以地理坐標(biāo)明確描述，或以街道地址、郵政
2024年02月16日
瀏覽(19)
Redis Geo：掌握地理空間數(shù)據(jù)的藝術(shù)
歡迎來到我的博客，代碼的世界里，每一行都是一個故事在移動互聯(lián)網(wǎng)和物聯(lián)網(wǎng)的時代，地理位置數(shù)據(jù)無處不在。從導(dǎo)航和配送到社交網(wǎng)絡(luò)和廣告，地理位置信息正在重新定義我們與世界的互動方式。但是，處理和分析這些大量的地理空間數(shù)據(jù)絕非易事。這時，Redis Geo應(yīng)運而
2024年02月02日
瀏覽(21)
地理信息系統(tǒng)原理-空間數(shù)據(jù)結(jié)構(gòu)（7）
?四叉樹編碼 1.四叉樹編碼定義四叉樹數(shù)據(jù)結(jié)構(gòu)是一種對柵格數(shù)據(jù)的壓縮編碼方法，其基本思想是將一幅柵格數(shù)據(jù)層或圖像等分為四部分，逐塊檢查其格網(wǎng)屬性值（或灰度）；如果某個子區(qū)的所有格網(wǎng)值都具有相同的值，則這個子區(qū)就不再分割，否則還要把這個子區(qū)再分割成
2024年02月06日
瀏覽(28)
HarmonyOS鴻蒙開發(fā)指南:基于ArkTS的聲明式開發(fā)范式聲明式UI開發(fā)實例繪圖與動畫
目錄繪制圖形繪制基本幾何圖形繪制自定義幾何圖形 animateTo實現(xiàn)閃屏動畫頁面轉(zhuǎn)場動畫繪制能力主要是通過框架提供的繪制組件來支撐，支持svg標(biāo)準(zhǔn)繪制命令。本節(jié)主要學(xué)習(xí)如何使用繪制組件，繪制詳情頁食物成分標(biāo)簽（基本幾何圖形）和應(yīng)用Logo（自定義圖形）。
2024年01月17日
瀏覽(26)

<strike id="zvsmc"><thead id="zvsmc"><li id="zvsmc"></li></thead></strike>

<menu id="zvsmc"></menu><dl id="zvsmc"><nav id="zvsmc"></nav></dl>