Файл: Коваль Л.А. Автоматизированная система обработки данных магниторазведки с применением ЭВМ (АСОМ-АМ).pdf

3. В случае, если ДОПИНФ, приготовленный "Матрицей", не под­ ходит (не хватает градаций, нужен другой номер шрифта, необходимо изменить масштаб, жестко задаваемый в системе и т.п.), нетрудно приготовить новый ДОПИНФ вручную.

На рис. 33, 34, 35 и 36 приведены образцы карт графиков и карт изодвнам лТ , вычерченных с использованием программы Ш на АТЛАСЕ. Сеть значений для каждого из образцов в масштабе изобра­ жений равна 2 х 10 мм или I х 10 мм, т.е. отличаются большой ани­ зотропией. Исходные значения поля не подвергались какому-либо сглаживанию, кроме того, который здесь именуется первым уровнем (сглаживания). Между тем карты имеют удовлеторительный рабочий вид, хотя и возможно, нуждаются в некоторой редакции.

Участки съемки обычно имеют сложную конфигурацию и,даже "раз­ бивая" площадь на части, как правило, не удается обеспечить стро­ го прямоугольную форму фрагментов. Заполнение строк прямоугольной матрицы в местах, где не обеспечена интерполяция исходных наблю­ дений, ведется линейной интерполяцией по строке или по краям про­ филя, присвоением ближайших верных значений, т.е., например, если первые (последние) 10 точек строки матрицы находились вне участка съемки, то им будет присвоено значение одиннадцатой (с конца) точки строки. Области пропусков съемки имеют своеобразный "поло­ сатый" вид, их легко выделить и заклеить.

Е. МИС-ГЗ. Тригонометрическое сглаживание-анализ (ячейки 3500-4170)

Программа выполняет гармонический анализ заданного профиля с выбором параметра сглаживания и наглядным изображением резуль­ татов анализа (описание алгоритма см. в § 8, раздел "В"). В про­ грамме используется ключ 0001 (включен,если анализ ведется с ри­ сованием; выключен, если рисовать не надо), исходный числовой ма­ териал на ЛПМ I . Номер профиля дается с пульта на клавишном набо ре. Перед обращением к МИС-13 ШАПКУ следует пополнить уровнем сглаживания

I . 3500-3513. Вызов с ЛПМ I ДОПИНФА-4, снятие о клавишного набора номера анализируемого профиля. Поиск в ДОПИНФЕ-4 информа­ ции о заданном профиле. Предусмотрен останов СчАК 3513 (с замком), если набран номер, отсутствующий в КАТАЛОГЕ: (СМ) = (клавишный набор). Используется СП 5 - обмен с НМД.

11Б

2.3514-3566. Определенно со (8-17) числа палеток г, умещаю­ щихся в профиле, определение по (8-18) и (8-19) местоположения и числа точек в каждой палетке. Используется СП 10 - перевод целого двоичного числа в 2 с.с. с плавающей запятой.

3.3567-3625. Начало цикла по г (число палеток). Вызов (во­ зобновление) профиля в ГОЗУ, забивание кодов ошибок интерполиро­ ванными значениями. Формирование и выдача на печать значений дТ,

относящихся к очередной палетке. Используются СП: 5 - обмен с НМЛ; 12 - ш1терполяция на ошибочные точки; 7 - перевод 2—10 с плавающей запятой массива чисел.

4.3626-3651. Вычисление по формулам (8-8) и (8-7) линейного фона и д. Последняя величина выдается на печать. Используются СП: 10 - перевод целого двоичного числа в 2 с.с. с плавающей за­ пятой; 7 - перевод 2—10 с плавающей запятой массива чисел.

5.3652-3703. Вычисление по (8-10) значений Ь(К) и выдача их на печать в 10 с.с. с плавающей запятой. Используются СП:

10 - перевод целого двоичного числа в 2 с.с. с плавающей запятой; 16 - stnx(cosx) ; 7 - перевод 2—10 с плавающей запятой массива чисел.

6.3704-3753. Определение из условия минимума (8-16) номера последней оставляемой при синтезе гармоники m , величины R ;

выдача на печать п" (число точек в палетке), m и R . Исполь­ зуются СП: 10 - перевод целого двоичного числа в 2 с.с. с плаваю­ щей запятой; 7 - перевод 2—10 с плавающей запятой массива чисел.

7.3754-4002. Вычисление по ( 8 - I I ) и (8-12) сглаженной кри­ вой. Выдача на печать сглаженной кривой. Используются СП: 16 - sinх(cosх) ; 7 - перевод 2—10 с плавающей запятой массива чисел.

8.4003-4005. Переход по ключу к рисованию на АТЛАСЕ. Цикл

по г (количество палеток). Останов СчАК 4006.

9. 4010-4170. Рисование результатов гармонического анализа на АТЛАСЕ. Используются СП: 14, 15 - "Линия", "Надпись"; 9 - пе­ ревод 2—10 с кодовой запятой и заданной точностью; 3 - 'График". Предусмотрен останов СчАК 4011 - в сумматоре порядковый номер па­ летки: включить АТЛАС, перописец в левый нижний угол планшета, пускI

