Autor Tema: Povijest i GIS  (Posjeta: 11629 )

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
Povijest i GIS
« : 22 Listopad 2009, 16:03:22 poslijepodne »

Pretenciozan naslov inspiriran temom "Rimski puti i naselja".

Pokušaj precrtavanja trase jedne rimske ceste u programu Google Earth je zapeo.

Izvornik je karta u članku "Rimska cesta od Prologa preko Donjeg Unca i Petrovca u dolinu Sane" dostupnom na adresi http://www.idoconline.info/digitalarchive/public/serve/index.cfm?fuseaction=form&ElementId=404565 u .pdf formatu.

Greška koja se pojavljuje u preklapanju slike u Google Eartu je puno veća od mogućeg zbroja ulaznih grešaka i sinergije. Trasa rimske ceste ucrtana je na kartu mjerila 1:300.000, složena za tisak potkraj 19. stoljeća, otisnuta u mjerilu 1:300.000, skenirana nedavno.
Papirnata, nedigitalna karta odmah uključuje grešku od 300 metara. Jedan milimetar predstavlja tih 300 metara.
Google Earth pokriva naše područje u dvije tipične nazivne rezolucije (0.5 i 1.0 m) što znači da jedan piksel predstavlja pola metra ili metar ovisno o rezoluciji satelitske slike za koju se točni podaci mogu naći kod Google-ovog dobavljača mapa (Digital Globe).

Postupak "digitaliziranja" dijelova povijesnih (i drugih) karata uključuje problem kartografskih projekcija, mreže i datuma.
Prva greška s kartom iz primjera bila je u zanemarivanju suženja meridijana. Pošto otisnuta karta nema podatke o zemljopisnoj duljini i širini napravljen je slijedeći pogrešan pokušaj georeferenciranja karte.

Karta je izrezana programom za obradu slika (GIMP) tako da pravokutnik za isijecanje prolazi kroz Banjaluku, Knin i Split tako da imamo rubne točke pravokutnika za unos koordinata u Google Earth. Grube ulazne podatke o koordinatama tih mjesta uzmemo sa http://www.geonames.org i kasnije fino podešavamo.
U finom podešavanju, između ostalog, zaboravljeno je sužavanje meridijana. Pokušaj da se nađu podaci o kordinatama za "izvornu" austrijsku mapu (1:300.000) na http://lazarus.elte.hu/hun/digkonyv/topo/3felmeres.htm i slični kartografskim web mjestima nisu urodili plodom.

Greške u našim digitaliziranim kartama se jasno vide u Google Earth i na maps.google.com gdje digitalizirani vektorski sadržaji odstupaju od satelitske rasterske slike.
Postoji stari hoax o pomaku od 400 metara na zapad. Podaci potrebni za ispravnu višeparametarsku transformaciju su kriptični i teško dostupni.
Tipičan izvor grešaka opisan je u mail-u: http://lists.maptools.org/pipermail/proj/2009-March/004485.html
...

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
Odg: Povijest i GIS
« Odgovori #1 : 22 Listopad 2009, 18:05:27 poslijepodne »
Za rekonstrukciju starih komunikacija zanimljive su stare karte, bez obzira na nepreciznost.
The Maps of Heinrich Kiepert | Ottoman Empire, 1867, sheet 4
http://www.lib.uchicago.edu/e/su/maps/kiepert/G7430-1867-K547-sheet4.html

http://www.davidrumsey.com ima gotove .kmz datoteke za Google Earth.
Preklapanje karte Italy 1800 daje neke zanimljive ideje.
Npr. o ubifikaciji Bilubia u Lunbio. Nema veze što je preklopljeno tako loše da je stvarni Sinj skoro na mjestu Salone.

Molim belevarca ako zaluta na ovu temu da mi odgonetne toponim "Schich Seria" na karti
http://www.davidrumsey.com/luna/servlet/detail/RUMSEY~8~1~33566~1171058:General-Karte-der-Europaischen-Turk?qvq=w4s:/where/Bosnia+and+Hercegovina;lc:RUMSEY~8~1&mi=0&trs=2
(malo sporo se učitava).
« Zadnja izmjena: 22 Listopad 2009, 19:01:29 poslijepodne glonga »

Offline Slovak

  • Gožđeno pero foruma
  • ***
  • Postova: 3.393
  • Bez cezure.
