| Повече ▼

Поща: WGS84 Lat Lon Point към Gauss-Kruger 3

Поща: WGS84 Lat Lon Point към Gauss-Kruger 3


Как да конвертирам WGS84 lat lon point в PostGIS в Guass Kruger 3 (31467)?

Някак не изглежда да съм разбрал правилно. Опитах обичайния начин за конвертиране на проекции:

st_transform (st_setsrid (st_makepoint (50, 8), 2346), 31467)

Това обаче дава грешка извън границите. Предполагам, че PostGIS очаква стойностите да бъдат декартови координати.


Може би само печатна грешка? WGS84 е 4326.


Поща: WGS84 Lat Lon Point към Gauss -Kruger 3 - Географски информационни системи

Това е част 1 от поредица от 3 части.

Ако погледнете интерактивната уеб карта на Cape Town Marathon, ще видите километрични маркери по състезателните маршрути. Създаването им не е толкова лесно, колкото си мислите. За екипа, планиращ маратона и другите маршрути, той включваше старателен ръчен процес в Google Земя.

  • Маршрутите се редактират доста често, водещи до събитието
  • Те не искат ръчно да създават точки всеки път, когато редактират маршрута
  • Те се нуждаят от разстоянията, за да бъдат възможно най -точни. В края на краищата маратонът е златно събитие на IAAF. Въпреки това дори най -добрата карта няма да има 100% точно разстояние поради грешки в картографирането и промени в надморската височина. За да установи официалната позиция на финалната линия, колоездач измерва маршрута по пътя, който бегачът би поел. Ето защо маркерите за крайно разстояние не са точно в края на маршрута.

Бях хвърлил око на решение, включващо линейни препратки от самото начало и опитах няколко подхода, преди да се спра на този, който подробно описвам по -долу. Намерих решения, използвайки v.segment от GRASS и приставката QGIS LRS, но смятах, че те са твърде сложни за специалист, който не е ГИС, и включват също няколко стъпки.

Последното ми решение беше изглед в PostGIS, който просто седи там без сървър на сървъра и генерира точките автоматично всеки път, когато някой прегледа уеб картата. Планиращите могат да редактират маршрутите, без да се притесняват за маркерите. Когато записват актуализация на маршрут, маркерите просто работят.

Първо създадох таблица със събития с идентификаторите на маршрутите и разстоянията, на които исках маркери. Това е пример:

След това, тъй като работим в местния CRS за Кейптаун, който не е стандартен в GIS софтуера, трябваше да добавя WGS19 CRS (наричан „HBK_NO19“ в QGIS) към PostGIS, като изпълнявам тази заявка, давайки й произволен id на 40019:

След това създадох изглед като този за всеки маршрут.

Подзаявката WITH подготвя таблица, в която маршрутът във всеки ред е проектиран към локалния CRS, заедно с разстоянието на маркера от таблицата на събитията и дължината на маршрута (много точно изчислено на сфероида и като се вземат предвид промените на котата, тъй като геометрията има Z стойности).

След това основната част от заявката интерполира точка по маршрута на всяко разстояние на маркера и прави известно форматиране на разстоянието, така че да може да се използва като етикет в картата. Той също така филтрира всеки маркер, който може да надхвърля дължината на маршрута (което би било безсмислено да се интерполира и би довело до неуспех на заявката).

Трябваше да проектирам маршрутите, тъй като те се съхраняват непроектирани (EPSG: 4326) и исках да минимизирам изкривяванията за линейните препратки, които работят само върху геометрии, а не върху географии. Ще забележите, че ST_Length () е направено с география, която е тип PostGIS, която принуждава функциите, приложени към нея, да работят на сфероида WGS84, а не на равнина. Дължината, изчислена на сфероида, е по -точна от дължината, изчислена във всяка равнинна проекция, и се увеличава допълнително, ако има стойности Z в координатите (ST_Length () в 3D география връща 3D дължина).

Ако имате елегантно решение на същия проблем, бих искал да знам за него. Използването на generator_series () по този начин вместо присъединяване към таблица на събития е доста елегантно.

