Autor Tema: Povijest i GIS  (Posjeta: 11630 )

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
GRASS
« Odgovori #15 : 06 Veljača 2012, 15:25:28 poslijepodne »

Gdje se još sukobljavaju proj4 i gdal, zapravo EPSG?
GRASS daje ideju o tome gdje treba tražit problem. Ovdje.
proj-proj4 je pravio Američanin.

»The parameters are commonly referred to defining the transformation ""from source coordinate reference system to target coordinate reference system"", whereby (Xs, Ys, Zs) are the coordinates of the point in the source geocentric coordinate reference system and (Xt, Yt, Zt) are the coordinates of the point in the target geocentric coordinate reference system. But that does not define the parameters uniquely; neither is the definition of the parameters implied in the formula, as is often believed. However, the following definition, which is consistent with the "Position Vector Transformation" convention is common E&P survey practice,
+
+(dX, dY, dZ): Translation vector, to be added to the point's position vector in the source coordinate reference system in order to transform from source system to target system; also: the coordinates of the origin of the source coordinate reference system in the target coordinate reference system.
+
+(Rx, Ry, Rz): Rotations to be applied to the point's vector. The sign convention is such that a positive rotation about an axis is defined as a clockwise rotation of the position vector when viewed from the origin of the Cartesian coordinate reference system in the positive direction of that axis; e.g. a positive rotation about the Z-axis only from source system to target system will result in a larger longitude value for the point in the target system. Although rotation angles may be quoted in any angular unit of measure, the formula as given here requires the angles to be provided in radians.
+
+M: The scale correction to be made to the position vector in the source coordinate reference system in order to obtain the correct scale in the target coordinate reference system. M = (1 + dS*10^-6), where dS is the scale correction expressed in parts per million.«


"Position Vector Transformation" vs. "Coordinate Frame Rotation"


» Although being common practice particularly in the European E&P industry, the Position Vector Transformation sign convention is not universally accepted.  A variation on this formula is also used, particularly in the USA E&P industry.  That formula is based on the same definition of translation and scale parameters, but a different definition of the rotation parameters. The associated convention is known as the "Coordinate Frame Rotation" convention (EPSG coordinate operation method code 9607).
 […]
+(dX, dY, dZ): Translation vector, to be added to the point's position vector in the source coordinate reference system in order to transform from source coordinate reference system to target coordinate reference system; also: the coordinates of the origin of source coordinate reference system in the target frame.
+
+(Rx, Ry, Rz): Rotations to be applied to the coordinate reference frame. The sign convention is such that a positive rotation of the frame about an axis is defined as a clockwise rotation of the coordinate reference frame when viewed from the origin of the Cartesian coordinate reference system in the positive direction of that axis, that is a positive rotation about the Z-axis only from source coordinate reference system to target coordinate reference system will result in a smaller longitude value for the point in the target coordinate reference system. Although rotation angles may be quoted in any angular unit of measure, the formula as given here requires the angles to be provided in radians.
+
+M: The scale factor to be applied to the position vector in the source coordinate reference system in order to obtain the correct scale of the target coordinate reference system. M = (1+dS*10^-6), where dS is the scale correction expressed in parts per million.
+
+In the absence of rotations the two formulas are identical; the difference is solely in the rotations. The name of the second method reflects this.
+
+Note that the same rotation that is defined as positive in the first method is consequently negative in the second and vice versa. It is therefore crucial that the convention underlying the definition of the rotation parameters is clearly understood and is communicated when exchanging datum transformation parameters, so that the parameters may be associated with the correct coordinate transformation method (algorithm).","The same example as for the Position Vector Transformation (coordinate operation method 9606) can be calculated, however the following transformation parameters have to be applied to achieve the same input and output in terms of coordinate values: […] «
« Zadnja izmjena: 06 Veljača 2012, 15:34:52 poslijepodne glonga »

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
Helmert nije Bursa-Wolf(e)
« Odgovori #16 : 06 Veljača 2012, 16:01:04 poslijepodne »

NATO/STANAG 2211 jasno kaže što je Helmertova 7-parametraska transformacija.
Objavljene specifikacije (službene i neslužbene) CRS-ova i njihovih transformacija navode Helmertovu metodu, što nije točno.
Uz sve zablude oko elipsoida, datuma… još i to.
»The seven-parameter Bursa-Wolf transformation, which is widely applied in Europe, uses the same formula, but the meaning of the rotation parameters is different.«*
Nikako u osnovnoj proj i dopunskoj proj3 (i dalje) dokumentaciji doć do ovog, a lino provlačit račun.

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
0.001"
« Odgovori #17 : 06 Veljača 2012, 16:56:57 poslijepodne »
Prema slijedećem Z i elipoidalne visine ne bi trebale bit uzrokom grešaka.
Nu "Expected results:" i Viduka narečeni. Tisućinka lučne sekunde
Ostaje testiranje postojećih dok se ne mogne ponovo do fizičkih točaka glede sniga.

