def METHOD_NAME():
    """Create map of DWD weather stations_result in Germany"""
    stations = DwdObservationRequest("climate_summary", "daily", "historical")
    stations_df = stations.all().df
    fig, ax = plt.subplots()
    # Rainbow colormap
    cmap = colors.LinearSegmentedColormap.from_list(
        "",
        [
            "#86007D",
            "#0000F9",
            "#008018",
            "#FFFF41",
            "#FFA52C",
            "#FF0018",
        ],
    )
    bounds = (
        stations_df.select(pl.col("height"))
        .to_series()
        .quantile([0, 0.16666667, 0.33333333, 0.5, 0.66666667, 0.83333333, 1])
        .values
    )
    norm = colors.BoundaryNorm(bounds, cmap.N)
    stations_df = stations_df.rename(columns={"height": "Height [m]"})
    stations_df.plot.scatter(x="longitude", y="latitude", c="Height [m]", cmap=cmap, norm=norm, ax=ax)
    ax.set_xlabel("Longitude [°]")
    ax.set_ylabel("Latitude [°]")
    ax.set_title("German weather stations_result")
    ax.text(0.3, 0.05, "Source: Deutscher Wetterdienst", ha="center", va="center", transform=ax.transAxes)
    plt.savefig("german_weather_stations.png")
    return