Как така има само няколко маркера, когато се намали и след това по един на всеки километър, когато се увеличи? Направихме това с рендиране, базирано на правила в QGIS с зависимост от мащаба във всяко правило и този филтър в „правилото за мащабиране“: dist % 2 = 0, което филтрира всичко, което има остатък, разделено на две (или пет в случай на маратона).


Изтеглянето на уебсайта на PostGIS съдържа информация за инсталиране на други платформи, включително macOS и други дистрибуции на Linux

След като PostGIS е инсталиран, ще трябва да конфигурирате базата данни да използва разширенията. Ако сте инсталирали PostGIS версия & gt 2.0, това е толкова просто, колкото да издадете следната команда с psql, използвайки адресната база данни от предишното ни упражнение.

В зависимост от вашата версия, бихте могли да намерите повече инструкции как да активирате пространствено база данни на адрес https://postgis.net/docs/postgis_administration.html#create_spatial_db.


16.1.5. Гледайки инсталираните функции на PostGIS¶

PostGIS може да се мисли като съвкупност от функции в базата данни, които разширяват основните възможности на PostgreSQL, така че да може да се справя с пространствени данни. Под ‘ справяне с ’ имаме предвид съхранение, извличане, запитване и манипулиране. За да направите това, в базата данни са инсталирани редица функции.

Нашият PostgreSQL адрес базата данни вече е геопазарно активирана, благодарение на PostGIS. В следващите раздели ще се задълбочим много по -дълбоко в това, но нека „#8217“ ви даде бърз дегустатор. Нека ’s да кажем, че искаме да създадем точка от текст. Първо използваме командата psql, за да намерим функции, свързани с точка. Ако все още не сте свързани с адрес база данни, направете го сега. След това стартирайте:

Това е командата, която търсим: st_pointfromtext . За да прелиствате списъка, използвайте стрелката надолу, след което натиснете q за да излезете обратно към черупката psql.

Essayez d ’exécuter cette commandnde:

  • Дефинирахме точка в позиция 1,1 (EPSG: 4326 се приема), използвайки ТОЧКА (1 1) ,
  • Изпълнихме sql изявление, но не на която и да е таблица, само на данни, въведени от SQL подканата,
  • Полученият ред няма особен смисъл.

Полученият ред е във формат OGC, наречен ‘Well Known Binary ’ (WKB). Ще разгледаме този формат подробно в следващия раздел.

За да върнем резултатите като текст, можем да направим бързо сканиране през списъка с функции за нещо, което връща текст:

Заявката, която ние ’ търсим сега, е st_astext . Нека ’s го комбинират с предишната заявка:

Тук въведохме низа ТОЧКА (1,1) , го превърна в точка с помощта st_pointfromtext () , и го превърна обратно в читава за човека форма с st_astext () , което ни върна оригиналния низ.

Един последен пример, преди наистина да влезем в детайлите на използването на PostGIS:

Какво направи това? Той създаде буфер от 1 градус около нашата точка и върна резултата като текст.


Основни елементи на нашето решение за бази данни

В този раздел ние показваме с малък пример как сме изградили нашите бази данни върху NetCDF4 файлове, които са монтирани към локална файлова система. Следният код далеч не е готов за производство, но трябва да подчертае някои от дизайнерските решения, които направихме.

Следният код умерено експлоатира NumPy, Scipy и Pandas. Ако не сте запознати с тези цифрови инструменти, вероятно ще трябва да се консултирате с тяхната документация, за да разберете подробностите.

В този пример използваме само един файл NetCDF4 за температурата на въздуха. Ако искате да стартирате кода, трябва да включите този файл в текущата си работна директория. Този конкретен файл с данни не се използва в нашите продукти, но структурата прилича на тази на нашите действителни данни.

Започваме примера с импортиране на цифровите инструменти и клас, който може да отваря NetCDF4 файлове за четене:

