Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Программы и алгоритмы для обработки данных дистанционного зондирования: ERDAS, ENVI и другие.
MBG
Участник
Сообщения: 83
Зарегистрирован: 02 июн 2012, 12:02
Репутация: 18

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение MBG » 12 июн 2019, 15:23

gamm писал(а):
12 июн 2019, 12:41
У вас вся корреляцият "наведенная", за счет использования спектрального представления.
Это что за ерунда, простите? Вы снятием порчи подрабатываете, видимо. Не бывает "наведенной" корреляции, а спектральное представление это известный и стандартный метод, гугл вам в помощь.

Update: выделение локальной и региональной составляющей полей это тоже пространственная фильтрация high-pass и low-pass фильтрами, технически может реализовываться и в пространственной и в частотной плоскостях и эти подходы эквивалентны.

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 1828
Зарегистрирован: 04 янв 2011, 22:00
Статьи: 12
Проекты: 1
Репутация: 1141
Откуда: Казань

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение Игорь Белов » 15 июн 2019, 20:37

MBG писал(а):
12 июн 2019, 13:39
Вот и подумайте, что именно вы измеряете вашими гравиметрами.
Земная кора неоднородна, поэтому гравиразведка работает. Откройте для себя геологию и геофизику.

gamm
Гуру
Сообщения: 2962
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 687
Ваше звание: программист
Откуда: Казань

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение gamm » 16 июн 2019, 07:05

MBG писал(а):
12 июн 2019, 15:23
Это что за ерунда, простите?
гистограмму невязок покажите, потом можно будет что-то обсуждать.

У вас же спектр разлагается глобально, насколько я понимаю, и к локальным вариациям высокие частоты имеют мало отношения (если только локальные вариации не периодические). У вас гармоники от глобальной составляющей преобладают, если так понятнее - а они одинаковые, коррелятор вы наш.

MBG
Участник
Сообщения: 83
Зарегистрирован: 02 июн 2012, 12:02
Репутация: 18

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение MBG » 16 июн 2019, 07:20

Игорь Белов писал(а):
15 июн 2019, 20:37
Вы даже с цифрами в руках не понимаете, что в измеренных значениях гравиметрии влияние рельефа превышает влияние погребенных неоднородностей?

Окей, смотрите учебник "Г Р А В И Р А З В Е Д К А: СПОСОБЫ УЧЕТА ВЛИЯНИЯ РЕЛЬЕФА МЕСТНОСТИ"
http://www.ifz.ru/fileadmin/user_upload ... sticyn.pdf
Необходимый элемент обработки гравиметрических материалов – это вычисление поправок за влияние рельефа местности, которое обычно выполняется после завершения полевого сезона. В соответствии с требованиями «Инструкции по гравиразведке» [24] поправки за влияние рельефа местности вводятся в наблюденные значения силы тяжести в тех пунктах, где они превышают 50% величины проектной точности съемки.
В данном топике речь идет о технологии для "полная автоматизация вычислений поправки во всей области учитываемого влияния рельефа" из вот этого пайплайна (см. указанный учебник):
Во второй части учебного пособия дается описание современной компьютерной технологии высокоточного определения поправок за влияние рельефа местности при гравиметрической съемке. Ее основными элементами являются: использование в качестве исходных данных векторизованных крупномасштабных топографических карт и цифровых моделей рельефа (GTOPO30, SRTM); построение аналитической модели рельефа путем двухмерного преобразования Фурье; полная автоматизация вычислений поправки во всей области учитываемого влияния рельефа; стохастическая оценка точности получаемых результатов.
Поправка эта представляет собой 1) полосу частот спектров Фурье гравики и топографии, в которой компоненты спектров совпадают и 2) коэффициент масштаба для этой полосы частот. Очевидно, эта поправка может быть использована как для вычитания поправки за рельеф из гравиметрических измерений, так и для синтеза данных гравиметрии более высокого разрешения.

