Python 用底图绘制小河

我试图用Python中的Basemap模块显示特拉华河。然而，我能找到的最近的河流是延伸到特拉华湾的河流的较低部分。我试过了

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

m = Basemap(width=120000,height=90000,projection='stere',resolution='f',lat_1=39.82,lat_2=39.88,lat_0=39.86,lon_0=-75.23694)
m.drawcoastlines()
m.drawrivers()
plt.show()

---

但这只提供了河流的较低部分和部分河流的中心线。有没有办法从Basemap获得特拉华河大部分地区的详细地图，或者我需要导入该河的形状文件？

当您需要地图上的特定数据要素，而地图数据库没有提供这些要素时，您需要从某处获取这些要素并投入使用。在这里，我通过在

GG My Maps

上绘图获得数据（长，纬度）。在地图上绘制河流路径沿线后，我可以将线数据导出为

KML

。从该KML文件中，我只获取（long，lat，h）列表，并进行一些编辑，以获取以下代码中

lons

和

lats

的数据

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

lons = np.array([-74.8739737, -74.8633163, -74.8585042, -74.8564394, -74.8524035,
       -74.8459655, -74.830347 , -74.7874564, -74.7661703, -74.7593039,
       -74.7531241, -74.7345847, -74.7235983, -74.724285 , -74.7380179,
       -74.7469443, -74.7641104, -74.7826498, -74.8005026, -74.8142355,
       -74.8259085, -74.832775 , -74.8355215, -74.852001 , -74.8623007,
       -74.8698538, -74.9014395, -74.9158591, -74.9302786, -74.9453848,
       -74.9701041, -74.9982565, -75.0140494, -75.0270956, -75.047695 ,
       -75.0552481, -75.0655478, -75.0731009, -75.0916403, -75.111553 ,
       -75.1294058, -75.1369589, -75.1383322, -75.1300925, -75.1300925,
       -75.1355856, -75.1403921, -75.1695615, -75.1971884, -75.2130589,
       -75.2385357, -75.2697054, -75.3127467, -75.3521119, -75.3813978,
       -75.3987804, -75.4230237, -75.4458246])

lats = np.array([40.3023524, 40.2900427, 40.2829709, 40.2758991, 40.2654202,
       40.2567749, 40.2447233, 40.2248056, 40.2127464, 40.1959649,
       40.1828514, 40.1739328, 40.1597655, 40.1471699, 40.1392965,
       40.1356219, 40.134047 , 40.1219717, 40.1219717, 40.1282721,
       40.1277471, 40.1203965, 40.0988651, 40.0941378, 40.0878342,
       40.0825808, 40.0725981, 40.0704963, 40.0715472, 40.062614 ,
       40.050526 , 40.0310756, 40.0221371, 40.0179303, 40.0079382,
       40.0005746, 39.9884756, 39.9795315, 39.9742697, 39.9711125,
       39.9616398, 39.9542712, 39.9321608, 39.9179432, 39.9031959,
       39.8937139, 39.8847574, 39.8809411, 39.8785011, 39.8662113,
       39.8560259, 39.8505758, 39.8494219, 39.8388161, 39.8197959,
       39.8092274, 39.8007649, 39.7931383])

fig = plt.figure(figsize=[8,8])
ax = plt.gca()

m = Basemap(width=120000, height=90000, projection='stere', resolution='f', \
            lat_1=39.82, lat_2=39.88, lat_0=39.86, lon_0=-75.23694)

m.drawcoastlines()
# m.drawrivers()  # no needed now

# plot Delaware river using arrays of data
m.plot( *m(lons, lats), linewidth=2.5, color='blue', zorder=18)
plt.show()

结果图：

---

当您需要地图上的特定数据要素，而地图数据库没有提供这些要素时，您需要从某处获取这些要素并投入使用。在这里，我通过在

GG My Maps

上绘图获得数据（长，纬度）。在地图上绘制河流路径沿线后，我可以将线数据导出为

KML

。从该KML文件中，我只获取（long，lat，h）列表，并进行一些编辑，以获取以下代码中

lons

和

lats

的数据

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

lons = np.array([-74.8739737, -74.8633163, -74.8585042, -74.8564394, -74.8524035,
       -74.8459655, -74.830347 , -74.7874564, -74.7661703, -74.7593039,
       -74.7531241, -74.7345847, -74.7235983, -74.724285 , -74.7380179,
       -74.7469443, -74.7641104, -74.7826498, -74.8005026, -74.8142355,
       -74.8259085, -74.832775 , -74.8355215, -74.852001 , -74.8623007,
       -74.8698538, -74.9014395, -74.9158591, -74.9302786, -74.9453848,
       -74.9701041, -74.9982565, -75.0140494, -75.0270956, -75.047695 ,
       -75.0552481, -75.0655478, -75.0731009, -75.0916403, -75.111553 ,
       -75.1294058, -75.1369589, -75.1383322, -75.1300925, -75.1300925,
       -75.1355856, -75.1403921, -75.1695615, -75.1971884, -75.2130589,
       -75.2385357, -75.2697054, -75.3127467, -75.3521119, -75.3813978,
       -75.3987804, -75.4230237, -75.4458246])

lats = np.array([40.3023524, 40.2900427, 40.2829709, 40.2758991, 40.2654202,
       40.2567749, 40.2447233, 40.2248056, 40.2127464, 40.1959649,
       40.1828514, 40.1739328, 40.1597655, 40.1471699, 40.1392965,
       40.1356219, 40.134047 , 40.1219717, 40.1219717, 40.1282721,
       40.1277471, 40.1203965, 40.0988651, 40.0941378, 40.0878342,
       40.0825808, 40.0725981, 40.0704963, 40.0715472, 40.062614 ,
       40.050526 , 40.0310756, 40.0221371, 40.0179303, 40.0079382,
       40.0005746, 39.9884756, 39.9795315, 39.9742697, 39.9711125,
       39.9616398, 39.9542712, 39.9321608, 39.9179432, 39.9031959,
       39.8937139, 39.8847574, 39.8809411, 39.8785011, 39.8662113,
       39.8560259, 39.8505758, 39.8494219, 39.8388161, 39.8197959,
       39.8092274, 39.8007649, 39.7931383])

fig = plt.figure(figsize=[8,8])
ax = plt.gca()

m = Basemap(width=120000, height=90000, projection='stere', resolution='f', \
            lat_1=39.82, lat_2=39.88, lat_0=39.86, lon_0=-75.23694)

m.drawcoastlines()
# m.drawrivers()  # no needed now

# plot Delaware river using arrays of data
m.plot( *m(lons, lats), linewidth=2.5, color='blue', zorder=18)
plt.show()

结果图：

如果使用全分辨率（

“f”

）但未看到某个功能，则该功能根本不是数据的一部分。您可以使用自己的数据来绘制该特征。如果您使用全分辨率（

“f”

）并且未看到某个特征，则该特征根本不是数据的一部分。您可以使用自己的数据绘制该特征。