Тогава нека ’s дефинират какъв ще бъде крайният изход. Искаме да създадем клас база данни, от който можем да поискаме прогнозни данни за маршрут на кораб#8217s. Бихме искали базата данни да ни предоставя редовно решетъчни данни, за да можем след това лесно да използваме интерполаторните прогнозни стойности на SciPy ’s за маршрута. Изразено в код:

По -горе изрязването на данните означава, че ние извличаме “cube ” от растерни точки, така че всички триизмерни маршрутни точки се съдържат в най -външните клетки на куба. ос_точки тогава е кортеж от три масива, които съдържат стойности за време, ширина и дължина, които съответстват на клетките вътре в куба.

Сега можем да продължим да прилагаме действителния клас база данни:

Нищо сложно там. Идеята в clip_data е, че екземплярът от базата данни първо разглежда исканите часове и след това отваря съответните NetCDF4 файлове. След това векторът за заявените времена се разделя според отворените файлове.

Сега индексът на файловете на базата данни напълно липсва, както и логиката, която решава кои файлове да се отворят. Това са малки подробности за внедряването, които избрахме да пропуснем тук. Вместо това просто отваряме един файл в текущата работна директория. В действителност различните изрязани данни също трябва да бъдат свързани заедно, но това е друга подробност, която пропускаме за простота.

Функцията, която все още не е дефинирана, _clip_one_dataset, осъществява достъп само до необходимите стойности във файловете на NetCDF4 и ги връща като първата им възвращаема стойност. Вторите връщащи стойности съдържат стойности на оси, които съответстват на върнатите стойности на точки. В този пример функцията е дефинирана, както следва:

Първо функцията дефинира по-малка функция, която намира подмножества от оси с тези съраунд заявени стойности на оси. След като се получат подмножествата, съответните прогнозни стойности се извличат с помощта на API за набор от данни. Тук има допълнителна ос променливи, височина, която също се съдържа в този пример данни.

И това е#8217 това е! Има нашата елементарна реализация на база данни.


Съдържание

Концепцията за дължина е разработена за първи път от древногръцки астрономи. Хипарх (2 век пр.н.е.) използва координатна система, която приема сферична земя и я разделя на 360 °, както правим и до днес. Първият му меридиан премина през Александрия. [3]: 31 Той също така предлага метод за определяне на географската дължина, като сравнява местното време на лунно затъмнение на две различни места, като по този начин демонстрира разбиране за връзката между дължина и време. [3]: 11. [4] Клавдий Птолемей (2 век от н.е.) разработва система за картографиране, използвайки извити паралели, които намаляват изкривяванията. Той също така събира данни за много места, от Великобритания до Близкия изток. Той използва основен меридиан през Канарските острови, така че всички стойности на географската дължина да са положителни. Докато системата на Птолемей беше здрава, използваните от него данни често бяха лоши, което водеше до груба надценка (с около 70%) от дължината на Средиземноморието. [5] [6]: 551–553 [7]

След падането на Римската империя интересът към географията силно намалява в Европа. [8]: 65 индуски и мюсюлмански астрономи продължават да развиват тези идеи, добавяйки много нови места и често подобрявайки данните на Птолемей. [9] [10] Например ал-Батани използва едновременни наблюдения на две лунни затъмнения, за да определи разликата в географската дължина между Антакия и Ракка с грешка по-малка от 1 °. Това се счита за най -доброто, което може да бъде постигнато с наличните тогава методи - наблюдение на затъмнението с невъоръжено око и определяне на местното време с помощта на астролаба за измерване на височината на подходяща „звезда с часовници“. [11] [12]

През по -късното Средновековие интересът към географията се възражда на запад, тъй като пътуванията се увеличават и арабската наука започва да се познава чрез контакти с Испания и Северна Африка. През XII век астрономическите таблици са изготвени за редица европейски градове въз основа на работата на ал-Заркали в Толедо. Лунното затъмнение от 12 септември 1178 г. е използвано за установяване на разликите в географската дължина между Толедо, Марсилия и Херефорд. [13]: 85

