Potrzebuję obliczyć pole powierzchni akwenów (obszarów) będących wielokątami (30-50 wierzchołków). Akweny mam w pliku GPX (ale mogę to przerobić na TXT, XLS czy coś innego). Wierzchołki mają współrzędne geograficzne w układzie WGS-84.
Znacie jakiś program (bezpłatny) który liczyłby powierzchnię czegoś takiego?
Może być on line, może być excel...
Albo ma ktoś inne pomysły?
Przykład w załączeniu.

_________________
wolne publikacje nautyczne

***********************************************
* Gdyby wszyscy byli bogaci nikt nie chciałby wiosłować *
***********************************************

Kiedyś pomógłbym szybko, bo stara (pod DOS) wersja programu CAD, na którym pracowałem, łykała pliki TXT ze współrzędnymi X-Y a nawet X-Y-Z i rysowała, co tam byś chciał, np. właśnie wielokąty, i liczyła też pole powierzchni, o ile liczba punktów wsadowych była mniejsza od 200, ale można było taki większy skleić z mniejszych...

Dopiero we czwartek mógłbym zajrzeć do wersji pod Windows, którą mam, ale z tej funkcji wczytywania współrzędnych z pliku TXT nie korzystam, więc nawet nie wiem, czy jest.

_________________
Robert Hoffman
***
https://stareaparatyimojefotografowanie.blogspot.com/

Mam CADa i korzystam z obliczania powierzchni czy objętości (choć wierzchołków z pliku nie wczytywałem), ale... co ze zboczeniem nawigacyjnym?

_________________
wolne publikacje nautyczne

***********************************************
* Gdyby wszyscy byli bogaci nikt nie chciałby wiosłować *
***********************************************

Da się tę figurę podzielić na prostokąty? Pole prostokąta łatwe, x= a x h
Jaka dokładność?


Wzór na odległość między dwoma punktami na mapie:

_________________
noone

Teoretycznie - wszystko się da, ale
Obawiam się że wiedza o odległości poszczególnych punktów - niekoniecznie mi się przyda...

Zresztą - odległości pomiędzy punktami (i obwód) mogę odczytać przekształcając akwen w drogę i odczytując odległości poszczególnych odcinków (w OpenCPNie).

_________________
wolne publikacje nautyczne

***********************************************
* Gdyby wszyscy byli bogaci nikt nie chciałby wiosłować *
***********************************************

Gdybym dostał współrzędne X-Y wierzchołków, to nawet wklepałbym, do mojego CADa z palca, na klawiaturze, taki wielokąt w parę minut!

_________________
Robert Hoffman
***
https://stareaparatyimojefotografowanie.blogspot.com/

Jeżeli będziesz miał dane o długości poszczególnych odcinków to możesz akwen przekształcić w kwadrat lub prostokąt lub trójkąt lub koło, na które wzory powierzchni są dostępne.

_________________
Wojciech1968 Port Rybnik
W życiu piękne są tylko chwile...

Ja z matmy miałem poprawkę, ale pamiętam, że figury o tej samej sumie długości boków nie muszą mieć tej samej powierzchni!!!
Ale może źle pamiętam...
Kwadrat i romb, o boku "a", nie mają tej samej powierzchni...

_________________
Robert Hoffman
***
https://stareaparatyimojefotografowanie.blogspot.com/

OK. wypukłość geoidy przy takim akwenie można pominąć. Ale musiałbym kupę nowych punktów wyliczyć:
- kwadrat, prostokąt - przecięcia odcinków łączących punkty z południkami i/lub równoleżnikami innych punktów.
- trójkąt - potrzebowałbym wysokości (więc wyliczanie przecięcia prostopadłych do jednego odcinka z innym) lub kątów...
- koło? - - zaszalałeś! U mnie jpg widać, a u Was?
To już wołałbym pokratkować drobno i liczyć kratki...
Zwracam uwagę na użycie słowa "akwen" w pierwszym poście w liczbie mnogiej!

_________________
wolne publikacje nautyczne

***********************************************
* Gdyby wszyscy byli bogaci nikt nie chciałby wiosłować *
***********************************************

Zaraz dostaniesz, o ile .jpg jest w skali. W jakich jednostkach? m2, mile2, sążnie2....

km2 (może być w m2, już sobie "przeliczę").
A czym to robisz? Drugi akwen załączam.
A JPG! Ale to z odwzorowania wiernokątnego (Merkator), które wiernopowierzchniowym nie jest... Nie potrafię "z głowy" oszacować błędu...
to JPG drugiego akwenu też wrzucam:

_________________
wolne publikacje nautyczne

***********************************************
* Gdyby wszyscy byli bogaci nikt nie chciałby wiosłować *
***********************************************