Все это отлично работает практически, см.
"Study of the Earth׳s short-scale gravity field using the ERTM2160 gravity model", https://doi.org/10.1016/j.cageo.2014.09.001
Residual gravity model ERTM2160 computed from the SRTM topography at 250 m resolution.
convert the high-pass filtered SRTM topography – complemented with bathymetric information in coastal zones – to implied short-scale gravity effects.
Фурье преобразование имеет свои проблемы, поэтому сейчас чаще пользуются дискретными преобразованиями (вейвлетами), в опубликованном мною коде используется преобразование, аналогичное следующей публикации:
"High-resolution regional gravity field recovery from Poisson wavelets using heterogeneous observational techniques", https://doi.org/10.1186/s40623-017-0618-2

Изображение

MBG
Участник
Сообщения: 83
Зарегистрирован: 02 июн 2012, 12:02
Репутация: 18

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение MBG » 16 июн 2019, 07:24

gamm писал(а):
16 июн 2019, 07:05
У вас же спектр разлагается глобально, насколько я понимаю
Абсолютно неверно. В опубликованном мною коде используется свертка с двумерным гауссианом. Это локальная операция, очевидно.

gamm
Гуру
Сообщения: 2962
Зарегистрирован: 15 окт 2010, 08:33
Репутация: 687
Ваше звание: программист
Откуда: Казань

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение gamm » 16 июн 2019, 08:01

MBG писал(а):
16 июн 2019, 07:24
В опубликованном мною коде используется свертка с двумерным гауссианом.
откуда гауссиан, если на картинке написано Пуассон? Если там вейвлеты, все равно гармоники получите, хоть и меньше, разложите меандр низкочастотный, и посмотрите.

Пока нет гистограммы локальных невязок модели с фактическими измерениями, говорить все равно не о чем - это вопросы веры, и красивых картинок. Это ровно те же методы, которые применяют в космосъемке при раскрашивании панхрома по спектральным каналам низкого разрешения, картинка красивая, к реальности отношения имеет мало.

MBG
Участник
Сообщения: 83
Зарегистрирован: 02 июн 2012, 12:02
Репутация: 18

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение MBG » 16 июн 2019, 11:09

gamm писал(а):
16 июн 2019, 08:01
откуда гауссиан
Открываем код по ссылке из статьи:
https://github.com/mobigroup/gis-snippe ... tion.ipynb

Код: Выделить всё

from scipy.ndimage.filters import gaussian_filter
# low-pass filter
def raster_filter(src, gamma):
    ...
    dst.values = gaussian_filter(...)
    ...
Вам на самом деле не понятно, что gaussian_filter это фильтр Гаусса?

MBG
Участник
Сообщения: 83
Зарегистрирован: 02 июн 2012, 12:02
Репутация: 18

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение MBG » 16 июн 2019, 11:15

gamm писал(а):
16 июн 2019, 08:01
Если там вейвлеты, все равно гармоники получите
И словами написано, и код опубликован, а вы все гадаете. Есть много способов спектрального разложения, теоретически они эквивалентны, на практике есть определенные проблемы у каждого из них. Притом вы уже сами поняли, почему нам не подходят и Фурье и вейвлеты. Использованная свертка с гауссианом требует больше вычислений, зато избавляет от вышеназванных проблем.

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 1828
Зарегистрирован: 04 янв 2011, 22:00
Статьи: 12
Проекты: 1
Репутация: 1141
Откуда: Казань

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение Игорь Белов » 16 июн 2019, 11:28

MBG писал(а):
16 июн 2019, 07:20
Вы даже с цифрами в руках не понимаете, что в измеренных значениях гравиметрии влияние рельефа превышает влияние погребенных неоднородностей? …
Всё это очень интересно. Вот только не содержит ответа на цитату. Плюс очевидно, что смысл введения поправок Вам неясен. Между тем после их введения как раз остаётся поле аномалий силы тяжести, демонстрирующее неоднородности плотности земных пород.

Диагноз: отсутствуют базовые знания в геонауках, без чего крайне самонадеянно считать, что понимаешь смысл современных методов.

