Рельеф Земли - отклонения геоида

Следствия:

Влияние аномалий гравитационного поля Земли на динамику Мирового океана - Геоид - Гравитация - Карта Земли после сдвига полюса на 17 градусов - Категория:Погрешности модели рельефа - Решение геометрического искажения - Тяпкин, 1998 - Тяпкин, Довбнич, 2009 - Тяпкин 2012, Идентичность структурообразования континентов и океанов - Тяпкин 2012, Изменение положения оси вращения в теле Земли: причина, механизм и использование

^{комментарии в LiveJournal} Нарисуем форму Земли математически точно.

Известно, что Земля - это эллипсоид с искажениями. Искажения называются "отклонениями геоида" (по английски "geoid undulations").

Беру формулу эллипса в полярных координатах из википедии - . Подставляю в данные Earth2014 - и тут меня ожидало разочарование: диапазон отклонений геоида оказался на 50 метров больше ожидаемого. Вместо 192м я получил 250м отклонений.

Ожидаемые отклонения геоида беру из английской и русской википедии от -107 м до +85,4 м; диапазон 192,4 м.

Проверил параметры эллипсоида: Earth2014 использует GRS80, википедия - WGS84, но радиусы у них идентичны (радиус оси экватора a = 6 378 137,0 м; радиус полярной оси b = 6 356 752,3141 м).

Оказалось, что не ту формулу эллипса я взял. Нужно брать формулу от автора данных Earth2014 - Кристиана Херта. Находится она в файле Hirt_Rexer2015_Earth2014.pdf.

Использование GIS

Пытался нарисовать карту в геоинформационной системе (называю в дальнейшем GIS). Использовал бесплатную GIS GRASS. Наткнулся на ограничение: медленное преобразование карт при изменениях датума. GIS вроде как и поддерживают расчет наклона оси вращения оси, но в GRASS этих настроек не нашел. Также выяснилось, что при интенсивном изменении датума скорость GIS существенно падает. А моделировать изменение скорости и угла вращения Земли нужно обязательно.

Поэтому рисую карты самодельными алгоритмами.

Нарисовал Землю Файл:Earth.gif

Мерцания на береговых линиях показывает погрешности преобразования из HEALPix в равнопромежуточную проекцию. Палитры цветов получены из GRASS: water и модифицированный gyr.

Алгоритм

Алгоритм расчета отклонений геоида выглядит так:

        protected override double GetAltitude(HealCoor pixel)
        {
            var a = 6378137d; // big axis, angle is measured from here
            var b = 6356752.3141;
            var e_2 = 1 - (b * b) / (a * a);

            var altitude = Relief.GetAltitude(pixel);
            var altitudeShape = ReliefBed.GetAltitude(pixel) + 6371000;

            var angle = Math.PI / 2 - pixel.Theta.Value;
            var sin_2 = Math.Sin(angle) * Math.Sin(angle);
            var ellipsoid = a * Math.Sqrt((1 - e_2 * (2 - e_2) * sin_2) / (1 - e_2 * sin_2));
            var undulation = altitudeShape - ellipsoid - altitude;
            return undulation;
        }

При его выполнении нашел закладку от Кристиана Херта. Я почти отправил ему письмо про как бы ошибку, но подумав, уразумел, что это просто "защита от дурака".

Смысл в том, представляя данные Earth2014 в двух вариантах точности: в точности 1 минута дуги и [5 минут дуги, Christian Hirt публикует файлы в равнопромежуточной проекции. А известно, что отклонения геоида версии EGM96 сделаны в сферических функциях. Преобразования между ними нетривиальные, поэтому предполагаю, что для неподготовленной публики он преобразовал EGM96 в равнопромежуточную проекцию красиво Но реально полезные данные из такого изображения в проекции HEALPix будут искажены

Я уже начал сомневаться в выборе проекции HEALPix для расчетов.

Но потом глянул данные Earth2014 в точности 1 минута дуги.

Оказалось, что там EGM96 преобразована правильно для расчетов: в равнопромежуточной проекции отклонения геоида визуально искажены , а в проекции HEALPix наоборот красивы для расчетов