Iztržci iz:
http://www.ogp.org.uk/pubs/430-3.pdf

Test procedure: GIGS 5201. Geographic/geocentric conversions.
Test purpose: To verify that conversions for the geographic/geocentric method are consistent with EPSG method 9602 (Geographic/geocentric conversions).

Results should agree to within 0.01m/0.0003” of the test data. See file GIGS_tfm_52xx_output_ [version_date].xls, worksheet 5201 GeogGeocen.
For the round trip calculation the initial coordinates of the point should change by less than 6mm/0.0002” after 1000 iterations. Test result will be pass or fail. If fail, details of failure should be recorded.

Expected results: Results should agree to within 0.01m/0.0003” of the test data. See file GIGS_tfm_52xx_output_ [version_date].xls, worksheet 5201 GeogGeocen.
For the round trip calculation the initial coordinates of the point should change by less than 6mm/0.0002” after 1000 iterations. Test result will be pass or fail. If fail, details of failure should be recorded.


Test procedure: GIGS 5203. Position Vector (geographic domain) transformations

Test purpose: To verify that transformations for the Position Vector 7-parameter transformation method which operate from and to geographic coordinates are consistent with EPSG methods. [Either the Position Vector transformation (geog2D domain), method 9606 or Position Vector transformation (geog3D domain) method 1037 is acceptable.]

Expected results: Results should agree to within 0.001” of the test data. See file GIGS_tfm_52xx_output_[version_date].xls, worksheet 5203 PosVec.
For the round trip calculation the initial coordinates of the point should change by less than 0.0002” after 1000 iterations. Test result will be pass or fail. If fail, details of failure should be recorded.
Some geoscience software ignores ellipsoidal heights during these transformations (or set the input ellipsoidal heights to zero). If that is the case, results will match the results as shown for the geog2D case, EPSG method 9606. Horizontal coordinates obtained for those points with ellipsoidal heights significantly different from zero will be incorrect, whereas correct results may be generated for points with zero (or near zero) ellipsoidal heights.
For large ellipsoidal heights (either positive or negative), the correct results are given by the geog3D EPSG method 1037
Which results are matched should be clearly documented in the report on the test results.

Ctrl+F 9606
Ctrl+F 9607

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
grid
« Odgovori #18 : 06 Veljača 2012, 17:43:55 poslijepodne »

Ako zboksaš grid datoteku, uz nju neš dobit dokumentaciju potrebnu (i dovoljnu) za uporabu u proj4/gdal.
Njemačka.
Brezplatni njemački grid.
http://crs.bkg.bund.de/crseu/crs/descrtrans/BeTA/de_dhdn2etrs_beta.php
Icrpna njemačka dokumentacija. Deutsche Baustelle.
http://crs.bkg.bund.de/crseu/crs/descrtrans/BeTA/BETA2007dokumentationV14.pdf
Ne znači da će bit primjenjiva. Kriva srastanja. Namjerna?

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
1901.
« Odgovori #19 : 06 Veljača 2012, 19:46:52 poslijepodne »
Uglavnom, glede stoljetnog šlamperaja sve podatke iz razdoblja 1901. (1918.) do pojave subcentimetarskog GPS-a moš bacit kako sugeriraju.

Lapaine, Tutić, Tutić. Numeričke vrijednosti geometrijskih konstanti elipsoida GRS 80. sažetak*, članak, .pdf.*:

»Koliko je poznato autorima ovoga rada, u Hrvatskoj se upotrebljavao i još se upotrebljava Besselov elipsoid kao osnova za geodetska računanja, ali ne postoji službeni dokument u kojem bi bile objavljene službene numeričke vrijednosti osnovnih parametara toga elipsoida. Po svemu sudeći do takve objave neće ni doći jer je nedavno službeno objavljen prijelaz na nove geodetske datume i kartografske projekcije (NN, 2004a).«

Ili "Jovo nanovo" pazit na nadnevak nastanka podataka i nadat se da je nakon pada austrougarske netko ipak poštovao neke pisane odrebe i preporuke. 1. do 1924. 2. 1924.-1951. 3. 1951.-1992. 4. 1992.-1996. (2004.)
U softveru modelirat 4 nova elipsoida za različite "epohe" (1.,2.,3.,4.,…)

