Skript pro výpočet souřadnice středu cesty?

Dobrý den, existuje nějaká (nejlépe Pythonní) knihovna pro programatický přístup k OSM? Chtěl bych být schopen rychle vytvářet Geo-URI vesnice, abych ji mohl zmínit třeba na Twitteru. Když někdo zmíni, že Ukrajinská armáda dobyla Novomykilske, tak si rychle najdu, že se jedná asi o https://www.openstreetmap.org/way/30722274, ale chtěl bych něco, co by mi z tohoto URL vygenerovalo něco jako 49.2943095,38.1804967 (přibližně, ne nutně tento bod). Jedná se tedy o dva problémy: 1. Jak vygenerovat seznam všech souřadnic dané cesty. 2. Jak vypočítat něco jako ?střed? obecného geometrického útvaru. Doufám, že na 2. někde existuje nějaký více méně standardní algoritmus (existuje?) a existuje nějaká knihovna na 1.? Děkuji za jakékoli nakopnutí správným směrem. Matěj -- https://matej.ceplovi.cz/blog/, @mcepl na floss.social GPG Finger: 3C76 A027 CA45 AD70 98B5 BC1D 7920 5802 880B C9D8 You either die a hero or you live long enough to see yourself become the villain. -- Harvey Dent in The Dark Knight _______________________________________________ talk-cz mailing list talk-cz na openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz https://openstreetmap.cz/talkcz

Rozkliknu, posunu/zazoomuju mapu dle potřeby, dostanu něco jako https://www.openstreetmap.org/way/30722274#map=9/49.2705/38.2104 , hotovo. Netřeba skriptů, pokud to nepotřebuješ dělat automatizovaně.

Rozkliknu, posunu/zazoomuju mapu dle potřeby, dostanu něco jako https://www.openstreetmap.org/way/30722274#map=9/49.2705/38.2104 , hotovo. Netřeba skriptů, pokud to nepotřebuješ dělat automatizovaně.

Doufám, že na 2. někde existuje nějaký více méně standardní algoritmus (existuje?) a existuje nějaká knihovna na 1.? Děkuji za jakékoli nakopnutí správným směrem.

Doufám, že na 2. někde existuje nějaký více méně standardní algoritmus (existuje?) a existuje nějaká knihovna na 1.? Děkuji za jakékoli nakopnutí správným směrem.

(...) Jedná se tedy o dva problémy: 1. Jak vygenerovat seznam všech souřadnic dané cesty. 2. Jak vypočítat něco jako ?střed? obecného geometrického útvaru. Doufám, že na 2. někde existuje nějaký více méně standardní algoritmus (existuje?) a existuje nějaká knihovna na 1.?

(...) Jedná se tedy o dva problémy: 1. Jak vygenerovat seznam všech souřadnic dané cesty. 2. Jak vypočítat něco jako ?střed? obecného geometrického útvaru. Doufám, že na 2. někde existuje nějaký více méně standardní algoritmus (existuje?) a existuje nějaká knihovna na 1.?

Ahoj On Mon, 14 Nov 2022 at 18:48, Matěj Cepl <mcepl na cepl.eu> wrote:

(...) Jedná se tedy o dva problémy: 1. Jak vygenerovat seznam všech souřadnic dané cesty. 2. Jak vypočítat něco jako ?střed? obecného geometrického útvaru. Doufám, že na 2. někde existuje nějaký více méně standardní algoritmus (existuje?) a existuje nějaká knihovna na 1.?