Но нет повода так переживать, всё можно наверстать. Начните с азбуки, с базовых учебников. Правда, прежде чем изучать геофизику, настоятельно рекомендую проштудировать теорию вероятностей и мат. статистику, иначе не будете понимать, что такое измерения. Удачи!

MBG
Участник
Сообщения: 83
Зарегистрирован: 02 июн 2012, 12:02
Репутация: 18

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение MBG » 16 июн 2019, 11:31

gamm писал(а):
16 июн 2019, 08:01
Пока нет гистограммы локальных невязок модели с фактическими измерениями
Еще раз повторяю - покажите ваш датасет с фактическими измерениями, и проверим. Иначе получается как с вашим товарищем по секретности из этого топика, который никак не поймет, зачем нужна коррекция измеренных значений гравиметрии на рельеф - немудрено, что результаты его работы лучше засекретить навсегда :)
gamm писал(а):
16 июн 2019, 08:01
Это ровно те же методы
Вы никак не можете понять, что у нас за метод обсуждается, зато готовы утверждать про "ровно те же методы". Чуть ранее вы уж про порченный спектр писали, и про наведенную на него корреляцию... Впрочем, если вы готовы показать ссылки на работы, где спектры/вейвлеты "применяют в космосъемке при раскрашивании панхрома по спектральным каналам низкого разрешения" покажите. Хотя даже этот же самый форум вам объяснит, что это делается намного проще:
viewtopic.php?t=4973

MBG
Участник
Сообщения: 83
Зарегистрирован: 02 июн 2012, 12:02
Репутация: 18

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение MBG » 16 июн 2019, 11:44

Игорь Белов писал(а):
16 июн 2019, 11:28
Итак, несмотря на приведенные ссылки и объяснения, вы снова утверждаете про отсутствие корреляции между высокочастотными компонентами измеренных значений гравиметрии и рельефа? Как я писал выше, для использования детальной гравиметрической съемки для поиска погребенных неоднородностей необходимо выполнить коррекцию за рельеф с точностью значений высот порядка миллиметра для измеренных гравиметрических значений порядка микрогала. Это легко посчитать, благо значение первой производной силы тяжести известно - и выше я уже привел пример вычислений.

Также смотрим публикации, например,
ОПЫТ ИСПОЛЬЗОВАНИЯ ВЫСОКОТОЧНОЙ ГРАВИРАЗВЕДКИ ДЛЯ КАРТИРОВАНИЯ ЗОН
РАЗУПЛОТНЕНИЯ ПРИ РЕШЕНИИ ИНЖЕНЕРНО-ГЕОЛОГИЧЕСКИХ ЗАДАЧ.

Точность измерений поля силы тяжести в условиях городских помех среднего уровня (отсутствие в
окрестности объекта предприятий, относительно слабая загруженность улиц, линия метро не ближе
500м) составила ±2-3 мкГал, точность определения высот точек наблюдения ±2-3 мм.
точность гравиметрической съемки ±5 мкГал, точность определения высот точек наблюдения ±1 см
Вы же, видимо, используете измеренные данные гравиметрии для поиска малых неоднородностей без коррекции за рельеф, да еще утверждаете об отстутствии связи между этими данными.

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 1828
Зарегистрирован: 04 янв 2011, 22:00
Статьи: 12
Проекты: 1
Репутация: 1141
Откуда: Казань

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение Игорь Белов » 16 июн 2019, 11:55

MBG писал(а):
16 июн 2019, 11:44
вы снова утверждаете про отсутствие корреляции между высокочастотными компонентами измеренных значений гравиметрии и рельефа?
Вот покажите, где я это утверждал! Речь вовсе не об отсутствии корреляции, а о том, что её высокая величина (средняя по больнице) сама по себе не позволяет получить карту аномалий со сколько-нибудь удовлетворительной точностью. Это факт, подтверждённый полевыми измерениями, хотя и без них тут всё было очевидно.
MBG писал(а):
16 июн 2019, 11:44
Вы же, видимо, используете измеренные данные гравиметрии для поиска малых неоднородностей без коррекции за рельеф, да еще утверждаете об отстутствии связи между этими данными.
Да что Вы говорите :lol:
Вот что должно быть в голове у человека, только получившего диплом, чтобы поучать специалистов, отработавших десятки лет, как им надо обрабатывать измерения!?