»Postavlja se pitanje izbora početnih vrijednosti za poluosi. Dva bi pristupa mogla biti ispravna. Prvi je da se u skladu s Pravilnikom za državni premer (1951) kao početne vrijednosti parametara Besselova elipsoida prihvate poluosi u logaritamskom obliku (1.2). [ log a = 6.80464 34637 log b = 6.80318 92839 (1.2) ] Druga je mogućnost, prema naredbi ministra financija iz 1924. o prihvaćanju Gauss-Krügerove projekcije (Osnovni geodetski radovi, 1953), prihvaćanje numeričkih vrijednosti poluosi u metrima (1.1). [ a = 6 377 397.15500 b = 6 356 078.96325 (1.1) ] Lapaine i dr. (1992) predložili su da se točnija geodetska računanja temelje na preporuci iz priručnika Jordan/ Eggert/Kneissl (1958), prema kojoj Helmertove logaritme za poluosi log a = 6.80464 34637 00… log b = 6.80318 92839 00… (1.4) treba uzeti kao osnovne i bespogrešne veličine. Na temelju njih, svi se ostali parametri, kao što su primjerice polumjer zakrivljenosti u polu, prvi i drugi numerički ekscentricitet, prva, druga i treća spljoštenost, mogu izračunati s proizvoljnom točnošću.«
« Zadnja izmjena: 06 Veljača 2012, 19:50:53 poslijepodne glonga »

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
cs2cs
« Odgovori #20 : 07 Veljača 2012, 16:43:32 poslijepodne »
Gomila bilješki. Nastavak.

Daklem, nejasno je (ili nije dokumentirano) iz dokumentacije kako QGIS barata argumentima za transformaciju i kako ih prosljeđuje u proj4.
Ne pomaže (zasad) ni pretraga izvornog kôda proj4 (dostupan ovdje*.
Manpage proj* za cs2cs spominje argumente nu potraga za formatom i konzistencijom argumenata ne da rezultat.

cs2cs [ -eEfIlrstvwW [ args ] ] [ +opts[=arg] ]

      [ +to [+opts[=arg]] ] file[s]  

Ništa o formatu argumenata ni na manpages ni u izvornom kôdu zasad.
[+opts[=arg]

GRASS je jasniji glede argumenta dfPPM* dS i M_BF i formata argumenta kojeg očekuje kao ulazni podatak. Ctrl+F na tom linku zadnjem. "M = (1 + dS*10^-6)".
Tako je u GRASS-u, i u ILWIS-u, otud Henglu možda ideja da je njegov M_BF < 1 (0.99999445824).

Toliko za M_BF i neusklađenosti QGIS-a i njegove implementacije proj4 i cs2cs.

Ostaju i rotacije i libraryji definicije elipsoida. proj.4 ima to definitivno(?) na xxxxxxxx.155 koliko je vidit.

Austrijske rotacije (MGI) u EPSG su pozitivne. Promijenimo im svima konzekventno znak. Ostavimo M_BF većim od 1 (=2.4232) i pustimo u jednom prolazu.
Očekivani rezultat. Negativne rotacije će pomaknut točke na zapad i na jug i smanjit im visinu? Pozitivne translacije će to fino poštelat malo nazad.
Napravimo drugi prolaz s M_BF manjim od 1. M_BF = ( 1 - 2.4232 * 10^-6 ).
Očekivani rezultat. Točke bi trebale bit južnije za neki iznos. Vidi ćemo je li to značajno.

Pogledajmo nakon toga kakve probleme su bilježe Austrijanci na tom svom setu parametara.
Po (ne)vjerodostojnim austrijskim izvorima geocache-ri imali su problema s predznakom rotacija.
I ne samo oni. Tehničke specifikacije na službenom izvoru su (ovih dana) nedostupne.

Viduka. Lapaine. Broj decimala. Preciznost. Točnost.
Negdje nekad na webu je bio jedan stari .pdf dokument koji je točno opisivao matematički Gauss-Krügerov postupak. Triba ga nać. Bilježim.
Nadalje, u izvornom kôdu treba nać koju EPSG metodu koristi proj4 za izračun x i y koordinata za Transverse Mercator.
G7-2.pdf* JHS ili USGS formule? Str. 48-50 aktualne verzije G7-2 dokumenta.
Vjerujem da proj4 (američki je) koristi USGS formule. I onda opet dolazimo do nečeg što spominje Viduka. Do kojeg reda su razvijani polinomi za fi1? O tomu ovisi preciznost preračunavanja x i y u geocentrične X i Y. U izvornom kôdu i dokumentaciji proj4 treba pokušat to nać.
Nadalje. Vrijedi oprobati i preciznije definiranje parametara o kojem govori jedan ne osobito vjerodostojan izvor, ali u nedostatku i nemogućnosti pristupa tehničkoj specifikacji možda vrijedi pogledat i oprobat ovo* i zabilježit rezultate.
MGI
"AT_MGI_to_ETRS89",
"AT_MGI_to_ETRS89",
"Bessel 1841",
577.32600000000002000,
90.12900000000000500,
463.91899999999998000,
0.00002490487879860,
0.00000714615365955,
0.00002568058068837,
0.00000242320000000
Tri minusa i radijani(?) na rotacijama, ppm (part per million) na mjerilu. Prepravit ono za što je proj4 jasno napisao ulazni format (lučne sekunde za rotacije)? Više prolaza. Global Mapper vjerojatno očekuje takve ulazne podatke, no to je legacy i propietary, kad je ovako otvoreni kôd dokumentiran a neprohodan moš mislit na što bi ličio reverse engineering u takvim programima. I što ovdje* zapravo podrazumijeva pojam "scale correction" dS, M_BF…?

Te stranice su jasnije, nu od drveća se ne vidi uvik šuma:

»******************* Attention *******************
In some transformation applications a different formula
than the formula above could be in use. The signs and order
of the rotation parameters are then define differently and
rotation unit can be arcseconds. For test purposes you can
find in this information system verification data to check
your application to get correct results.
 
Distinguish the following cases:
a) change of signs of rotation (Coordinate Frame Rotation)
   Rx = -Rx
   Ry = -Ry
   Rz = -Rz

b) change of signs and order of rotation
   Rx = -Rz
   Ry = -Ry
   Rz = -Rx
