# Genera docs/calle.svg: calle pintada que se disuelve, con «presencias» que,
# juntas, forman una figura grande con el puño en alto (a la manera de Suso33).
import sys, random, math
random.seed(26)
OUT = sys.argv[1]
W, H = 1600, 900
VX, VY = 800, 70                 # punto de fuga
INK, ORA = "#111", "#ff7a1a"
o = []
add = o.append

def esc(y):                      # tamaño según la distancia
    return max(0.0, (y - VY) / (H - VY))

def edge(xb, y):                 # x de una recta que va del punto de fuga a (xb, H)
    return VX + (xb - VX) * esc(y)

Y0 = 210                         # lo de arriba queda en blanco: se recorta
add(f'<svg xmlns="http://www.w3.org/2000/svg" viewBox="0 {Y0} {W} {H - Y0}" width="{W}" height="{H - Y0}">')
add('''<defs>
<filter id="brocha" x="-10%" y="-10%" width="120%" height="120%"><feTurbulence type="fractalNoise" baseFrequency="0.015 0.15" numOctaves="3" seed="4" result="t"/><feDisplacementMap in="SourceGraphic" in2="t" scale="6"/></filter>
<filter id="pintura" x="-20%" y="-20%" width="140%" height="140%"><feTurbulence type="fractalNoise" baseFrequency="0.6" numOctaves="2" seed="9" result="t"/><feDisplacementMap in="SourceGraphic" in2="t" scale="2.2"/></filter>
<filter id="sombra" x="-20%" y="-20%" width="140%" height="140%"><feGaussianBlur stdDeviation="1.2"/></filter>
<filter id="grano"><feTurbulence type="fractalNoise" baseFrequency="0.012 0.05" numOctaves="4" seed="2"/><feColorMatrix values="0 0 0 0 0.07  0 0 0 0 0.07  0 0 0 0 0.07  0 0 0 -2.2 1.25"/></filter>
<linearGradient id="fv" x1="0" y1="0" x2="0" y2="1"><stop offset="0" stop-color="#000"/><stop offset=".3" stop-color="#000"/><stop offset=".62" stop-color="#fff"/></linearGradient>
<radialGradient id="fh" cx=".5" cy=".95" r=".75"><stop offset=".55" stop-color="#fff"/><stop offset="1" stop-color="#000"/></radialGradient>
<linearGradient id="gl"><stop offset="0" stop-color="#000"/><stop offset=".12" stop-color="#fff"/><stop offset=".88" stop-color="#fff"/><stop offset="1" stop-color="#000"/></linearGradient>
<mask id="lados" maskUnits="userSpaceOnUse" x="0" y="0" width="1600" height="900"><rect width="1600" height="900" fill="url(#gl)"/></mask>
<mask id="borde" maskUnits="userSpaceOnUse" x="0" y="0" width="1600" height="900"><rect width="1600" height="900" fill="url(#fh)"/></mask>
<mask id="disipa" maskUnits="userSpaceOnUse" x="0" y="0" width="1600" height="900"><g mask="url(#borde)"><rect width="1600" height="900" fill="url(#fv)"/></g></mask>
</defs>''')

# --- calle pintada: todo dentro de la máscara que la disipa ---
add('<g mask="url(#lados)"><g mask="url(#disipa)">')
road = [(VX, VY), (edge(-150, H), H), (edge(1750, H), H)]
add(f'<path d="M{VX},{VY} L-150,{H} L1750,{H}Z" fill="#f1efe9"/>')
def brush(x0, y0, x1, y1, w, n, col=INK, op=.85):
    for _ in range(n):
        j = lambda: random.uniform(-w, w) * .5
        add(f'<line x1="{x0+j():.1f}" y1="{y0+j():.1f}" x2="{x1+j():.1f}" y2="{y1+j():.1f}" '
            f'stroke="{col}" stroke-width="{random.uniform(w*.25, w*.6):.1f}" stroke-linecap="round" opacity="{random.uniform(op*.5, op):.2f}"/>')
add('<g filter="url(#brocha)">')
for xb, w in ((-150, 26), (1750, 26)):                         # bordillos
    y0 = VY + 40
    brush(edge(xb, y0), y0, xb, H + 30, w, 7)
add('</g></g></g>')

# --- figura grande: siluetas sencillas en coordenadas de la imagen ---
def cap(px, py, ax, ay, bx, by, r):              # cápsula (segmento con radio)
    dx, dy = bx - ax, by - ay
    t = max(0, min(1, ((px - ax) * dx + (py - ay) * dy) / (dx * dx + dy * dy)))
    return math.hypot(px - (ax + t * dx), py - (ay + t * dy)) < r
