Как мы показали в предыдущих статьях, поверхностные геологические нарушения (линеаменты) могут быть выделены с помощью линейного преобразования Радона (Хафа), в то время как геологическая плотность может быть восстановлена с помощью кольцевого преобразования Радона (Хафа). При использовании пространственной фильтрации исходных данных можно сопоставить длину волны с глубиной и получить трехмерную модель. При этом, кольцевое преобразование имеет понятные ограничения — кольцевые статистики радиусом меньше 3х пикселов недостоверны, так что, скажем, при использовании открытого рельефа всей планеты разрешением 30 м верхние 40 метров (2*30/sqrt(2)) трехмерной модели являются результатом интерполяции. В случае же линейного преобразования мы избавлены от этого ограничения, поскольку анализируем именно выходы глубинных структур на поверхность.
Сегодня мы построим объемную линеаментную модель и сравним ее с плотностной моделью и интерферограмммой сейсмического события из предыдущих статей. Как будет показано, эти модели согласуются между собой и дополняют друг друга.