******************* Attention *******************«

I kad pogledaš Z i HRVS71, imati u vidu s iste stranice:

»The three-dimensional coordinates of (X,Y,Z) of AT_MGI
were derived under using ellipsoidal heights, which are
computed from levelling heights related to Molo Sartorio
(Trieste) and a Geoid related to AT_MGI Datum Hermannskogel
and Josefstadt in Bohemia

REMARK: for transformation the (ellip) coord. must refer
all to Prime Meridian Greenwich «
« Zadnja izmjena: 07 Veljača 2012, 16:58:22 poslijepodne glonga »

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
egm2008
« Odgovori #21 : 07 Veljača 2012, 23:36:05 poslijepodne »

Undulacije geoida u raličitim razlučivostima su odavno dostupne. EGM2008*
I online za 3 modela istodobno - GeoidEval*

Nu treba li u proj4, gdal i sl. zarotirati parametre rotacije osim mijenjanja samog predznaka?
Valda su ovi QGIS i mnogi slični kopirali GRASS, a tamo opet ispliva još nešto ponovno:
»M = (1+dS*10^-6), where dS is the scale correction expressed in parts per million«*
M_BM = M? M_BM=ds?
Na "službenim" (crs.bkg.bund.de) stranicama doslovno piše:
»Operation parameter name | correction of scale
Operation parameter value | +2.42 ppm«
U BEV dokumentu projektionen_transf.pdf stoji "Scale difference" "2,4232".
Znači li "Scale difference" dS?

Aha, triba (i) li zarotirat rotacije?
Na crs.bkg.bund.de stoji:

7 Parameter Helmert Transformation,
 |X|   |X|    |Tx|   |  0  -Rz   Ry |   |X|       |X|
 |Y| = |Y| +  |Ty| + | Rz   0   -Rx | * |Y| + D * |Y|
 |Z|   |Z|    |Tz|   |-Ry   Rx    0 |   |Z|       |Z|
    T     S                                S         S

T          ... Target Datum
S          ... Source Datum
Tx, Ty, Tz ... geocentric X/Y/Z translations  [m]
Rx Ry, Rz  ... rotations around X/Y/Z axis    [radian]
D          ... correction of scale            [ppm]


GRASS devel stranica iz 2007* za Bursa-Wolf (Position Vector Transformation) navodi:

(Xt)              (  1       -Rz    +Ry)     (Xs)     (dX)
(Yt)  =      M *  ( +Rz      1      -Rx)  *  (Ys)  +  (dY)
(Zt)              ( -Ry     +Rx      1 )     (Zs)     (dZ)


a za Helmerta (Coordinate Frame Rotation)