Z pliku JPG będzie spory błąd wynikający z rozdzielczości pliku. Tak to ja też mógłbym zrobić - mój CAD też łyka bitmapy. Trzeba znać jeszcze długość przynajmniej jednego boku w celu wyskalowania obrazu JPG w pamięci CAD. Ale chodziło chyba o dokładne współrzędne...

_________________
Robert Hoffman
***
https://stareaparatyimojefotografowanie.blogspot.com/

Wojtek, ten akwen ma wymiary ok. 3x2 Mm.

Zboczeniem jest mówienie o zboczeniu nawigacyjnym

_________________
Pozdrawiam,
Marek Grzywa

Użyj strony Geoportal2. Wejdź na mapę , przesuń się na interesujący cię obszar, kliknij na górnej zakładce "pomiar powierzchni". Przesuwając kursor do twoich kolejnych współrzędnych robisz klik, jak objedziesz wszystkie punkty pojawi ci się powierzchnia.
Jak potrzebujesz wsparcia dzwoń. 501 077 854
ZR

Z błędem skali ten obszar ma 7,67649731 km2

* ~ * ~ * ~ * ~ * ~ * ~ * ~ * ~ * ~ *



Brak podziałki.

Chcesz olać, że akwen leży na szerokości 54°17' czy fakt, że cosinus tej wartości = 0,58 (przyjąć = 1 czy jak?)? Czy może inaczej definiujesz "zboczenie nawigacyjne"?

* ~ * ~ * ~ * ~ * ~ * ~ * ~ * ~ * ~ *

_________________
wolne publikacje nautyczne

***********************************************
* Gdyby wszyscy byli bogaci nikt nie chciałby wiosłować *
***********************************************

Niestety. Nijak nie mogę określić wierzchołków na wodzie. Nie ma współrzędnych geograficznych.

_________________
wolne publikacje nautyczne

***********************************************
* Gdyby wszyscy byli bogaci nikt nie chciałby wiosłować *
***********************************************

Za późno już - nie wiem dlaczego myślalem o deformacji tych potencjalnych prostokątów na które chciał dzielić Adaś.

_________________
Pozdrawiam,
Marek Grzywa

6152660,08 m2

Wyjaśniając o co chodzi:

Decyzja środowiskowa o budowie kanału przez Mierzeję Wiślaną z grudnia 2018 wyznaczała akwen z zakazem żeglugi, który ktoś radośnie namalował tak: , nie przejmując się ani torami wodnymi, ani kotwicowiskiem...
Urząd Morski w Gdyni wydał (13 lipca 2020) zarządzenie porządkowe nr 2020/8, określające dwa akweny zabronione - przez podanie współrzędnych ich wierzchołków. Odpowiednia wiadomość (428/20) ukazała się w WŻ 29-30/20. Akweny wyglądały tak:

W styczniu 2021 - pojęcia nie mam dlaczego - dyrektor UM wydał zarządzenie 2021/2 które zmieniło poprzednie zarządzenie w tej sprawie. Z zarządzenia wyleciały współrzędne wierzchołków akwenów... W ślad za nim w WŻ 03-04/21 ukazała się wiadomość 47/21.

Jak się wydaje de iure akweny są obecnie określone enigmatycznie:
- obszar I zlokalizowany pomiędzy ujściem rzeki Szkarpawa i rzeki Nogat a ujściem cieku Wisła Królewiecka, o powierzchni ok. 9 km²;
- obszar II zlokalizowany pomiędzy ujściem rzeki Szkarpawa i rzeki Nogat a ujściem cieku Cieplicówka, o powierzchni 6 km².

Zastanawiam się kto (BHMW?, UM?), jak (jakim "algorytmem") i na jakiej podstawie określił ukazane na mapie aweny zwane "obszarami"?

"Obszary" na mapie ENC wygkądają obecnie tak:
(indagowany na okoliczność: "skądś ty?" obszar - twierdzi że z zarządzenia DUM)

W międzyczasie doszły żółte pławy - były odpowiednie wiadomości "wrysowujące" je na mapy. Nie wiem w jaki sposób "akwen II" ma obecnie ścięty SW narożnik (bo wcześniej kończył się sporo za ujściem Cieplicówki...).

I stąd dociekanie powierzchni.
Obszar II się ~zgadza - około 6 km², obszar I - nie (9 km² > 7,7km²) - mniejszy o ponad 15%.

_________________
wolne publikacje nautyczne

***********************************************
* Gdyby wszyscy byli bogaci nikt nie chciałby wiosłować *
***********************************************

Da się to zrobić posługując się wzorem.

1. Najpierw musisz przygotować wstępnie dane w postaci zestawu punktów o współrzędnych (xi, yi), gdzie i = 0, 1, 2, ..., n; gdzie n to liczba boków wielokąta (lub wierzchołków), który opisują te punkty.

