Ich hab den Troll(o) deswegen gemeldet, Nazi-Polizei-Vergleiche gehen hier gar nicht, finde ich.Ptolemaios hat geschrieben: ↑05.08.2026, 14:45 Bezogen auf meine Person, bin zumindest ich hier kein Gestapospitzel
Dann würde ich vielleicht bald gegen mich selbst sein. Und Ausländer bin ich auch - nämlich in über 150 anderen Staaten... 
Seine Überlegung „bei vier Stücken \(4×R\), also bei 16 entsprechend …“ hat keine geometrische Grundlage; der Zusammenhang ist nicht linear. 
kurz...Code: Alles auswählen
#!/usr/bin/env python3
"""
Erzeugt Weltkarten mit verändertem Meeresspiegel aus ETOPO 2022.
Benötigt:
pip install rasterio numpy pillow
Datei:
etopo.tif
Ausgabe:
flood_10m_klein.png
flood_20m.png
flood_50m.png
"""
import numpy as np
import rasterio
from rasterio.enums import Resampling
from PIL import Image
# ---------------------------------------------------
# Einstellungen
# ---------------------------------------------------
ETOPO = "etopo.tif"
OUTPUT_WIDTH = 3000
OUTPUT_HEIGHT = 1500
SEA_LEVELS = [0, 1, 2, 5, 10, 15, 20, 25, 30, 40, 50, 70, 100, 150, 200, 250, 300, 350, 400, 500, 1000]
# Farben
DEEP_WATER = np.array([25, 70, 170], dtype=np.uint8)
SHALLOW_WATER = np.array([90, 160, 255], dtype=np.uint8)
LOWLAND = np.array([170, 185, 120], dtype=np.uint8)
HIGHLAND = np.array([120, 105, 85], dtype=np.uint8)
MOUNTAIN = np.array([245, 245, 245], dtype=np.uint8)
# ---------------------------------------------------
# Hilfsfunktion
# ---------------------------------------------------
def blend(c1, c2, t):
return (c1 * (1 - t) + c2 * t).astype(np.uint8)
# ---------------------------------------------------
# ETOPO laden
# ---------------------------------------------------
print("Lade ETOPO...")
with rasterio.open(ETOPO) as src:
elev = src.read(
1,
out_shape=(OUTPUT_HEIGHT, OUTPUT_WIDTH),
resampling=Resampling.bilinear,
)
print("Größe:", elev.shape)
# ---------------------------------------------------
# Karten erzeugen
# ---------------------------------------------------
for sea in SEA_LEVELS:
print(f"Erzeuge +{sea} m")
rgb = np.zeros((OUTPUT_HEIGHT, OUTPUT_WIDTH, 3), dtype=np.uint8)
water = elev <= sea
land = ~water
# ---------------------------
# Wasser
# ---------------------------
depth = np.clip((sea - elev) / 6000.0, 0, 1)
rgb[water] = blend(
SHALLOW_WATER,
DEEP_WATER,
depth[water][:, None]
)
# ---------------------------
# Land
# ---------------------------
h = np.clip((elev - sea) / 4000.0, 0, 1)
low = h < 0.5
high = ~low
land_index = np.where(land)
hh = h[land]
colours = np.empty((hh.size, 3), dtype=np.uint8)
a = hh < 0.5
colours[a] = blend(
LOWLAND,
HIGHLAND,
(hh[a] * 2)[:, None]
)
colours[~a] = blend(
HIGHLAND,
MOUNTAIN,
((hh[~a] - 0.5) * 2)[:, None]
)
rgb[land_index] = colours
outfile = f"flood_{sea}m_klein.png"
Image.fromarray(rgb).save(outfile)
print("Gespeichert:", outfile)
print("Fertig.")