Христофор Колумб направи два опита да използва лунните затъмнения, за да открие своята дължина, първият на остров Саона на 14 септември 1494 г. (второ плаване), а вторият в Ямайка на 29 февруари 1504 г. (четвърто пътешествие). Предполага се, че той е използвал астрономически таблици за справка. Определянията му по географска дължина показват големи грешки съответно 13 и 38 ° W. [14] Randles (1985) документира измерване на географска дължина от португалците и испанците между 1514 и 1627 г. както в Америка, така и в Азия. Грешките варират от 2-25 °. [15]

Телескопът е изобретен в началото на 17-ти век. Първоначално устройство за наблюдение, развитието през следващия половин век го превърна в точен инструмент за измерване. [16] [17] Часовникът с махало е патентован от Christiaan Huygens през 1657 г. [18] и дава увеличение на точността с около 30 пъти спрямо предишните механични часовници. [19] Тези две изобретения биха революционизирали наблюдателната астрономия и картография. [20]

На сушата в периода от развитието на телескопи и махални часовници до средата на 18-ти век се наблюдава непрекъснато увеличаване на броя на местата, чиято дължина е определена с разумна точност, често с грешки по-малки от градуса и почти винаги в рамките на 2-3 °. Към 1720 -те грешките са били постоянно по -малко от 1 °. [21] В морето през същия период ситуацията беше много различна. Два проблема се оказаха неразрешими. Първият беше необходимостта от навигатор за незабавни резултати. Второто беше морската среда. Правенето на точни наблюдения при океанско подуване е много по -трудно, отколкото на сушата, а часовниците с махало не работят добре при тези условия.

В отговор на проблемите с корабоплаването редица европейски морски сили предложиха награди за метод за определяне на географската дължина в морето. Най-известният от тях е Законът за географската дължина, приет от британския парламент през 1714 г. [22]: 8 Той предлага две нива на награди за решения в рамките на 1 ° и 0,5 °. Награди бяха дадени за две решения: лунни разстояния, направени на практика чрез таблиците на Тобиас Майер [23], разработени в морски алманах от астронома Крал Невил Маскелин и за хронометрите, разработени от йоркширския дърводелец и часовникар Джон Харисън. Харисън построи пет хронометра за повече от три десетилетия. Тази работа е подкрепена и възнаградена с хиляди лири от борда на Longitude [24], но той се бори да получи пари до най -високата награда от £ 20,000, като накрая получава допълнително плащане през 1773 г. след намесата на парламента [22]: 26. Измина известно време преди двата метода да се използват широко в навигацията. В първите години хронометрите бяха много скъпи, а изчисленията, необходими за лунните разстояния, все още бяха сложни и отнемаха много време. Лунните разстояния влизат в обща употреба след 1790 г. [25] Хронометрите имат предимствата, че както наблюденията, така и изчисленията са по-опростени и тъй като поевтиняват в началото на 19-ти век, те започват да заменят луни, които рядко се използват след 1850 г. [26]

Първите работещи телеграфи са създадени във Великобритания от Уитстоун и Кук през 1839 г., а в САЩ от Морс през 1844 г. Бързо се осъзнава, че телеграфът може да се използва за предаване на времеви сигнал за определяне на географска дължина. [27] Методът скоро беше на практика използван за определяне на географска дължина, особено в Северна Америка, и на все по -големи разстояния с разширяването на телеграфната мрежа, включително Западна Европа с завършването на трансатлантически кабели. Проучването на брега на САЩ беше особено активно в това развитие, и не само в Съединените щати. Проучването установи вериги от картографирани местоположения през Централна и Южна Америка и Западна Индия и чак до Япония и Китай през годините 1874–90. Това допринесе значително за точното картографиране на тези области. [28] [29]

Докато морските пехотинци се възползваха от точните карти, те не можеха да приемат телеграфни сигнали, докато са в ход, и затова не можеха да използват метода за навигация. Това се промени, когато безжичната телеграфия стана достъпна в началото на 20-ти век. [30] Безжичните сигнали за времето за използване на кораби се предават от Халифакс, Нова Скотия, започвайки от 1907 г. [31] и от Айфеловата кула в Париж от 1910 г. [32] Тези сигнали позволяват на навигаторите да проверяват и настройват хронометрите си на честа основа. [33]

