calculatorNávrh GCP + kontrolasnímky + body → gcp_list s pixely, pak kontrola proklikání
Návrh: nahraj složku snímků z mise a CSV s GCP body (EPSG:5514; Z nepovinné,
chybějící výšky dopočtu z ČÚZK DMR 5G). Vyfiltruju nadir + oblique snímky, spočítám skutečnou projekci bodů do
obrazu (GPS + gimbal úhly z XMP) a vygeneruju gcp_list_ODHAD.txt — 5 nadir + 5 oblique
snímků na bod s předvyplněnými pixely. V GCP Editor Pro ho načteš přes Resume work on GCP
file a jen doladíš kříže. Kontrola: porovná odhad s ručně proklikaným souborem
z editoru a najde řádky, které zůstaly neopravené. Vše běží v prohlížeči, nic se nikam neposílá.
Oblique pitch od do
°
(nadir = pitch ≤ −85°, všechno ostatní se zahodí)
Korekce roviny DSM (z DTM checku): GX
GY
C m
— odečtu GX·(X−X̄)+GY·(Y−Ȳ)+C, X̄/Ȳ = průměr GCP bodů
bod
kandidáti nadir
kandidáti oblique
occlusion −
vybráno
⚠️ Předvyplněné pixely jsou odhad, typicky 200–400 px od skutečné polohy.
Slouží k nalezení snímku a přibližného místa, ne k přesnému určení. Každý řádek musí projít
ruční kontrolou v GCP Editor Pro (Resume work on GCP file). Pixely schválně na 4 desetinná
místa — v editoru je okem odlišíš od ručně kliknutých (ty mají plnou float přesnost).
Checkpointy pojmenuj s prefixem CHK-, WebODM je vynechá z rekonstrukce a použije
na nezávislé reziduum.
bod
snímek
posun [px]
kategorie
Páruje se přes dvojici (bod, snímek). NEOPRAVENO = posun < 2 px — nikdo
netrefí ruční klik do 2 px od odhadu, který se v mediánu mýlí o stovky px; řádek zůstal
nezkontrolovaný a vyhazuje se. PODEZŘELÝ = posun > 800 px — klik na jiný objekt, nebo
špatně odečtený bod; zůstává, rozhodni sám. Limity: bod < 5 pozorování ⚠️, < 3 ❌ (celý bod
se vyhazuje), < 4 control body v setu 🛑 STOP — tři body určují rovinu bez redundance a chyba
se tiše rozpustí do řešení, checkpointy (CHK-) se do počtu nepočítají.
GCP výška z DMR 5GCSV bodů → Z (Bpv) z ČÚZK DMR 5G
Nahraj CSV s body v S-JTSK (Krovak EPSG:5514). Dopočítám výšku terénu (Z, Bpv)
z ČÚZK DMR 5G a vrátím doplněné CSV pro GCP Editor Pro. Vše běží v prohlížeči.
name
X
Y
Z (Bpv)
Zdroj: ČÚZK DMR 5G ImageServer (identify). Přesnost terénu ~0,18 m v otevřeném terénu.
DMR = reliéf/terén — body ber na holé zemi (na střeše dostaneš terén pod ní). Informativní podklad.
Polygon → WebODM boundaryGeoJSON z QGISu → hranice pro WebODM
Nahraj GeoJSON z QGISu (Save Features As → GeoJSON). Vezmu první polygon,
jeho vnější prstenec, a vrátím jediný Feature ve WGS84 (EPSG:4326) jako boundary pro WebODM.
S-JTSK metry (EPSG:5514) reprojektuju automaticky. Vše běží v prohlížeči.
DTM checktransekt + DTM → náklon modelu vůči DMR 5G
Nahraj CSV transekt bodů v S-JTSK (EPSG:5514) a raw DTM GeoTIFF z WebODM
(float32 výšky, export s EPSG:5514). Referenční výšku Z dotáhnu z ČÚZK DMR 5G, DTM navzorkuju
mediánem 0,30 m okna a z trendu delta vs. vzdálenost podél transektu vyhodnotím
náklon modelu: PASS / WARN / FAIL. Alternativa bez TIFFu: sloupec z_dtm v CSV
(ruční sampling v QGIS). Vše běží v prohlížeči.
name
d [m]
Z DMR
Z DTM
delta
delta = Z_DTM − Z_DMR; d = vzdálenost podél hlavní osy transektu. Trend = lineární fit
delta vs. d; náklon přes AOI = trend × rozpětí transektu. PASS < 0,10 m, WARN 0,10–0,20 m,
FAIL > 0,20 m nebo náklon > 2× rozptyl reziduí. Trend s R² < 0,2 je nespolehlivý
(šum/vegetace, ne náklon) — verdikt se pak dává z |mean Δ| stejnými prahy. Body
s Δ > +0,30 m mají hint podezření na vegetaci. DMR 5G ~±0,2 m — informativní kontrola
deformace modelu, ne geodetické měření.
Diff dvou epoch (DTM)dva DTM → plošný posun mezi scany
Nahraj dva raw DTM GeoTIFFy téže lokality z různých epoch (EPSG:5514, float32).
Odečtu je od sebe na společné mřížce a spočítám plošný posun — medián rozdílu výšek.
Dvě epochy můžou obě projít DTM checkem vůči DMR a přesto mezi sebou nést systematický step,
který zabíjí odečitatelnost časové řady. Práh: PASS < 5 cm, WARN 5–10 cm,
FAIL > 10 cm. Heatmapa ukáže, kde posun je. Vše běží v prohlížeči.
posun = Z_novější − Z_starší na společné mřížce (průnik obou rastrů, převzorkováno bilineárně
na hrubší rozlišení z obou). Kladná hodnota = terén stoupl, záporná = klesl. Verdikt z |mediánu|
plošně: PASS < 0,05 m, WARN 0,05–0,10 m, FAIL > 0,10 m. Medián je robustní vůči lokálním změnám
(stavba, vegetace) — chytá systematický plošný offset, ne bodové rozdíly. Heatmapa: modrá = pokles,
červená = růst, kolem nuly bílá; prázdná místa = nodata nebo mimo průnik. Škála ±95. percentil.