Кстати, ещё раз о Вашем непонимании смысла поправок: они нужны именно для выделения неоднородностей, которые могут маскироваться влиянием рельефа. Но! Выраженные структуры отлично видны даже на карте сырых измерений.

MBG
Участник
Сообщения: 83
Зарегистрирован: 02 июн 2012, 12:02
Репутация: 18

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение MBG » 16 июн 2019, 12:28

Игорь Белов писал(а):
16 июн 2019, 11:55
Вот покажите, где я это утверждал!
Да легко:
Игорь Белов писал(а):
11 июн 2019, 11:07
На картинках разности между значениями аномалии в свободном воздухе из полевых измерений и моделью GGMPlus. Размах разностей составляет 40% от размаха значений аномалий, и это до снятия тренда.
Модель GGMplus использует рельеф SRTM (точность по высоте ±16 m). Поскольку вы сравниваете свои данные с данными GGMplus, необходимо, чтобы ваши данные тоже были скорректированы по этому же рельефу - разномасштабные данне сравнивать нельзя, разумеется, могу ссылки дать, если сомневаетесь. Это соответствует точности гравиметрических измерений порядка десятка миллигал. Никаких микронеоднородностей плотности вы там не найдете, пока не откорректируете по рельефу с вертикальной точностью порядка миллиметра для измеренных значениях порядка микрогала.

P.S. Если же вы, вдруг, сами не зная о том, корректируете по детальному рельефу с вертикальной точностью в миллиметры - это полностью бессмысленно, сравнивать результат с данными GGMplus.

Аватара пользователя
Игорь Белов
Гуру
Сообщения: 1828
Зарегистрирован: 04 янв 2011, 22:00
Статьи: 12
Проекты: 1
Репутация: 1141
Откуда: Казань

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение Игорь Белов » 16 июн 2019, 12:39

MBG писал(а):
16 июн 2019, 12:28
Да легко:

Игорь Белов писал(а): ↑
11 июн 2019, 12:07
На картинках разности между значениями аномалии в свободном воздухе из полевых измерений и моделью GGMPlus. Размах разностей составляет 40% от размаха значений аномалий, и это до снятия тренда.
И? Коэффициенты корреляции там практически единичные, но расхождения между моделями говорят о неприменимости GGMplus как замены реальных измерений.
MBG писал(а):
16 июн 2019, 12:28
Если же вы, вдруг, сами не зная о том, корректируете по детальному рельефу с вертикальной точностью в миллиметры - это полностью бессмысленно, сравнивать результат с данными GGMplus.
Поздравляю, вот мы и пришли к тому, с чего я начал: модель неприменима для практических целей.
Последний раз редактировалось Игорь Белов 16 июн 2019, 12:41, всего редактировалось 1 раз.

MBG
Участник
Сообщения: 83
Зарегистрирован: 02 июн 2012, 12:02
Репутация: 18

Re: Build Super-resolution Gravity from GGMplus Free-Air Gravity Anomaly (200m) enhanced by SRTM topography (30m)

Сообщение MBG » 16 июн 2019, 12:40

Игорь Белов писал(а):
16 июн 2019, 11:55
Вот что должно быть в голове у человека, только получившего диплом
Это вы про себя? Откуда нам знать, что у вас в голове :) Моя магистерская диссертация по созданию голограмм в нелинейных средах была признана лучшей работой по нелинейной оптике по России 15 лет назад. На линкедине в профиле есть пруф. Так что не волнуйтесь - и студентам объяснял, что такое пространственный спектр, и даже вам объясню, независимо от вашего желания :)

Ответить

Вернуться в «Обработка ДДЗ»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 1 гость