====== Precerpavaci elektrarna Cechy-stred ======

Prehradi-li se Labe na vytoku z CR, muze po osazeni telesa hraze turbinami vzniknout nejvetsi precerpavaci elektrarna v Ceske republice [citation needed]. Dopady stavby na zivotni prostredi jsou zanedbatelne, nebot moznost nezatopena mista osadit dle nadmorske vysky bud vetrniky, nebo repkou, o 666% vylepsi zivotni prostredi v CR. 

[[http://brmlab.s0c4.net/brmdam.tgz|Zdrojaky]] k mergeWENS, hgt2pngn, damnflood


==== Data o nadmorske vysce ====

Nadmorskou vysku pevniny zmapoval a zcenzurovanou verzi zverejnil spojenostatoamericky projekt [[https://en.wikipedia.org/wiki/Shuttle_Radar_Topography_Mission|SRTM]]. Data lze stahnout ve forme matic o rozmerech 1201x1201 16b integeru o kroku 1m; s vertikalni presnosti ~10m (hodnota skalaru) a plosnym rozlisenim 3" (tri uhlove sekundy; ~90m).

Stazeni dat pro CR (zabere ~100MB):
<code>
#!/bin/bash

for i in {48..51..1}
do
  for e in {011..019..1}
  do
    wget http://dds.cr.usgs.gov/srtm/version2_1/SRTM3/Eurasia/N"$i"E"$e".hgt.zip
    unzip ./N"$i"E"$e".hgt.zip -d ./
    rm ./N"$i"E"$e".hgt.zip
  done
done
</code>

==== Zpracovani SRTM dat ====

Existuje projekt gdal, ale je zabugovany k nepouzitelnosti, proto jsem napsal vlastni C programy "mergeWENS.c" (slouci SRTM matice do jedne) a "hgt2pngn.c" (prevede matici na png obrazek). Stazeni: [[http://brmlab.s0c4.net/brmdam.tgz]]

=== Slouceni .hgt matic do jedne velke matice ===

SRTM matice nelze trivialne napojit na sebe, protoze prvni/posledni radek/sloupec se prekryva se sousedni matici. Je tedy potreba pouzit muj mergeWENS.c:

<note tip>
<code>
clang ./mergeWENS.c -o ./mergeWENS
mergeWENS #WE #NS input.hgt [...] output.hgt
</code>

  * #WE = pocet west-to-east matic
  * #NS = pocet north-to-south matic
  * input.hgt ... = vstupni hgt matice v poradi od severozapadu, pres severovychod a jihozapad k jihovychodu.
  * output.hgt = vystupni matice 
</note>

Pro CR a zaokrouhlene okoli (48-51 SS; 11-19 VD) je ten prikaz nasledujici:

<code>
./mergeWENS 9 4 ./N51E011.hgt ./N51E012.hgt ./N51E013.hgt ./N51E014.hgt ./N51E015.hgt ./N51E016.hgt ./N51E017.hgt ./N51E018.hgt ./N51E019.hgt ./N50E011.hgt ./N50E012.hgt ./N50E013.hgt ./N50E014.hgt ./N50E015.hgt ./N50E016.hgt ./N50E017.hgt ./N50E018.hgt ./N50E019.hgt ./N49E011.hgt ./N49E012.hgt ./N49E013.hgt ./N49E014.hgt ./N49E015.hgt ./N49E016.hgt ./N49E017.hgt ./N49E018.hgt ./N49E019.hgt ./N48E011.hgt ./N48E012.hgt ./N48E013.hgt ./N48E014.hgt ./N48E015.hgt ./N48E016.hgt ./N48E017.hgt ./N48E018.hgt ./N48E019.hgt ./NallEall.hgt
</code>

=== Konverze do png ===

Matici je pro potreby zobrazeni prevest do png, opet se opakuje situace se zabugovanymi nastroji v rezozitari/na internetech. Takze pouziju svuj hgt2pngn.c (potrebuje opencv-dev):

<note tip>
<code>
clang $(pkg-config --cflags --libs opencv) hgt2pngn.c -o hgt2pngn
./hgt2pngn width height mode input.hgt output.png
</code>

  * width = pocet sloupcu vstupni matice
  * height = pocet radku vstupni matice
  * mode = rezim konverze vstupniho 16b skalaru:
    * 8 = grayscale 8b png (ztrati se presnost)
    * 16 = grayscale 16b png (bez ztraty presnosti, ale vetsina sw neumi zpracovat)
    * 3 = RGB 8b png.
      * Cervena = oblasti kde nejsou dostupna SRTM data
      * Zelena = (vyska/256)*32
      * modra = vyska%256
  * input.hgt = vstupni matice
  * output.png = vystupni png
</note>

MergeWENS vypise rozmery vystupni matice, coz NENI prosty soucet rozmeru vstupu (kvuli presahum matic do sousednich).

<code>
./hgt2pngn 10793 4798 8 ./NallEall.hgt ./NallEall.png
./hgt2pngn 10793 4798 16 ./NallEall.hgt ./NallEall.png
./hgt2pngn 10793 4798 3 ./NallEall.hgt ./NallEall.png
</code>

==== Vypocet a vykresleni zatopene plochy ====

K vypoctu zatopene plochy je potreba 16bit grayscale PNG z hgt2pngn. O vypocet se stara program damnflood, ktery krom grafickeho vystupu vypocte i nejnizsi bod hraze, jeji vysku, objem zadrzene vody v m^3 a zatopenou plochu v m^2.

<note tip>
<code>
clang $(pkg-config --cflags --libs opencv) -lm damnflood.c -o damnflood
damnflood Ax Ay Bx By h input.png output.png
</code>
  * Ax, Ay = souradnice pocatku hraze
  * Bx, By = souradnice konce hraze
  * h = jake vysky ma hladina dosahnout (v metrech nad morem)
    * 0 znamena vykresleni pouze hraze (pro debugovani umisteni)
  * input.png = 16b grayscale png
  * output.png = 16b RGB png
</note>

Barvy:
  * cervena = teleso hraze
  * modra = zatopena oblast
  * bily bod = chyba mereni SRTM
  * odstiny sede (v podstate cerna) = nethnute informace o terenu z SRTM

Priklad:
<code>
./damnflood 3800 1400 4200 1450 423 ./NallEall16.png ./cechy16.png
Root : Top : Height = 123 : 501 : 378
Filling west/south side!
Stats:
	Volume: 3523475351700 m^3
	Area: 30450289500 m^2
</code>

Vystup:

[[http://brmlab.s0c4.net/flood8.png|8b PNG]] zatopene oblasti (ztrata presnosti 16b->8b skalovanim)

[[http://brmlab.s0c4.net/flood16.png|16b PNG]] zatopene oblasti (puvodni vyska terenu v m.n.m.)


===== PE Morava-stred =====

Projekt Precerpavaci elektrarna Morava-stred je megalomanska odpoved na PECS. Zajimavost: do tela hraze byly zakomponovany vsechny bojove saliny.

<code>
./damnflood 6892 4090 7364 4418 236 ./NallEall16.png ./morava16.png
Root : Top : Height = 142 : 246 : 104
Stats:
	Volume: 349313128788 m^3
	Area: 8896651200 m^2
</code>

Vystup:

[[http://brmlab.s0c4.net/morava8.png|8b PNG]] zatopene oblasti (ztrata presnosti 16b->8b skalovanim)

[[http://brmlab.s0c4.net/morava16.png|16b PNG]] zatopene oblasti (puvodni vyska terenu v m.n.m.)


===== Analýza proveditelnosti =====

Doporučujeme provést přehrazení Labe na území Německa, těsně za českou hranicí. Díky tomu bude stačit v EIA (Environmental Impact Analysis) vyřešit krátký úsek údolí a stavební povolení tak bude snadné získat. Případné negativní externality zatěžující okolní země (v tomto případě zatopení třetiny území) není třeba řešit, jak již ukázala [[http://dfens-cz.com/o-nas-bez-nas-po-druhe-a-opet-v-mnichove/|kauza Energiewende]].

==== Výpočet uložené energie ====

~20000 km^2 s průměrnou hloubkou 100 metrů (nejspíš více) → 2000 km^3 vody ve výšce 100 metrů (nejspíš více) → (2000e12*10*100)/3.6e3 = 550 TWh = **80 Temelínů** (při akumulaci půl roku)

FIXME, navrhujeme numerickou integraci hmotnosti vody * tíhové zrychlení podle výšky hladiny.

==== Další příznivé vlivy ====

Stavba zjednoduší realizaci kanálu Labe-Odra-Dunaj. Kanál bude moct být mnohem kratší.

==== Pilotní projekt ====

V případě, že se nepodaří investory přesvědčit hned ze začátku ke vstupu do celého projektu, doporučujeme jako technology preview postavit přehradu Suchdol-Bohnice.

==== Budoucí rozvoj ====

Pokud se projekt osvědčí, lze pokračovat [[https://en.wikipedia.org/wiki/Atlantropa|vypouštěním Středozemního moře]].