Построение стереограмм трещиноватости является повседневной задачей в структурной геологии. Для её решения уже разработано большое количество программных средств под Windows, но под Linux специальной утилиты пока не существует. Есть лишь небольшое количество малофункциональных программ.
Но так как строить такие стереограммы мне необходимо, да еще и в больших количествах, я написал скрипт, осуществляющий такие построения, используя функционал GMT.
#!/bin/bash
fig=test.ps
title=»Проверка»
grid=test.grd
cpt=test.cpt
gmtset CHAR_ENCODING ISO-8859-5
echo ‘0, -90’ | psxy -R0/360/-90/0 \
-JS0/-90/16c -Sx0.4c -Ba90f90/g90:.$title: \
—BASEMAP_TYPE=plain -K -V > $fig
#рисуем нормали к плоскостям
cat test_n.txt | awk -F»\t» ‘{print $1, $2-90}’ \
| psxy -R -J -Sc0.2c -N -Gblack -Wwhite -O -K >> $fig
#рисуем направления смещений по плоскостям
cat test_v.txt | awk -F»\t» ‘{print $1, $2-90, -$3+90, «1»}’ \
| psxy -R -J -SV0.08/0.3c/0.12c -N -Gblack -O -K >> $fig
echo «50 5 12 0 0 CT Тестовая стереограмма» | pstext -R0/50/0/50 -JX16c -N -O >> $fig
Результатом выполнения скрипта будет вот такая стереограмма: