|
|
| Řádek 1: |
Řádek 1: |
| {{Geoinformatika}} | | {{Geoinformatika}} |
| {{toc|right}} | | {{toc|right}} |
|
| |
| == Výpočet výměr==
| |
|
| |
| K dispozici je textový soubor seznamu souřadnic [http://geo102.fsv.cvut.cz/vyuka/155GIT1/data/body.txt body.txt] ve formátu:
| |
| cb y x
| |
| Vypočítejte výměru parcely ohraničené body seznamu. Postup výpočtu zapište do souboru ''vymera.m''.
| |
|
| |
| Pozn.: Pro výpočet lze použít L´Huillierovy vzorce, viz {{wikipedia| Výpočet plochy pomocí L´Huillierových vzorců}}
| |
|
| |
| <math>
| |
| P = \sum\limits_{i=1}^N \frac{x_i(y_{i+1} - y_{i-1})}{2}
| |
| </math>
| |
|
| |
| nebo:
| |
|
| |
| <math>
| |
| P = \sum\limits_{i=1}^N \frac{y_i(x_{i-1} - x_{i+1})}{2}
| |
| </math>
| |
|
| |
| kde <math>N</math> je počet vrcholů plochy (pro <math>i=1</math> je <math>y_{i-1} = y_N</math> a pro <math>i=N</math> je <math>y_{i+1} = y_1</math>; analogicky pro <math>x</math>).
| |
|
| |
| Pozn.: Pro výpočet výměry musejí být lomové body parcely seřazeny za sebou po obvodu parcely (jinak vyjde výměra špatně), a to správným směrem (jinak vyjde výměra záporně). V podkladových souborech ''body.txt'' a ''body2.txt'' jsou body seřazeny správně.
| |
|
| |
| [[Image:Body.svg|center]]
| |
|
| |
| Dále proveďte:
| |
|
| |
| # Vypočtenou výměru ve skriptu zaokrouhlete na celé metry čtvereční a vypište na obrazovku.
| |
| # Ověřte správnost postupu výpočtu výměr - vytvořte nový datový soubor s několika málo body, které budou ohraničovat oblast, jejíž výměra je evidentní (např. čtverec: 0 0, 0 1, 1 1, 1 0).
| |
| # Doplňte do skriptu test na nedostatečný počet vstupních bodů (pokud počet bodů bude menší než 3, skript vypíše příslušné upozornění a skončí).
| |
|
| |
| Jako rozšíření úlohy lze volitelně provést další výpočet (postup uložte do souboru ''vymery.m'').
| |
|
| |
| K dispozici je textový soubor [http://geo102.fsv.cvut.cz/vyuka/155GIT1/data/body2.txt body2.txt], který obsahuje parcelní čísla a souřadnice lomových bodů několika parcel.
| |
|
| |
| Formát souboru je následující:
| |
|
| |
| <pre>
| |
| čp1 % parcelní číslo
| |
| čb1 y1 x1 % číslo bodu, souřadnice y, x 1.bodu
| |
| čb2 y2 x2
| |
| ...
| |
| čbN yN xN
| |
| prázdný řádek
| |
| další parcela, nebo konec souboru
| |
| </pre>
| |
|
| |
| Pro představu je k dispozici [http://geo102.fsv.cvut.cz/vyuka/155GIT1/data/nacrt.pdf náčrt] (ve formátu pdf, orientace os S-JTSK).
| |
|
| |
| Vypočítejte výměry všech parcel a zapište je do textového souboru ''vymery.txt'' ve formátu:
| |
|
| |
| <pre>
| |
| čp1 výměra1
| |
| čp2 výměra2
| |
| ...
| |
| celková výměra
| |
| </pre>
| |
|
| |
| Pozn.: Výměry zapište do souboru zaokrouhlené na celé metry čtvereční.
| |
|
| |
| Pro kontrolu výpočtu je k dispozici protokol o hromadném výpočtu výměr vytvořený v programu Kokeš: [http://geo102.fsv.cvut.cz/vyuka/155GIT1/data/vymeryKokes.txt vymeryKokes.txt].
| |
