У меня есть любимый проект по созданию изображений карт, где я рисую дороги и другие объекты на контурной диаграмме рельефа местности. Он предназначен для планирования маршрутов на горных велосипедах (в прошлом я делал несколько векторных рисунков вручную, и они отлично работают для визуализации).

В настоящее время я загружаю Цифровую модель рельефа в 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 создайте эти более легкие строки, либо через установку параметров, либо через хак (изменение исходного кода), аналогично предложенному здесь:

Как отформатировать контурные линии из Matplotlib

Также, если кто-нибудь знает, как создать такую карту более элегантным способом, используя другие библиотеки, я был бы очень признателен за совет!

Спасибо за чтение.