Радионавигационните системи влизат в обща употреба след Втората световна война. Всички системи зависят от предавания от фиксирани навигационни маяци. Приемник на борда на кораб изчислява позицията на кораба от тези предавания. [34] Те позволяват точна навигация, когато лошата видимост възпрепятства астрономическите наблюдения, и се превръщат в утвърдения метод за търговско корабоплаване, докато не бъдат заменени с GPS в началото на 90 -те години.

Основните методи за определяне на географската дължина са изброени по -долу. С едно изключение (магнитно отклонение) всички те зависят от общ принцип, който беше да се определи абсолютно време от събитие или измерване и да се сравни съответното местно време на две различни места.

    . В орбитата си около Земята Луната се движи спрямо звездите със скорост малко над 0,5 °/час. Ъгълът между Луната и подходяща звезда се измерва със секстант и (след консултации с таблици и дълги изчисления) дава стойност за абсолютно време.
  • Сателити на Юпитер. Галилео предложи, че с достатъчно точно познаване на орбитите на спътниците, техните позиции могат да осигурят мярка за абсолютно време. Методът изисква телескоп, тъй като луните не се виждат с просто око.
  • Устреми, потъмнения и затъмнения. Пулсът е най -малкото видимо разстояние между два обекта (Луната, звезда или планета), затъмнение се случва, когато звезда или планета минава зад Луната - по същество е вид затъмнение. Лунните затъмнения продължават да се използват. Времената на всяко от тези събития могат да се използват като мярка за абсолютно време. . Часовникът е настроен към местното време на начална точка, чиято дължина е известна, а дължината на всяко друго място може да бъде определена чрез сравняване на местното време с времето на часовника.
  • Магнитно отклонение. Иглата на компас като цяло не сочи точно на север. Разликата от истинския север варира в зависимост от местоположението и се предполага, че това може да осигури основа за определяне на географската дължина.

С изключение на магнитното отклонение, всички доказани практически методи. Развитието на сушата и морето обаче беше много различно.

Няма друг физически принцип, определящ географската дължина директно, но с времето. Географската дължина в дадена точка може да бъде определена чрез изчисляване на времевата разлика между тази на нейното местоположение и координираното универсално време (UTC). Тъй като има 24 часа на ден и 360 градуса в кръг, слънцето се движи по небето със скорост 15 градуса на час (360 ° ÷ 24 часа = 15 ° на час). Така че, ако часовата зона, в която се намира човек, е три часа преди UTC, това лице е близо до 45 ° географска дължина (3 часа × 15 ° на час = 45 °). Думата близо до се използва, тъй като точката може да не е в центъра на часовата зона, а часовите зони са определени политически, така че техните центрове и граници често не лежат на меридиани, кратни на 15 °. За да извърши това изчисление обаче, човек трябва да има хронометър (часовник), настроен на UTC и трябва да определи местното време чрез слънчево или астрономическо наблюдение. Детайлите са по -сложни от описаните тук: вижте статиите за Универсалното време и за уравнението на времето за повече подробности.

Географската дължина се дава като ъглово измерване, вариращо от 0 ° на началния меридиан до +180 ° на изток и −180 ° на запад. Гръцката буква λ (ламбда), [35] [36] се използва за обозначаване на местоположението на място на Земята източно или западно от Първия меридиан.

Всяка степен на географска дължина е разделена на 60 минути, всяка от които е разделена на 60 секунди. По този начин географската дължина се посочва в шестнадесетична нотация като 23 ° 27 ′ 30 ″ E. За по -голяма точност секундите се посочват с десетична дроб. Алтернативно представяне използва градуси и минути, където части от минута се изразяват в десетична нотация с дроб, като по този начин: 23 ° 27,5 ′ Д. Степените могат да бъдат изразени и като десетична дроб: 23,45833 ° Д. За изчисления ъгловата мярка могат да бъдат преобразувани в радиани, така че географската дължина може също да бъде изразена по този начин като подписана част от π (pi) или беззначна част от 2 π.