Ad 2) jakože.. vždycky si to můžeš vypočítat ručně ( https://math.stackexchange.com/questions/90463/how-can-i-calculate-the-centroid-of-polygon). Ale šel bych tou cestou Shapely ;-) Ad 1) nerozumím moc problému. Možná přesněji formuluj, co by mělo být vstupem toho Tvého procesu a co očekáváš za výstup. Cesta je definována svými souřadnicemi. Resp. v OSM modelu je definována seznamem svých uzlů, ale není problém exportovat si data to nějakého standardního formátu, kde samozřejmě budou souřadnice, třeba GeoJSON, GPX, GML, KML, ... H. _______________________________________________ talk-cz mailing list talk-cz na openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz https://openstreetmap.cz/talkcz

Dobrý den, existuje nějaká (nejlépe Pythonní) knihovna pro programatický přístup k OSM? Chtěl bych být schopen rychle vytvářet Geo-URI vesnice, abych ji mohl zmínit třeba na Twitteru. Když někdo zmíni, že Ukrajinská armáda dobyla Novomykilske, tak si rychle najdu, že se jedná asi o https://www.openstreetmap.org/way/30722274, ale chtěl bych něco, co by mi z tohoto URL vygenerovalo něco jako 49.2943095,38.1804967 (přibližně, ne nutně tento bod). Jedná se tedy o dva problémy: 1. Jak vygenerovat seznam všech souřadnic dané cesty. 2. Jak vypočítat něco jako ?střed? obecného geometrického útvaru. Doufám, že na 2. někde existuje nějaký více méně standardní algoritmus (existuje?) a existuje nějaká knihovna na 1.? Děkuji za jakékoli nakopnutí správným směrem. Matěj -- https://matej.ceplovi.cz/blog/, @mcepl na floss.social GPG Finger: 3C76 A027 CA45 AD70 98B5 BC1D 7920 5802 880B C9D8 You either die a hero or you live long enough to see yourself become the villain. -- Harvey Dent in The Dark Knight _______________________________________________ talk-cz mailing list talk-cz na openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz https://openstreetmap.cz/talkcz

LMGTFY (tj. nepouzil jsem to) https://stackoverflow.com/questions/23583649/how-to-compute-the-centroid-of-a-polygon-in-python#23591228

LMGTFY (tj. nepouzil jsem to) https://stackoverflow.com/questions/23583649/how-to-compute-the-centroid-of-a-polygon-in-python#23591228

On 2022-11-14, 18:21 GMT, Miroslav Suchý wrote:

LMGTFY (tj. nepouzil jsem to)

https://stackoverflow.com/questions/23583649/how-to-compute-the-centroid-of-a-polygon-in-python#23591228 Zajímavé. Matěj -- https://matej.ceplovi.cz/blog/, @mcepl na floss.social GPG Finger: 3C76 A027 CA45 AD70 98B5 BC1D 7920 5802 880B C9D8 At the day of judgement it shall not be asked of us what we have read, but what we have done: nor how well we have said, but how religiously we have lived. -- Thomas á Kempis The Imitation of Christ, III. _______________________________________________ talk-cz mailing list talk-cz na openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz https://openstreetmap.cz/talkcz

Ale šel bych tou cestou Shapely ;-)

Ad 1) nerozumím moc problému. Možná přesněji formuluj, co by mělo být vstupem toho Tvého procesu a co očekáváš za výstup.

Ale šel bych tou cestou Shapely ;-)

Ad 1) nerozumím moc problému. Možná přesněji formuluj, co by mělo být vstupem toho Tvého procesu a co očekáváš za výstup.

Je ale psaná v C++ a Python binding nemáme, takže nevím jak moc je pro tebe vnodná. Data pro danou zemi je buď potřeba importovat, nebo stáhnout ze serveru pomocí OSMScout2 UI.

Je ale psaná v C++ a Python binding nemáme, takže nevím jak moc je pro tebe vnodná. Data pro danou zemi je buď potřeba importovat, nebo stáhnout ze serveru pomocí OSMScout2 UI.

Pokud to má být dostatečně mrňavý polygon (jednotky km a méně), lze předstírat placatou Zemi, a prostě udělat průměr z. šířky a průměr délky. Je to ošklivé, je to nepřesné, je to kacířské, a je to triviální a mnohdy dostačující ;)

Pokud to má být dostatečně mrňavý polygon (jednotky km a méně), lze předstírat placatou Zemi, a prostě udělat průměr z. šířky a průměr délky. Je to ošklivé, je to nepřesné, je to kacířské, a je to triviální a mnohdy dostačující ;)

Pokud to má být dostatečně mrňavý polygon (jednotky km a méně), lze předstírat placatou Zemi, a prostě udělat průměr z. šířky a průměr délky. Je to ošklivé, je to nepřesné, je to kacířské, a je to triviální a mnohdy dostačující ;)

Pokud to má být dostatečně mrňavý polygon (jednotky km a méně), lze předstírat placatou Zemi, a prostě udělat průměr z. šířky a průměr délky. Je to ošklivé, je to nepřesné, je to kacířské, a je to triviální a mnohdy dostačující ;)

