不使用 Matplotlib 绘制的线/轨迹答案

【问题标题】：Lines/Tracks Not Plotting with Matplotlib不使用 Matplotlib 绘制的线/轨迹
【发布时间】：2017-08-09 19:39:56
【问题描述】：

我正在尝试使用matplotlib.在地图上绘制 CSV 文件中的线条和标记

数据：

AL99,2017080912,SHIP,0,17.1,-55.6,25,0
AL99,2017080912,SHIP,12,18.1,-57.6,27,0
AL99,2017080912,SHIP,24,19.0,-59.2,29,0
AL99,2017080912,SHIP,36,20.1,-60.2,34,0
AL99,2017080912,SHIP,48,21.5,-61.6,39,0
AL99,2017080912,SHIP,60,23.3,-63.0,47,0
AL99,2017080912,SHIP,72,25.4,-65.2,54,0
AL99,2017080912,SHIP,84,27.9,-68.1,61,0
AL99,2017080912,TABD,0,17.1,-55.7,0,0
AL99,2017080912,TABD,6,17.5,-56.7,0,0
AL99,2017080912,TABD,12,17.8,-57.3,0,0
AL99,2017080912,TABD,18,18.1,-57.9,0,0
AL99,2017080912,TABD,24,18.5,-58.3,0,0
AL99,2017080912,TABD,30,19.0,-58.6,0,0
AL99,2017080912,TABD,36,19.6,-58.8,0,0

Python 代码：

tc = np.recfromcsv(csv_file, unpack=True, names=['stormid', 'initdate', 'mems', 'times', 'tclat', 'tclon', 'tcwind', 'tcpres'], dtype=None)
for j in range(len(tc.times)):
lon, lat = tc.tclon[j], tc.tclat[j]
xpt, ypt = m(lon, lat)
lonpt, latpt = m(xpt, ypt, inverse=True)

if tc.mems[j] == 'TABD':
    tccolor = '--bo'
elif tc.mems[j] == 'AEMN':
    tccolor = '-ro'
else:
    tccolor = '-k'

m.plot(xpt, ypt, tccolor)

结果：

我正在让标记正确绘制颜色，但是线条不存在。

【问题讨论】：

提醒一下，您的代码不完整（m 是什么？）并且格式不正确（for 循环后没有缩进）。此外，您的数据似乎与您的图不匹配（例如，您的数据中的lat 仅介于 15 到 20 之间；但您的图看起来像在大西洋上，纬度至少高达 40）。因此，其他人可能难以解决您的问题并提供帮助。
不要指望在不提供完整代码时得到答案。除非你非常幸运。通常，愿意回答的人会付出比必要更多的努力。

标签： python csv matplotlib matplotlib-basemap

【解决方案1】：

要启用线图，plot() 函数必须提供多点 (x, y) 数组，而不是像您所做的那样提供单个值。这是一个工作代码和生成的输出图（基于您提供的数据）。

from matplotlib import pyplot as plt
from mpl_toolkits.basemap import Basemap
import numpy as np

fig = plt.gcf()
fig.set_size_inches([8,8])

m = Basemap(projection='cyl', \
        llcrnrlat= 10, urcrnrlat= 30, \
        llcrnrlon= -70, urcrnrlon= -50, \
        resolution='l')

csv_file = "storm_data.csv"
tc = np.recfromcsv(csv_file, unpack=True, names=['stormid', 'initdate', 'mems', 'times', 'tclat', 'tclon', 'tcwind', 'tcpres'], dtype=None)

# for each line segment, (xs, ys) is initialized here
# two of them are done here for demo purposes 
xs1 = []
ys1 = []
xs2 = []
ys2 = []

for j in range(len(tc.times)):
    lon, lat = tc.tclon[j], tc.tclat[j]
    xpt, ypt = m(lon, lat)
    lonpt, latpt = m(xpt, ypt, inverse=True)

    if tc.mems[j] == 'TABD':
        xs1.append(lon)
        ys1.append(lat)

    elif tc.mems[j] == 'SHIP':
        xs2.append(lon)
        ys2.append(lat)

    # *** elif for other line segmems ***

    else:
        pass


# plot the collected line segments
m.plot( xs1, ys1, '--bo', xs2, ys2, '--ro' )

# draw coastline
m.drawcoastlines()

plt.show()

结果图：

【讨论】：