За изчисления суфиксът Запад/Изток се заменя с отрицателен знак в западното полукълбо. Международната стандартна конвенция (ISO 6709)-че Изтокът е положителен-е в съответствие с десницата на декартовата координатна система със Северния полюс нагоре. След това определена географска дължина може да се комбинира със специфична географска ширина (положителна в северното полукълбо), за да се даде точна позиция на земната повърхност. Объркващо, понякога се вижда и конвенцията за негатив за Изтока, най -често в Съединените щати Лабораторията за изследване на земната система я използва на по -стара версия на една от техните страници, за да „направи координатното влизане по -малко неудобно“ за приложения, ограничени до западното полукълбо. Оттогава те преминаха към стандартния подход. [37]

Обърнете внимание, че географската дължина е единична при полюсите и изчисленията, които са достатъчно точни за други позиции, могат да бъдат неточни при или близо до полюсите. Също така с прекъсванията на меридиана ± 180 ° трябва да се работи внимателно при изчисленията. Пример е изчисляването на изместване на изток чрез изваждане на две дължини, което дава грешен отговор, ако двете позиции са от двете страни на този меридиан. За да избегнете тези сложности, помислете за замяна на географската ширина и дължина с друго представяне на хоризонтално положение при изчисление.

Дължината на степен на дължина (разстояние изток -запад) зависи само от радиуса на кръг с географска ширина. За сфера с радиус a този радиус на географска ширина φ е а cos φ , а дължината на една градусова (или π / 180 радиан) дъга по кръг от географска ширина е

Δ l o n g 1 = π 180 ∘ a cos ⁡ ϕ < displaystyle Delta _ < rm >^<1> = < frac < pi> <180^< circ >>> a cos phi>

φ Δ 1
лат
Δ 1
дълго
110.574 км 111.320 км
15° 110.649 км 107.551 км
30° 110.852 км 96.486 км
45° 111.133 км 78.847 км
60° 111.412 км 55.800 км
75° 111.618 км 28.902 км
90° 111.694 км 0,000 км

Когато Земята е моделирана от елипсоид, тази дължина на дъгата става [38] [39]

където е, ексцентрицитетът на елипсоида, е свързан с голямата и втората оси (съответно екваториалният и полярният радиус) чрез

Алтернативна формула е

Cos φ намалява от 1 на екватора до 0 на полюсите, което измерва как кръговете на географската ширина се свиват от екватора до точка на полюса, така че дължината на степен на дължина намалява по същия начин. Това контрастира с малкото (1%) увеличение на дължината на степен на географска ширина (разстояние север -юг), екватор до полюс. Таблицата показва както за елипсоида WGS84 с a = 6 378 137 .0 m и b = 6 356 752 .3142 m. Обърнете внимание, че разстоянието между две точки на 1 градус една от друга в същия кръг на географската ширина, измерено по този кръг от географската ширина, е малко повече от най -краткото (геодезическо) разстояние между тези точки (освен ако на екватора, където те са равни) разликата е по -малко от 0,6 м (2 фута).

Географска миля се определя като дължина на една дъгова дъга по екватора (една екваториална минута дължина), следователно степента на дължина по екватора е точно 60 географски мили или 111,3 километра, тъй като има 60 минути в градус . Дължината на 1 минута дължина по екватора е 1 географска миля или 1,855 км или 1,153 мили, докато дължината на 1 секунда от нея е 0,016 географска миля или 30,916 м или 101,43 фута.


5 Въведение в библиотеката PROJ.4

The PROJ.4 Понастоящем библиотеката има три функционалности за преобразуване на координати и геодезически изчисления:

proj/invproj - директни и обратни картографски проекции,

cs2cs - преобразуване между координатни системи (географски или проектирани),

geod/invgeod - директни и обратни геодезически изчисления.

