python绘制随机地形地图

在Python中绘制随机地形地图通常涉及使用Perlin噪声或者其他噪声函数生成高度数据，然后使用可视化库如matplotlib或专业的地理信息系统库如geopandas、cartopy来绘制地图。以下是一个简单的例子，使用numpy生成Perlin噪声和matplotlib绘制二维地形图：

import numpy as np
import matplotlib.pyplot as plt
from scipy.ndimage import gaussian_filter

# 生成Perlin噪声
def perlin_noise(width, height, octaves=4, persistence=0.5, lacunarity=2.0):
    x, y = np.mgrid[0:width:1, 0:height:1]
    noise = np.zeros((height, width))
    freq = 1.0
    ampl = 1.0
    for i in range(octaves):
        noise += ampl * np.reshape(np.random.permutation(256), (256, 1)) * (x / freq + y / freq)
        noise *= 2.0 - 2.0 * (np.floor(noise) % 2.0)
        noise = (noise - noise.min()) / (noise.max() - noise.min())
        freq *= lacunarity
        ampl *= persistence
    return noise

# 参数设置
width, height = 800, 800
octaves = 4
persistence = 0.5
lacunarity = 2.0

# 生成噪声数据
noisy_heights = perlin_noise(width, height, octaves, persistence, lacunarity)

# 应用高斯滤波平滑噪声
smooth_heights = gaussian_filter(noisy_heights, sigma=2)

# 将高度值映射到颜色
colors = plt.cm.get_cmap('terrain')(np.linspace(0, 1, smooth_heights.shape[0] * smooth_heights.shape[1])).reshape(smooth_heights.shape[0], smooth_heights.shape[1], 4)

# 绘制地图
plt.imshow(smooth_heights, extent=[0, width, 0, height], cmap='terrain', interpolation='bicubic')
plt.colorbar(label='Height')
plt.title('Random Terrain Map')
plt.show()

这个例子创建了一个Perlin噪声地形图，并使用了高斯滤波器来平滑噪声。terrain颜色映射给地形添加了色彩，使其看起来更像真实的地形。你可以调整octaves、persistence和lacunarity参数来改变地形的细节和复杂性。

请注意，这只是一个基本的二维视图，如果你想要更复杂的地形图，例如带有海岸线、河流和地标的地图，可能需要结合使用其他库和数据源，如geopandas和地理信息数据。