Ad 1) nerozumím moc problému. Možná přesněji formuluj, co by mělo být vstupem toho Tvého procesu a co očekáváš za výstup.

Potřebuji přístup do OSM, abych z toho URL dostal seznam souřadnic jednotlivých bodů polygonu. Dává smysl?

Ad 1) nerozumím moc problému. Možná přesněji formuluj, co by mělo být vstupem toho Tvého procesu a co očekáváš za výstup.

Potřebuji přístup do OSM, abych z toho URL dostal seznam souřadnic jednotlivých bodů polygonu. Dává smysl?

Pokud to má být dostatečně mrňavý polygon (jednotky km a méně), lze předstírat placatou Zemi, a prostě udělat průměr z. šířky a průměr délky. Je to ošklivé, je to nepřesné, je to kacířské, a je to triviální a mnohdy dostačující ;)

Good point.

On 2022-11-15, 19:11 GMT, Jan Martinec wrote:

Pokud to má být dostatečně mrňavý polygon (jednotky km a méně), lze předstírat placatou Zemi, a prostě udělat průměr z. šířky a průměr délky. Je to ošklivé, je to nepřesné, je to kacířské, a je to triviální a mnohdy dostačující ;)

Good point.

Ani moc ne. Takový postup ti nedá střed polygonu, jak jej obvykle chápeme

(tj. jeho těžiště), ale získáš střed jeho minimálního ohraničujícího obdélníku (bounding box / envelope). Tzn. dost snadno dostaneš třeba bod který leží zcela mimo ten vstupní polygon. V Tvém případě by to mohlo mít ten důsledek, že třeba pro polygon města Kdyně <https://www.openstreetmap.org/relation/439157#map=13/49.3785/13.0482> (sry, ukrajinský vesnice tak dobře neznám jako ty český) bys dostal bod někde v Brnířově. Protože ano, i takhle může vypadat území obce... Ale ani těžiště nemusí nutně ležet na území polygonu. Jsou obce s

Ani moc ne. Takový postup ti nedá střed polygonu, jak jej obvykle chápeme

(tj. jeho těžiště), ale získáš střed jeho minimálního ohraničujícího obdélníku (bounding box / envelope). Tzn. dost snadno dostaneš třeba bod který leží zcela mimo ten vstupní polygon. V Tvém případě by to mohlo mít ten důsledek, že třeba pro polygon města Kdyně <https://www.openstreetmap.org/relation/439157#map=13/49.3785/13.0482> (sry, ukrajinský vesnice tak dobře neznám jako ty český) bys dostal bod někde v Brnířově. Protože ano, i takhle může vypadat území obce... Ale ani těžiště nemusí nutně ležet na území polygonu. Jsou obce s

út 15. 11. 2022 v 21:39 odesílatel Jan Macura <macurajan na gmail.com> napsal:

Ani moc ne. Takový postup ti nedá střed polygonu, jak jej obvykle

chápeme (tj. jeho těžiště), ale získáš střed jeho minimálního ohraničujícího obdélníku (bounding box / envelope). Tzn. dost snadno dostaneš třeba bod který leží zcela mimo ten vstupní polygon. V Tvém případě by to mohlo mít ten důsledek, že třeba pro polygon města Kdyně <https://www.openstreetmap.org/relation/439157#map=13/49.3785/13.0482> (sry, ukrajinský vesnice tak dobře neznám jako ty český) bys dostal bod někde v Brnířově. Protože ano, i takhle může vypadat území obce... Ale ani těžiště nemusí nutně ležet na území polygonu. Jsou obce s

exklávami i obce s enklávami. --- JAnD _______________________________________________ talk-cz mailing list talk-cz na openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz https://openstreetmap.cz/talkcz

út 15. 11. 2022 v 21:39 odesílatel Jan Macura <macurajan na gmail.com>

Ani moc ne. Takový postup ti nedá střed polygonu, jak jej obvykle chápeme (tj. jeho těžiště), ale získáš střed jeho minimálního ohraničujícího obdélníku (bounding box / envelope). Tzn. dost snadno dostaneš třeba bod který leží zcela mimo ten vstupní polygon. V Tvém případě by to mohlo mít ten důsledek, že třeba pro polygon města Kdyně <https://www.openstreetmap.org/relation/439157#map=13/49.3785/13.0482> (sry, ukrajinský vesnice tak dobře neznám jako ty český) bys dostal bod někde v Brnířově. Protože ano, i takhle může vypadat území obce... Ale ani těžiště nemusí nutně ležet na území polygonu. Jsou obce s