Выдачу на печать МИС-13 открывают два десятичных числа, на­ печатанных в форме "с запятой". Первое число - номер профиля, второе - номер палетки. Далее через интервал печатается п* де­ сятеричных чисел с плавающей запятой, составляющих палетку

117

анализируемой величины.

Графический вывод на АТЛАС (рис. 27) открывается заголовком, в котором после номера профиля указывается номер палетки. Во вто­ рой строке заголовка указываются значения N, М, R , масштаб графиков.

Графики - исходный, сглаженный и коэффициентов Ь(К) вычер­ чиваются с шагом 0,5 см между точками. На первых двух графиках оцифровка означает номера точек (в палетке), на последнем - но­ мера гармоник. На сглаженном графике первая и последняя точки, не входящие в сглаживание, не показываются. Аналогично, нулевая и последняя гармоники, не имеющие смысла, не показываются на графике коэффициентов Ь(К) .

Ж. МИС-14. Тригонометрическое сглаживание числового материала (ячейки 5500-6010)

Программа обеспечивает тригонометрическое сглаживание профи­ лей в пределах участка съемки или фрагмента (по включенному ключу 0001). При сглаживании по фрагменту профили, не имеющие общих то­ чек с полосой y m l n , y m m t , подвергаются первому уровню сгла­ живания (линейная интерполяция на ошибочные точки). Исходный чис­ ловой материал на ЛПМ I,сглаженные профили вместе с Д0ПИНФ0М-4бис переписываются на ЛПМ 2. Перед обращением к МИС-14 ШАПКУ следует пополнить параметром сглаживания R и, если расчеты будут вес­ тись по фрагменту, уточнить координаты его угловых точек.

1. 5500-5505, 5750-5771. Вызов с ЛПМ I ДОПИНФА-4, определе­ ние полосы ymi.„, ymox,включающую фрагмент. Используется СП-5

-обмен с НМЛ.

2.5506-5514, 5772-6010. Начало цикла по профилям. Вызов с М03У профиля с ЛПМ I . Первый уровень сглаживания и определение (при включенном ключе 0001) имеет ли профиль общие точки с поло­

сой y m l n , y m e x . Если имеет - переход к следующему блоку; ина­

116

и числа точек в каждой палетке. Используется СП 10 - перевод це­ лого двоичного числа в двоичную систему счисления с плавающей запятой.

4. 5577-5716, 5745-5747. Цикл по г (число палеток) внутри профиля. Выделение палеток, вычисления по формулам (8-8) и (8-7)-

и (8-12) сглаженной кривой. Композиция с учетом перекрытий сгла­ женного проф;ля. Используется СП: 10 - перевод целого двоичного числа в 2 с.с. с плавающей запятой; 16 - sinx (cosx).

5.5717-5736. Циклическое суммирование сглаженного профиля. Запись профиля на ЛПМ 2. Исправление КАТАЛОГА для Д0ПИНФА-4-бис. Цикл по профилям. Используется СП 5 - обмен с НМЛ.

6.5733-5744. Заключительный блок в программе: циклическое суммирование Д0ПИНФА-4-бис, переписывание его на ЛПМ 2. Исполь­ зуется СП 5 - обмен с НМЛ. СчАК 5745 - останов.

3. МИС-17, 18. Итерационное сглаживание (ячейки 6160-6776, 7671-7776)

Две программы МИС-17 (анализ сглаживания по Савинскому) и МИС-18 (сглаживание числового материала по Савинскому) объединя­ ются здесь в одном описании в силу того, что имеют ряд общих блоков, выполненных в стиле СП, в отличие от более автономных

-31 00 6640 0017 в ячейку 0040 заносится код +п-1 0000 0000, где п - число точек в сглаживаемой кривой, а сама кривая распо­ лагается в ячейках МОЗУ, начиная с 1500 (действует обычное в си­ стеме ограничение - не больше I 02 IJQ точки). На выходе коэффици­

т о ч к а м (ячейки 7671-7712). Блок осуществляет однократное сглаживание по формуле (8-28) числового материала, располагающе­

)

Сглаженная кривая располагается на месте исходной. Блок исполь­

индексные ячейки 0015-0017, СП 10 - перевод целого двоичного чис­ ла в двоичную систему счисления с плавающей запятой и СП 6 - ко­ рень квадратный.

циальная перезапись исходной сглаживающей матрицы в ячейках 1100-

ячейки 6340-6637). В программе сохранен блок, позволяющий анали­ зировать сглаживание по Савинскому для кривых, специально подго­ товленных и введенных в память машины. Поэтому программа имеет два начала.

6340: перед обращением к этой команде в МОЗУ с ячейки 1500 •должна находиться испытуемая кривая, а в ячейках 6000-6004, соот­ ветственно, номер профиля, п - число точек, к - число итера­

обращением к МИС-17 с началом 6360 необходимо пополнить ШАПКУ

120