The PROJ.4 библиотеката е разработена в началото на 1980 г. като програма RATFOR (език за програмиране RATional FORtran), като по -голямата част от кода произхожда от GCTP (Общият пакет за картографска трансформация на геоложката служба). The PROJ.4 програмата е прекодирана в езика за програмиране C, когато пакетът MAPGEN (от който PROJ.4 е основен компонент) е въведен в операционната система UNIX. Към програмата бяха добавени много нови прогнози.

Библиотеката се основава на работата на Джералд Евенден, по -късно се поддържа от Франк Уормърдам, създателят на библиотеката GDAL, а сега и от програмната общност. The PROJ.4 библиотеката се предлага с лиценз MIT (Масачузетски технологичен институт), позволяващ безплатно използване, копиране, редактиране, подлицензиране и комерсиално използване, ако оригиналният лиценз е кредитиран.

Инсталирането на PROJ.4 библиотеката зависи от операционната система по -долу, ще бъде дадена кратка илюстрация как да я използвате. Самата библиотека съдържа параметри на много елипсоиди, както и предварително дефинирани геодезически точки със стойности на параметрите.

Списъци с налични елипсоиди с параметри могат да бъдат получени с помощта на командата cs2cs -le, където -le е функционалният аргумент за показване на списъка с елипсоиди.

По същия начин може да се получи списък на основните меридиани, спрямо меридиана на Гринуич.

Списък с предварително дефинирани дати се получава с помощта на аргумента -ld.

Използвайки PROJ.4 нотация, може да се дефинира всяка координатна референтна система, а не само предварително дефинираните дати, съдържащи се в библиотеката.

За Сърбия е дадена координатна референтна система в проекцията на Гаус-Крюгер със следните параметри:

+proj = tmerc - проекция на Гаус Кригер (напречен меркатор),
+lat_0 = 0 +lon_0 = 21 - начало на координатите на елипсоида,
+k = 0,9999 - мащаб по средния меридиан,
+x_0 = 7500000 +y_0 = 0 - начало на координатите в равнината на картата
+ellps = bessel - елипсоид
+тегличи84 = 574.027,170.175,401.545,4.88786, -0.66524, -13.24673,6.89
три превода, три завъртания, скала спрямо WGS84, изразена в ppm (части на милион)
+= m - единици, метри.

Освен това, координатната референтна система в PROJ.4 нотация може да бъде дефинирана с помощта на EPSG кода. Наборът от данни за геодезически параметри на EPSG е съвкупност от дефиниции на координатна референтна система и параметри, които ги описват.

При http://spatialreference.org/ или http://epsg.io, дефиниции на координатни референтни системи могат да бъдат намерени в няколко означения: EPSG, proj4, WKT, GML, JSON ... На фигура 5.1 е даден пример за координатна система WGS84 с код EPSG 4326.

Фигура 5.1: Пример за използване http://spatialreference.org с възможност за изтегляне на параметрите на координатната референтна система в множество нотации.

Пример за директно изчисление е даден за проекцията на Гаус-Крюгер:

Когато изпълнява командата proj, потребителят се очаква да въведе елипсоидни координати чрез терминал или от текстов файл, за да получи правоъгълни координати. Използвайки функционалността на Linux, командата може да се изпълни директно от командния ред.

При обратни картографски изчисления е необходимо да се въведе командата -I, която следва да бъде последвана от въвеждане на правоъгълни координати и в резултат на това се получават елипсоидни координати (дължина и ширина).

Трансформацията на координатите от една система в друга е илюстрирана от примерното преобразуване от старата държавна система на Гаус-Крюгер (елипсоид на Бесел и датата на Херманскогел) в новата система за координати на UTM за Сърбия (елипсоид GRS80 с дата за ETRS за Сърбия):

Даден е и пример за преобразуване от системата Gauss-Krüger за Сърбия в WGS84, в който WGS84 е дефиниран с помощта на EPSG кода, а подробности за функцията могат да бъдат намерени на: http://proj4.org/apps/cs2cs.html.

The PROJ.4 библиотеката е внедрена в множество ГИС софтуер и може лесно да се използва чрез графичен интерфейс.


1 отговор 1

