[Talk-cz] Skript pro výpočet souřadnice středu cesty?
Vlákno 14.11. - 16.11.2022, počet zpráv: 22
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
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>
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.
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
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>
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>
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>
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>
On 2022-11-14, 18:21 GMT, Miroslav Suchý wrote:
zobrazit citaci
> 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.
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>
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
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
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
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>
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>
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>
ú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>
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>
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>
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>
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)
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