2. Punkty muszą być ułożone w porządku zgodnym kierunkiem ruchu wskazówek zegara lub odwrotnym, tzn. nie może być tak, że będą w kolejności dowolnej. Bryła musi wychodzić wypukła, w tym sensie, że żaden bok wieloboku nie powinien przecinać żadnego innego boku. Jeżeli punkty są ułożone zgodnie z kierunkiem ruchu wskazówek zegara, suma wyjdzie ujemna. Przeciwnie do kierunku ruchu wskazówek zegara: dodatnia.

3. Współrzędne będziesz miał w postaci liczby zmiennoprzecinkowej (zapomnij o postaci "stopnie minuty, sekundy"), ale te liczby będą wyrażały współrzędne wyrażające wartość kąta. W sumie można to policzyć, ale w jakich jednostkach wyjdzie ci pole powierzchni, stopniach kwadratowych? Żeby uzyskać pole w kilometrach kwadratowych, albo nm kwadratowych, musisz najpierw przeliczyć wartości współrzędnych w stopniach, na wartości odpowiadające tym współrzędnym, ale wyrażone w nm lub km odległości od punktu (0,0) wsp geograficznych. Pewnie zrobisz to w MS Exelu, lub w jakimś Calcu z pakietu Libre lub Open Office, zatem będziesz miał w dwóch kolumnach współrzędne x i y kolejnych punktów.

4. Wzór wygląda w sposób następujący:

Ale... ten wzór zdaje się nie uwzględniać zboczenia nawigacyjnego - nigdzie nie widzę cosinusa szerokości geograficznej... Tak jak pisałem wcześniej - odwzorowanie Merkatora jest wprawdzie wiernokątne, ale bynajmniej nie wiernopowierzchniowe...

_________________
wolne publikacje nautyczne

***********************************************
* Gdyby wszyscy byli bogaci nikt nie chciałby wiosłować *
***********************************************

Tak, ten wzór jest po prostu matematyczną formułą na układ kartezjański z założeniem że jest to po prostu powierzchnia płaska i nie uwzględnia zboczenia nawigacyjnego. Jeżeli ci zależy to poszukam jakiegoś wzoru na układ sferyczny.

Zapodałem ten wzór, bo nie kojarzę tak na szybko wzoru na wycinek powierzchni sfery liczony ze współrzędnych punktów na tej sferze. Postaram się czegoś poszukać.

Od bidy jeżeli obszar nie jest zbyt duży w porównaniu do promienia ziemi, to da się to troszkę nagiąć i oszacować z minimalnym błędem. Mianowicie przeliczyć wartości współrzędnych punktów, ustalając jeden z nich za (0,0), a współrzędne pozostałych przeliczyć do tego nowego układu biorac odległości z uwzględnieniem zboczenia nawigacyjnego dla danej szerokości geograficznej. Dopiero później wrzucić to w ten wzór. Nadal wynik będzie zawierał błąd, ale jeżeli stosunek znaczącego wymiaru liczonego obszaru do promienia ziemi będzie bardzo mały, to błąd też będzie niewielki.

Tymczasem poszukam czegoś bardziej bezpośredniego.

Pozdrawiam.

Zaproponuję dość łatwy i szybki sposób na policzenie pola powierzchni i obwodu obszaru bez liczenia na piechotę w excelu.

Polega on na stworzeniu wieloboku w "Twoich mapach" w google.maps.

1. Wchodzisz w google maps i wybierasz moje miejsca:


Następnie "Mapy":


Tutaj wybierasz jakąś mapkę, którą miałeś utworzoną wcześniej, albo tworzysz nową. Zakładam że tworzysz nową:


Jezeli wybrałeś istniejącą, to musisz Kliknąć "Otwórz w moich mapach":


Jeżeli utworzyłeś nową, to będziesz od razu mógł edytować warstwy. Dodajesz warstwę:


Jak już masz warstwę, to możesz do niej zaimportować instrukcje utworzenia twojego obszaru, ale żeby to zrobić musisz mieć przygotowany plik csv. To przykładowa zawartość takiego pliku:



Z tego co masz w tym przykładzie, wymieniasz tylko:
1. współrzędne punktów (muszą być w takim formacie jak w zamieszczonym powyżej przykładzie, zwróć uwagę na kolejność lon i lat),
2. ostatni punk powinien być powtórzeniem pierwszego,
3. ostatni człon, będzie nazwą kolekcji (tutaj "S-6")

Zapisujesz ten plik pod jakąś nazwą z rozszerzeniem .csv, na przykład "Obszar S-6.csv". Nazwa pliku po zaimportowaniu go, stanie się nazwą warstwy, do której został zaimportowany.

CD w następnym poście

_________________
Adam O.

Importowanie pliku csv:


Klikasz "Wybierz plik z urzadzenia":


Zaznaczasz WKT oraz "Dalej":


nazwa i "Zakończ":


Powinieneś zobaczyć na mapie twój polygon. Klikasz w niego i w infowindow na dole masz pole powierzchni oraz obwód:

_________________
Adam O.