Тогава prasters = readandproject (c ("file1.nc", "file2.nc")) трябва да върне списък с два проектирани растра от имената на файловете. сюжет (прастери [[1]]), за да видите първия.

Също така заредих един от растерите в QGIS и замених CRS с WGS84 (epsg: 4326) и го сравних с базовия слой на OpenStreetMap и корелацията изглеждаше доста добре до разделителната способност на растъра.

Това е елегантно - странно spTransform дава грешка "не може да намери наследен метод за функция", но използването на projectRaster тук му позволява да се изпълнява. Имам профилен файл за границата на страната, който препроектирах и добавих към сюжета и той съвпада добре - достатъчно ли е това, за да сме сигурни, че е 4326? (Наистина оценявам вашата помощ, все още щях да остана в момента!)


2 отговора 2

За да конвертирате от long/lat в Albers равна площ, от тази статия в Wikipedia за проекцията на Albers:


Координатите от сферична точка могат да се трансформират в конични проекции с равни площи на Albers със следните формули, където λ е географската дължина, ╎ референтната дължина, φ географската ширина, ▼ референтната ширина и & #9661 и ▾ стандартните паралели:


Според О.С. Общата теория на Адамс за еквивалентни проекции (1945) (стр. 37, от статията, а не PDF):


Освен представянето на двете току -що описани полукълба, понякога е така
иска да представи целия свят на една карта. Ако желаете да имате равна площ
карта от този вид, би било необходимо да се използва конична проекция с минимум
деформации между северния полюс и 50 º южна ширина. Деформацията отвъд
паралелът на 50 º на юг не би бил обезпокоителен, тъй като отвъд не се намира важна земя
тази точка, тъй като само върхът на Южна Америка се простира по -на юг. Северният полюс
трябва да се вземе като център и разделянето да се извърши на 170 и#186 западна дължина
който преминава през Беринговия проток и не среща никаква земна площ. Тази проекция
съответства на м= 0,432, той не предизвиква никаква деформация по паралела на
18 䣁 'южно на северния полюс, особена точка на проекцията, 2 δ възлиза на
118 䣂 '. Най -голямата стойност на 2 δ освен тази точка е 58 䣓 'от а², 1.710 и от а², 2.924.



Координатна система

Man merke sich (für Anwendungen mit Schweiz-bezug) folgendes:

Für übriges Europa/Deutschland EPSG: 25832 (ETRS89 / UTM зона 32N). Wenn der Liegenschaftskataster von EPSG: 31467 (DHDN / 3-степенна Гаус-Крюгер зона 3) auf epsg: 25832 umgestellt wird, erfolgt of auch eine Bereinigung (Трансформация, Beseitigung historischer Fehler). Es reicht daher nicht, die Standard-Formeln für diese Koordinatensysteme anzuwenden. Es bleibt ein Restfehler von 1/2 bis 3 Metern е nach Gebiet. Man muss dann zusätzlich eine Passpunkt-Transformation nach NTv2 durchführen. Als Passpunkt-Datei kann z.B. Beta2007.gsb verwendet werden: https://duckduckgo.com/?q=NTv2+BETA2007.gsb

  • Siehe inkl. Formeln und Software (u.a. Applets): http://www.ottmarlabonde.de/
  • lat/lon или lon/lat ?:
    • Die geographische Breite (географска ширина, "Norting") hat mit dem Äquator einen natürlichen Nullpunkt mit +/- 90 in nördlicher und südlicher Richtung.
    • Die geographische Länge (дължина, "Easting") е ein Winkel ausgehend vom Nullmeridian (0 °) до 180 ° в östlicher и 180 ° в Westlicher Richtung.
    • Bei Webkarten wird meist lon/lat - също X/Y, Easting/Northing (z.B. 47.22 8.82 für Rappi) - erwartet.
    • 50% от Tools verwenden die eine oder andere Reihenfolge. (вижте http://www.macwright.org/lonlat/)

    Гледай видеото: #LatLongtoUTMConverter Latitude and Longitude to UTM Converter Free. Tutorial by RK Production