Ispanijos geodezinė reforma: Real Decreto 1071/2007 ir perėjimas iš ED50 į ETRS89
Daugiau nei pusę amžiaus visa Ispanijos valstybinė kartografija, kadastro registrai (Dirección General del Catastro) ir civilinės inžinerijos projektai rėmėsi European Datum 1950 (ED50) koordinačių sistema. Sukurta po Antrojo pasaulinio karo Jungtinių Valstijų kariuomenės žemėlapių tarnybos (AMS), siekiant suderinti Vakarų Europos geodezinius tinklus, ED50 naudojo istorinį International 1924 (Hayford) elipsoidą su pradiniu atskaitos tašku Potsdamo observatorijoje Vokietijoje. Nors ši sistema dešimtmečius sėkmingai tarnavo nacionaliniams topografiniams žemėlapiams (MTN50 ir MTN25), palydovinės navigacijos (GNSS, GPS) ir šiuolaikinės kosminės geodezijos atsiradimas atskleidė esminius regioninius netolygumus bei poslinkius.
Siekdama užtikrinti visišką suderinamumą su Europos Sąjungos INSPIRE direktyva ir moderniomis palydovinėmis sistemomis, Ispanijos Karalystės Vyriausybė 2007 m. liepos 27 d. priėmė Real Decreto 1071/2007. Šis teisės aktas nustatė privalomą ir neatšaukiamą perėjimą nuo ED50 prie ETRS89 (European Terrestrial Reference System 1989) žemyninėje Ispanijos dalyje bei Balearų salose, o Kanarų salose – prie REGCAN95 sistemos (visiškai suderinamos su WGS84 ir ITRF93). Pasibaigus pereinamajam laikotarpiui 2015 m. sausio 1 d., visi teisiniai kadastro planai, statybos projektai ir GIS sluoksniai Ispanijoje privalo būti pateikiami išskirtinai ETRS89 sistemoje.
Elipsoidų geometriniai skirtumai ir fizinis poslinkio vektorius
Matematinis poreikis perskaičiuoti koordinates kyla iš pamatinių atskaitos elipsoidų geometrijos ir jų padėties Žemės atžvilgiu skirtumų:
- ED50 (Hayford 1924 / International 1924): Didysis pusašis a = 6,378,388.0 m, suplokštėjimo atvirkštinis dydis 1/f = 297.0, pirmasis ekscentriciteto kvadratas e2 = 0.00672267002233. Šis elipsoidas buvo orientuotas taip, kad geriausiai atitiktų vietinį Vakarų Europos geoidą, o ne Žemės masės centrą.
- ETRS89 / GRS80 (Suderinamas su WGS84): Didysis pusašis a = 6,378,137.0 m, suplokštėjimo atvirkštinis dydis 1/f = 298.257222101, pirmasis ekscentriciteto kvadratas e2 = 0.00669438002290. GRS80 yra geocentrinis elipsoidas, kurio centras sutampa su tikruoju Žemės masės centru.
Kadangi Hayford 1924 ir GRS80 elipsoidų centrai 3D erdvėje yra pasislinkę vienas kito atžvilgiu per šimtus metrų, atidėjus senąsias ED50 koordinates tiesiai ant modernaus ETRS89 / GPS žemėlapio, Ispanijos teritorijoje susidaro 130–210 metrų fizinis geodezinis poslinkis. Madride ir centrinėje Ispanijoje šis poslinkis siekia ~166.7 metro pietvakarių kryptimi (azimutas ~217.9° SW), Barselonoje – ~154.8 m, o Sevilijoje – ~160.2 m. Kadastro ir žemėvaldos srityje toks neperskaičiuotas poslinkis reikštų, kad sklypo riba persikelia į kaimyno kiemą arba pastatą.
Oficiali IGN Ispanijos 7 parametrų Position Vector Helmert transformacija
Tiksliam trimačiui perskaičiavimui tarp ED50 ir ETRS89 be vietinių iškraipymų Ispanijos nacionalinis geografijos institutas (Instituto Geográfico Nacional – IGN) patvirtino oficialius 7 parametrų Helmert (Position Vector R2007) transformacijos parametrus žemyninei Ispanijai ir Balearams:
• Poslinkiai (Translacija): ΔX = -129.549 m, ΔY = -82.688 m, ΔZ = -163.641 m
• Pasukimo kampai (Rotacija): Rx = -0.426″, Ry = -0.355″, Rz = +0.705″
• Mastelio pokytis: Δs = +1.036 ppm (+1.036 × 10-6)
Matematinis perskaičiavimo algoritmas atliekamas keturiais etapais:
- Atvirkštinė UTM projekcija: Pradinės plokštuminės UTM koordinatės (X, Y) perskaičiuojamos į elipsoidines geodezines koordinates (φ, λ) naudojant Redfearn / Krüger formules.
- Geodezinių koordinačių vertimas į 3D Dekarto koordinates: (φ, λ, h) paverčiamos į (X, Y, Z) per elipsoido pirmojo vertikalo kreivumo spindulį N(φ) = a / √(1 − e2 sin2φ).
- Helmert 3D erdvinė matrica (Position Vector): Koordinačių vektorius padauginamas iš posūkio matricos, mastelio faktoriaus ir sudedamas su linijiniais poslinkiais:
XETRS89 = ΔX + (1 + Δs) · ( XED50 − Rz·YED50 + Ry·ZED50 )
YETRS89 = ΔY + (1 + Δs) · ( Rz·XED50 + YED50 − Rx·ZED50 )
ZETRS89 = ΔZ + (1 + Δs) · ( −Ry·XED50 + Rx·YED50 + ZED50 ) - Dekarto koordinačių vertimas atgal į elipsoidines ir tiesioginė UTM projekcija: Naujas (X, Y, Z) vektorius Bowring algoritmu paverčiamas į naująją platumą ir ilgumą GRS80 elipsoide, o vėliau suprojektuojamas į tikslinę UTM zoną (k0 = 0.9996, X0 = 500,000 m).
Keturios Ispanijos UTM zonos (Huso 28, 29, 30, 31) ir regionai
Ispanijos Karalystės teritorija driekiasi per keturias universalias skersines Merkatoriaus (UTM) zonas:
- UTM Zona 28 (Huso 28, Centrinis meridianas λ0 = -15° W): Apima Kanarų salas (Tenerifę, Gran Kanariją, Lanzarotę, Fuerteventūrą, La Palmą, La Gomerą, El Hierro). Oficiali sistema – REGCAN95 (EPSG:4083), visiškai suderinama su WGS84/ETRS89.
- UTM Zona 29 (Huso 29, Centrinis meridianas λ0 = -9° W): Apima Galisiją (A Koruniją, Lugą, Ourensę, Pontevedrą), Vakarų Astūriją, Vakarų Kastiliją ir Leoną bei Vakarų Estremadūrą (ETRS89: EPSG:25829 / ED50: EPSG:23029).
- UTM Zona 30 (Huso 30, Centrinis meridianas λ0 = -3° W): Pagrindinė Ispanijos ašis, apimanti Madrido bendruomenę, beveik visą Andalūziją, Kastiliją-La Mančą, Rytų Kastiliją ir Leoną, Baskų kraštą, Kantabriją, Navarą bei La Riochą (ETRS89: EPSG:25830 / ED50: EPSG:23030).
- UTM Zona 31 (Huso 31, Centrinis meridianas λ0 = +3° E): Apima Kataloniją (Barseloną, Žironą, Lleidą, Taragoną), Valensijos bendruomenę, Mursijos regioną, Rytų Aragoną bei Balearų salas (Maljorką, Menorką, Ibisą, Formenterą) (ETRS89: EPSG:25831 / ED50: EPSG:23031).
Kartografinio poslinkio metafizinė komedija: kaip ispaniška vila nukeliavo 167 metrus į kaimyno alyvmedžių giraitę
Žmogaus įsitikinime, jog žemė po kojomis yra absoliuti ir nepajudinama, slypi gili biurokratinė tragikomedija. Ištisus dešimtmečius Ispanijos žemvaldžiai, notarai ir orūs matininkai ramiai miegojo, šventai tikėdami, kad jų protėvių finca ar saulėtas andalūziškas kiemas amžiams įkalti į granitinę Žemės tiesą. Realybėje visa Pirėnų pusiasalio žemė plūduriavo ant Johno Fillmore'o Hayfordo – 1909 metų amerikiečių geodezijos genijaus – sukurto elipsoido, kurio centras dėl to meto technologijų ribotumo prasilenkė su tikruoju Žemės masės centru per šimtus metrų. Kai Real Decreto 1071/2007 teisiškai panaikino ED50, akmeniniai riboženkliai (mojones) nepajudėjo nė milimetro, tačiau skaitmeniniuose Kadastro registruose kiekviena sodyba, baseinas ir gotikinė varpinė akimirksniu teleportavosi maždaug 167 metrus į pietvakarius – tiesiai į gretimo sklypo pomidorų lysves.
Ši vaiduokliška migracija yra blaivus memento mori kiekvienam, painiojančiam skaitmeninę koordinačių eilutę su fizine tikrove. Žemei visiškai nerūpi mūsų elipsoidinės abstrakcijos, o šimtametis Chaeno alyvmedis neprivalo konsultuotis su Bursa-Wolf rotacijos matrica prieš leisdamas šaknis gilyn. Tačiau steriliame statybos leidimų ir nekilnojamojo turto registrų panteone neįvertintas geodezinis poslinkis reiškia, kad jūsų milijoninės vertės saulės jėgainė teisiškai atsidurs melioracijos griovyje, o autostrados viadukas pasibaigs tuštumoje. Geodezija – tai preciziškas, šiek tiek neurotiškas menas sutaikyti žmogaus geometrinę puikybę su netobulu, bulvės formos besisukančios Žemės geoidu.
Praktinis taikymas: Ispanijos kadastras (Catastro), GML failai ir IGN NTv2 tinklelis
Šiuolaikiniai topografai, architektai ir geodezijos inžinieriai nuolat susiduria su senosiomis ED50 koordinatėmis archyviniuose savivaldybių planuose, nuosavybės registruose (Registro de la Propiedad) bei kasybos koncesijose. Oficialiam teisiniam Ispanijos kadastro sklypų ribų įteisinimui bei netikslumų šalinimui (Archivos GML de subsanación de discrepancias pagal Ley 13/2015 įstatymą) ir alternatyvaus grafinio patvirtinimo ataskaitoms (Informes de Validación Gráfica Alternativa IVGA / FXCC) teikti į Sede Electrónica del Catastro, privalomas valstybinis Instituto Geográfico Nacional (IGN) standartas yra NTv2 distorsijos tinklelis (PENR2009.gsb / BALR2009.gsb).
Šis konverteris apjungia oficialią 7 parametrų Position Vector Helmert matricą (IGN R2007) su regionine konformine NTv2 tinklelio interpoliacija, garantuojančia aukštą geodezinį tikslumą (±0.1 – 0.5 m) visoje Ispanijoje. Jis palaiko 3 patogius režimus (geografinių koordinačių su DMS, plokštuminių UTM koordinačių ir masinio taškų sąrašo su formatų parinkikliu), atlieka momentinį abipusį perskaičiavimą, automatiškai nustato UTM zoną iš ilgumos, rodo tašką Google Maps ir OpenStreetMap bei eksportuoja duomenis į CSV ir TSV (tinkamus QGIS, ArcGIS ir AutoCAD).