Процедура
FILTER
Назначение
процедуры: фильтрация поля по размерам и простиранию аномалий.
При
работе процедура запрашивает:
1. Filter Type - Тип фильтра ( 1, 2 или 3),
где
1 - полосовой по простиранию,
2 - прожекторный по простиранию,
3 - полосовой размера.
2. Параметры
Freq1,
Freq2,
Freq3
и Freq4 - узловые точки фильтра ( все
целые). Для фильтра по простиранию ( тип 1 или 2 )
параметры задаются углами с осью Х в градусах. Для полосового фильтра размера ( тип 3 ) значения параметров определяются в (1/км)*1000.
Примечания:
1.
Для типов фильтра 1, 2 значения углов должны
заключаться в пределах (-90, 90) градусов.
2.
В процедуре FILTER
при определении параметров допускается F1=F2
и F3=F4
(
идеальный фильтр ).
Можно
также задавать F1 = F2 = 0, в этом случае мы получаем
низкочастотный фильтр.
Чтобы
задать фильтр высоких частот, необходимо задать
очень большие значения
параметрам F3 и F4.
Процедура присвоит им значения частоты Найквиста
F3=F4=0.5/DX,где DX - шаг сетки карты.
На рис. 2
показаны примеры операторов фильтра размера.
1) низкочастотный фильтр
1.I*********
I .*
F1=F2=0
I . *
I . *
I . *
I . *
0.I -------------*----------------> KX
F3 F4
2) высокочастотный фильтр
1.I *******************
I
*. F3=10000000
I * . F4=10000000
I * .
I * .
0.I------*-------------------------> KX
F1 F2
3) полосовой
фильтр
I **************
I *.
.*
I * .
. *
I * .
. *
I * .
. *
0.I------*----------------------*--------> KX
F1 F2
F3 F4
Пусть,
например, необходимо выделить аномалии с размерами от 1000 м до 2000 м, шаг
сетки равен 250 м, тогда F2=250, F3=500. Чтобы уменьшить явления Гиббса,
определяем с помощью F1 и F4 "переходные области" :
F1=200, F4=600.