Odg: Povijest i GIS
« Odgovori #2 : 22 Listopad 2009, 19:36:16 poslijepodne »
Molim belevarca ako zaluta na ovu temu da mi odgonetne toponim "Schich Seria" na karti


Ćićarija, daleko etimološki, ili raščlanjenjo - serija šiša. šiša je tavan. gomila brežuljaka. ;)
jer na karti nisam našao pojam, možda i cincar-planina.
Nisam umislio da sam belevarac, vidit ćemo šta čovik znade. samo nek mi ne govori da sam polupa lončiće jer unada oden tamo u kujinju i sve polupan, koliko to mene zaboli.
« Zadnja izmjena: 22 Listopad 2009, 19:39:39 poslijepodne Slovak »
Dosta parola.

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
Odg: Povijest i GIS
« Odgovori #3 : 22 Listopad 2009, 21:11:38 poslijepodne »

Još jedno web mjesto ima dobre stare karte za Google Earth u .kml/.kmz formatu:
Citat:
With funding from the Institute of Museum and Library Services (IMLS), the team scanned and geo-referenced the map series "Spezialkarte der Oesterreichisch-ungarischen Monarchie im Masse 1:75,000 der Natur" (Special Map of the Austro-Hungarian Empire at the scale of 1:75,000). This set was published in Vienna over the period, 1877-1914. The set contains 776 sheets, with all editions of all sheets numbering over 3,665. Through a search engine, the maps can be viewed through a web browser or through Google Earth.
http://www4.uwm.edu/libraries/AGSL/auhu.cfm
http://imlsmap.lib.uconn.edu/search_engine.html

Slovače, ona mapa ima zoom (povećaj) i pan (pomiči) opcije, pa možeš doć do Imackog i malo na istok. Ne mogu ovdje umetnut sliku sa zaslona zbog limita forumskog, a postavljanje linka na tu sliku zaslona na nekom online servisu za postavljanje fotografija nema smisla jer je slika već dostupna, kartu treba samo pomaknuti da prikazuje naše područje i uvećati.
Ma evo je na ovom linku http://img406.yfrog.com/img406/654/scich.jpg
Schich Seria je bagavo zaokružena plavozelenim.

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
Odg: Povijest i GIS
« Odgovori #4 : 22 Listopad 2009, 22:35:20 poslijepodne »

Još jedno korisno web mjesto.
http://english.mapywig.org/Austro-Hungarian.html
a trebalo je naći skelet od Generalkarte von Central-Europa 1:300 000

Offline likarija

  • Gožđeno pero foruma
  • ***
  • Postova: 2.599
  • Nu me!
    • Oldtimer klub Zagreb
Odg: Povijest i GIS
« Odgovori #5 : 23 Listopad 2009, 09:19:21 prijepodne »
Pretenciozan naslov inspiriran temom "Rimski puti i naselja".
Pokušaj precrtavanja trase jedne rimske ceste u programu Google Earth je zapeo.
...
Greška koja se pojavljuje u preklapanju slike u Google Eartu je puno veća od mogućeg zbroja ulaznih grešaka i sinergije. Trasa rimske ceste ucrtana je na kartu mjerila 1:300.000, složena za tisak potkraj 19. stoljeća, otisnuta u mjerilu 1:300.000, skenirana nedavno.
Papirnata, nedigitalna karta odmah uključuje grešku od 300 metara.

Imaš problem odmah u startu.
Uz to što je, kako si i sam naveo:
Postupak "digitaliziranja" dijelova povijesnih (i drugih) karata uključuje problem kartografskih projekcija, mreže i datuma.
Prva greška s kartom iz primjera bila je u zanemarivanju suženja meridijana. Pošto otisnuta karta nema podatke o zemljopisnoj duljini i širini napravljen je slijedeći pogrešan pokušaj georeferenciranja karte.


