14. prosinca 2012.

Kreiranje NTv2 grid datoteke za transformaciju između HDKS i HTRS96 sustava


Kao što je spomenuto u prethodnom postu NTv2 format grid datoteke koristi se u velikom broju državama kao što su Kanada, Australija, Novi Zeland, Velika Britanija, Švicarska, Njemačka, Brazil, Francuska itd. NTv2 format je razvijen u Kanadi, zatim je prihvaćen i prilagođen u Australiji, a nakon toga u ostalim državama.


Grid polja različitih intervala

NTv2 format grid datoteke može sadržavati različita grid polja ( sub-grid ), tj gridove s različitim razmacima (intervalima) između točaka grida. Također se može po potrebi nadopunjavati s novim grid poljima. NTv2 datoteka može biti u ASCII i u binarnom formatu.
Treba napomenuti da NTv2 datoteka služi za 2D transformaciju, tj za transformaciju geodetske dužine i geodetske širine iz jednog sustava u drugi. NTv2 format grid datoteke je možda najbolje dokumentiran u Australiji gdje su razvijeni programi GDAit  i GDAy  koji služe za transformaciju između starog Australian Geodetic Datum (AGD) i novog Geocentric Datum of Australia (GDA) australskog datuma, a koriste NTv2 format grid datoteke za transformaciju. Nešto više o NTv2 formatu možete pronaći u GDAitUsers Guide  i GDAyUsers Guide . NTv2 datoteka se sastoji od zaglavlja i redova u kojima se nalaze vrijednosti korekcija za geodetsku dužinu i geodetsku širinu u svim točkama grida. U točkama grida tj. u korekcijama je modelirana distorzija.

Zaglavlje NTv2 datoteke:

NUM_OREC 11
NUM_SREC 11
NUM_FILE  1
GS_TYPE SECONDS
VERSION NTv2.0 
DATUM_F HDKS 
DATUM_T HTRS96
MAJOR_F  6377397.155
MINOR_F  6356078.963
MAJOR_T  6378137.0
MINOR_T  6356752.3141
SUB_NAME HRGRID
PARENT  NONE   
CREATED 12-10-15
UPDATED 12-10-15
S_LAT     151200.000000
N_LAT     167760.000000
E_LONG    -70200.000000
W_LONG    -46800.000000
LAT_INC       60.000000
LONG_INC      90.000000
GS_COUNT 72297

Dio NTv2 datoteke:

  0.801326 18.827497  0.000000  0.000000
  0.798218 18.816583  0.000000  0.000000
  0.795113 18.805667  0.000000  0.000000
  0.792010 18.794749  0.000000  0.000000
  0.788909 18.783831  0.000000  0.000000

gdje je u prvom redu korekcija geodetske širine, zatim korekcija geodetske dužine i točnost geodetske širine i dužine za prvu točku grida. Podaci o području obuhvata grida, intervalu unutar točaka grida, kao i o korištenim elipsoidima nalaze se u zaglavlju datoteke.



Grafički prikaz grid polja

Na slici je grafički prikaz grid polja gdje je  n = broj kolona   i   m = broj redova
Za dohvat podataka grida za traženu točku P koriste se sljedeće formule:

n = int(W-E / dLon)+1   
m = int(N-S / dLat)+1
i = int(LatP - S / dLat) + 1
j = int(LonP - E / dLon) + 1
A = n * (i - 1) + j
B = A + 1
C = A + n
D = C + 1

Nakon dohvata podataka koristi se metoda bilinearne interpolacije za računanje parametara transformacije tražene točke.

Za potrebe transformacije iz HDKS sustava u HTRS96 sustav i obratno kreirao sam NTv2 grid datoteku za područje Hrvatske. Vrijednosti  korekcija tj. parametara transformacije u točkama grida su izračunate programom T7D za područje od 42° do 46,6°  u smjeru jug - sjever s korakom od 60" i od 19,5° do  do 13 u smjeru istok - zapad s korakom od 90".


Područje obuhvata kreirane NTv2 datoteke

Datoteka s vrijednostima geodetske širine i geodetske dužine u HDKS sustavu za sve točke grida je kreirana programski (VB6). Nakon toga su izračunate vrijednosti tih koordinata u sustavu HTRS96 programom T7D u kojem je ugrađen model distorzije. Tako je dobiven drugi set koordinata ( geodetska širina i dužina ) za sve točke grida u HTRS96 sustavu. Razlike tih koordinata su ugrađene u NTv2 datoteku nakon kreiranog zaglavlja također programski. Time je kreirana NTV2 datoteka u ASCII formatu.
Linkovi za skidanje NTv2 datoteka u ASCII i binarnom formatu:

HRNTv2.gsb – binarna NTv2 datoteka - http://www.mediafire.com/?76kkkqpmhy2d75s
HRNTv2 .asc – ASCII NTv2 datoteka  - http://www.mediafire.com/?m7p4foehw4d6qcv