Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки

На основе анализа основных соотношений полиномиального преобразования изображений космической съемки наземных объектов, используемого в процессе их топографической привязки, показана возможность оценок ее точности. Предложена методология соответствующих расчетов и моделирования, позволяющая ка...

Повний опис

Збережено в:
Бібліографічні деталі
Опубліковано в: :Штучний інтелект
Дата:2011
Автори: Белозерский, Л.А., Мурашко, Н.И.
Формат: Стаття
Мова:Російська
Опубліковано: Інститут проблем штучного інтелекту МОН України та НАН України 2011
Теми:
Онлайн доступ:https://nasplib.isofts.kiev.ua/handle/123456789/59849
Теги: Додати тег
Немає тегів, Будьте першим, хто поставить тег для цього запису!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Цитувати:Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки / Л.А. Белозерский, Н.И. Мурашко // Штучний інтелект. — 2011. — № 3. — С. 207-219. — Бібліогр.: 9 назв. — рос.

Репозитарії

Digital Library of Periodicals of National Academy of Sciences of Ukraine
_version_ 1859675476702789632
author Белозерский, Л.А.
Мурашко, Н.И.
author_facet Белозерский, Л.А.
Мурашко, Н.И.
citation_txt Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки / Л.А. Белозерский, Н.И. Мурашко // Штучний інтелект. — 2011. — № 3. — С. 207-219. — Бібліогр.: 9 назв. — рос.
collection DSpace DC
container_title Штучний інтелект
description На основе анализа основных соотношений полиномиального преобразования изображений космической съемки наземных объектов, используемого в процессе их топографической привязки, показана возможность оценок ее точности. Предложена методология соответствующих расчетов и моделирования, позволяющая как выполнять расчеты в конкретных задачах, так и вскрывать пути снижения ошибок, а также получить потенциальные характеристики точности, сняв имеющие место ограничения в распространении достигнутых результатов на любой практически возможный случай. На основі аналізу основних співвідношень поліноміального перетворення зображень космічної зйомки наземних об’єктів, використовуваних в процесі їх топографічної прив’язки, показана можливість оцінок її точності. Запропонована методологія відповідних розрахунків і моделювання, що дозволяє як виконувати розрахунки в конкретних завданнях, так і розкривати шляхи зниження помилок, а також отримувати потенційні характеристики точності, знявши обмеження, що мають місце, в поширенні досягнутих результатів на будь-який практично можливий випадок. On the basis of analysis of the basic proportion of polynomial transformation of images of space survey of the land objects, which is used during the process of the topographic reference of these objects, possibility of accuracy estimations is shown. The methodology of special calculations and modeling for practical calculation, for error reduction and also for getting potential characteristics of accuracy is offered. Thus the restriction in dissemination of the achieved results in any possible situation is removed.
first_indexed 2025-11-30T15:52:24Z
format Article
fulltext «Штучний інтелект» 3’2011 207 4Б УДК 004.93’1; 004.932; 528.85/.87(15) Л.А. Белозерский, Н.И. Мурашко Объединенный институт проблем информатики НАН Беларуси, г. Минск blzrsk@newman.bas-net.by, mni@newman.bas-net.by Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки На основе анализа основных соотношений полиномиального преобразования изображений космической съемки наземных объектов, используемого в процессе их топографической привязки, показана возможность оценок ее точности. Предложена методология соответствующих расчетов и моделирования, позволяющая как выполнять расчеты в конкретных задачах, так и вскрывать пути снижения ошибок, а также получить потенциальные характеристики точности, сняв имеющие место ограничения в распространении достигну- тых результатов на любой практически возможный случай. Введение Тематическое использование изображений космической съемки, как правило, имеет целью выяснение положений тех или иных природных или антропогенных образований на поверхности Земли. Для каждого из таких применений допустимые отклонения соответствующих координат от истинных значений могут существенно отличаться. При этом очевидно, что существует целый ряд специфических областей использования результатов дистанционного зондирования Земли, где востребованы наиболее высокие точности знания топографического положения объекта. К ним, на- пример, могут быть отнесены задачи предотвращения или устранения последствий чрезвычайных ситуаций, практика так называемого «точного земледелия», а также такие наиболее критичные к ошибкам топопривязки проблемы, как сопоставление изображений разновременной съемки. В целом возникновение ошибок привязки изображений спутниковой съемки может быть обусловлено: – особенностями применяемой аппаратуры съемки земной поверхности; – способом устранения отличий в пропорциях изображения дистанционной съемки и соответствующего участка земной поверхности. Ограничения здесь связаны с невозможностью безошибочных знаний параметров внешнего и внутреннего ориентирования аппаратуры в момент съемки. В то же время моделирование полета носителя этой аппаратуры (космического аппарата) и условий съемки с некоторым объемом сопутствующих измерений дает возможность получать координаты привязки с достаточно низкими ошибками [1]. Согласно имеющимся упоминаниям [2-5] «минимальные уровни обработки» кос- мических снимков (например, Ikonos Pro и QuickBird Orthorectified), предлагаемые ком- паниями-поставщиками, предполагают достижение среднеквадратической ошибки (СКО) их привязки не хуже 5 – 7 м. Однако пользователю придется заплатить за такие снимки цену, превышающую стоимость «сырых» изображений в 2 – 3 раза. Если же вести речь Белозерский Л.А., Мурашко Н.И. «Искусственный интеллект» 3’2011 208 4Б о снимках, имеющих СКО привязки до 2 – 4 м (уровни обработки Precision, Precision stereo), то соответствующее увеличение достигает уже 6 раз. Причем в ряде случаев для обеспечения указанных точностей ортотрансформированного изображения компанией- поставщиком от пользователя требуется предварительное предоставление реперных (опорных) точек, что наряду с высокой стоимостью еще и существенно удлиняет сроки получения материалов съемки [3]. Поэтому рассмотренные особенности позволяют считать [3], что пользователи могут заказать изображения первого уровня (в 5,5 раза дешевле продукта Precision), а дальнейшую обработку произвести самостоятельно. При этом наиболее распростра- ненным в тематической обработке изображений космической съемки является при- менение полиномиальной геометрической коррекции. Возможности ее достаточно широки в выборе параметров осуществления, что в конечном итоге отражается на степени повышения интересующей точности плоского геометрического сведения со- поставляемых изображений. В то же время, приступая к такому виду обработки, исполнитель должен иметь не просто представления о характере зависимости ошибок привязки изображений косми- ческой съемки от тех или иных факторов, а знания и реализуемые возможности оценок и контроля точностей пространственного преобразования, достигаемых в процессе вы- полнения этого этапа геометрической коррекции. На создание соответствующих условий и ориентирована настоящая статья. 1 Аналитические возможности оценки ошибок полиномиальной геометрической коррекции В основе полиномиальной геометрической коррекции лежит очевидное преобразо- вание растровых координат изображения в плоскую прямоугольную систему топографи- ческих координат: , ~~ ; ~~ vW vQ   T T T T y x (1) где T v 2233221 pppppppppppp yxyxyxyxyxyx  – вектор, построенный на растровых координатах пикселя изображения, пересчет которых в топографиче- ские координаты указанным способом осуществляется (зависит от порядка регрессии); pp yx , – растровые координаты пикселя, топографические координаты ко- торого TT yx ~,~ являются конечной целью преобразования; T Q 110 ~...~~~  nqqq ; T W 110 ~...~~~  nwww – статистические (в общем случае) оценки коэффициентов полиномиального преобразования. Следует сразу обратить внимание, что компоненты векторов v представляют собой детерминированные значения. В то же время оценки векторов преобразований WQ ~ , ~ , в основе определения которых лежат регрессионные зависимости между раст- ровыми и топографическими координатами реперных точек, определенных на местности с ошибками, являются случайными величинами. Следствие этого – ошибки оценок прямоугольных координат каждого из пикселей изображения TT yx ~,~ . Отсюда согласно (1) дисперсия возникающих ошибок для координаты Tx~ (для координаты Ty~ аналогична) может быть определена так:            vKvvQQvvQQvvQvQ  TTTTTTTT T MMMxD ~~~~~~~ , (2) Закономерности поведения оценок точности полиномиального преобразования… «Штучний інтелект» 3’2011 209 4Б где Q ~  – означает разность случайной величины вектора преобразования Q ~ и его математического ожидания (центрированный вектор значений); K – ковариационная матрица коэффициентов полиномиального преобразо- вания. Присутствующая здесь (2) ковариационная матрица полиномиального преобразова- ния согласно [6], [7] выражается следующей зависимостью: 1 CK  2 PT . (3) Здесь PT – среднеквадратический разброс ошибок (СКО) прямоугольных ко- ординат привязки реперных точек в топографической системе; При этом сама матрица C в данной задаче строится на основе использовании матрицы предикторных переменных V , представляемой аргументами полиномиального преобразования: ,VVС  Т (4) где T nmpppp   T m T 3 T 2 T 1 vvvvV  ; T i v 2233221 ppippippippippippippippippippippippip yxyxyxyxyxyx  – компо- нента вектора, построенная на полиномиально представленных предикторах растро- вых координат i-й реперной точки ( _____ ,1 mi  ). В итоге, с учетом (3), дисперсия координат любого пикселя преобразуемого изображения (2) принимает окончательную запись:   vCv 12~  T PTTxD  . (5) Таким образом, последнее выражение позволяет уже более конкретно судить о том, что дисперсии оценок ошибок топопривязки координат геометрически преобразо- ванного изображения  TxD ~ ,  TyD ~ зависят явно от дисперсии реперных точек ( 2 PT )1 и неявно – от порядка полиномиального преобразования, соответствующего степени ис- кажений изображения при съемке (через растровые координаты расставленных реперных точек); – от числа реперных точек, формирующих матрицу предикторных переменных V и через нее – матрицу C ; – от значений растровых координат, определяющих матрицу C , а следовательно, от расстановки реперных точек. Отсюда очевидна возможность получения соответствующих количественных данных при любых сочетаниях параметров явных и неявных связей. 2 Представления о системе исходных данных, действиях, моделировании и решаемых задачах Создание представлений о величинах ошибок и определяющих их факторах обеспечивает построение системы задания данных, взаимосвязанных расчетов и мо- делирования, основанное на варьировании условий и параметров реализации выра- жения (5). Структурная схема такой системы может быть представлена, как показано на рис. 1. 1 Здесь без ущерба для общности считаем 2 PT = 2 xPT = 2 yPT . Белозерский Л.А., Мурашко Н.И. «Искусственный интеллект» 3’2011 210 4Б Рисунок 1 – Структурная схема системы задания, расчетов и моделирования в оценке ошибок топопривязки пространственно преобразованного изображения Рассмотренная структурная схема определяет полный перечень исходных данных в комбинации с взаимосвязанными действиями, обеспечивающими оценку точности привязки любой точки реального изображения космической съемки. Получаемые здесь ошибки – это результат полиномиальной пространственного преобразования для за- данной расстановки реперных (опорных) точек, их числа и точностей топопривязки. Причем, располагая точно известными характеристиками перечисленных данных, модель оказывается полностью идентифицированной для достижения указанной цели. Соответ- ствие натуре при этом обеспечивает, как использование одних и тех же данных не- посредственно для геометрической коррекции и ее моделирования, так и применение единой для них аналитической основы в преобразованиях координат и в определении ошибок. Естественно, что указанное соответствие (в том числе и при статистическом задании ошибок) сохраняется в случаях необходимости решения задач не только кон- кретного расчета, но и определения путей улучшения результата и поиска закономер- ностей при изменении: – размеров изображения; – пространственной разрешающей способности съемки; – растровых координат расстановки реперных точек; – точностных характеристик топографического позиционирования реперных точек; – случайных ошибок расстановки реперных точек при их интерактивном раз- мещении на изображении относительно мест точного положения на местности. То есть, изменяя в интересующих пределах перечисленные параметры, можно оценивать реакцию на них ошибок привязки геометрически корректируемого изобра- жения космической съемки во всем факторном пространстве внешних воздействий, определяемых вышеприведенным перечнем. Остается понять, насколько индивидуаль- ны получаемые реакции в каждой конкретной комбинации приведенных параметров, что и позволяют выполнить имитационное моделирование любой практически име- ющей место ситуации согласно приведенной схеме (рис. 1). Задание (выбор) числа реперных точек PTn Расстановка реперных точек (определение растровых координат при расстановке) Задание размеров изображения и пространственной разрешающей способности съемки Формирование, расчет и обращение матрицы C Расчет СКО ошибок  TxD ~ ,  TyD ~ Задание вектора растровых координат v преобразуемого изображения (выбор из числа возможных) Задание СКО топопривязки реперных точек PT Графическое оформление заданий и расчетов Формирование ошибок расстановки реперных точек (датчик случайных чисел) Закономерности поведения оценок точности полиномиального преобразования… «Штучний інтелект» 3’2011 211 4Б 3 Расчет параметров, его ограничения и общие закономерности Оценка ошибок геометрической коррекции изображения космической съемки тре- бует некоторой начальной определенности параметров. Для конкретизации последующего анализа поведения определяемых характеристик используем следующий пример: – размеры изображения 20002000 HG NN пикселей; – пространственная разрешающая способность съемки 5 м; – расстановка реперных точек – геометрически равномерная; – ошибки расстановки отсутствуют; – порядок регрессии топографических координат по предикторам растровых координат. На рис. 2 в условиях указанных ограничений приведено семейство зависимостей СКО ( )~( TXoo xD или )~( TYoo yD ) топографической привязки северо-западного угла геометрически откорректированного изображения 0000 ,YX для первого и второго порядков регрессии. При этом количественно отличающиеся наборы реперных точек ( ,10,5PTn 30,25,20,15 ) и среднеквадратические точности позиционирования каждой из них (1, 2 и 5 м) выбраны произвольно1. 0 1 2 3 4 5 1 2 3 4 5 6 Число реперных точек/5 С К О т оп оп ри вя зк и из об ра ж ен ия [ м ] СКО реперной точки = 1 м СКО реперной точки = 2 м СКО реперной точки =5 м 0 1 2 3 4 5 1 2 3 4 5 Число реперных точек/5 - 1 [м] С К О т оп оп ри вя зк и из об ра ж ен ия [ м ] СКО реперной точки = 1 м СКО реперной точки = 2 м СКО реперной точки = 5 м а) б) Рисунок 2 – СКО топопривязки преобразованного изображения для отличающегося числа и точностей определения координат реперных точек: а) для линейной регрессии, б) для квадратичной регрессии Приведенные результаты, прежде всего, демонстрируют, что для линейной ре- грессии среднеквадратический разброс рассматриваемой ошибки снижается практически вдвое в принятом диапазоне изменений числа реперных точек при любой из точностей их топографического позиционирования. В то же время для регрессии второго порядка такого влияния числа реперных точек на ошибку привязки изображения не наблюдается при более высоком ее относительном уровне. На первый взгляд, полученные оценки кажутся ограниченными, характерными, прежде всего, для приведенного случая в силу единственности размера изображения, заданного для моделирования. 1 Произвольный характер выбора объема первоначальных оценок обусловливался необходимостью наблюдения имеющих место закономерностей при всегда имеющей здесь место возможности дополнить состав реперных точек и задать интересующие дополнительно ошибки их топопривязки на местности . Белозерский Л.А., Мурашко Н.И. «Искусственный интеллект» 3’2011 212 4Б Однако вопреки этому, повторяя расчеты, но при увеличении изображения (с со- ответствующим изменением координат реперных точек), например, в 10 раз (то есть до 20000 × 20000 пикселей), достаточно просто обнаружить, что точности привязки остаются теми же, что и приведены на рис. 2. Причем к такому же результату приходим и с уменьшением изображения. Одновременно с равным успехом к нему же ведут отли- чающиеся изменения размеров его по каждой из координат. Объяснение получаемого результата следует из пропорциональности определи- теля (D) матрицы C (4) и ее алгебраического дополнения 11A размеру изображения. В результате компонента обращенной матрицы C не изменяет своего значения: constDAс  /11 1 11 . При этом вектор растровых координат точки привязки изображения (северо-за- падный угол изображения) выглядит для линейной и квадратичной регрессии доста- точно просто: T 001v и T 00001v . Соответственно среднеквадратическая ошибка топопривязки согласно (5) ока- жется одинаковой для любого увеличения изображения и равной: 11 1 00 CPT T PTX    vCv . Наряду с установленной особенностью преобразования координат легко убедить- ся, что увеличение или уменьшение пространственной разрешающей способности съемки при постоянном размере изображения не оказывает никакого влияния на относительное расположение расставленных реперных точек, а значит, и на проведенные оценки. Эти изменения в анализируемом преобразовании оказываются равносильны изменениям раз- меров изображения. В результате можно считать, что первое и второе из введенных начальных огра- ничений моделирования, связанные с выбором геометрических параметров изображения и разрешающей способности, могут быть сняты, а полученные зависимости точностей привязки пространственно преобразуемого изображения (рис. 2) следует относить к не- изменным характеристикам, потенциально достигаемым (оценка снизу) для любого по размерам изображения, для любой разрешающей способности. Пропорционально увели- чивая или уменьшая изображение (до реально полученных размеров в космической съемке), в любом практически имеющем место случае заданных величин1 PT при со- ответствующем изменении координат идеально расставленных реперных точек приходим к тем же зависимостям потенциальной точности привязки. Несмотря на такие преимущества, рассмотренные характеристики все-таки ка- саются оценки только нижней границы ошибок. Само появление приведенного выше определения «потенциальная точность» обусловлено здесь тем, что «чистая» равно- мерность расстановки реперных точек может быть обеспечена только при моделиро- вании. Практическое достижение такого же качества подготовки геометрической кор- рекции реальных изображений космической съемки оказывается невозможным. Здесь при стремлении к желаемой равномерности начинает играть роль сложность выбора достаточного числа характерных мест расстановки реперных точек, не говоря уже о невозможности топопривязки их в каждом из заранее выбранных мест. 1 В силу незначительности принятого здесь для моделирования диапазона изменений СКО топопривязки реперных точек возможно также ошибочное отнесение результатов, полученных на рис. 2, к трудностям рас- пространения их на любой практически имеющий место случай противоречит выражению (5). Пропорционально увеличивая или уменьшая приведенные значений )~( TXoo xD в число раз, в которое отличается реально интересующее значение PT от моделируемых, легко удовлетворить любые потребности. Закономерности поведения оценок точности полиномиального преобразования… «Штучний інтелект» 3’2011 213 4Б Являясь ориентиром при организации обработки реальных изображений косми- ческой съемки и пределом в их привязке, потенциальная точность не исключает оценку реально достигаемых значений ошибок. Для практически имеющих место расстановок реперных точек ничто не препятствует их определению на основе использования той же, уже рассмотренной методологии. Таким образом, сняв предшествующие ограничения, пока остается необходимым считаться с очередным из перечисленных в начале раздела ограничений: моделирование выполнялось в предположении, что ошибки расстановки реперных точек отсутствуют. Насколько серьезно такое допущение, позволяет установить теперь уже статисти- ческое моделирование, предусмотренное согласно структурной схеме рис. 1. Для этого достаточно задействовать датчик случайной расстановки реперных точек оператором. Причем в качестве закона распределения соответствующих ошибок оператора можно принять равномерное распределение в диапазоне его способностей. Выбирая для исследо- ваний два крайне возможных диапазона этих ошибок ( 100 пикселей – недостаточно внимательный оператор и 3 пикселя – аккуратный оператор), удается проследить ин- тересующие закономерности. Согласно изложенным условиям на рис. 3 приведена статистически получен- ная выборка среднеквадратических ошибок (СКО) привязки изображения космиче- ской съемки при расстановке 10 реперных точек с точностью знания их положения на местности мPT 1 . 0,663 0,665 0,667 0,669 0,671 1 3 5 7 9 11 13 15 17 19 21 23 25 Номер эксперимента С КО п ри вя зк и из об ра ж ен ия [м ] 0,665 0,6655 0,666 0,6665 0,667 1 3 5 7 9 11 13 15 17 19 21 23 25 Номер эксперимента С КО п ри вя зк и из об ра ж ен ия [м ] а) б) Рисунок 3 – Реализации СКО привязки изображения при растровой расстановке 10 реперных точек с ошибками: а) 100 пикселей; б) 3 пикселя) Качественно оценивая один и другой результат, можно отметить, что относительно ошибки привязки изображения 0,667, достигаемой при точной расстановке 10 реперных точек с мPT 1 (рис. 2а), разброс, обусловливаемый неточностями операторов в их пра- ктической работе по ручной расстановке, незначителен. В табл. 1 для тех же исходных данных по разрешающей способности и точности топографического позиционирования реперных точек сведены статистические оценки привязки геометрически преобразованного изображения в диапазонах ошибок расстанов- ки, более приемлемых для практики. Как можно увидеть по приведенным данным (рис. 2, табл. 1), неточность растро- вой расстановки заданного числа реперных точек с ростом ошибок оператора растет не существенно. Отклонения от математического ожидания при возрастании диапазона ошибок расстановки от 3 до 30 пикселей составляют 0,2% – 2%. Белозерский Л.А., Мурашко Н.И. «Искусственный интеллект» 3’2011 214 4Б Таблица 1 Диапазон ошибок расстановки реперных точек [пикселей] Математическое ожидание СКО привязки [м] Среднеквадратический разброс СКО привязки [м] Интервал практически возможных значений точности привязки (по выборке) [м] 3030  0,667 0,0044 675,0656,0  2020  0,667 0,0028 672,0661,0  1010  0,667 0,0013 668,0663,0  33  0,667 0,0004 667,0665,0  Здесь же следует обратить внимание на то, что, как показывает моделирование, изменения разрешающей способности для того же диапазона операторских ошибок расстановки реперных точек не сказывается на ошибках топопривязки изображения космической съемки. Таким образом, и начальное ограничение в рассмотренном моделировании, со- стоящее в отказе от учета ошибок расстановки реперных точек, оказывается допустимым, если не предполагать их необоснованно большими. Например, имея дело со случаем, когда ширина проезжей части пересекающихся дорог составляет 8 м, при желании уста- новить реперную точку в центре их перекрестка, хорошо отображаемого космическим снимком с пространственной разрешающей способностью 1 м, ошибка установки практи- чески никогда не превысит 2 – 3 пикселя. Но даже необоснованно выставляя эту точку на краю одной из дорог, хотя и в районе пересечения вместо естественного стремления к центру, придется иметь дело с максимально возможной ошибкой всего в 4 пикселя. При равной тщательности расстановки других реперных точек это не приведет к ухудшению ошибки привязки геометрически скорректированного изображения (рис. 2) более, чем на 0,2%. Аналогично обстоит дело и при ухудшении разрешающей способности. Если в отличие от предшествующего случая интересующее пересечение дорог выглядит всего лишь как пересечение линий, то в качестве места установки реперной точки не может быть принят другой пиксель, кроме того, в пределы которого это пересечение попадает. 4 Точность геометрической коррекции и степень равномерности расстановки реперных точек Возвращаясь к вопросу равномерности расстановки реперных точек, придется под- черкнуть, что в качестве действий по преодолению ограниченности моделирования в соответствующем случае предлагалось выполнение отдельного расчета согласно рас- смотренной модели для реально имеющих место условий геометрической коррекции. При этом, как следует из рис. 2, точность топографической привязки изображений может быть повышена увеличением числа реперных точек и снижением ошибок определения их положений на местности. Такой путь, как можно было понять, часто более сложен, чем стремление к равномерной расстановке. Здесь приходится обратить внимание, что рассмотренная методология подготовки данных, расчетов и моделирования наряду с перечисленными способами приближения к желаемой точности геометрической коррекции открывает одновременно и возмож- Закономерности поведения оценок точности полиномиального преобразования… «Штучний інтелект» 3’2011 215 4Б ности для изменения расстановки реперных точек. К открытию ее приводят попытки варьирования расположением реперных точек в пределах изображения, а соответст- венно на отображаемой части земной поверхности, с оценками точности линейной регрессии. В процессе исследования различных вариантов моделирования может оказаться, что и реально неравномерное расположение реперных точек в некоторых случаях обеспе- чивает ту же потенциальную точность привязки северо-западного угла изображения. В частности этому соответствует сосредоточение реперных точек в непосредственной близости от указанного угла, в то время как в остальной части изображения они вообще не расставлялись. То есть графики зависимостей ),( 00 РТРТx n в таком случае остаются теми же, что и приведены на рис. 2. Для расширения представлений о таких ситуациях следует не ограничивать инте- ресы оценки точности топографической привязки единственной рассмотренной точкой изображения ( 0000 ,YX ), играющей важную роль в его применении, а распространить расчет на все пространство изображения. В результате можно обнаружить закономерные изменения точностей привязки изображения космической съемки по всей его площади. Так, на рис. 4 приведены ошибки топографической привязки сетки пикселей изображе- ния, полученные моделированием линейной регрессии для 10 равномерно расставленных реперных точек и мPT 1 . Получаемая поверхность ошибок в физической аналогии представляет собой «платок, подвешенный закреплением его на равных высотах за четыре угла». Причем вертикально вниз от каждой точки закрепления располагаются углы рассматриваемого изображения, занимающего горизонтальную часть рисунка. Здесь можно увидеть, что максимальная точность привязки, обеспечиваемая гео- метрической коррекцией, достигается в точке максимального провисания «платка», сов- падающей с центром изображения (0,33 м), в то время как значение 0,67 м, приведенное на рис. 2 для одного указанного выше угла, характерно и для всех его угловых положений («высоты закрепления»). Причем, увеличивая число реперных точек, принимающих участие в такой оцен- ке, придется увидеть, что представленная на рис. 4 поверхность перемещается вниз по оси Z (уменьшаются высоты точек закрепления «платка» и точки его максимального провисания, соответствующие ошибкам топопривязки). Тем самым отражается повышение точностей для каждого из угловых положений изображения и, в частности, центра его. 1 2 3 4 5 6 Р1 Р3 Р5 0,25 0,3 0,35 0,4 0,45 0,5 0,55 0,6 0,65 0,7 СКО топопривязки пикселей изображения [м] Координаты точек по оси Х изображения =(x-1)*2000[м] Координаты точек по оси Y изображения Py=y*2000 [м] Рисунок 4 – Среднеквадратические ошибки топопривязки сетки пикселей изображения при равномерной расстановке 10 реперных точек ( 1PT м) Белозерский Л.А., Мурашко Н.И. «Искусственный интеллект» 3’2011 216 4Б Закономерность такого поведения оценок на сетке пикселей легко проверяется теоретически, если попытаться определить экстремум среднеквадратической ошибки топографических координат как функции растровых координат, представляемой вы- ражениями (2) или (4), в плоскости изображения:     .0~;0~       T P T P xD y xD x Тогда с учетом особенностей обращенных матриц C для первого порядка ре- грессии в соответствии с приведенными экспериментальными данными получим, что экстремум соответствует минимуму СКО и достигается в точке с растровыми ко- ординатами, определяемыми через коэффициенты корреляции, как составляющие ко- вариационной матрицы K (2): 22 12min )~(minarg K K xDx T x p p  ; 33 13min )~(minarg K K xDy T y p p  . При конкретных значениях коэффициентов корреляции линейной регрессии, получаемых в процессе моделирования, и равномерной расстановке реперных точек используемого изображения 10000 × 10000 м получим: .5000;5000 minmin мyмx pp  Теперь можно воспользоваться моделированием для объяснения упомянутой выше ситуации сохранения потенциальной точности Xoo , Yoo при группировании реперных точек в пределах некоторой части изображения, примыкающей к северо-западному углу его. В частности, для этого при моделировании достаточно уменьшить, например, в 2 раза значения заданных координат реперных точек (равномерного размещения) при сохране- нии размеров моделируемого изображения. В результате такого действия область расстановки реперных точек окажется за- нимающей левую верхнюю часть пространства изображения, составляя четверть его, в пределах которой те же 10 реперных точек размещаются по-прежнему равномерно. Соответственно этому на рис. 5 приведены ошибки топографической привязки сетки пикселей всего изображения. 1 2 3 4 5 6 Р1 Р3 Р5 0 0,2 0,4 0,6 0,8 1 1,2 1,4 1,6 1,8 СКО топопривязки пикселей изображения [м] Координаты по оси X изображения =(х-1)*2000 [м] Координаты по оси Y изображения Py=y*2000 [м] Рисунок 5 – Среднеквадратические ошибки топопривязки сетки пикселей при сосредоточенной равномерной расстановке 10 реперных точек в левой верхней четверти изображения ( 1PT м) Закономерности поведения оценок точности полиномиального преобразования… «Штучний інтелект» 3’2011 217 4Б Приведенные результаты в отличие от предшествующего случая позволяют обнару- жить, что «платок» оказывается подвешенным за углы на отличающихся высотах. То есть – среднеквадратическая ошибка топопривязки северо-западного угла изображения (левый передний угол поверхности на рис. 5) сохранила свое значение ( Xoo = 0,67 м), соответствующее равномерной расстановке реперных точек в пределах всего изобра- жения (рис. 4); – среднеквадратические ошибки других угловых положений изображения су- щественно возросли и отличаются между собой; – минимальная ошибка топопривязки сохранила свое значение (0,33 м); – местоположению минимальной ошибки соответствуют растровые координаты изображения .2500;2500 minmin мyмx pp  Более тщательное изучение рассматриваемого результата дает возможность увидеть: – минимальная ошибка топопривязки (0,33 м) соответствует не центру всего изо- бражения, а центру области сосредоточения равномерно распределенных реперных точек (5000 × 5000 м); – на углах полученной здесь области сосредоточения реперных точек ошибки топографической привязки пикселей равны между собой и их СКО составляет 0,67 м, как это было для всего изображения при равномерной расстановке реперных точек в его пределах. На основе полученных результатов придется вспомнить, что полиномиальная ап- проксимация зависимости между растровыми и топографическими координатами для регрессии первого порядка строится на данных области расстановки реперных точек. Полиномиальное определение координат каждой точки изображения по восстановленной указанным образом регрессии представляет собой интерполяцию. При ее выполнении достигается минимум ошибки в центре этой области и равные значения в равноудален- ных и симметричных относительно центра точках. Чем более удалены такие точки, тем больше СКО их привязки. Центр, определяемый по минимуму ошибки интерполяции, всегда соответствует центру области равномерной расстановки. В случае, если область равномерной расстановки реперных точек по занимаемой площади меньше площади изображения, то для определения ошибок привязки точек изображения, лежащих вне указанной области, полиномиальная интерполяция пере- ходит в экстраполяцию с естественным ростом ошибок. Таким образом, с какой бы расстановкой реперных точек ни пришлось иметь дело, определение соответствующего центра (по минимуму ошибки привязки сетки пикселей изображения) и его смещения относительно истинного центра интересующего изобра- жения для регрессии первого порядка позволяет: – установить степень равномерности имеющей место расстановки в пределах всего изображения; – дополнить имеющийся состав реперных точек для достижения большей равно- мерности; – сопоставляя ошибки привязки углов изображения с потенциально возможными, определить степень удаления области расстановки от соответствующих углов изобра- жения. Заключение Предпринятый в статье поиск закономерностей поведения ошибок полино- миальной коррекции изображения космической съемки основывался на расчетах и моделировании таких оценок во всем факторном пространстве изменений соответ- Белозерский Л.А., Мурашко Н.И. «Искусственный интеллект» 3’2011 218 4Б ствующих внешних воздействий. Предложенная методика оценки ошибок полиномиаль- ного преобразования координат открыла возможности углубленного анализа причин их возникновения, наглядно и количественно демонстрируя их связи с такими факто- рами, как – число реперных точек, участвующих в восстановлении регрессионной зави- симости между растровой и топографической системами координат; – точности привязки координат реперных точек в топографической системе; – порядок регрессии связи топографических координат точек изображения с их растровыми координатами; – степень равномерности распределения реперных точек на изображении, соответ- ствующем отображаемому участку земной поверхности. В процессе исследований, результаты которых приводятся в статье, для каждого порядка регрессии установлена независимость потенциальной точности топографи- ческой привязки координат изображений космической съемки от – размеров изображения; – разрешающей способности космической съемки; – ошибок растровой расстановки реперных точек на изображении космической съемки. На основе исследования результатов моделирования предлагается методология, позволяющая пользователю в процессе тематической обработки данных космической съемки: – располагать заранее детерминированными графическими или табличными дан- ными для оценки возможного нижнего предела ошибок (потенциальная точность) при- вязки полиномиально преобразуемого изображения применительно к его параметрам; – получать оценки ошибок топопривязки геометрически трансформированных изображений космической съемки для реального числа, заданной расстановки репер- ных точек и знаний ошибок их положений на местности; – определять требуемые ошибки позиционирования на местности и количество используемых реперных точек при необходимости повышения точности привязки изображений космической съемки; – оценивать равномерность расстановки реперных точек и при необходимости повышения точности привязки изображения корректировать ее на основе оценки ошибок топографической привязки сетки пикселей изображения в целом. Литература 1. Болсуновский М.А. Геометрическая коррекция данных со спутника QUICKBIRD / М.А. Болсуновский // Компания «СОВЗОНД» 2006-03-12 [Электронный ресурс]. – Режим доступа : http://www.sovzond.ru/dzz/publications/542/3072.html. 2. Фрейзер К. Перспективы использования данных со спутников IKONOS и QuickBird для создания картографической продукции / К. Фрейзер// Австралия. Университет Мельбурна [Электронный ресурс]. – Режим доступа:http://www. sovzond.net/ articles/ article cartog.html. 3. Бехдиниан Б. Создание ортоисправленных изображений со спутника IKONOS на основе ЦМР со спутника SPOT. (Перевод компании «СОВЗОНД») / Б. Бехдиниан // Иранский центр дистанционного зондирования Земли [Электронный ресурс]. – Режим доступа : http://www.sovzond.net/articles/article_orto.html. 4. Болсуновский М.А. Перспективные направления развития дистанционного зондирования Земли из космоса / М.А.Болсуновский // Компания «СОВЗОНД» 2010-28-12 [Электронный ресурс]. – Режим доступа : http://www. sovzond. ru/dzz/publications/540/3239.html. 5. Лавров В.В. Космические съемочные системы сверхвысокого разрешения / В.В. Лавров // ГИА «Инно- Тер» [Электронный ресурс]. – Режим доступа : http://www.innoter.com/articles/imaging_systems. 6. Линник Ю.В. Метод наименьших квадратов и основы теории обработки наблюдений / Ю.В. Линник. – М. : ГИФМЛ, 1962. – 350 с. Закономерности поведения оценок точности полиномиального преобразования… «Штучний інтелект» 3’2011 219 4Б 7. Дрейпер Норманн Р. Прикладной регрессионный анализ / Н.Р. Дрейпер, Г. Смит ; пер. с англ. – [3-е изд.]. – М. : Издательский дом «Вильямс», 2007. – 912 с. 8. Гонсалес Р. Цифровая обработка изображений / Р. Гонсалес, Р. Вудс ; [пер. с англ. ; под ред. П.А. Чочиа]. – М. : ТЕХНОСФЕРА, 2006. – 1072 с. 9. Белозерский Л.А. Особенности полиномиальной геометрической коррекции применительно к задачам анализа изображений разновременной космической съемки / Л.А. Белозерский, Н.И. Мурашко, Д.С. Сущеня // Искусственный интеллект.– 2010. – № 3. – С. 299-311. Literatura 1. Bolsunovskij M.A. Geometricheskaja korrekcija dannyh so sputnika QUICKBIRD. http://www.sovzond.ru/dzz/publications/542/3072.html 2. Frejzer K. Perspektivy ispol'zovanija dannyh so sputnikov IKONOS i QuickBird dlja sozdanija kartograficheskoj produkcii. http://www. sovzond.net/ articles/ article cartog.html 3. Behdinian B. Sozdanie ortoispravlennyh izobrazhenij so sputnika IKONOS na osnove CMR so sputnika SPOT. http://www.sovzond.net/articles/article_orto.html 4. Bolsunovskij M.A. Perspektivnye napravlenija razvitija distancionnogo zondirovanija Zemli iz kosmosa. http://www. sovzond. ru/dzz/publications/540/3239.html 5. Lavrov V.V. Kosmicheskie s"emochnye sistemy sverhvysokogo razreshenija. http://www.innoter.com/articles/imaging_systems 6. Gonsales R. Cifrovaja obrabotka izobrazhenij. M. : TEHNOSFERA. 2006. 1072 s. 7. Belozerskij L.A. Iskusstvennyj intellect. № 3. 2010. S. 299-311. 8. Linnik Ju.V. Metod naimen'shih kvadratov i osnovy teorii obrabotki nabljudenij. M.: GIFML. 1962. 350 s. 9. Drejper N.R. Prikladnoj regressionnyj analizM.: Izdatel'skij dom “Vil'jams”. 2007. 912 s. Л.А. Білозерський, М.І. Мурашко Закономірність поведінки оцінок точності поліноміального перетворення геометричної корекції зображень космічної зйомки На основі аналізу основних співвідношень поліноміального перетворення зображень космічної зйомки наземних об’єктів, використовуваних в процесі їх топографічної прив’язки, показана можливість оцінок її точності. Запропонована методологія відповідних розрахунків і моделювання, що дозволяє як виконувати розрахунки в конкретних завданнях, так і розкривати шляхи зниження помилок, а також отримувати потен- ційні характеристики точності, знявши обмеження, що мають місце, в поширенні досягнутих результатів на будь-який практично можливий випадок. L.A. Belozerskij, N.I. Murashko Regularity of Behavior of Accuracy Estimation at Polynomial Transformation of Image Space Survey Geographic Correlation On the basis of analysis of the basic proportion of polynomial transformation of images of space survey of the land objects, which is used during the process of the topographic reference of these objects, possibility of accuracy estimations is shown. The methodology of special calculations and modeling for practical calculation, for error reduction and also for getting potential characteristics of accuracy is offered. Thus the restriction in dissemination of the achieved results in any possible situation is removed. Статья поступила в редакцию 31.05.2011.
id nasplib_isofts_kiev_ua-123456789-59849
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
issn 1561-5359
language Russian
last_indexed 2025-11-30T15:52:24Z
publishDate 2011
publisher Інститут проблем штучного інтелекту МОН України та НАН України
record_format dspace
spelling Белозерский, Л.А.
Мурашко, Н.И.
2014-04-10T12:26:54Z
2014-04-10T12:26:54Z
2011
Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки / Л.А. Белозерский, Н.И. Мурашко // Штучний інтелект. — 2011. — № 3. — С. 207-219. — Бібліогр.: 9 назв. — рос.
1561-5359
https://nasplib.isofts.kiev.ua/handle/123456789/59849
004.93’1; 004.932; 528.85/.87(15)
На основе анализа основных соотношений полиномиального преобразования изображений космической съемки наземных объектов, используемого в процессе их топографической привязки, показана возможность оценок ее точности. Предложена методология соответствующих расчетов и моделирования, позволяющая как выполнять расчеты в конкретных задачах, так и вскрывать пути снижения ошибок, а также получить потенциальные характеристики точности, сняв имеющие место ограничения в распространении достигнутых результатов на любой практически возможный случай.
На основі аналізу основних співвідношень поліноміального перетворення зображень космічної зйомки наземних об’єктів, використовуваних в процесі їх топографічної прив’язки, показана можливість оцінок її точності. Запропонована методологія відповідних розрахунків і моделювання, що дозволяє як виконувати розрахунки в конкретних завданнях, так і розкривати шляхи зниження помилок, а також отримувати потенційні характеристики точності, знявши обмеження, що мають місце, в поширенні досягнутих результатів на будь-який практично можливий випадок.
On the basis of analysis of the basic proportion of polynomial transformation of images of space survey of the land objects, which is used during the process of the topographic reference of these objects, possibility of accuracy estimations is shown. The methodology of special calculations and modeling for practical calculation, for error reduction and also for getting potential characteristics of accuracy is offered. Thus the restriction in dissemination of the achieved results in any possible situation is removed.
ru
Інститут проблем штучного інтелекту МОН України та НАН України
Штучний інтелект
Интеллектуальные интерфейсы и распознавание образов. Системы цифровой обработки изображений
Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки
Закономірність поведінки оцінок точності поліноміального перетворення геометричної корекції зображень космічної зйомки
Regularity of Behavior of Accuracy Estimation at Polynomial Transformation of Image Space Survey Geographic Correlation
Article
published earlier
spellingShingle Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки
Белозерский, Л.А.
Мурашко, Н.И.
Интеллектуальные интерфейсы и распознавание образов. Системы цифровой обработки изображений
title Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки
title_alt Закономірність поведінки оцінок точності поліноміального перетворення геометричної корекції зображень космічної зйомки
Regularity of Behavior of Accuracy Estimation at Polynomial Transformation of Image Space Survey Geographic Correlation
title_full Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки
title_fullStr Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки
title_full_unstemmed Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки
title_short Закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки
title_sort закономерности поведения оценок точности полиномиального преобразования геометрической коррекции изображений космической съемки
topic Интеллектуальные интерфейсы и распознавание образов. Системы цифровой обработки изображений
topic_facet Интеллектуальные интерфейсы и распознавание образов. Системы цифровой обработки изображений
url https://nasplib.isofts.kiev.ua/handle/123456789/59849
work_keys_str_mv AT belozerskiila zakonomernostipovedeniâocenoktočnostipolinomialʹnogopreobrazovaniâgeometričeskoikorrekciiizobraženiikosmičeskoisʺemki
AT muraškoni zakonomernostipovedeniâocenoktočnostipolinomialʹnogopreobrazovaniâgeometričeskoikorrekciiizobraženiikosmičeskoisʺemki
AT belozerskiila zakonomírnístʹpovedínkiocínoktočnostípolínomíalʹnogoperetvorennâgeometričnoíkorekcíízobraženʹkosmíčnoíziomki
AT muraškoni zakonomírnístʹpovedínkiocínoktočnostípolínomíalʹnogoperetvorennâgeometričnoíkorekcíízobraženʹkosmíčnoíziomki
AT belozerskiila regularityofbehaviorofaccuracyestimationatpolynomialtransformationofimagespacesurveygeographiccorrelation
AT muraškoni regularityofbehaviorofaccuracyestimationatpolynomialtransformationofimagespacesurveygeographiccorrelation