Способ геометрической коррекции космических сканерных снимков линейной развертки

(ц 1) С 01 С 11/00 он При налици сканирозания векторограмма рассчитана заВектора см выбранные оси х ара данных о борта ИСмещений подобнаожет быт нее.щений и икоординат проекции ногут быть ГОСУДАРСТВЕННЫЙ НОМИТЕТПО ИЗОБРЕТЕНИЯМ И ОТНРЫТИНМПРИ ГННТ ССО(71) Лаборатория азрометодов Производственного геологического объединения "Аэрогеология"(56) Сажин С.М. О способе геометрического преобразования многозональной спутниковой информации: Труды//ГосНИЦИПР, 1979, вып.5, с.46-49,(54) СПОСОБ ГЕОМЕТРИЧЕСКОЙ КОРРЕКЦИИКОСМИЧЕСКИХ СКАНЕРНЫХ СНИМКОВ ЛИНЕЙНОЙ РАЗВЕРТКИ(57) Изобретение может быть использовано в геологии, сельском хозяйстве,природоохранной тематике и т,пЦель Изобретение относится к фотограмметрии и может быть использовано в геологии, сельском хозяйстве, природоохранной тематике и других отраслях народного хозяйства при необходимости выполнения геометрической коррекции космических сканерных снимков линейной развертки, обладающих существенными геометрическими искажениями.. Цель изобретения - снижение стоимости работ за счет исключения прикладного программного обеспечения,На фиг. 1 представлено исходное изображение; на Фиг. 2 - векторограмма смещений точек снимка относительно идентичных точек карты; на фиг. 3 и 4 - контуры проекции Гаусса-Крюгера,ЛО 1553 5 А 1 2изобретения - повышение экономичности, С 1 оссб включает изготовлениеконтратипа исходного изображения,построение векторограммы смещенийтоцек снимка относительно идентичныхточек карты и усредненного векторасмещений, разделение контратипа начасти перпендикулярно усредненномувектору, расчет величины шага дискретизации для каждой из частейконтратипа. Затем контратип вводятв блок обработки путем установкиусредненного вектора смещений перпендикулярно строке считывания изображения, госле чего выводят преобразованное изображение., растягивая или сжимая каждую часть контратипа в соответствии с рассцитанным шагом записиизображения, 4 ил,и контуры геометрически откорректированного изображения.Способ осуществляют следующим об разом,Предварительно изготавливают на фотопленке исходное изображение (к тратип) для цифрового ввода его в ЭВМ.Строят векторограмму смещений точек снимка относительно идентицных тоцек карты.(2) определены также и с помощью карты. Для этого-переносят на две разные кальки опознаки со снимка и опознаки с карты. На кальках выбирают два опознака, располагающиеся на прямой примерно параллельно изображению трассы космического летательного аппарата (КЛА) и кальку с опознаками со снимка и ка.пьку с изображением опознаков, снятых с карты, редуцируют в один масштаб, т,е два выбранных опознака на карте и на снимке должны совпадать. Графицески определяют вектора смещений для остальных опознаков.По векторам смещений строят усредненный вектор смещений.Условно разбивают контра гип изображения на части перпендикулярно хсредненному вектору по формуле.Г(1 зх, - 1 ьнх) - сО где 0 - размер части, мм;г) - шаг дискретизации УВВИ, ммЬ минимальный шаг дискретизации УВВИ, мм,Лзаданная величина точностигеометрицеской коррекцииизображения в масштабе снимка, мм;1 - максимальный вектор смещенияданного участка мм;,номер обрабатываемого участкаДалее, используя возможность УВВИбарабанного типа в процессе вводавывода растягивать или сжимать изображение, рассчитывают требуемый шагдискретизации между центрами строкЧ для каждой части по формуле Г Я(1, - ,)Ч = еп -+ где +" при 1 о, - 1 и,11 11при 1; . - ,-. с О.Устанавливают контратип изображения на барабане УВВИ таким образом,чтобы усредненный вектор смещенийбыл перпендикулярен строкам сканирования на барабане УВВИ., и поочередно вводят все части контратипа вЭВМ и затем выводят из ЗВМ, растягИ. вая или сжимая каждую часть введенного изображения в соответствии срассчитанным шагом дискретизации цэг,П р и м е р, На фиг. 1 показано"исходное изображение космическогосканерного снимка линейной развертки,на котором представлен Крымский полуостров. Масштаб снимка вдоль трассы полета ИСЗ равен 1;2570090,Необходимо геометрически откорректировать это исходное изображениев масштабе снлмка в карту проекцииГаусса-Крюгера, с точностью +1 мм сшагом дискретизации Ц = 0,1 мм.Первой операцией геометрицескойкоррекции является изготовление контратипа, представляющего собой негатив лли позитив исходного изображения, показанного на фиг. 1.20Вторая операция " определение поля смещения точек снимка относительно идентичных точек касты,На фиг. 2 показано поле векторовсмещения, полученное с помощью карты,На фиг. 2 показаны контуры исходногоизооражения 1 и контуры выбраннойкартографицеской проекции 2, приведенные к масштабу исходного изображения по опорным точкам 3 и Й.По опорным точкам 5-7 определяютнаправление и размеры соответствующих векторов смещения 8-10,Осуществляют выбор преобладающего(генерального) направления необходимых смещений как среднее из направлений имеюшихся векторов смещения.На фиг. 2 это направление обозначенопозицией 11,,О Затем проводят условное разделение контратипа на части перпендикулярно преобладающему (генеральному)направлению 11 смещений,По фиг. 2 визуально определяют,,:.5 цто контуры исходного изображения 1на участке 12 размером 30 мм совпадают с выбранной картографическойпроекцией 2 с требуемой точностьюл= +1 мм и поэтому не требуют преобразования.Размеры вектора 9 превышают заданную точность преобразования и требуют геометрической коррекции,Рассцитывают размер части 13 дляопорной точки Ь, соответствующей5."вектору 9 смещений от предыдущегоучастка 12,Подставляют в выражение (1) следующие значения: э. = 1; 1, = 1 мм;= 0,0125 мм;= +1 мм.Получают, что размер участка должен находиться в интервале 16Р й( 32 мм.Рассчитывают размер части дляопорНой точки 7 по формуле (1) созначениями д = 1; 1 дх;-1 = 1 мм;С 3 = 0,100 мм; А = 0,0125 мм;= +1 мм.Получают, что 96Р й 112 мм.При таких расстояниях опорнаяточка 7 выходит за пределы контуракарты 2, т.е, необходима следующаячасть 14 обработки. Размеры первойобрабатываемой части 13 принимаютравными 30 мм, Рассчитывают размерчасти 14 для опорной точки 7 приследующих значениях по выражению (1):= ;.1 мм.Размер части 14 должен находиться в интервале 36Р б 44 мм. Принимают размер части 14 равным 44 мм.Всего исходное изображение (контратип) условно разделили,на три частиразмерами 30, 30 и 44 мм.Далее определяют шаг дискретиза"ции ц для каждой части по выражению(г) .Часть 12 не требует геометрической коррекции, так как удовлетворяетточности 6 = +1 мм. Поэтому Я х= Ц, мм =,0,100 мм.Для части 13 шаг дискретизациипреобразования находят по выражению(2) при значениях Я = 0,100 мм; й0,0125 мм; 1 ъах = 4 мм; 1 ах- 1 мм:, Р = 30 мм.Получают Ч= 0,1125 мм.Для участка 14 рассчитывают Опо формуле (2) при значениях= 0,100 мм; Е = 0,0125 мм; 1 ;- 14 мм; 1, = 4 мм; Р = 44 мм.Получают х = О, 1250 мм.Устанавливают контратип изображения на барабане УВВИ такЧм образом,чтобы строки УВВИ были перпендикулярны направлению, преобладающего смещения 11.Осуществляют ввод контратипа вЭВМ через УВВИ барабанного типа, например "фотомейшн".Части 12-14 .контратипа вводят поочередно в ЭВМ с шагом дискретизацииЧ = 0,100 мм,Вывод контратипа из ЭВМ производятчерез то же УВВИ. 6Части 12-14 контратипа выводят поочередно из ЭВМ с рассчитанными шагами дискретизации а часть 12 с Чх,- 0,100 мм; часть 13 с ц, = 0,1125 мм;часть 14 с сх = 0,1250 мм,В результате получают геометрически от корр екти рова нное и зображ ение,представленное на фиг. 3,На фиг, 4 приведены контуры проекции 2 карты Гаусса-Крюгера и контурыгеометрически откорректированногоизображения,Из фиг. 4 видно, что контуры геометрически откорректированного изображения совпадают с контурами картыс заданной точностью. формула изобретения Способ геометрической коррекции космических сканерных снимков линейной развертки, включающий изготовление на светочувствительном материале исходного изображения, построчное считывание и ввод его в блок обработки, обработку введенного изображения с учетом координат опорных точек карты и вывод из блока обработки преобразованного изображения, о т л и ч аю щ и й с я тем, что, с целью повышения экономичности, перед вводом изображения в блок обработки строят векторограмму смещений точек исходного изображения относительно опорных точек карты, определяют усредненный вектор смещений, по направлению которого разбивают исходное изображение на части, и рассчитывают для каждой части величину шага записи преобразованного изображения, при вводе в блок обработки предварительно ориентируют усредненный вектор смещений исходного изображения перпендикулярно строке считывания изображения, а при выводе преобразованного изображения производят последовательное изменение шага записи, при этом размер каждой части Р исходного изображения рассчитывают из выражениялгде- показатель точности геометрической коррекции изображения. в масштабе снимка, мм;- шаг записи изображения, мм;6 в .минимальный шаг записи изображения, определяемый конструктивными параметрамииспользованных средств ввода и вывода изображения, мм;гпе 1- максимальный вектор смещения "+" причасти изображения мм; "- при- номер обрабатываемой части, а величину шага записи части с х опре деляют из выражения1553825 ректор С.Шевкун еда Тираж 3 каз 450 одписн Государственного комитета по изобретениям и открытиям при ГКНТ СССР 113035,. Москва, Ж, Раушская наб., д. 4/5 Производственно-издательский комбинат "Патент", г, Ужгород, ул. Гагарина, 10 Составитель В.ЖивиА.Огар Техред Л.Олийнык Фи