(X')              (  1      +Rz      -Ry)     (X)      (dX)
(Y')  =      M *  ( -Rz       1      +Rx)  *  (Y)  +   (dY)
(Z')              ( +Ry    -Rx        1 )     (Z)      (dZ)


Indeksi u zapisima ne mijenjaju raspored (order), mijenjaju samo (pred)znakove u svakom od tih zapisa. U prvom zapisu su samo izvukli onaj jedan vektor pa su u 3x3 vektoru na dijagonali nule.
Izvukli su i D parametar. Izvukli su (D+1)=M? Odo ponovit matrični račun.

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
SRTM 3"
« Odgovori #22 : 08 Veljača 2012, 15:57:05 poslijepodne »

SRTM/DEM podaci se često znaju zagubit. Google će vas odvest na web stranice na kojima ćete se gnjavit s registracijom, logiranjem…
Izravan pristup podacima razlučivosti 3" (1") je ovdje* na USGS web mjestu.
Npr. N43E17 je donji lijevi kut za podatke na našem području. Datoteke pokrivaju po jedan stupanj dužine i širine, tako da narečena pokriva 43-44º širine i 17-18º dužine.

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
Krüger
« Odgovori #23 : 08 Veljača 2012, 18:34:10 poslijepodne »

Potraga za Krügerovim formulama postupcima i numeričkim primjerima online i po starim diskovima.

Wikipedijin članak* i njegovi vanjski linkovi.
Henrik Seidel: Die Mathematik der Gauss-Krüger-Abbildung.*
(razvijeno do e10)
Koordinator*. Testiranje koordinatora pokazuje robustnost i uvjetnu irelevanciju GK zona. Unos (6420000,4818000) u zadnji redak "Gauss-Krüger (DHDN/Bessel)" pokazat će vam (neočekivano) na karti (približno?) točnu lokaciju ponad Lovreća, desni prozori će (očekivano, centralni meridijan 18º) pokazivat gluposti.

Na starim diskovima ima svašta.
Ivan Lončarić. Prijelaz iz Gauss-Krügerovih u UTM koordinate putem Helmertove transformacije za područje grada Sarajeva hrčak. sažetak i prisup .pdf članku:*
Hiroshi Masaharu. Re-evaluation of the first formula given in Krueger (1912) for the Gauss-Krueger projection. .pdf*
google it (naslov Hiroshijevog rada), ima zanimljivih rezultata, npr.
Kazushige Kawase. A General Formula for Calculating Meridian Arc Length and its Application to Coordinate Conversion in the Gauss-Krüger Projection. .pdf*
( u formuli (2) išli su do C6 tj. do e10)

Konačno, digitalizirani izvornik (google "krueger2.pdf").
Johann Heinrich Louis Krüger. Konforme Abbildung des Erdellipsoids in der Ebene. Veroffentlichung Koniglich Preuszischen geodatischen Institutes. 1912. Neue Folge Nr. 52. 172p. .pdf*

Ako se opet zagubi u bespućima mrežne zbiljnosti tu je DOI*: 10.2312/GFZ.b103-krueger28

Karney unatoč bombaznu naslovu ne ide dalje od osmog stupnja. Ili ovo triba pozornije proštit?

Charles F. F. Karney. Transverse Mercator with an accuracy of a few nanometers. .pdf* DOI: 10.1007/s00190-011-0445-3
« Zadnja izmjena: 08 Veljača 2012, 19:01:05 poslijepodne glonga »

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
minus i plus
« Odgovori #24 : 08 Veljača 2012, 22:17:01 poslijepodne »

Odgovor na minus i plus glede rotacija je cilo vrime bijo prid nosom u proj4 dokumentaciji*.

x_out = M_BF*(       x[io] - Rz_BF*y[io] + Ry_BF*z[io]) + Dx_BF;
y_out = M_BF*( Rz_BF*x[io] +       y[io] - Rx_BF*z[io]) + Dy_BF;
z_out = M_BF*(-Ry_BF*x[io] + Rx_BF*y[io] +       z[io]) + Dz_BF;


Međutizin za M_BF i je li 1+dS… ili 1-dS…, zasad je oprobat modifikacije EPSG:31276 u proj4 i usporedit rezultate za slijedeća 3 proj4 stringa. Kakve su tu konvencije?

mgi6cft
+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

mgi6cp
+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,0.9999975768 +units=m +no_defs

mgi6cm
+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,1.0000024232 +units=m +no_defs

« Zadnja izmjena: 08 Veljača 2012, 22:22:08 poslijepodne glonga »

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
qgis bug #4583
« Odgovori #25 : 09 Veljača 2012, 01:34:17 prijepodne »

QGIS 1.7.3. Wroclaw na Mac-u i na Windowsima nije davao suvislu poruku o grešci, osim dojma bugovitosti u Mac slučaju. Tek je pod FreeBSD-om kvrknuo nešto o narečenom bugu
Error: non-convergent inverse meridinal dist
Svaki ti će java i na te.
Nu: http://hub.qgis.org/issues/4583

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
gdje je nestao prozor?
« Odgovori #26 : 13 Veljača 2012, 06:42:46 prijepodne »

http://hub.qgis.org/issues/3743
QGIS 1.7.3 Wroclaw je prilično neuporabiv na Mac-u. U ovom bagu je teško otkrit di je prozor za unos uopće, a prozor s canvasom se skloni kad ti treba i iskoči kad ti ne treba. Lakše ručno utipkavat točke nego pepat prozor gori-doli.
R za GRASS se pokazao problematičnim u studenta. Paket rgdal traži prethodnu instalaciju paketa sp. Preporučena je instalacija binarnog paketa sp. Nije uspjela s romanom poruka o greškama. Instalacija iz source-code-a nije uspjela jer ta nema svoju staru verziju gfortrana iz svog binarnog instalera* i to tamo gdje njoj paše.
$ which gfortran
/usr/local/bin/gfortran
$ gfortran --version
GNU Fortran (GCC) 4.6.1
Copyright (C) 2011 Free Software Foundation, Inc.

GNU Fortran comes with NO WARRANTY, to the extent permitted by law.
You may redistribute copies of GNU Fortran
under the terms of the GNU General Public License.
For more information about these matters, see the file named COPYING
$ gfortran -v
Using built-in specs.
COLLECT_GCC=gfortran
COLLECT_LTO_WRAPPER=/usr/local/gfortran/libexec/gcc/x86_64-apple-darwin10/4.6.1/lto-wrapper
Target: x86_64-apple-darwin10
Configured with: ../gcc-4.6.1/configure --prefix=/usr/local/gfortran --with-gmp=/Users/fx/devel/gcc/deps-static/x86_64 --enable-languages=c,c++,fortran,objc,obj-c++ --build=x86_64-apple-darwin10
Thread model: posix
gcc version 4.6.1 (GCC)

Posix ti je mater. Bez sačuvanog source koda od gfortrana koji je instaliran. 4.6.1. — odavde*. To treba maknut jer ovaj očekuje verziju 4.2.3. odavde* iz binarne instalacije. Ovi lažu:
»Note: more recent Fortran binaries for Xcode are available from*«
Tamo nema novih bajnariza, a da ih i ima rgdal će odbit suradnju s njima jer na njihovoj adresi /usr/local/bin očekuje gfortran ver. 4.2.3. konfiguriran po ko zna čemu što je zamotano u binarnoj instalaciji.
R paketima ne moš stat u kraj s depedencijama, ne moš pratit stalno što je 32-bitno, što 64-bitno i koju verziju depedencije će tražit i oće li ju nać, ko 32/64-bitna java i geotrans pa na entu.
Aha, zapravo R (rgdal) traži (i) 32-bitni gfortran, a onaj gore je (samo) 64-bitni (x86_64-apple-darwin10) piše u prozoru jasno samo što zbunjuje neuobičajen odgovor na upit o verziji.
Misliš da je gcc verzija 4.6.1. (gnu c compiler, što i jest slučaj kod tebe i što zbunjuje).
Kod tebe je (aha) verzija gfortrana točno iz ovog upisa*
Makni ga ko što tamo piše i probaj ovaj koji traži R (rgdal). Ako zapne, i ako bude vrimena, proba ću reproducirat grešku u virtualnom stroju.
« Zadnja izmjena: 13 Veljača 2012, 06:49:38 prijepodne glonga »

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
vraćaj tu točku
« Odgovori #27 : 13 Veljača 2012, 19:48:31 poslijepodne »

Bilješke o testungu između badiljanja. Greške su moguće. Prsti nisu baš sinkronizirani s ostalim čimbenicima (Ctrl+C i Ctrl+V greške).



Hengl* uzima obostrano određenu točku T49 i provjerava rezultate transformacije. Dok se ne mogne na točke od sniga ista točka je propuštena kroz qgis-ovu računaljku korisnički definiranog CRS-a koji nudi testnu transformaciju WGS84 ° ' " koordinata u korisničke metarske CRS-ove.

T49 WGS84: 45°33'46.3998" N, 18° 41' 47.80491" E (45.562888833333300, 45.56288883)
T49 GK6: (x) 5046738.0300 N, (y) 6554781.0700 E

Mijenjani su predznaci rotacija rx, ry, rz, te njihov poredak rz, ry, rx u sedamparametarskim transformacijama te se uvrštavala vrijednost dS, dM, u različitim kombinacijama predznaka i vrijednosti, ±dS, ±dM = ±(1±dS+10^6).

Halabuka oko preciznosti 7-parametarske transformacije bez datoteke s podacima o distorziji u mreži točaka nema nikakve osnove niti 7-parametarske transformacije vrijede truda za mnoge korisnike GPS-a. Pokazano je rezultatima troparametraskih transformacija.

Hengl:
+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
Greška koju je Hengl dobio na svojoj stranici uz svoje preporučene parametre je u metrima u odnosu na poznate y,x koordinate:
∆E = -0.5300, ∆N = -1.1700, ∆EN = 1.284445405
QGIS uz njegove parametre daje slijedeće rezultate:
∆E = -0.5376, ∆N = -1.1515, ∆EN = 1.270813129

EPSG:31276
+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
Za ovaj proj4 string nam je jasno iz očitanja GPS-a da znakovi ne mogu biti takvi (svi +).
∆E = 305.8957, ∆N = 3.4487, ∆EN = 305.9151399
To je tih tristotinjak metara koje se vuku od Austrougarske preko Austrije, ostalih stotinjak metara je od Jugoslavije i od WGS84, to će se pojavit kasnije.

EPSG:31276 "konzekventni minusi" na rotacijama daju još veću grešku.
+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
∆E = 310.8455, ∆N = -11.9781, ∆EN = 311.0761961

EPSG:31276, -,-,-,dM = (1-ds*10^6)
+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,0.9999975768 +units=m +no_defs
∆E = 310.8452, ∆N = -11.9478, ∆EN = 311.074731
3 centimetra razlike po sjeveru, a toliko priče o skali dM

EPSG:31276, -,-,-,dM = (1+ds*10^6)
+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,1.0000024232 +units=m +no_defs
∆E = 310.8452, ∆N = -11.9478, ∆EN = 311.074731
Isti rezultat (do četvrte decimale) u prethodna dva slučaja i neznatno različit u prethodna tri pokazuju da velike promjene na dM uzrokuju neznatne promjene na ∆N što ne ne bi činilo na prvi pogled u transformacijske jednadžbe. Ili proj4 ima ugrađen mehanizam prepoznavanja "prave" vrijednosti dM iz ostalih podataka (elipsoid, projekcija)? Ovo treba provjerit "pješke" uz dM = 1 i puno decimala (double precision).

EPSG:31276, -,-,-,-dM
+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
∆E = 310.8446, ∆N = -11.8748, ∆EN = 311.0713362

GPS (ana)logika navodi na +rx, +ry, -rz.

EPSG:31276, +,+,-,-dM
+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
∆E = 76.2010, ∆N = 1.5542, ∆EN = 76.21684813
Tu je negdje još onih stotinjak metara iz VGI i WGS84? 3.39" E
Krüger je jednadžbe objavio 1912. Ovamo ovi zovu sustav Hermanskogel 1901. Zagužvalo je 1914. Kad je Ferro napušten u korist Greenwicha?

EPSG:31276, -,-,+,dM
+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
∆E = 540.5392, ∆N = -9.9803, ∆EN = 540.6313283

EPSG:31276, +,-,+,dM
+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
∆E = 325.9038, ∆N = 89.8106, ∆EN = 338.0521124

EPSG:31276, +rz,+ry,-rx,-dM
+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.297,1.474,-5.137,2.4232 +units=m +no_defs
∆E = 76.3275, ∆N = 3.1384, ∆EN = 76.39199442

EPSG bi se još dao kombinirat, međutim nala razlika u apsolutnim vrijednostima zadnja 4 parametra mi govori da to ne može bitno primaknut.

Varijacije službenog GK6 (HDKS).
GK6, +,+,-,(1-ds*10^6)
+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
∆E = -0.4737, ∆N = -0.9911, ∆EN = 1.09848573

GK6, -,-,+,(1-ds*10^6)
+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
∆E = 773.5522, ∆N = 10.0911, ∆EN = 773.6180171

GK6, +,+,-,-|dM|
+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,-4.45886242 +units=m
∆E = -0.4747, ∆N = -0.8748, ∆EN = 0.995296503

GK6, +,+,-,|dM|
+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,4.45886242 +units=m
∆E = -0.4730, ∆N = -1.0648, ∆EN = 1.165130052

GK6, +,+,-,(1+ds*10^6)
+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,1.00000445886242 +units=m
∆E = -0.4737, ∆N = -0.9911, ∆EN = 1.09848573

GK6, +rz,+ry,-rx,-|dM|
+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,11.51062442,2.07314165,5.90560751,-4.45886242 +units=m
∆E = 3.9557, ∆N = 54.6206, ∆EN = 54.76365133

Županijske neslužbene.*

SDŽ, ++-,-|dM|
+proj=tmerc +pm=greenwich +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bessel +towgs84=561.8621,167.6479,460.5585,5.899175,1.460372,-11.744642,5.332886 +units=m
∆E = -0.9628, ∆N = -0.7513, ∆EN = 1.221243436

OBŽ, ++-,-|dM|
+proj=tmerc +pm=greenwich +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bessel +towgs84=519.3653,169.9679,454.9108,6.675248,2.541343,-11.300560,0.180196 +units=m
∆E = 0.7246, ∆N = -0.7513, ∆EN = 0.888681838. Očekivano je i bolje od ovog.


3 parametarske

GK6 Namibia (macgpspro)
+proj=tmerc +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bess_nam +towgs84=653,-212,449 +units=m +no_defs
∆E = -2.8490, ∆N = -68.2902, ∆EN = 68.3496029

NIMA/NGA Namibia
+proj=tmerc +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bess_nam +towgs84=682,-203,480 +units=m +no_defs
∆E = -3.6230, ∆N = -68.2669, ∆EN = 68.36297071

NIMA/NGA
+proj=tmerc +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bessel +towgs84=682,-203,480 +units=m +no_defs
∆E = -2.8807, ∆N = 0.6406, ∆EN = 2.951067747

"Fugawi" rec.
+proj=tmerc +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bessel +towgs84=634.12,-216.43,436.10 +units=m +no_defs
∆E = -0.2907, ∆N = 5.3846, ∆EN = 5.392441345

"Fugawi" rec. bess_nam
∆E = -1.0330, ∆N = -63.5229, ∆EN = 63.53129869
VGI-jevci su mogli (na)pravit onih sedamdesetak metara razlike na zamjeni x-y i E-N, isto ko i na zamjeni bessel i bess_nam, ili skupa.

"Garmin" "Ozi"
+proj=tmerc +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bessel +towgs84=653,-212,449 +units=m +no_defs
∆E = -2.1066, ∆N = 0.6173, ∆EN = 2.195181735

Biram "Garmin/Ozi". Ako je već u pitanju greška reda veličine metar (ili dva), biram lakšu metodu, koju je lakše rekonstruirat, pogotovo uz masovnost ručnih GPS uređaja koji su upravo na tim parametrima.

Ručni GPS novijeg datuma koji dopušta definiciju korisničkog datuma sa 7 parametara u slučaju EPSG:31276 pravi nešto čudno, i u MGI 6 (° ' ") i u WGS84 (° ' ") na istom mjestu daje ista očitanja.
To nas svrće na nepodržani MGI (deg) EPSG:63126405:
+proj=longlat +ellps=bessel +towgs84=426.9,142.6,460.1,4.91,4.49,-12.42,3.5271281868253483 +no_defs
mod
+proj=tmerc +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bessel +towgs84=426.9,142.6,460.1,4.91,4.49,-12.42,3.5271281868253483 +units=m +no_defs
∆E = 3.1925, ∆N = -2.5763, ∆EN = 4.102362483

I kakva je razlika EPSG:31276
+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
i EPSG:31266
+proj=tmerc +lat_0=0 +lon_0=18 +k=1 +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

k=1? Ne š ti.

Vidjeti cs2cs i ovakvo nešto za tmerc etc.:
cs2cs +proj=latlong +datum=WGS84 \
+to +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 \
+x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +no_defs

EPSG&proj4?
»How do I use EPSG coordinate system codes with PROJ.4?

There are known to be problems with some coordinate systems, and any coordinate systems with odd axes, a non-greenwich prime meridian or other quirkyness are unlikely to work properly. Caveat Emptor!«*

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
a pedja?
« Odgovori #28 : 13 Veljača 2012, 20:01:21 poslijepodne »
Neslužbeni VGI. Prema Pedji.*

VGI
+proj=tmerc +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bessel +towgs84=653,-212,431 +units=m +no_defs
∆E = -1.9964, ∆N = -11.9836, ∆NE = 12.14875639

VGI bess_nam
+proj=tmerc +lat_0=0 +lon_0=18 +k=0.9999 +x_0=6500000 +y_0=0 +ellps=bess_nam +towgs84=653,-212,431 +units=m +no_defs
∆E = -2.7387, ∆N = -80.8913, ∆NE = 80.93764818

Osamdeset. Ima nešto s tom Namibijom i zavrćanjem x-N, y-E.

Offline glonga

  • Forumska Al-tipkovnica
  • Gožđeno pero foruma
  • ***
  • Postova: 10.154
Ferro
« Odgovori #29 : 13 Veljača 2012, 23:03:30 poslijepodne »
»It was thought to be exactly 20 degrees west of the Paris meridian, so indeed the exact position of Ferro was never considered. Old maps (outside of Anglo-America) often have a common grid with Paris degrees at the top and Ferro degrees offset by 20 at the bottom. Louis Feuillée also worked on this problem in 1724.

It was later found that the actual island of El Hierro itself is in fact 20° 23' 9" west of Paris, but the Ferro meridian was still defined as 20 degrees west of Paris.

According to the European longitude adjustment of Carl Theodor Albrecht (ca. 1890) the Ferro meridian is 17° 39' 46.02" west of the Greenwich meridian. But for the geodetic networks of Austria, Germany and Czechoslovakia, the value 17° 40' 00" was adopted in the 1920s, not only for practical reasons but also as it was discovered that the longitude of the Berlin (Rauenberg) fundamental point was miscalculated by 13.39". For the geodetic networks of Hungary and Yugoslavia, the value of Albrecht was used prior to the switch to the Greenwich prime meridian.«*

Prati Mađarsku, ne Austriju.
∆E = 13.98"