[Talk-cz] Výšky budov
Od kud se dají získat výšky budov? Nyní jsou někde .27 což je naprostá hovadina asi si to někdo spletl s výměrou nebo co. Nebo tři vyznačené stejné budovy na sobě vykreslené pokaždé s jinou výškou. Doplněk pointInfo (ver.30470) mě to nezobrazí. S díky za odpověď Lukáš ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Výšky budov
IMHO RUIAN - počet pater J Dne 25. července 2014 19:02 Lukas Novotny lenoc...@tiscali.cz napsal(a): Od kud se dají získat výšky budov? Nyní jsou někde .27 což je naprostá hovadina asi si to někdo spletl s výměrou nebo co. Nebo tři vyznačené stejné budovy na sobě vykreslené pokaždé s jinou výškou. Doplněk pointInfo (ver.30470) mě to nezobrazí. S díky za odpověď Lukáš ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz -- Jachym Cepicky e-mail: jachym.cepicky gmail com URL: http://les-ejk.cz GPG: http://les-ejk.cz/pgp/JachymCepicky.pgp Give your code freedom with PyWPS - http://pywps.wald.intevation.org ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Výšky budov
0.27 pater? Nějaké divné ne. Nechápu tu poznámku k ruian. V něm žádná výška budovy není. Pouze počet pater a ty do tagu height=* nedáváme. Netuším, jestli někdo na základě počtu pater zkouší odhadnout výšku budovy, ale bude to jen hodně, hodně přibližné. Osobně bych to nedělal. Marián On 25. července 2014 19:29:55 CEST, Jachym Cepicky jachym.cepi...@gmail.com wrote: IMHO RUIAN - počet pater J Dne 25. července 2014 19:02 Lukas Novotny lenoc...@tiscali.cz napsal(a): Od kud se dají získat výšky budov? Nyní jsou někde .27 což je naprostá hovadina asi si to někdo spletl s výměrou nebo co. Nebo tři vyznačené stejné budovy na sobě vykreslené pokaždé s jinou výškou. Doplněk pointInfo (ver.30470) mě to nezobrazí. S díky za odpověď Lukáš ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz -- Odesláno z mého telefonu s Androidem pomocí pošty K-9 Mail. Omluvte prosím moji stručnost. ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
[Talk-cz] Skript na parcely RUIANu
Ahoj,
obrázky parcel od Petra Vejsady mi nějak nedaly spát, tak jsem včera v
noci ubastlil skript na konverzi parcel do OSM.
Skriptík nemá jiné ambice než vyzkoušet si natažení parcel do JOSM a
vůbec se seznámit s VFR formátem RUIANu, takže nečekejte zázraky. Mj.
nefungují multipolygony, beru jen polygony typu gml:LinearRing apod.
Tagování jsem narychlo opsal z
http://wiki.openstreetmap.org/wiki/Cs:RUIAN, knihovny na zpracování
GML apod. jsem zatím nehledal.
Vyžaduje perl, proj, Geo::Proj4
(http://search.cpan.org/dist/Geo-Proj4/), a grid jezek_czech08.llb
přejmenovaný na czech (vyšťouráno tady v archívu,
http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid).
Pár obrázků:
http://www.maatts.cz/obrazek/3/lomnice-ruian-parcely-png/
http://www.maatts.cz/obrazek/3/blansko-ruian-parcely-png/
http://www.maatts.cz/obrazek/3/brno-ruian-parcely-png/
Vyzkoušeno na uvedených třech obcích = na čtvrté klidně může
zhavarovat. ;-)
BACHA - ve velkých městech skript žere mraky paměti, optimalizacema
jsem se netrápil. :-)) Třeba Brno (kód obce 582786) má ve výsledku
1.122.529 uzlů a 251.328 parcel, a skriptík na to potřebuje aspoň 12GB
RAM a/nebo hodně swapu na SSD.
Martin
#!/usr/bin/perl -w
use strict;
use Geo::Proj4;
use XML::Simple;
use Data::Dumper;
#
# !!! NEPOUZIVAT PRO OSTRE MAPOVANI DO OSM, POKUSNA ALFA VERZE !!!
#
#
# Pouziti:
#
# (1) Stahnout Kompletni datovou sadu vybrane oblasti z http://vdp.cuzk.cz/vdp/ruian/vymennyformat/vyhledej, napr. podle kodu
# 581976 = Lomnice u Tisnova
# 581283 = Blansko
# 582786 = Brno (velke!)
#
# (2) gunzip 20140630_OB_581976_UKSH.xml.gz
#
# (3) par2osm.pl 20140630_OB_581976_UKSH.xml lomnice.osm
#
#
#
my $input_file = $ARGV[0];
die No input file specified unless (defined $input_file);
die Input file not found unless (-f $input_file);
my $epsg_srs_name = 'urn:ogc:def:crs:EPSG::5514';
my $proj_from = Geo::Proj4-new('+proj=krovak +ellps=bessel +nadgrids=czech');
my $proj_to = Geo::Proj4-new('+proj=longlat +datum=WGS84');
die Undefined source projection unless (defined $proj_to);
die Undefined destination projection unless (defined $proj_from);
my $nodemap = {};
my $idgen = -1;
my $druh_pozemku_map = {
# 1 = [] ???
2 = [ 'landuse', 'farmland' ],
3 = [ 'landuse', 'hop_garden' ],
4 = [ 'landuse', 'vineyard' ],
5 = [ 'leisure', 'garden' ],
6 = [ 'landuse', 'orchard' ],
7 = [ 'landuse', 'meadow' ],
8 = [ 'landuse', 'meadow' ],
10 = [ 'landuse', 'forest', 'wood', 'yes' ],
};
my $zpusoby_vyuziti_map = {
1 = [ 'landuse', 'greenhouse_horticulture' ],
2 = [ 'landuse', 'plant_nursery' ],
3 = [ 'landuse', 'plantation', 'wood', 'yes' ],
4 = [ 'natural', 'wood', 'wood', 'yes' ],
5 = [ 'landuse', 'forest', 'wood', 'yes' ],
6 = [ 'natural', 'water', 'water', 'pond' ],
7 = [ 'waterway', 'river' ],
8 = [ 'waterway', 'canal' ],
9 = [ 'natural', 'water', 'water', 'lake' ],
10 = [ 'natural', 'water', 'water', 'reservoir' ],
11 = [ 'natural', 'wetland' ],
#12 = [ ] ???
13 = [ 'landuse', 'brownfield' ],
14 = [ 'landuse', 'railway' ],
15 = [ 'landuse', 'highway' ],
16 = [ 'landuse', 'highway' ],
17 = [ 'landuse', 'highway' ],
18 = [ 'landuse', 'transport' ],
19 = [ 'landuse', 'grass' ],
20 = [ 'landuse', 'recreation_ground' ],
21 = [ 'landuse', 'cemetery' ],
#22 = [] ???
23 = [ 'highway', 'service', 'area', 'yes' ],
24 = [ 'landuse', 'quarry' ],
25 = [ 'landuse', 'landfill' ],
#26 = [] ???
27 = [ 'natural', 'scrub' ], # ... casty vyskyt
28 = [ 'natural', 'water' ],
29 = [ 'landuse', 'industrial' ],
};
sub warning
{
my ($text) = @_;
print STDERR $text\n;
return undef;
}
sub nodemap_get
{
my ($lon, $lat) = @_;
my $key = $lon/$lat;
my $node = $nodemap-{$key};
return $node if (defined $node);
my $krpos = [$lon, $lat];
my $osmpos = $proj_from-transform ($proj_to, $krpos);
$node = {
id = --$idgen,
krpos = $krpos,
osmpos = $osmpos,
};
$nodemap-{$key} = $node;
return $node;
}
sub parse_gml_pos
{
my ($pos) = @_;
return undef unless ($pos =~ /^\s*(\S+)\s+(\S+)\s*$/);
my $lon = $1;
my $lat = $2;
return nodemap_get ($lon, $lat);
}
sub parse_gml_poslist
{
my ($poslist) = @_;
$poslist =~ s/^\s*//;
$poslist =~ s/\s*$//;
my @list = split (/\s+/, $poslist);
die Odd number of poslist entries: $poslist
unless ((int(@list) % 2) == 0);
my $waynodes = [];
for (my $i = 0; $i int(@list); $i += 2)
{
my $lon = @list[$i];
my $lat = @list[$i + 1];
push (@$waynodes, nodemap_get ($lon, $lat));
}
my $way = {
id = --$idgen,
nodes = $waynodes
};
return $way;
}
sub parse_pai_defpoint
{
my ($xdefpoint) = @_;
my $gml_point = $xdefpoint-{'gml:Point'};
return undef unless (defined $gml_point);
my $gml_pos = $gml_point-{'gml:pos'};
return undef unless (defined $gml_pos);
my $srs_name = $gml_point-{'srsName'};
return undef unless (defined $srs_name);
die Unsupported EPSG: $srs_name unless ($srs_name eq $epsg_srs_name);
return parse_gml_pos ($gml_pos);
}
sub
Re: [Talk-cz] Skript na parcely RUIANu
Ahoj, vypadá to hezky, ale fakt chceme mít v OSM přesný obraz KM? Sem tam si zemědělci něco přiorají, občas zase kus nechají a ten brzy zaroste nějakým křovím nebo lesem. Někde jsem zase viděl velkou zahradu, o kterou se už majitel nechtěl starat, tak jí kus pronajal sousedovi a ten si o to rozšířil své pole. Jak už tu padlo, v KM je zanesen právní stav, ale v OSM potřebujeme skutečnost. Z tohohle pohledu mi přijde pLPIS lepší. Plus to, že je dostupný pro celou republiku. Marián On 25. července 2014 21:22:23 CEST, Martin Švec - OSM o...@maatts.cz wrote: Ahoj, obrázky parcel od Petra Vejsady mi nějak nedaly spát, tak jsem včera v noci ubastlil skript na konverzi parcel do OSM. Skriptík nemá jiné ambice než vyzkoušet si natažení parcel do JOSM a vůbec se seznámit s VFR formátem RUIANu, takže nečekejte zázraky. Mj. nefungují multipolygony, beru jen polygony typu gml:LinearRing apod. Tagování jsem narychlo opsal z http://wiki.openstreetmap.org/wiki/Cs:RUIAN, knihovny na zpracování GML apod. jsem zatím nehledal. Vyžaduje perl, proj, Geo::Proj4 (http://search.cpan.org/dist/Geo-Proj4/), a grid jezek_czech08.llb přejmenovaný na czech (vyšťouráno tady v archívu, http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid). Pár obrázků: http://www.maatts.cz/obrazek/3/lomnice-ruian-parcely-png/ http://www.maatts.cz/obrazek/3/blansko-ruian-parcely-png/ http://www.maatts.cz/obrazek/3/brno-ruian-parcely-png/ Vyzkoušeno na uvedených třech obcích = na čtvrté klidně může zhavarovat. ;-) BACHA - ve velkých městech skript žere mraky paměti, optimalizacema jsem se netrápil. :-)) Třeba Brno (kód obce 582786) má ve výsledku 1.122.529 uzlů a 251.328 parcel, a skriptík na to potřebuje aspoň 12GB RAM a/nebo hodně swapu na SSD. Martin ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz -- Odesláno z mého telefonu s Androidem pomocí pošty K-9 Mail. Omluvte prosím moji stručnost.___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Skript na parcely RUIANu
No, poté co jsem si to vyzkoušel jsem si prakticky jistý, že 1:1 hromadný převod všech parcel do OSM není průchozí :-) Ty miliony rozdrobených parcelek by nadělaly solidní binec v databázi i na mapě. A jak říkáš, zrovna LPIS u zemědělské půdy líp kopíruje realitu. Spíš je to o tom hledat způsoby, jak ty RUIAN data rozumně využít, když už jsou k dispozici. Martin Dne 25.7.2014 21:54, Marián Kyral napsal(a): Ahoj, vypadá to hezky, ale fakt chceme mít v OSM přesný obraz KM? Sem tam si zemědělci něco přiorají, občas zase kus nechají a ten brzy zaroste nějakým křovím nebo lesem. Někde jsem zase viděl velkou zahradu, o kterou se už majitel nechtěl starat, tak jí kus pronajal sousedovi a ten si o to rozšířil své pole. Jak už tu padlo, v KM je zanesen právní stav, ale v OSM potřebujeme skutečnost. Z tohohle pohledu mi přijde pLPIS lepší. Plus to, že je dostupný pro celou republiku. Marián On 25. července 2014 21:22:23 CEST, Martin Švec - OSM o...@maatts.cz wrote: Ahoj, obrázky parcel od Petra Vejsady mi nějak nedaly spát, tak jsem včera v noci ubastlil skript na konverzi parcel do OSM. Skriptík nemá jiné ambice než vyzkoušet si natažení parcel do JOSM a vůbec se seznámit s VFR formátem RUIANu, takže nečekejte zázraky. Mj. nefungují multipolygony, beru jen polygony typu gml:LinearRing apod. Tagování jsem narychlo opsal z http://wiki.openstreetmap.org/wiki/Cs:RUIAN, knihovny na zpracování GML apod. jsem zatím nehledal. Vyžaduje perl, proj, Geo::Proj4 (http://search.cpan.org/dist/Geo-Proj4/), a grid jezek_czech08.llb přejmenovaný na czech (vyšťouráno tady v archívu, http://freegis.fsv.cvut.cz/gwiki/S-JTSK_/_Grid). Pár obrázků: http://www.maatts.cz/obrazek/3/lomnice-ruian-parcely-png/ http://www.maatts.cz/obrazek/3/blansko-ruian-parcely-png/ http://www.maatts.cz/obrazek/3/brno-ruian-parcely-png/ Vyzkoušeno na uvedených třech obcích = na čtvrté klidně může zhavarovat. ;-) BACHA - ve velkých městech skript žere mraky paměti, optimalizacema jsem se netrápil. :-)) Třeba Brno (kód obce 582786) má ve výsledku 1.122.529 uzlů a 251.328 parcel, a skriptík na to potřebuje aspoň 12GB RAM a/nebo hodně swapu na SSD. Martin -- Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz -- Odesláno z mého telefonu s Androidem pomocí pošty K-9 Mail. Omluvte prosím moji stručnost. ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Skript na parcely RUIANu
Ahoj, ze 40 milionů drobných parcelek se dá jedním příkazem udělat 60 multipolygonů select druh_pozemku,zpusob_vyuziti_pozemku,st_union(hranice) from ruian.rn_parcela group by druh_pozemku,zpusob_vyuziti_pozemku a z nich dalším jedním či dvěma příkazy se dá udělat N polygonů, které budou sdružovat sousedící parcely se stejným druh_pozemku a zpusob_vyuziti_pozemku st_dump(hranice) ... atd. a toto zkombinovat s LPIS. LPIS by měl mát přednost, ano, a tam, kde LPIS není, protože se majitel půdy nezaregistroval či jde o landuse mimo zemědělství, tam použit RUIAN. ? Dne Pá 25. července 2014 22:04:50, Martin Švec - OSM napsal(a): No, poté co jsem si to vyzkoušel jsem si prakticky jistý, že 1:1 hromadný převod všech parcel do OSM není průchozí :-) Ty miliony rozdrobených parcelek by nadělaly solidní binec v databázi i na mapě. A jak říkáš, zrovna LPIS u zemědělské půdy líp kopíruje realitu. Spíš je to o tom hledat způsoby, jak ty RUIAN data rozumně využít, když už jsou k dispozici. -- Petr ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Skript na parcely RUIANu
Ahoj, ano, tak bych si představoval to rozumné využití RUIANu. Akorát mě nenapadlo, že je to až tak snadné. Asi je čas naučit se PostGIS ;-) Martin On 25.7.2014 22:25, Petr Vejsada wrote: Ahoj, ze 40 milionů drobných parcelek se dá jedním příkazem udělat 60 multipolygonů select druh_pozemku,zpusob_vyuziti_pozemku,st_union(hranice) from ruian.rn_parcela group by druh_pozemku,zpusob_vyuziti_pozemku a z nich dalším jedním či dvěma příkazy se dá udělat N polygonů, které budou sdružovat sousedící parcely se stejným druh_pozemku a zpusob_vyuziti_pozemku st_dump(hranice) ... atd. a toto zkombinovat s LPIS. LPIS by měl mát přednost, ano, a tam, kde LPIS není, protože se majitel půdy nezaregistroval či jde o landuse mimo zemědělství, tam použit RUIAN. ? Dne Pá 25. července 2014 22:04:50, Martin Švec - OSM napsal(a): No, poté co jsem si to vyzkoušel jsem si prakticky jistý, že 1:1 hromadný převod všech parcel do OSM není průchozí :-) Ty miliony rozdrobených parcelek by nadělaly solidní binec v databázi i na mapě. A jak říkáš, zrovna LPIS u zemědělské půdy líp kopíruje realitu. Spíš je to o tom hledat způsoby, jak ty RUIAN data rozumně využít, když už jsou k dispozici. -- Petr ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
Re: [Talk-cz] Import of farmland from LPIS
Ahoj! ani nic nequotuji, jen co se mi zdá, tak podle toho zákona, co sem dával odkaz hanoj, jsou v LPIS jen ty pozemky, u kterých to jejich uživatelé chtějí. Kdo tam nechce být, tak tam není. Chápu správně? Jistej si nejsem, ale taky mam ten pocit. Na druhou stranu -- je to kvuli dotacim, tak tam vetsina pozemku bude. Udělal jsem narychlo vizualizaci RUIAN, pouze druh pozemku. V tabulkách je ještě způsob využití pozemku, což je údaj upřesňující. Tím jsem se teď nezabýval. Obrázky jsou na http://pedro.poloha.net/osm/landuse vždy OSM odklad a pak RUIAN obarvené parcely podle druhu pozemku. Zdá se, že to nevypadá zle, na to, že je to jen to hrubé členění. Pekne, moc pekne. Kdyz to budu chtit zopakovat vlastnim qgisem -- jakou adresu musim zadat? Je tam i screenshot z qgisu, kde je vysvětlení, co která barva znamená. S qgisem moc neumím a tak nevím, jak vyexportovat legendu; proto ten screenshot. Jsou tam vidět zahrady i trávníky u domů, lesy i orná půda, koryto potoka atd. S tím upřesněním by to bylo samozřejmě ještě lepší, ale to je 30 kategorií a nevím, jaké barvičky použít :-) Kdyby se z toho podarilo dostat udaje o zahradkach a plotech, hned by mapa vypadala zas o kousek lip :-) Pavel -- (english) http://www.livejournal.com/~pavelmachek (cesky, pictures) http://atrey.karlin.mff.cuni.cz/~pavel/picture/horses/blog.html ___ Talk-cz mailing list Talk-cz@openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz
