28

我有一个宠物项目来创建地图图像，我在地形高程的等高线图上绘制道路和其他东西。它旨在规划山地自行车路线（过去我手工制作了一些矢量图，它们非常适合可视化）。

目前，我从这里下载 GeoTIFF 中的数字高程模型：http: //www.ecologia.ufrgs.br/labgeo/arquivos/downloads/dados/SRTM/geotiff/rs.rar

然后使用 GDAL 和 Matplotlibcontourf函数创建绘图：

from osgeo import gdal
import matplotlib
import matplotlib.pyplot as plt
from pylab import cm
import numpy

f = 'rs.tif'

elev = gdal.Open(f)

a = elev.GetRasterBand(1).ReadAsArray()

w = elev.RasterXSize
h = elev.RasterYSize
print w, h

altura  = (0.35, 0.42)
largura = (0.70, 0.82)

a = a[int(h*altura[0]):int(h*altura[1]),
      int(w*largura[0]):int(w*largura[1])]


cont = plt.contourf(a, origin='upper', cmap=cm.gist_earth, levels=numpy.arange(0,1000,20))
plt.title('Altitudes - max: %d m; min: %d m' % (numpy.amax(a), numpy.amin(a)))
plt.show()

这使：

在此处输入图像描述

问题是等高线是“白色的”，会产生一些视觉污染，这是不希望的，因为我想稍后绘制道路和河流。

所以，我正在尝试修改contourf创建这些较轻线的方式，或者通过参数设置，或者通过 hack（更改源代码），类似于这里提出的方法：

如何从 Matplotlib 格式化轮廓线

另外，如果有人知道如何使用其他库以更优雅的方式生成这样的地图，我将非常感谢您的提示！

谢谢阅读。

4

4 回答 4

53

我终于找到了解决这个长期存在的问题（目前在 Matplotlib 3 中）的适当解决方案，它不需要多次调用轮廓或光栅化图形。

请注意，问题中说明的问题仅出现在保存的出版质量图形格式（如 PDF）中，而不出现在质量较低的光栅文件（如 PNG）中。

我的解决方案受到这个答案的启发，与颜色条的类似问题有关。类似的解决方案也可以解决等高线图，如下所示：

import numpy as np
import matplotlib.pyplot as plt

np.random.seed(123)
x, y = np.random.uniform(size=(100, 2)).T
z = np.exp(-x**2 - y**2)
levels = np.linspace(0, 1, 100)

cnt = plt.tricontourf(x, y, z, levels=levels, cmap="ocean")

# This is the fix for the white lines between contour levels
for c in cnt.collections:
    c.set_edgecolor("face")

plt.savefig("test.pdf")

下面是修复前的轮廓示例

在此处输入图像描述

下面是上述修复后的相同图

在此处输入图像描述

于 2015-10-02T16:03:06.440 回答
7

尽管不是真正的解决方案，但更多的是一种解决方法，但该问题的解决方案很简单：只需重复相同的contourf命令，这将神奇地消除虚假轮廓。

正如 OP 所述，在进行轮廓填充contourf（我们可以通过设置大量间隔来重现这种行为，例如：

plt.contourf(plon,plat,ssh,np.arange(-1,1.001,0.001)) # 2001 intervals

这给了我们作为输出：

在此处输入图像描述

薄的虚假轮廓明显影响轮廓填充的净颜色。

如果您执行该命令两次：

plt.contourf(plon,plat,ssh,np.arange(-1,1.001,0.001)) # Not once,
plt.contourf(plon,plat,ssh,np.arange(-1,1.001,0.001)) # but twice!

给我：

在此处输入图像描述

现在好多了。这是最好的，有 3 个连续的contourf命令：

在此处输入图像描述

我再也看不到任何细细的轮廓了！不幸的是，这可能会显着减慢您的脚本速度，具体取决于数组大小和等高线间隔的数量。如果使用更多的等高线间隔，则虚假等高线会更加突出。通常最适合我的是使用 50 到 100 个等高线间隔，然后重复contourf两次。

请注意，我使用的 matplotlib 版本不是最新版本。此问题可能已在版本中得到解决1.1.0。如果有，请告诉我。

Python 2.7.1 |EPD 7.0-2 (32-bit)| (r271:86832, Nov 29 2010, 13:52:51)
In [1]: matplotlib.__version__
Out[1]: '1.0.1'
于 2012-05-17T19:59:15.690 回答
6

发布的解决方案不适用于我将 alpha 设置为 <1，但我在此解决方案中找到了答案：Matplotlib Contourf Plots Unwanted Outlines when Alpha < 1

添加论点antialiased=True为我解决了这个问题。

于 2018-03-21T18:17:21.067 回答
4

尝试将 kw 参数添加到plt.contourf(...)call:lw=0或ls=None. http://matplotlib.sourceforge.net/api/pyplot_api.html#matplotlib.pyplot.contourf

于 2011-11-25T10:48:01.063 回答