« zpět na výpis měsíce |

[Talk-cz] Skript pro výpočet souřadnice středu cesty?

Vlákno 14.11. - 16.11.2022, počet zpráv: 22


14.11.2022 06:45:54 (#1)
gravatar

Matěj Cepl

<mcepl at cepl.eu>
161
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

14.11.2022 06:52:10 (#2)
gravatar

Jan Martinec

<jan at martinec.name>
471 3565
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ě. Zdar, HPM Dne po 14. 11. 2022 18:48 uživatel Matěj Cepl <mcepl na cepl.eu> napsal: zobrazit citaci
> 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 >
------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221114/b1d2001a/attachment.htm>

14.11.2022 07:02:51 (#3)
gravatar

Matěj Cepl

<mcepl at cepl.eu>
161
On 2022-11-14, 17:52 GMT, Jan Martinec wrote: zobrazit citaci
> 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ě.
Automatizace byla pochopitelně celý point a chtěl jsem komunikovat v platform-independent stylu (takže Geo-URL raději nežli normální https do osm.org). Představoval jsem si něco jako: $ osm_stred 'https://www.openstreetmap.org/way/30722274#map=13/49.2925/38.1816' geo:49.2943095,38.1804967 $ (asi spíše tak, nežli hledání podle názvu, protože potřebuji prohledávat možné varianty). Matěj -- https://matej.ceplovi.cz/blog/, @mcepl na floss.social GPG Finger: 3C76 A027 CA45 AD70 98B5 BC1D 7920 5802 880B C9D8 They that can give up essential liberty to obtain a little temporary safety deserve neither liberty nor safety. -- Benjamin Franklin, Historical Review of Pennsylvania, 1759.

14.11.2022 07:21:39 (#4)
gravatar

Miroslav Suchý

<mirek at rodina-sucha.cz>
112 5668
Dne 14. 11. 22 v 18:45 Matěj Cepl napsal(a): zobrazit citaci
> 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.
LMGTFY (tj. nepouzil jsem to) https://stackoverflow.com/questions/23583649/how-to-compute-the-centroid-of-a-polygon-in-python#23591228 Mirek

14.11.2022 10:06:42 (#5)
gravatar

Jan Macura

<macurajan at gmail.com>
661 2731
Ahoj On Mon, 14 Nov 2022 at 18:48, Matěj Cepl <mcepl na cepl.eu> wrote: zobrazit citaci
> (...) > > 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. ------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221114/26c8009c/attachment.htm>

14.11.2022 10:20:40 (#6)
gravatar

Jiří Sedláček

<jirisedlacek at gmail.com>
43 539
Ahoj, se shapely jsem pracoval, ale bacha ? pokud si to pamatuju dobře, tak pro větší (více bodů) polygony ty výpočty nebyly úplně fofrem. Spíš to bylo místama dost pomalý. Ale zas se s tím pracuje snadno, jen tam musíš nějak dostat ta data ... (geojson to myslím umí). On Mon, Nov 14, 2022 at 10:08 PM Jan Macura <macurajan na gmail.com> wrote: zobrazit citaci
> 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 >
-- S pozdravem, Jirka Sedláček --- jirisedlacek na gmail.com ------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221114/6c14671b/attachment-0001.htm>

15.11.2022 09:17:15 (#7)
gravatar

ramirezhav na seznam.cz

<ramirezhav at seznam.cz>
103
Tohle téma mě taky zajímá, ale jen z hlediska kachlí. Jelikož vím, že země je kulatá, dá se postupovat v případě kachle tak, že se jen spočítá rozdíl souřadnic (x2-x1)/256 * poloha, která mě zajímá, nebo se musí počítat se zaoblením země jako s nějakou konstantou, která se má započítat?
---------- Původní e-mail ---------- Od: Jan Macura <macurajan na gmail.com> Komu: OpenStreetMap Czech Republic <talk-cz na openstreetmap.org> Datum: 14. 11. 2022 22:13:50 Předmět: Re: [talk-cz] Skript pro výpočet souřadnice středu cesty? " Ahoj On Mon, 14 Nov 2022 at 18:48, Matěj Cepl <mcepl na cepl.eu (mailto: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 (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 " ------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221115/25aebed1/attachment.htm>

15.11.2022 10:08:20 (#8)
gravatar

Lukáš Karas

<lukas.karas at centrum.cz>
52
Ahoj, pokud by jsi měl zájem o čistě offline řešení, mohu odkázat na OSMScout knihovnu: https://github.com/Framstag/libosmscout/ 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. Hledání místa podle názvu pak použít LocationLookup demo: ./Demos/LocationLookup ~/Maps/europe-ukraine-20221022-121415 "Novomykilske" = Region (????????? ??????? ???????) * ????????? ??????? ???????/?????????????? ?????/????????? ???????/??????? ... - Area 304695217 boundary_administrative Poté pomocí DumpData toolu zjistit střed polygonu 304695217... ./DumpData/DumpData ~/Maps/europe-ukraine-20221022-121415 -ao 304695217 Area { OSM id: 13631713 fileOffset: 304695217 type: boundary_administrative boundingBox: [49,18544 N 39,76960 E - 49,61710 N 40,22758 E] center: 49,40127 N 39,99859 E cell level: 8 Name: ????????? ??????? ??????? AdminLevel: 7 ... S pozdravem, Lukáš Dne pondělí 14. listopadu 2022 18:45:54 CET Matěj Cepl napsal(a): zobrazit citaci
> 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
------------- další část --------------- A non-text attachment was scrubbed... Name: signature.asc Type: application/pgp-signature Size: 488 bytes Desc: This is a digitally signed message part. URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221115/13cd613d/attachment.sig>

15.11.2022 07:12:40 (#9)
gravatar

Matěj Cepl

<mcepl at cepl.eu>
161
On 2022-11-14, 18:21 GMT, Miroslav Suchý wrote: zobrazit citaci 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.

15.11.2022 08:11:06 (#10)
gravatar

Jan Martinec

<jan at martinec.name>
471 3565
Ahoj, 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í ;) Zdar, HPM Dne út 15. 11. 2022 20:08 uživatel Matěj Cepl <mcepl na cepl.eu> napsal: zobrazit citaci
> 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 >
------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221115/6bf6b2d1/attachment.htm>

15.11.2022 08:36:01 (#11)
gravatar

Matěj Cepl

<mcepl at cepl.eu>
161
On 2022-11-14, 21:06 GMT, Jan Macura wrote: zobrazit citaci
> Ale šel bych tou cestou Shapely ;-)
To vypadá užitečně. zobrazit citaci
> 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? Matěj -- https://matej.ceplovi.cz/blog/, @mcepl na floss.social GPG Finger: 3C76 A027 CA45 AD70 98B5 BC1D 7920 5802 880B C9D8 Be kind, for everyone you meet is fighting a hard battle. -- Ian MacLaren

15.11.2022 08:38:12 (#12)
gravatar

Matěj Cepl

<mcepl at cepl.eu>
161
On 2022-11-15, 09:08 GMT, Lukáš Karas wrote: zobrazit citaci
> 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.
Tak to by šlo nějak obejít pomocí Cythonu, ale nemáme tu knihovnu balenou pro openSUSE, takže si ještě počkám jestli nevymyslím něco jednoduššího. Děkuji za radu. Matěj -- https://matej.ceplovi.cz/blog/, @mcepl na floss.social GPG Finger: 3C76 A027 CA45 AD70 98B5 BC1D 7920 5802 880B C9D8 Science is meaningless because it gives no answer to our question, the only question important to us: ?What shall we do and how shall we live?? -- Lev Nikolaevich Tolstoy

15.11.2022 08:43:50 (#13)
gravatar

Matěj Cepl

<mcepl at cepl.eu>
161
On 2022-11-15, 19:11 GMT, Jan Martinec wrote: zobrazit citaci
> 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. Matěj -- https://matej.ceplovi.cz/blog/, @mcepl na floss.social GPG Finger: 3C76 A027 CA45 AD70 98B5 BC1D 7920 5802 880B C9D8 Reality is merely an illusion, albeit a very persistent one. -- Albert Einstein

15.11.2022 08:47:31 (#14)
gravatar

Jáchym

<jachym at cepicky.ch>
2
Pardon, že se vkládám do debaty, asi trochu mimo libovolná  knihovna pro práci s geometriemi dává všem objektům atribut `centroid` bod nemusí nutně ležet na linii/v polygonu ale třeba to pomůže? ---- dne Tue, 15 Nov 2022 20:43:50 +0100, Matěj Cepl <mcepl na cepl.eu> napsal --- On 2022-11-15, 19:11 GMT, Jan Martinec wrote: zobrazit citaci
> 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. Matěj -- https://matej.ceplovi.cz/blog/, @mailto:mcepl na floss.social GPG Finger: 3C76 A027 CA45 AD70 98B5 BC1D 7920 5802 880B C9D8 Reality is merely an illusion, albeit a very persistent one. -- Albert Einstein _______________________________________________ talk-cz mailing list mailto:talk-cz na openstreetmap.org https://lists.openstreetmap.org/listinfo/talk-cz https://openstreetmap.cz/talkcz ------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221115/8b3fc07b/attachment.htm>

15.11.2022 08:50:46 (#15)
gravatar

Jan Macura

<macurajan at gmail.com>
661 2731
On Tue, 15 Nov 2022 at 20:44, Matěj Cepl <mcepl na cepl.eu> wrote: zobrazit citaci
> > 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? >
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..? H. ------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221115/c6e465ac/attachment.htm>

15.11.2022 09:37:42 (#16)
gravatar

Jan Macura

<macurajan at gmail.com>
661 2731
On Tue, 15 Nov 2022 at 20:45, Matěj Cepl <mcepl na cepl.eu> wrote: zobrazit citaci
> 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... Ergo možná to může být někdy dostačující, ale pro zmiňovaný use-case bych to jako vhodné řešení neviděl. Obzvlášť, když zmiňovaná knihovna Shapely ti umí vyplivnout očekávaný výsledek bez velké námahy. H. ------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221115/c1400b62/attachment-0001.htm>

16.11.2022 07:34:57 (#17)
gravatar

Jan Dudík

<jan.dudik at gmail.com>
326 645
út 15. 11. 2022 v 21:39 odesílatel Jan Macura <macurajan na gmail.com> napsal: zobrazit citaci
> > >> >> 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 ------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221116/447d10d0/attachment.htm>

16.11.2022 08:18:43 (#18)
gravatar

Mirek Dlask

<dlask.m at gmail.com>
152
Centroid v postgisu má ještě sourozence https://postgis.net/docs/ST_Centroid.html Střed plochy mi nezaručí, že zjištěný bod bude blízko obydlené části města, vesnice. Část bodů (odkazů) skončí v lese, nebo na poli. V OSM jsou nejen hranice, ale každá obec má i bod v zastavěné části obce, nebo přinejhorším hodně blízko. place=* Nebylo by to lepší řešení? Bez ohledu na výše uvedené bude muset řešit konflikty názvů. Naše Babice jsou v okresech: Prachatice, Praha-východ, Hradec Králové, Olomouc, Třebíč, Uherské Hradiště? Takže něco na ekvivalent našeho okresu, kraje. https://postgis.net/docs/ST_Contains.html Mir st 16. 11. 2022 v 7:37 odesílatel Jan Dudík <jan.dudik na gmail.com> napsal: zobrazit citaci
> > > ú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 >
------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221116/a7b2a35d/attachment.htm>

16.11.2022 08:20:07 (#19)
gravatar

Jan Macura

<macurajan at gmail.com>
661 2731
On Wed, 16 Nov 2022 at 07:37, Jan Dudík <jan.dudik na gmail.com> wrote: zobrazit citaci
> > ú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. ------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221116/8938b050/attachment.htm>

16.11.2022 08:36:33 (#20)
gravatar

Jan Martinec

<jan at martinec.name>
471 3565
Typicky Středočeský kraj s těžištěm v Praze;) Ale napadá mě, že obec může mít ve své OSM relaci admin_centre - teda pokud je to zmapovaný, no. Zdar, HPM Dne st 16. 11. 2022 8:24 uživatel Jan Macura <macurajan na gmail.com> napsal: zobrazit citaci
> > 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 >
------------- další část --------------- HTML příloha byla odstraněna... URL: <http://lists.openstreetmap.org/pipermail/talk-cz/attachments/20221116/b1c4e94b/attachment-0001.htm>

16.11.2022 02:07:03 (#21)
gravatar

Matěj Cepl

<mcepl at cepl.eu>
161
On 2022-11-15, 19:50 GMT, Jan Macura wrote: zobrazit citaci
> 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..?
Přesně tak, https://git.sr.ht/~mcepl/osm_where D?oh! Matěj -- https://matej.ceplovi.cz/blog/, @mcepl na floss.social GPG Finger: 3C76 A027 CA45 AD70 98B5 BC1D 7920 5802 880B C9D8 Extremism in the defense of liberty is no vice; moderation in the pursuit of justice is no virtue. -- Barry Goldwater (actually written by Karl Hess)

16.11.2022 02:11:24 (#22)
gravatar

Matěj Cepl

<mcepl at cepl.eu>
161
On 2022-11-15, 20:37 GMT, Jan Macura wrote: zobrazit citaci
> 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...
Tak jsem se obešel bez shapely, OSMPythonTools to dává taky, takže stitny~$ osm_where -l cs -a CZ Kdyně geo:49.3970014,13.0507001 stitny~$ což je docela slušné a je to uvnitř Kdyně. Matěj -- https://matej.ceplovi.cz/blog/, @mcepl na floss.social GPG Finger: 3C76 A027 CA45 AD70 98B5 BC1D 7920 5802 880B C9D8 Pain is inevitable, but misery is optional. We cannot avoid pain, but we can avoid joy. -- Tim Hansel

« zpět na výpis měsíce