exklávami i obce s enklávami.

út 15. 11. 2022 v 21:39 odesílatel Jan Macura <macurajan na gmail.com> napsal:

Ani moc ne. Takový postup ti nedá střed polygonu, jak jej obvykle chápeme (tj. jeho těžiště), ale získáš střed jeho minimálního ohraničujícího obdélníku (bounding box / envelope). Tzn. dost snadno dostaneš třeba bod který leží zcela mimo ten vstupní polygon. V Tvém případě by to mohlo mít ten důsledek, že třeba pro polygon města Kdyně <https://www.openstreetmap.org/relation/439157#map=13/49.3785/13.0482> (sry, ukrajinský vesnice tak dobře neznám jako ty český) bys dostal bod někde v Brnířově. Protože ano, i takhle může vypadat území obce... Ale ani těžiště nemusí nutně ležet na území polygonu. Jsou obce s

exklávami i obce s enklávami.

On Wed, 16 Nov 2022 at 07:37, Jan Dudík <jan.dudik na gmail.com> wrote:

út 15. 11. 2022 v 21:39 odesílatel Jan Macura <macurajan na gmail.com> napsal:

Ani moc ne. Takový postup ti nedá střed polygonu, jak jej obvykle chápeme (tj. jeho těžiště), ale získáš střed jeho minimálního ohraničujícího obdélníku (bounding box / envelope). Tzn. dost snadno dostaneš třeba bod který leží zcela mimo ten vstupní polygon. V Tvém případě by to mohlo mít ten důsledek, že třeba pro polygon města Kdyně <https://www.openstreetmap.org/relation/439157#map=13/49.3785/13.0482> (sry, ukrajinský vesnice tak dobře neznám jako ty český) bys dostal bod někde v Brnířově. Protože ano, i takhle může vypadat území obce... Ale ani těžiště nemusí nutně ležet na území polygonu. Jsou obce s

exklávami i obce s enklávami.

Good point. H. _______________________________________________ talk-cz mailing list talk-cz na openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz https://openstreetmap.cz/talkcz

Mě moc ne, ale možná si sedím na vedení :-) Možná hledáš něco jako Overpass (https://overpass-turbo.eu/)? Nebo Query API ( https://www.openstreetmap.org/query?lat=49.77675&lon=18.45250)? Nebo Nominatim? Nebo něco úplně jiného..?

Mě moc ne, ale možná si sedím na vedení :-) Možná hledáš něco jako Overpass (https://overpass-turbo.eu/)? Nebo Query API ( https://www.openstreetmap.org/query?lat=49.77675&lon=18.45250)? Nebo Nominatim? Nebo něco úplně jiného..?

Ani moc ne. Takový postup ti nedá střed polygonu, jak jej obvykle chápeme (tj. jeho těžiště), ale získáš střed jeho minimálního ohraničujícího obdélníku (bounding box / envelope). Tzn. dost snadno dostaneš třeba bod který leží zcela mimo ten vstupní polygon. V Tvém případě by to mohlo mít ten důsledek, že třeba pro polygon města Kdyně <https://www.openstreetmap.org/relation/439157#map=13/49.3785/13.0482> (sry, ukrajinský vesnice tak dobře neznám jako ty český) bys dostal bod někde v Brnířově. Protože ano, i takhle může vypadat území obce...

Ani moc ne. Takový postup ti nedá střed polygonu, jak jej obvykle chápeme (tj. jeho těžiště), ale získáš střed jeho minimálního ohraničujícího obdélníku (bounding box / envelope). Tzn. dost snadno dostaneš třeba bod který leží zcela mimo ten vstupní polygon. V Tvém případě by to mohlo mít ten důsledek, že třeba pro polygon města Kdyně <https://www.openstreetmap.org/relation/439157#map=13/49.3785/13.0482> (sry, ukrajinský vesnice tak dobře neznám jako ty český) bys dostal bod někde v Brnířově. Protože ano, i takhle může vypadat území obce...