Problem je i u tisku ("printanju") i u scan-iranju.
Oba u startu "lažu" i 10 % (fotokopiranje i više"). Tu nisam računao kopiju kopije.
Da bi se bar to eliminiralo, treba scan-er s georeferenciranjem,
odnosno da točno prenese bez iskrivljavanja, pa kad se riješi ta greška u ulazu nastaviti kako si sam objasnio.

:starac

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
Odg: Povijest i GIS
« Odgovori #6 : 23 Listopad 2009, 23:59:06 poslijepodne »
Što prigovoriti starim kartama?

maps.ggogle.com osim što se utopio u Jadranu na primjeru pomorskih ruta s naših obalnih i otočnih luka do luka na poluotoku Garagano, zagubio se i na rijeci Sani.
Ako pođete iz Sanskog Mosta na jug jedinom digitaliziranom cestom bez broja doći ćete opet u još jedan Google-ov Sanski Most koji je zapravo Ključ na digitaliziranoj cesti označenoj brojem 11-01 pučanima poznatoj kao put AVNOJ-a koja iz pravca Bihaća i Petrovca prolazi kroz Ključ, Mrkonjić i trebala bi do Jajca međutim digitalizirana cesta 11-01 ničim izazvana kod Plivskog Jezera ode na vrh brda i taj njezin kraj povezuje se makadamskim nedigitaliziranim serpentinama na cestu E761 za Titovo Lijevo Jajce.

A predmetna karta sigurno nije tkz. Scheda-Karte (1873−1876 „Generalkarte von Central-Europa“, 47 Blätter im vergrößerten Maßstab 1:300.000, bekannt als „Scheda-Karte“) - ne uklapa se u MGI mrežu za te karte.

-
Češka zbirka austrougarskih digitaliziranih karata: http://mapy.vkol.cz
« Zadnja izmjena: 24 Listopad 2009, 00:30:37 prijepodne glonga »

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
Odg: Povijest i GIS
« Odgovori #7 : 26 Listopad 2009, 01:28:13 prijepodne »

Štivo i alati


Netočnost konverzije
http://spatial-analyst.net/wiki/index.php?title=MGI_/_Balkans_coordinate_systems

Formule za konverziju Transverse Mercatora s primjerom
http://www.scribd.com/doc/4533298/GPS-Basics-GPSX02007

Australska datoteka tabličnog kalkulatora za konverziju bez makro-naredbi
http://www.icsm.gov.au/gda/gdatm/transxyz.xls

Formule za konverziju Transverse Mercatora s primjerom
http://posc.org/Epicentre.2_2/DataModel/ExamplesofUsage/eu_cs34h.html
Nova verzija je u dokumentu G7-2.pdf

GRASS GIS
http://grass.bologna.enea.it/

Quantum GIS
http://qgis.org/index.php

Parametri za copy-paste u onoj australskoj datoteci.

Bessel 1841
6377397,155
299,15281285

ETRS89
6378137
298,257222101

B>E
551,70
162,90
467,90
-6,04
-1,96
11,38
-4,82

E>B
-464.2453
-116.2612
-438.2892
-4.455840
-3.146930
11.613150
-9.151447

"Novi" B>E
-550.4985
-164.1161
-475.1416
-5.80967190
-2.07901633
11.62385702
5.54176398
(ponovo provjeriti +/-, ovaj put i za translaciju)

Z7 B>E
574.027
170.175
401.545
-4.88786
0.66524
13.24673
6.88933

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
Simple Block Shift & Google Earth
« Odgovori #8 : 10 Prosinac 2010, 22:40:52 poslijepodne »
Vut?
« Zadnja izmjena: 10 Prosinac 2010, 23:07:08 poslijepodne glonga »

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
Odg: Povijest i GIS
« Odgovori #9 : 10 Prosinac 2010, 23:47:26 poslijepodne »
Vut? TF?

Error.
Nemam više volje ni živaca, ni vrimena.
0.47693
-17.73249
0.48
-17.73
« Zadnja izmjena: 10 Prosinac 2010, 23:59:21 poslijepodne glonga »

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
BOFH-u (ptičar)
« Odgovori #10 : 04 Veljača 2012, 22:13:24 poslijepodne »

I da je vidimo komu se podoba.

Vraćaju se Suljo i MuLJo iz lova s ptičarom. Prazni ruku.
— Ovaj ćuko ili nije tičar, ja li smo ga nisko bacali.

Nu na stvar. QGIS-u 1.7.3 (Wroclaw) ne virujen na Mac-u. Slutim po ponašanju da ne radi kako triba. Nije kompajliran iz izvornog kôda već po kyngchaos proceduri.
Sumnjam u python koji se negdje zavadio s sistemskim python-om, ko u nekim slučajevima perl i 32-bitna java.
Nu, ne virujem da GDAL i PROJ4 u red rade i na Ubuntu-u (10.4 LTS) s lege artis instaliranim QGIS-om 1.7.3. GRASS mi je ovih dana dok se ne očiste teste overkill.
Osnovano sumnjam da ti je problem u nedefiniranom datumu što popriko moš riješit instalacijom ILWIS-a ili uzduž kroz loše dokumentiranu sintaksu PROJ4 (i GDAL).
ILWIS koliko znam imaš samo za windoze. Ne vidim da imaju dostupan source code, samo .exe binary install. Pogubim se u ovim licencama, znači li no GNU GPL u licenci i Open u imenu - otvoren i dostupan izvorni kôd?
Otkud sumnja u izvor greške? Tomo Hengl!
Oprobaj na http://twcc.free.fr/
Pa ako ti ne valja vidi što veli Hengl:
»Note that the definition of the MGI / Balkans coordinate systems (EPSG:31275, EPSG:31276, EPSG:31277) available via http://spatialreference.org is incomplete because this database does not contain any definition of the geodetic datum. If the geodetic datum is missing, the difference from the true location will be even few hundreds of meters.
 […] even if you correctly set-up the (global) coordinate system parameters, you will not be able to achieve precision better than 0.5 meters«*
Uzput, parametri T7D* su službeni od 31. 3. 2011.*

Za moje potribe básanja po šumama i gorama glede gljiva i diviji voćaka i druge prirode brez društva, godinama na starinskom ručnom robustnom uređaju koristim obični pomak (shift) od 17.7" (sedamnaest točka sedam lučnih sekundi) na zapad jerbo isti nema "naprednu" maglu od datuma. To ti je zbog suženja meridijana na području na kojem básam pomak od do 400 metara na zapad.
Overkill mi je namišćat druge pomake pa i onaj nautički koji vridi južnije, a vridijo bi i za moje područje i za moju potribnu razinu točnosti. Taj pomak je u smjeru sjevera 0.01' (jedna stotinka lučne minute na sjever tj. 0.6" (nula točka šest lučnih sekundi) i 0.28' - dvadesetosam stotinki lučne minute na zapad tj. 16.8" (šesnaest točka osam lučnih sekundi) tako da mu to vektorski zbrojeno uz magnetni odklon sjevera daju što se mene tiče isto (i oprobano po srednjem Jadranu uz normalan oprez u zoni cca. 400 metara od problema).
Uzput, mafiji devetičnog ostatka je na ruku išlo (ili obratno) to što je Garmin pogriješio datum, zapravo mu je datum u redu ali je primijenjen na krivi elipsoid. Usporedi parametre elipsoida u Garmina i na NIMA (.mil) s tro(pet)parametarskom transformacijom. U šumu dostupnih netočnih parametara ubačena je i nedosljedna nomenklatura elipsoidâ, Bessel, Bessel 1841, Bessel Namibia, Bessel Etiopia… u softveru ručnih uređaja je posebna magla definicija Hermannskogel datuma - ko bombonjera - nikad ne znaš šta š dobit. Stvar se zove Hermannskogel, a koje su joj sastavnice to ti ne će nijedan proizvođač reć.

Imaju li smisla desetinke lučnih sekundi u slijedećoj tablici, vidi sâm.

#GK6 x [m]GK6 y [m]MGI lon E (I) º ' " MGI lat N (S) º ' " WGS84 E (I) lon º ' " WGS84 lat N (S) º ' "
0016420000481800017 00 38.143 30 07.317 00 20.643 30 07.5
0026429000481800017 07 18.843 30 10.517 07 01.243 30 10.8
0036429000483100017 07 12.743 37 11.817 06 55.143 37 12.0
0046420000483100017 00 31.243 37 08.517 00 13.743 37 08.7
0056419141481778617 00 00.043 30 00.016 59 42.543 30 00.3
0066429248481767217 07 30.043 30 00.017 07 12.443 30 00.3
0076429394483155717 07 30.043 37 30.017 07 12.443 37 30.2
0086419308483167117 00 00.043 37 30.016 59 42.543 37 30.2

Meni je i na tom području bio dovoljan shift od -17.7" lon. Razlike moš vidit u svakom redku.
Ono što je meni na ovoj temi problem još uvijek - greške koje dobijem za karte s Ferro meridijanom unatoč dostupnim podatcima. Ne vjerujem da su kartografi toliko griješili u vrijeme uporabe Ferro meridijana. Ili jesu?
« Zadnja izmjena: 04 Veljača 2012, 22:21:06 poslijepodne glonga »

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
puno teksta ni za što
« Odgovori #11 : 05 Veljača 2012, 06:08:05 prijepodne »

Sva ta buka oko Helmertove 7-parametarske transformacije mi odavno smrdi. Još od transxyz.xls modifikacije. Nikad se njome nije postigla preciznost veća od tkz. troparametarske.
Program MacGPSpro (i niz GPS uređaja) koriste "Croatia" datum i troparametarsku transformaciju i meni je to dovoljno precizno za moje potrebe. Podatci o tome su slijedeći (obrati pozornost na naziv elipsoida: Bessel 1841, Bessel 1841 Namibia, a da doista jest Namibia - jest a=a(BN), f= f(BN) BN=Bessel_Namibia, PROJ4 stringovi imaju bessel i bess_nam i onaj treći norveški koji se (još) nije upetljao u naške transformacije ili ja ne znam za to):
Datum name: Croatia
Ellipsoid: Bessel 1841 / a=6377483.865 m / f=1/299.152812800
Geocentric displacement (m):
dx=653
dy=-212
dz=449

Navigo je svojedobno mantao s elipsoidom (modificirana Etiopija) za:
za eTrex DA=+740m*
za GPS III isto*
za GPSMAP 76/76s DA=745m *
etc.
Međutim, ne moš im implicirat (veliko) neznanje i zlu nakanu jerbo na pathu /Library/Frameworks/GDAL.framework/Resources/elipssoid.csv piše i:

»7004,Bessel 1841,6377397.155,9001,299.1528128,,1,"Original Bessel definition is a=3272077.14 and b=3261139.33 toise. This used a weighted mean of values from several authors but did not account for differences in the length of the various toise: the ""Bessel toise"" is therefore of uncertain length.","US Army Map Service Technical Manual; 1943.",OGP,1999/04/22,1998.321 1998.340,0
7005,Bessel Modified,6377492.018,9001,299.1528128,,1,Used in Norway and also in Sweden with a 1mm increase in semi-major axis.,,OGP,1999/04/22,1998.321,0
7006,Bessel Namibia,6377483.865,9001,299.1528128,,1,a = 6377397.155 German legal metres. This is the same value as the Bessel 1841 figure (code 7004) but in different units.  Used in Namibia.,"Chief Directorate: Surveys and Mapping, Mowbray, South Africa.",OGP,1999/04/22,1997.160,1«

NIMA, sad NGA ima slijedeće:
Hermannkogel*
Ellipsoid: Bessel 1841
da=739.845m
df x 104=0.10037483
dx=682
dy=-203
dz=480

WGS84 elipoid ima a=6378137 m pa se u ova dva potonja slučaja radi o BR (Etiopia) a ne BN (Namibia) Besselovu elipsoidu.
da=a(WGS84)-a(BR)=6378137-6377397.155=739.845m (Bessel 1841 / BR/Ethiopia)
da=a(WGS84)-a(BN)=6378137-6377483.865=653.135 m (Bessel 1841 / BN/Namibia)
skoro pa kao dx iz troparametarskog seta za isti elipsoid.

Što su ti veliki brojevi? Na dva elipsoida skoro pa isti rezultati kad se gađa određena točnost transformacije.

Sedamparametarska transformacije i navodna točnost o kojoj trube naši državni i paradržavni geodeti mi je meni overkill, a ne mi daje ništa bolju točnost od troparametarske bez dodatka koji se zove T7D ili GRID datoteka. To se plaća.
Čini mi se da jedino još Kanada naplaćuje grid ((malim slovom) projekcijsku mrežu).
Meni je provjera overkil zbog geocentrične Z koordinate, ne da mi se gnjavit s nadmorskom visinom (koju mi GPS ni ne prikazuje, već onu drugu), geoidom, elipsoidom, normalnom ortometrijskom… ali i u uvjetima kad se imao Z i njegove sastavnice - sedamparametarska je uredno promašivala službeno obećanu točnost.
Još dok je na linkovima web mjesta crs-geo.eu / crs.bkg.bund.de* bilo i točaka za provjeru uz ondašnjih 7 parametara na temelju ondašnjih 120 identičnih točaka.
dx=+551.7 m
dy=+162.9 m
dz=+467.9 m
rx=+6.04" (!to be in agreement with formulas the rotation parameter has to be converted to Radians)
ry=+1.96"
rz=-11.38"
scale -4.82 ppm

U slučaju PROJ4 rotacije su u lučnim sekundama a skala je posebna priča. Skala za gornji slučaj -4.82 ppm za unos u PROJ4 se dobije;
M_BF=(1000000-4,82)/1000000=0.99999518
pazi to nije k parametar u PROJ4 već M_BF.
»The seven parameter case uses delta_x, delta_y, delta_z, Rx - rotation X, Ry - rotation Y, Rz - rotation Z, M_BF - Scaling. The three translation parameters are in meters as in the three parameter case. The rotational parameters are in seconds of arc. The scaling is apparently the scale change in parts per million.«*

Službeni T7D* parametri su:
tX = 546.70776
tY = 162.37348
tZ = 469.53683
rX = -5.90560751
rY = -2.07314165
rZ = 11.51062442
dM = -4.45886242

Da. Ali za T7D. Za PROJ4 rotacijama treba promijeniti predznak i presložiti dM.
Tako da je PROJ4 korisnički string za transformaciju i opis korisničkog CRS-a npr. zone GK6 slijedeći:

+proj=tmerc +pm=greenwich +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bessel +towgs84=546.70776,162.37348,469.53683,5.90560751,2.07314165,-11.51062442,0.99999554113758 +units=m
QGIS i drugi ti za to koriste EPSG:31276 MGI / Balkans zone 6* čiji je PROJ4 string* smrljan:
+proj=tmerc +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bessel +towgs84=577.326,90.129,463.919,5.137,1.474,5.297,2.4232 +units=m +no_defs
Hengl* veli:
The MGI / Balkans Zone 6 in the PROJ.4 should be set as:
+proj=tmerc +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bessel +towgs84=550.499,164.116,475.142,5.80967,2.07902,-11.62386,0.99999445824 +units=m
http://twcc.free.fr/ također koristi EPSG:31276 i druge stringove sa spatialreference.org, a meni se ne da ić tamo postavljat "ispravne", "svoje" user stringove.
Otkud im samo 2.4232 na mjestu parametra M_BF?
Uoči gdje je (sve) maknut problematični parametar +no_defs.

Dok se snig ne raščisti nemam točaka osim one T249 sa spomenute Henglove web stranice i ne da mi se radit dvosmjerna provjera.

Na našem području Google Earth i Google Maps imaju pomaknutu satelitsku sliku, na Google Earthu su čak i povijesne satelitske slike (tamo gje ih ima) međusobno smaknute po 30-50 metara.

S onih tvojih podloga za točke dobijem skoro iste rezultate prikaza na Google Earthu i  Google Maps točaka konvertiranih troparametarskom i sedamparametarskom transformacijom.
Troparametarska je narečena (prva iz ovog upisa) iz MacGPSPro (kao i OziExplorer, Fugawi i sl.) i to su "neozbiljni" GPS programi. Prikaz na Googleovim online satelitskim slikama za takve točke mi je na našem području oko trideset metara na azimutu 200-300 stupnjeva od slike (mjereno u Google Earthu, npr. do crkvenih vrata od prikaza točke, još ako je zvonik uz vrata pa se vidi više satelita mwahahaha).

Sedamparametarska ne da ništa bolje rezultate ni s ovim najnovijim i službenim parametrima, tako da je to za me i dalje overkill i gnjavaža s mafijom devetičnog ostatka. T7D program ni GRID datoteku s distorzijama ne kanim kupovat ni da imam viška novaca ni da mi ih nude za šankom u birtiji. Ne znam što bih s tim. Ni da sam na tvom mistu.

QGIS jest malo bugovit. Ponekad ga triba pričekat. Nekad mu triba ponovit. Često mu treba restart. Postupak. Georeferencirat u gk6 (u metrima). Postavit na glavnim postavkama "transformaciju CRS u letu". Izvoz georeferencirano u gk6 (u metrima). Tu u postavkama izvoza Helmert nema veze s imenom 7-parametarske transformacije. Ima samo je semantički pomaknuto. Uglavnom nije to transformacija CRS-a nego slike. Gotovo georeferenciranje. (Ludi QGIS mi je u prvom pokretanju po polinomu 3 stupnja uz kubično resampliranje dva dana pravio datoteku veliku skoro terabajt. Mo š mislit skočne obavijeti "Destinatiom disk allmost full." ko da je 1998. Poslije mu uz iste postavke nije takvo što padalo na pamet.) Dodat rasterki sloj - tiff datoteku u CRS gk6. Pokrenut dodatak "Čitanje koordinata" i na njegovu panelu škljocnut "Započni čitanje" pa gađat točke. Ako na panelu ne prikazuje i stupnje i metre u odgovarajućim prozorima panela (u statusnoj traci je uvijek samo jedna vrsta prikaza), pokušat škljocat odgovarajuće ikone u donjem desnom kraju statusne trake (Prikaži EPSG:4326 ukl/isklj, Prikaži USER:100000 uklj./isklj., škljoc na ikonu "CRS status" u donjem desnom kutu). Nekad fercera, nekad nije.

Mislim da vam je bilo do krivog PROJ4 stringa iz EPSG:31276 i drugih, nakon toga i do problematičnog no_defs parametra. Uostalom, mislim da je sve to overkill.
OziExplorer ruljz. Za moje potrebe.
Samo jednokratno napišeš redak teksta i spremiš ga u datoteku users.dat. Zapravo datums.dat. User.fd je me š čini u Fugawija. Khmda, Billov rikverc slash path \Program Files\Common Files\Fugawi\user.fd*
Croatia, Bessel 1841, 739.845, 0.10037483, 634.12, -216.43, 436.10 ne š ti razlike do NIMA/NGA
Croatia, Bessel 1841, 739.845, 0.10037483, 682, -203, 480
Slično Ozi-ju fercera i na Androidu androzic*.

MGI, 4, 682, -203, 480
(4 - redni broj elipoida Bessel 1841 u programu androzic, u OziExploreru je 3)

MacGPSpro, Garmin, Magellan… ovaj datum za 3-param transformaciju zovu "Croatia", "Croatian", ali ko Tadiću za Boškovića, to Peđi Supuroviću ne smeta da na više mjesta koristi naziv "Hermannskogel Serbia"* pa i u datoteci datums.dat*
Hermannskogel Serbia, 3, 653, -212, 449za koju ne piše da je nema u korijenskoj Ozi-jevoj mapi i da ju winblowz user može napravit čuj sad Notepad-om …\OziExplorer\datums.dat.
Pedja, Pedja, praviš hrkljuš od elipsoida i parametara, ne veruješ američkoj vojsci ni australskim programerima (OziExplorer, androzic, Fugawi):
MGI, 3, 682, -203, 480Druga je stvar što čarobni slijed brojeva i preciznost (ručnih) uređaja daje skoro iste rezultate na tri različita Besselova elipsoida. Npr. u androzicu probaj s istog mjesta očitat poziciju na tri različite postavke elipsoida (a sva tri Bessel 1841, dušu dalo za zabune, toise):
User, 4, 653, -212, 449User, 22, 653, -212, 449User, 28, 653, -212, 449
Bessel napravi elipsoid i definira ga u toisama. Svaka država uzme sebi svoju kartografsku toise i dobije drukčiji elipsoid. Odraz razlike u veličini toise je taman takav da topnik jedne vojske ne može pogodit cilj ni sa svoje ni s tuđe karte na tlu koje nisu kartirali i premjerili njegovi kartografi (najčešće ni na tom). Pa ti danas radi digitalne 3D modele povijesnih zgrada na takvim geospacijalnim podatcima. Kani se ćorava športa.

I što? Hengl veli ovdje*, a i ja mislim da je tebi najbolje ovo napravit ako ti već treba preciznost za tu spacijalnu arheologiju:
»Hence the difference is in fact significant (but can be ignored for more general scales e.g. <1:25000). The only way to increase the precision of the datum shift parameters would be to re-estimate them locally using local referent points. This can be done for example by using 3-4 local reference points and the find datum transformation parameters function in ILWIS.«

Svejedno će ti 3D prikaz povijesne građevine bit 30-50 metara dalje od satelitske slike mjesta u Google Earthu, barem dok Google ne poštela svoje satelitske slike na prava mjesta u ovom dijelu Europe. I znaš šta? Google ti ne će objavit koliki je pomak bio kad to poštelaju. Pa š ondak ponovo modelirat 3D prikaz, ili pogađat tri ili sedam novih parametara za pomak 3D prikaza, samo š sad bit u okviru od 40 a ne više u okviru od 400 metara, što je ipak - red veličine. He. Hehe.
« Zadnja izmjena: 05 Veljača 2012, 06:51:50 prijepodne glonga »

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
viduka
« Odgovori #12 : 06 Veljača 2012, 02:01:23 prijepodne »

BOFH. De kakvom svom đomeštru daj da pregleda malo ovaj blog:
http://vb-iostalebudalastine.com/
Ma kako se lik izražavao, ne mogu doć sebi od zbunjole da se svagdi u svitu stare projekcije i WGS84 slažu u milimetar u transformaciji, a ovuda digitalizirane ceste i granice parcela nisu u skladu ni š čim.

Kad pogledam digitalni katastar na arkod.hr počmem vjerovat da su Google Earth i Google Maps stvarno točno složili satelitske slikarije na WGS84 koordinate na našem području.
http://www.geo-portal.hr/Portal/ptk
http://preglednik.arkod.hr

Sad ti najozbiljnije preporučujem ono što veli Hengl - uzmi tri četiri poznate točke u blizini i sračunaj svoje transformacijske parametre.

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
grčka
« Odgovori #13 : 06 Veljača 2012, 12:49:14 poslijepodne »
Krčma Kunara je fizički nedostupna. A Grci?
Nešto nije u redu s rotacijama, mjerilom i njihovom obliku koji očekuju proj4 i gdal.
Ta četiri parametra iz fizičkog opisa CRS-a treba (izgleda) nekako preračunat.
Zasad RTFM za proj4 i gdal nisu od neke pomoći.
Otud valjda onaj 2.4232 sedmi parametar koji ostaje "sirov".
Zabilješka:
»The 7-parameter case uses delta_x, delta_y, delta_z, Rx-rotation X, Ry-rotation Y, Rz- rotation Z, M_BF - Scaling. The three translation parameters are in meters as in the 3- parameter case. The rotational parameters are not in physical units. They are something like the sine of the rotational angle times the ellipsoid axis length, but the exact details of how to derive these from a physical description of the rotation are not covered in this manual. The internal use of these parameters match the parameters used for transformation method 9606 in the EPSG database. The scaling is apparently the scale change in parts per million.
A more complete discussion of the 3 and 7 parameter transformations can be found in the EPSG database (trf_method's 9603 and 9606).«*pdf
Dostupna verzija baze u mySQL-u je smrljana, u M$Accesssu još gora. Prebrat još jednom G7-1 i G7-2.
Rotacijama je problematičan i (pred)znak. Nekonzistentni ulazni podatci.

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
proj4 & gdal nesporazum
« Odgovori #14 : 06 Veljača 2012, 14:10:04 poslijepodne »
Bilježim dalji mogući uzrok problema. Nedefiniran (krivo definiran) datum. Bessel. I što još? ništa. A Bessela ko u priči.
Parametar +wktext glede nesporazuma proj4 i gdal/ogr.

»Some parameters, such as grids, recognised by PROJ.4 may not be well understood and translated into the OGRSpatialReference model. It is possible to add the +wktext parameter which is a special keyword that OGR recognises as meaning "embed the entire PROJ.4 string in the WKT and use it literally when converting back to PROJ.4 format".«*

Dodatno:
Parameters:
dfDX - X child in meters.
dfDY  - Y child in meters.
dfDZ - Z child in meters.
dfEX    - X rotation in arc seconds (optional, defaults to zero).
dfEY    - Y rotation in arc seconds (optional, defaults to zero).
dfEZ    - Z rotation in arc seconds (optional, defaults to zero).
dfPPM - scaling factor (parts per million).
Prema ovom.
« Zadnja izmjena: 06 Veljača 2012, 14:27:07 poslijepodne glonga »