K, YT = .8, 250                                  # escala y altura de la figura
def dentro(x, y):
    x, y = 800 + (x - 800) / K, 60 + (y - YT) / K
    return (math.hypot(x - 772, (y - 205) * .9) < 56 or                  # cabeza
            cap(x, y, 790, 372, 815, 560, 88) or                         # tronco
            cap(x, y, 735, 360, 590, 270, 34) or cap(x, y, 590, 270, 560, 150, 32) or  # brazo en alto
            math.hypot(x - 556, y - 128) < 46 or                         # puño
            cap(x, y, 880, 360, 960, 470, 32) or cap(x, y, 960, 470, 935, 585, 30) or  # brazo bajo
            cap(x, y, 770, 590, 700, 735, 44) or cap(x, y, 700, 735, 650, 870, 40) or  # pierna
            cap(x, y, 860, 590, 940, 720, 44) or cap(x, y, 940, 720, 1010, 860, 40))   # pierna

# --- presencias ---
def presencia(x, y, s, col):
    h = 46 * s                                   # alto de la figura
    r = h * random.uniform(.10, .13)             # cabeza
    lean = random.uniform(-.12, .12) * h
    w = h * random.uniform(.13, .19)             # hombros
    fx = x + lean                                 # pies
    # sombra hacia abajo a la derecha
    L = h * random.uniform(.9, 1.3)
    add(f'<path d="M{fx-w*.25:.1f},{y:.1f} L{fx+L*.85:.1f},{y+L*.42:.1f} L{fx+L*.85+w*.25:.1f},{y+L*.42-w*.2:.1f} L{fx+w*.3:.1f},{y-w*.15:.1f}Z" fill="{INK}" opacity=".22" filter="url(#sombra)"/>')
    sx, sy = x, y - h + 2 * r                     # hombros
    split = random.random() < .35                # piernas abiertas al andar
    if split:
        d = (f'M{sx-w/2:.1f},{sy:.1f} Q{sx-w*.55:.1f},{sy+h*.35:.1f} {fx-w*.32:.1f},{y:.1f} '
             f'L{fx-w*.12:.1f},{y:.1f} L{fx:.1f},{y-h*.3:.1f} L{fx+w*.12:.1f},{y:.1f} L{fx+w*.32:.1f},{y:.1f} '
             f'Q{sx+w*.55:.1f},{sy+h*.35:.1f} {sx+w/2:.1f},{sy:.1f} Z')
    else:
        d = (f'M{sx-w/2:.1f},{sy:.1f} Q{sx-w*.45:.1f},{sy+h*.45:.1f} {fx-w*.1:.1f},{y:.1f} '
             f'L{fx+w*.1:.1f},{y:.1f} Q{sx+w*.45:.1f},{sy+h*.45:.1f} {sx+w/2:.1f},{sy:.1f} Z')
    add(f'<path d="{d}" fill="{col}"/>')
    add(f'<ellipse cx="{sx+lean*.1:.1f}" cy="{sy-r*.9:.1f}" rx="{r:.1f}" ry="{r*1.12:.1f}" fill="{col}"/>')

pts = []
def libre(x, y, s):
    return all(math.hypot(x - a, (y - b) * 1.6) > 17 * max(s, t) for a, b, t in pts)
for _ in range(200000):                            # la figura grande, densa
    x, y = random.uniform(400, 1150), random.uniform(240, 890)
    if dentro(x, y):
        s = .62 + .38 * esc(y)
        if libre(x, y, s * 1.35):
            pts.append((x, y, s))
grande = len(pts)
for _ in range(900):                              # poca gente suelta, por la calle
    y = random.uniform(300, 880)
    x = random.uniform(edge(-150, y), edge(1750, y))
    s = .62 + .38 * esc(y)
    if not dentro(x, y) and libre(x, y, s * 5) and random.random() < .35:
        pts.append((x, y, s))
pts.sort(key=lambda p: p[1])                      # de lejos a cerca
add('<g filter="url(#pintura)" mask="url(#lados)">')
for x, y, s in pts:
    presencia(x, y, s, ORA if random.random() < .05 else INK)
add('</g></svg>')
open(OUT, "w").write("\n".join(o))
print(f"{len(pts)} presencias ({grande} en la figura grande)", file=sys.stderr)