Posts

Visualization on The Globe

Image
Important note: we use only cross-platform Open Source Software and Open Datasets for our publications. You can freely download all the used software and data. These are our commonly used datasets in the same order as on the picture above: EMAG2 Earth Magnetic Anomaly Grid , Sandwell and Smith Free Air Gravity Grid , GEBCO Gridded Bathymetry Data , World Stress Map (vector fields ). That’s very helpful to produce a high-quality visualization for each of them to present our projects results.   ParaView   is the right tool to do it and we need to prepare the source data. Prepare datasets 1.1. Prepare georasters The original data grids are too large and that’s better to reduce the sizes. Also, only NetCDF files can be easy processed in   ParaView   by this way and so we convert the datasets to this format too.   GDAL tool   is suitable for the conversion: gdal_translate -of NetCDF -r average -tr .25 0.25 EMAG2_V2.tif EMAG2_V2.15min.nc gdal_translate -of NetCDF

Why gravity corrections for latitude/free-air/Bouguer/terrain are obsolete and the spatial filtered observed gravity must be used instead

Link:  https://www.linkedin.com/pulse/why-gravity-correction-latitudefree-airbouguerterrain-pechnikov/ Traditional geophysics uses a set of corrections (reduction) for the observed gravity measurements (https://en.wikipedia.org/wiki/Bouguer_anomaly). But what they mean? In terms of common physics these are just products of spatial bandpass filtering. However geophysics still based on these rough coefficients instead of using spatial bandpass filters. That’s really nonsense when these common tables of coefficients applied to high-resolution data observations! As example we have subcentimeter accuracy and ~100m spatial resolution of altimeter satellite measurements to calculate gravity on the sea surface level. In fact the table-based corrections are only tribute to traditions and these are absolutely useless. Let’s use spatial filters to extract the anomalies of all spatial scales.

Pearson correlation between spatial harmonics of free-air gravity anomaly on the Earth's surface and the Earth's topography is almost 1

Image
Link:  https://www.linkedin.com/pulse/pearson-correlation-between-spatial-harmonics-gravity-pechnikov/ We consider Free-air gravity anomaly and ETOPO1 topography in WGM2012 Earth's gravity anomalies datasets (see references below). There is high correlation (0.80) between topography and free air gravity for small areas only: For relatively large areas the correlation is low (0.41) but the scatter plot looks strange. Let's decompose in spatial domain the topography and the free air gravity and check the spatial harmonics. Pixel resolution of these datasets is about 3.7km. For spatial band pass pixel-wise filters 1-4px, 4-8px, 8-12px, 12-16px, 16-20px we have spatial scales 3.7-14.8km, and so on. We see on the pictures above that the correlation between the same spatial harmonics of topography and free air gravity is very high for as small as large areas. In fact we have the same results by all spatial harmonics based methods on topography and free air gra

Кольцевые структуры в геофизике

Опубликовал теоретическое обоснование и результаты численной модели для появления кольцевых структур в гравитационном поле. Поскольку найденные структуры существуют после применения фильтра высоких частот к полю, из этого следуют разные интересные вещи. Например, эти же структуры проявляются в цифровой модели рельефа и на космических снимках. Таким образом, подтверждаются практические работы Ю.И. Фивенского о существовании разномасштабных кольцевых структур, их видимости на рельефе и на снимках и пропорциональности их радиусов с глубиной залегания вызывающих их источников таких, как ловушки углеводородов, с коэффициентом около 0.7. Также можно предположить, что снежный покров работает как фильтр высокой частоты, что и позволяет лучше проявиться паттернам на поверхности - для тех конфигураций источников, которые иначе без фильтрации гравитационного поля не различимы. Ссылка на гитхаб-репозиторий со статьями: https://github.com/mobigroup/articles/  Первые две работы вводные, интересны

geomed3d-tclshp package for Shapefile creation from Tcl scripts

Понадобилось мне создавать шейпфайлы, попробовал делать это через gdal (типы полей в DBF, похоже, всегда строковые, не подходит, да и не понравилось - нужно создавать XML файл, в котором прописывать имена файлов, и проч. непотребства совершать) и Spatialite (типы полей могут быть и числовые, но все равно не устраивает - сложно и не быстро). Опробовав два вышеназванных пути, сделал свой пакет GeoMed3dSHP для языка Tcl (получается вдвое быстрее, чем делать то же через Spatialite). Пример использования: rm -f test.* && tclsh8.5 package require GeoMed3dSHP set id [SHPCreateFiles test] SHPWritePoint $id 374868.671871 2828378.81973 100.0 777 SHPWritePoint $id 374868.671871 2828378.81973 100.0 888 SHPCloseFiles $id exit Ну и заодно расширение GeoMed3dDBF для создания только DBF файлов: rm -f test.* && tclsh8.5 lappend auto_path . package require GeoMed3dDBF set id [DBFCreateFile test] DBFWrite $id 374868.671871 2828378.81973 100.0 777 DBFWrite $id 374868.671871 282837

Сканер CanoScan LiDE 210 в debian

Взял себе на днях CanoScan LiDE 210 взамен CanoScan LiDE 90, для которого поддержку так и не сделали в линуксе, а танцы с бубном мне надоели. Итак, подключил CanoScan LiDE 210, запустил xsane - все работает. Не настраивал вообще ничего, только установил xsane. Сканер шустрый и достаточно тихий. С поддержкой кнопок (в терминологии разработчиков sane - sensors) все хреново, как обычно, но есть финт ушами: $ scanimage -d "genesys:`sane-find-scanner|grep CanoScan|cut -d' ' -f 10`" -A|grep '\[hardware\]'|grep "\[yes\]"|wc -l 0 а теперь нажимаем любую кнопку: $ scanimage -d "genesys:`sane-find-scanner|grep CanoScan|cut -d' ' -f 10`" -A|grep '\[hardware\]'|grep "\[yes\]"|wc -l 1 Одна из кнопок сканера не вызывает никакой реакции, остальные работают (хотя их названия в выводе scanimage -A перепутаны). Так что несложно сделать скрипт для пакетного сканирования, реагирующий на кнопки: scan.sh #!/bin/bash SCANNER=

Открытый софт для научных расчетов

Задачки у меня бывают разные, кратенькое описание того, что может потребоваться я приводил здесь: Софт для спектрального анализа/численного моделирования/etc. Дело происходит в Linux Debian, соответственно, интересует софт, который устанавливается из стандартного репозитория (что подразумевает опенсорсность).  Для начала я протестировал двумерное Фурье преобразование (прямое и обратное) на широко известном тестовом изображении Lenna. Заодно в процессе поиска нужных функций увидел заметки об истории этого изображения  Просто Лена  и  Обработка изображений: кто такая Лена Выбор софта для такой типичной задачки вполне себе наличествует: Scientific Computing with Free software on GNU/Linux HOWTO Пакеты прикладных программ Инструменты визуализации данных для Linux На практике же оказалось не так радужно, как хотелось бы. Из того, что есть в дебиане, работает "из коробки" только пакет octave . Зато для него можно установить множество расширений - нап