RU2475852C1 - Method of constructing calculated shadow zones by controlled search of relief height matrix elements - Google Patents
Method of constructing calculated shadow zones by controlled search of relief height matrix elements Download PDFInfo
- Publication number
- RU2475852C1 RU2475852C1 RU2011126018/08A RU2011126018A RU2475852C1 RU 2475852 C1 RU2475852 C1 RU 2475852C1 RU 2011126018/08 A RU2011126018/08 A RU 2011126018/08A RU 2011126018 A RU2011126018 A RU 2011126018A RU 2475852 C1 RU2475852 C1 RU 2475852C1
- Authority
- RU
- Russia
- Prior art keywords
- matrix
- elevation
- shading
- zones
- digital
- Prior art date
Links
Images
Landscapes
- Processing Or Creating Images (AREA)
Abstract
Description
Изобретение относится к вычислительной технике и может быть использовано при построении и совершенствовании систем управления средствами мониторинга окружающей среды, получающими информацию различной физической природы об объектах или среде, находящихся на некотором удалении, а также активные информационные сигналы от удаленных технических средств, выполняющих различные задачи.The invention relates to computer technology and can be used in the construction and improvement of control systems for environmental monitoring tools that receive information of various physical nature about objects or the environment located at some distance, as well as active information signals from remote technical means that perform various tasks.
Предметом воздействия при разработке способа принята цифровая карта местности. Существуют различные способы построения цифровой карты местности, среди которых наиболее распространенным считается способ, когда цифровая карта местности составляется из атомарных элементов (площадей), имеющих свои координаты. Ввиду того, что атомарные элементы могут содержать множество характеристик местности, цифровая карта разбивается на двумерные массивы однотипных характеристик, часто их называют слоями. Индексами массивов выступают координаты, а значениями характеристики, например, для характеристик высоты рельефа каждому элементу соответствует усредненная в пределах периметра элемента высота.The subject of the impact in the development of the method adopted a digital map of the area. There are various ways to build a digital map of the area, among which the most common is the method when a digital map of the area is composed of atomic elements (areas) that have their coordinates. Due to the fact that atomic elements can contain many characteristics of the terrain, a digital map is divided into two-dimensional arrays of the same characteristics, often called layers. The coordinates are the indices of the arrays, and the values of the characteristics, for example, for the characteristics of the relief height, correspond to each element with the height averaged within the perimeter of the element.
Цифровые карты местности представляют собой информационные изделия. Особенностью данных изделий является то, что при их применении имеющиеся информационные возможности в ходе эксплуатации расширяются за счет добавления новых слоев пользователей. Данные возможности, как правило, допускаются изготовителями цифровых карт местности. В итоге пользователи в зависимости от характера воздействий на приобретенную цифровую карту местности в итоге всегда имеют свои оригинальным образом усовершенствованные изделия, которые применяют в соответствии со своими интересами.Digital terrain maps are information products. A feature of these products is that, when applied, the available information capabilities during operation are expanded by adding new layers of users. These features are generally permitted by manufacturers of digital terrain maps. As a result, users, depending on the nature of the impact on the acquired digital map of the area, as a result always always have their own original improved products that they use in accordance with their interests.
Одним из полезных усовершенствований являются слои, представляющие собой расчетные зоны затенения. Существуют разнообразные способы обработки информации цифровой карты местности с целью построения расчетных зон затенения.One useful improvement is the layers, which are the calculated shading zones. There are various ways of processing information of a digital map of the area in order to build settlement shading zones.
Известен способ построения расчетных зон затенения, который позволяет из точки стояния объекта оценить возможность прямой видимости удаленных объектов, находящихся на поверхности Земли или на некоторой высоте над Землей [2].There is a method of constructing settlement shading zones, which allows you to evaluate the possibility of direct visibility of distant objects located on the Earth’s surface or at a certain height above the Earth from an object’s standing point [2].
Исходными данными являются: цифровая карта высот рельефа местности, координаты точки стояния объекта, превышение объекта наблюдения над уровнем рельефа и превышение наблюдаемых объектов над уровнем рельефа, которые могут находиться в любой точке рассматриваемой области пространства. В качестве математического описания рельефа применяется матрица, содержащая значения высот, усредненных на площади квадрата, выбранного в качестве минимальной дискретной единицы измерения перепада высот. Матрица привязывается к координатной сетке применяемой системы координат и позволяет по координатам точки на карте выбирать конкретное значение усредненной высоты из двумерного массива матрицы высот.The initial data are: a digital map of the elevation of the terrain, the coordinates of the object’s standing point, the excess of the observation object over the relief level and the excess of the observed objects over the relief level, which can be located at any point in the considered area of space. As a mathematical description of the relief, a matrix is used that contains the values of the heights averaged over the area of the square selected as the minimum discrete unit of measurement for the height difference. The matrix is attached to the coordinate grid of the applied coordinate system and allows you to select the specific value of the averaged height from the two-dimensional array of the height matrix by the coordinates of the point on the map.
Для этого, Фиг.1, все пространство вокруг точки стояния объекта наблюдения О разбивается на равные сектора в горизонтальной плоскости Q, причем точка О∈Q. Угол сектора α в зависимости от требуемой точности может увеличиваться или уменьшаться. Так как обработка информации в секторах по всей окружности производится по одному и тому же алгоритму, на Фиг.1 показаны в виде равномерной прямоугольной сетки (50×50 квадратов) участок матрицы расположения дискретных квадратов усредненных высот, а также два смежных сектора: АОВ и BОС.For this, FIG. 1, the entire space around the standing point of the observation object O is divided into equal sectors in the horizontal plane Q, and the point O∈Q. The sector angle α, depending on the required accuracy, can increase or decrease. Since information processing in sectors around the circumference is carried out according to the same algorithm, Fig. 1 shows a plot of a matrix of discrete squares of averaged heights as a uniform rectangular grid (50 × 50 squares), as well as two adjacent sectors: AOB and BOC .
В секторах из точки О проводятся биссектрисы OF1 и OF2. Согласно методике на биссектрисах через равные дискретные расстояния Δd фиксируются точки. По координатам точек определяются индексы двумерного массива данных матрицы высот, а затем и конкретные числовые значения высоты фиксированных точек.In the sectors from point O, the bisectors OF 1 and OF 2 are drawn. According to the technique, points are fixed on bisectors through equal discrete distances Δd. The coordinates of the points determine the indices of the two-dimensional array of data of the DEM, and then the specific numerical values of the height of the fixed points.
Затем для каждой биссектрисы строится срез рельефа местности в вертикальной плоскости, проведенной так, что биссектриса принадлежит этой вертикальной плоскости. На Фиг.2 показана вертикальная плоскость U для случая, когда OF2∈U. Проекция кривой, соединяющей точки, размещенные на соответствующих высотах, на вертикальную плоскость U представляет собой линию "среза" рельефа местности, изменяющего значение высот h через дискретные расстояния Δd выбора высоты.Then, for each bisector, a section of the terrain is constructed in a vertical plane, drawn so that the bisector belongs to this vertical plane. Figure 2 shows the vertical plane U for the case when OF 2 ∈U. The projection of the curve connecting the points located at the corresponding heights onto the vertical plane U represents the “cut” line of the terrain, changing the value of the heights h through discrete distances Δd of the choice of height.
Зоны затенения представляют собой совокупность элементарных площадей усеченных секторов вида NLPR, Фиг.1. Принадлежность секторов зоне затенения или зоне видимости определяется видимостью из точки наблюдения центральной точки усеченных секторов, для случая, рассмотренного на Фиг.1, точка r4∈NLPR.Shading zones are a set of elementary areas of truncated sectors of the NLPR type, FIG. 1. The sectors belonging to the shading zone or visibility zone is determined by the visibility from the observation point of the center point of the truncated sectors, for the case considered in Figure 1, the point r 4 ∈NLPR.
Видимость точки определяется следующим образом. На "срезе" рельефа местности, Фиг.2, из точки наблюдения Тн проводятся лучи в точки, характеризующие дискретные высоты рельефа. Для случая, изображенного на Фиг.2, это точки r1, r2, r3, r4, r5.The visibility of a point is defined as follows. On the "slice" of the terrain, Figure 2, from the observation point T n rays are drawn to points characterizing discrete elevation of the relief. For the case depicted in FIG. 2, these are points r 1 , r 2 , r 3 , r 4 , r 5 .
Затем последовательно от ближней точки к самой дальней проводится анализ углов, образованных поворотом луча от Тнh вокруг точки наблюдения Тн по часовой стрелке до направлений в точки r1, r2, r3, r4, r5. Для случая на Фиг.2 углы α1, α2, α3, α4, α5.Then, successively from the closest point to the farthest, an analysis is made of the angles formed by the rotation of the beam from T n h around the observation point T n clockwise to directions to the points r 1 , r 2 , r 3 , r 4 , r 5 . For the case of FIG. 2, the angles α 1 , α 2 , α 3 , α 4 , α 5 .
На каждом k-м шаге определяется соотношение очередного угла αk с текущим углом закрытия βзакр, а затем проводится уточнение величины угла закрытия.At each k-th step, the ratio of the next angle α k with the current angle of closure β clos is determined, and then the closure angle is refined.
Условие видимости точки rk из точки наблюдения Тн состоит в том, чтоThe condition for the visibility of the point r k from the observation point T n is that
Во всех остальных случаях считается, что очередная точка rk невидима из точки наблюдения Tн, а соответствующий усеченный сектор принадлежит зоне затенения.In all other cases, it is believed that the next point r k is invisible from the observation point T n , and the corresponding truncated sector belongs to the shading zone.
На каждом шаге уточнение величины угла закрытия выполняется следующим образом. Начальное значение угла закрытия βзакр 0 принимается равным 180°. На каждом k-м шаге анализа закрытия после проверки видимости точки (1) выполняется операция сравнения:At each step, the refinement of the closing angle is as follows. The initial value of the angle of closure β clos 0 is taken equal to 180 °. At each k-th step of the analysis of closure, after checking the visibility of point (1), the comparison operation is performed:
и определяется очередное значение угла закрытия.and the next value of the closing angle is determined.
Недостаток данного способа состоит в том, что с увеличением расстояния элементарные площади, составляющие зоны затенения, увеличиваются, а усредненная высота в усеченном секторе при наличии пересеченной местности все меньше соответствует реальным высотам на направлениях. Попытки уменьшения углов секторов приводят к значительному количеству избыточных вычислений, так как в этом случае соседние биссектрисы многократно проходят по одним и тем же дискретным единицам измерения перепада высот цифровой карты рельефа местности, Фиг.3.The disadvantage of this method is that with increasing distance, the elementary areas that make up the shading zones increase, and the average height in the truncated sector in the presence of rough terrain less and less corresponds to the real heights in the directions. Attempts to reduce the angles of the sectors lead to a significant amount of redundant calculations, since in this case, the adjacent bisectors repeatedly pass through the same discrete units of measurement of the height difference of the digital map of the terrain, Fig.3.
Существенным недостатком данного способа является необходимость применения тригонометрических функций при определении дискретных квадратов, соответствующих точкам, находящимся на биссектрисах секторов через равные промежутки Δd.A significant drawback of this method is the need to use trigonometric functions in determining discrete squares corresponding to points located on the bisectors of sectors at equal intervals Δd.
С точки зрения повышения адекватности построения зон затенения более совершенным считается способ, когда прямая видимость определяется отдельно в каждой элементарной дискретной единице цифровой матрицы высот рельефа.From the point of view of increasing the adequacy of the construction of shading zones, a more perfect way is considered when the direct visibility is determined separately in each elementary discrete unit of the digital matrix of elevation of the relief.
Особенность данного способа состоит в том, что из точки наблюдения проводятся лучи в центр каждого квадрата, являющегося элементарной дискретной единицей цифровой матрицы высот рельефа, Фиг 4. Оценка возможности прямой видимости производится не по выбранной дискретной дальности Δd, а при каждом переходе из одной дискретной единицы перепада высот в очередную по направлению луча, направленного из точки стояния О в направлении выбранного для анализа квадрата, т.е. точки 1, 2, 3, 4, 5, Фиг.4.A feature of this method is that rays are drawn from the observation point to the center of each square, which is an elementary discrete unit of the digital elevation matrix of the relief, Fig 4. The line of sight visibility is estimated not by the selected discrete range Δd, but by each transition from one discrete unit the vertical difference in the next direction of the beam directed from the point of standing O in the direction of the square selected for analysis, i.e.
Применение данного способа позволяет строить зоны затенения не из усеченных секторов, обработка формы которых представляет затруднения, а из совокупности площадей элементарных квадратов, представляющих собой элементарные дискретные единицы матрицы перепада высот. Такое построение является предельным при заданном размере дискретной единицы перепада высот и, соответственно, обладает большей адекватностью формы зоны затенения рельефу, представленному матрицей высот, используемой в расчетах. Кроме того, вычисление угла закрытия определяется с учетом высоты каждой элементарной единицы матрицы высоты на пути продвижения луча.The application of this method allows you to build a zone of shading not from truncated sectors, the processing of the form of which is difficult, but from the set of areas of elementary squares, which are elementary discrete units of the elevation matrix. Such a construction is limiting for a given size of a discrete unit of height difference and, accordingly, has a greater adequacy of the shape of the shading zone to the relief represented by the matrix of heights used in the calculations. In addition, the calculation of the closing angle is determined taking into account the height of each elementary unit of the height matrix along the path of the beam.
Продвижение луча в очередной квадрат, выбранный для оценки прямой видимости, связано не с угловым, а с линейным перемещением от одного смежного квадрата к другому, что позволяет определять очередной квадрат без применения тригонометрических функций и тем самым значительно сократить процесс определения квадратов за счет применения алгоритма Брезенхэма [1], позволяющего проводить минимальное количество вычислений для выбора очередного смежного квадрата на луче.The advance of the beam into the next square, chosen for assessing line of sight, is associated not with angular, but with linear displacement from one adjacent square to another, which allows us to determine the next square without using trigonometric functions and thereby significantly reduce the process of determining squares by using the Bresenham algorithm [1], which allows the minimum number of calculations to select the next adjacent square on the beam.
Вместе с тем, как видно на Фиг.4, и при этом способе не удается избавиться от повторений, т.к. соседние лучи также проходят большое количество одних и тех же элементарных квадратов цифровой матрицы рельефа.However, as can be seen in Figure 4, and with this method it is not possible to get rid of repetitions, because neighboring rays also pass through a large number of the same elementary squares of the digital relief matrix.
Таким образом, оба способа при достижении максимальной адекватности содержат большое количество избыточных вычислений, что приводит к большой продолжительности решения задачи построения расчетной зоны затенения. В зависимости от величины рассматриваемого участка местности и уровня детализации перепада высот в цифровой карте продолжительность может составлять часы или даже сутки.Thus, both methods, upon reaching maximum adequacy, contain a large number of redundant calculations, which leads to a long duration of the solution to the problem of constructing the estimated shading zone. Depending on the size of the area under consideration and the level of detail of the elevation in a digital map, the duration can be hours or even a day.
В тех случаях, когда предполагается, что точка наблюдения будет долгое время находиться на одном месте, например вышка радиоретранслятора, увеличение продолжительности определения расчетных зон затенения не сказывается на процессе организации связи. Когда же точка наблюдения перемещается, например, движение самолета или лидарной установки, задержка вычислений даже на насколько минут приводит к тому, что полученные результаты будут бесполезными для выполнения задач наблюдения заданных участков местности или воздушного пространства.In cases where it is assumed that the observation point will be in one place for a long time, for example, a tower of a radio relay, an increase in the duration of the determination of the estimated shading zones does not affect the communication process. When the observation point moves, for example, the movement of an airplane or a lidar installation, a delay in calculations, even for how many minutes, leads to the fact that the results will be useless for performing the tasks of observing predetermined sections of the terrain or airspace.
В качестве прототипа данного изобретения выбран способ, сходный с заявляемым по сущности воздействия на информационное изделие, которое изменяет свое состояние за счет наращивания информационных возможностей введением слоев зон затенения для определенных условий. Прямая видимость определяется отдельно в каждой элементарной дискретной единице цифровой матрицы высот рельефа в пределах заданного участка просмотра.As a prototype of the present invention, a method similar to the claimed in essence impact on an information product is selected, which changes its state by increasing information capabilities by introducing layers of shading zones for certain conditions. Direct visibility is determined separately in each elementary discrete unit of the digital matrix of elevation reliefs within a given viewing area.
Указанный способ основывается на следующих положениях и исходных данных.The specified method is based on the following provisions and source data.
Под прямой видимостью понимается беспрепятственное прохождение луча из точки наблюдения в точку, для которой определяется прямая видимость. Зона затенения - зона, представляющая собой совокупность точек, для которых не выполняется условие прямой видимости.By direct visibility is meant the unhindered passage of a ray from an observation point to a point for which direct visibility is determined. Shading zone - a zone representing a set of points for which the line of sight condition is not fulfilled.
Исходные данные: Цифровая матрица высот рельефа, координаты точки стояния объекта наблюдения, значение высоты превышения объекта наблюдения над уровнем рельефа и значение превышения наблюдаемых объектов над уровнем рельефа.Initial data: Digital matrix of elevation of the relief, coordinates of the point of standing of the object of observation, the value of the height of the excess of the object of observation over the level of the relief and the value of the excess of the observed objects over the level of the relief.
Основой построения зон затенения является цифровая матрица высот рельефа.The basis for constructing shading zones is a digital matrix of elevation elevations.
Вся площадь местности, на которой производится построение зон затенения, разбивается на квадраты со стороной l. Каждый квадрат представляет собой минимальную дискретную единицу площади, на которой значение высоты принимается неизменным. Для определения этого значения на участке местности, ограничиваемом периметром квадрата, определяется hmax - максимальное и hmin - минимальное значение высоты местности (местных предметов), а затем определяется среднее значение, которое принимается за высоту, приписываемую дискретной единицеThe entire area of the area on which the construction of shading zones is performed is divided into squares with side l. Each square represents the minimum discrete unit of area over which the height value is assumed unchanged. To determine this value in the area of the area limited by the perimeter of the square, h max is determined as the maximum and h min is the minimum value of the terrain height (local objects), and then the average value is determined, which is taken as the height attributed to the discrete unit
Для корректного применения в математическом аппарате вычислений все значения дискретных единиц измерения высоты сводятся в прямоугольную матрицу H[m×n]. Таким образом, на основании матрицы при наличии координат квадрата i,j определяется высота рельефа на направлении распространения луча прямой видимости 
Координаты точки стояния объекта наблюдения позволяют определить исходный квадрат цифровой матрицы высот, из центра которого строятся лучи прямой видимости в центры квадратов, выбираемых для определения наличия прямой видимости.The coordinates of the point of observation of the object of observation make it possible to determine the initial square of the digital matrix of heights, from the center of which line of sight lines are built into the center of the squares selected to determine the presence of line of sight.
Значение высоты превышения объекта наблюдения над уровнем рельефа и значение превышения наблюдаемых объектов над уровнем рельефа позволяют определить расположение точек в пространстве. Высота представляет размер вертикальной линии из точки наблюдения (наблюдаемой точки) до пересечения с поверхностью, описывающей форму рельефа.The value of the height of the excess of the object of observation over the level of the relief and the value of the excess of the observed objects over the level of the relief make it possible to determine the location of points in space. The height represents the size of the vertical line from the observation point (observed point) to the intersection with the surface describing the shape of the relief.
Сущность способа состоит в том, что на горизонтальной плоскости фиксируются координаты точки стояния объекта наблюдения с точностью до квадрата цифровой матрицы высот. По координатам квадрата определяется высота рельефа. Сложением высоты рельефа в точке стояния объекта наблюдения и высоты превышения объекта наблюдения над уровнем рельефа определяется абсолютная высота объекта наблюдения. Затем из этой точки определяется совокупность квадратов, для которых выполняется условие прямой видимости объектов, находящихся в квадратах на фиксированной высоте превышения над рельефом, все остальные квадраты составляют зону затенения для заданных условий, Фиг.5.The essence of the method lies in the fact that on the horizontal plane the coordinates of the standing point of the observation object are fixed to the square of the digital DEM. The coordinates of the square determine the height of the relief. The addition of the height of the relief at the point of standing of the object of observation and the height of the excess of the object of observation over the level of the relief determines the absolute height of the object of observation. Then, from this point, a set of squares is determined for which the condition of direct visibility of objects that are in squares at a fixed height of elevation above the relief is satisfied, all other squares make up the shading zone for the given conditions, Figure 5.
Для определения наличия прямой видимости в квадратах на горизонтальной плоскости последовательно строятся прямые отрезки, соединяющие точку центра квадрата, в котором расположен объект наблюдения с точками, расположенными в центрах квадратов, рассматриваемых на возможность прямой видимости. Выборка для оценки возможности прямой видимости для квадратов с координатами (4,13); (4,14); (4,15); (4,16); (4,17); (7,5) показана на Фиг.6.To determine the presence of direct visibility in squares on a horizontal plane, straight segments are sequentially constructed connecting the point of the center of the square where the object of observation is located with points located in the centers of the squares considered for the possibility of direct visibility. A sample for assessing the possibility of line of sight for squares with coordinates (4.13); (4.14); (4.15); (4.16); (4.17); (7.5) is shown in FIG. 6.
Через каждый отрезок строится вертикальная плоскость U. На этой плоскости воспроизводится срез рельефа местности как линия, соединяющая точки на плоскости, которые получаются следующим образом, Фиг.7.A vertical plane U is constructed through each segment. A section of the terrain is reproduced on this plane as a line connecting the points on the plane, which are obtained as follows, Fig.7.
Если отрезок пересекает квадрат цифровой матрицы рельефа, Фиг.7, квадрат (2,1), из центра этого квадрата точка k21 проводится перпендикулярно к отрезку ОК. Из точки пересечения отрезка с перпендикуляром, точка a 1, строится перпендикуляр в вертикальной плоскости U, длина которого соответствует высоте рельефа в рассматриваемом квадрате. Искомая точка h1 находится на удаленном от отрезка конце перпендикуляра.If the segment intersects the square of the digital relief matrix, Fig. 7, the square (2,1), from the center of this square, the point k 21 is drawn perpendicular to the segment OK. From the point of intersection of the segment with the perpendicular, point a 1 , a perpendicular is constructed in the vertical plane U, the length of which corresponds to the height of the relief in the square under consideration. The desired point h 1 is located at the end of the perpendicular remote from the segment.
На Фиг.8 показано построение рельефа местности применительно для квадрата (7,5).On Fig shows the construction of the terrain in relation to the square (7.5).
Видимость точки очередного квадрата, оцениваемого на предмет прямой видимости, определяется следующим образом. На "срезе" рельефа местности, Фиг.9, из точки наблюдения Тн проводятся лучи в точки h1, h2,…,h9, характеризующие высоты рельефа квадратов, пересекаемых плоскостью U, а также луч в точку центра объекта наблюдения Тон, имеющую высоту предполагаемого объекта наблюдения в оцениваемом квадрате. Затем проводится анализ углов α, образованных вращением вертикали ТнН вокруг точки Tн по часовой стрелке до совпадения с проведенными лучами Тнh1, Тнh2, Тнh3 и т.д. до Tнh9 с углом φ, образованным вращением вертикали ТнН вокруг точки Tн по часовой стрелке до совпадения с лучом, проведенным из точки Тн в точку Тон, характеризующую высоту предполагаемого объекта наблюдения в оцениваемом квадрате.The visibility of the point of the next square, evaluated for line of sight, is determined as follows. On the "slice" of the terrain, Fig.9, from the observation point T n rays are drawn to the points h 1 , h 2 , ..., h 9 , characterizing the elevation of the relief of the squares intersected by the plane U, and also the beam to the point of the center of the object of observation T he having the height of the proposed object of observation in the estimated square. Then, an analysis is made of the angles α formed by the vertical rotation of T n H around the point T n clockwise until they coincide with the rays T n h 1 , T n h 2 , T n h 3 , etc. T h n to 9 with an angle φ, formed by rotating the vertical T n H around point T n clockwise into coincidence with the ray drawn from the point T H to the point T he characterizing the intended height of the object of observation in a ranked squared.
При данном построении условие прямой видимости точки Tон из Tн будет следующее:With this construction of the line of sight of the point T he T n will be the following:
Во всех остальных случаях считается, что очередная точка расположения объекта наблюдения Тон невидима из точки наблюдения Тн, и соответствующий квадрат принадлежит зоне затенения, определяемой при заданных исходных условиях. Значение δ=αmin из α∈(α1,α9) называется углом закрытия для очередного анализируемого квадрата. Данный угол всегда меньше или равен φ для квадратов зоны затенения, больше φ для квадратов зоны видимости.In all other cases, it is believed that the next point of the observation object T is invisible from the observation point T n , and the corresponding square belongs to the shading zone, which is determined under given initial conditions. The value δ = α min of α∈ (α 1 , α 9 ) is called the closing angle for the next analyzed square. This angle is always less than or equal to φ for the squares of the shadow zone, greater than φ for the squares of the visibility zone.
Анализ способа показал, что имеется большое количество избыточных вычислений, вызванных дублированием анализа углов α при вычислении видимости наблюдаемых объектов в соседних квадратах, что приводит к большой продолжительности решения задачи построения расчетной зоны затенения. В зависимости от величины рассматриваемого участка местности, заданного количества фиксированных высот превышения наблюдаемых объектов над рельефом, а также уровня детализации перепада высот в цифровой карте продолжительность может составлять часы и даже сутки. Это обстоятельство ограничивает эффективность применения аналитического анализа рельефа местности в системах, решающих задачи в условиях ограниченного времени.The analysis of the method showed that there is a large number of redundant calculations caused by duplication of the analysis of angles α when calculating the visibility of the observed objects in adjacent squares, which leads to a long duration of the solution to the problem of constructing the estimated shading zone. Depending on the size of the considered area, the specified number of fixed heights for the elevation of the observed objects over the terrain, as well as the level of detail of the elevation in a digital map, the duration can be hours or even a day. This circumstance limits the effectiveness of the application of analytical analysis of the terrain in systems that solve problems in a limited time.
Попытки исключения дублирования анализа углов α при вычислении видимости наблюдаемых объектов в соседних квадратах применением дополнительных вычислительных модулей позволило уменьшить количество полностью одинаковых вычислений, однако в рамках данного способа временные затраты на поиск и обработку результатов найденных повторяемых последовательностей на "срезе" рельефа местности соизмеримы с временными затратами на повторный анализ, а иногда и превосходят их.Attempts to exclude duplication of the analysis of angles α in calculating the visibility of observed objects in adjacent squares using additional computational modules made it possible to reduce the number of completely identical calculations, however, in the framework of this method, the time taken to search and process the results of the repeated sequences found on the “slice” of the terrain is comparable to the time taken to re-analyze, and sometimes surpass them.
Недостатком способа прототипа являются ограниченные возможности по исключению избыточных вычислений, вызванных дублированием анализа углов α при вычислении видимости наблюдаемых объектов в соседних квадратах цифровой карты местности.The disadvantage of the prototype method is the limited ability to eliminate redundant calculations caused by duplication of the analysis of angles α when calculating the visibility of the observed objects in adjacent squares of a digital terrain map.
Целью настоящего изобретения является сокращение времени на построение расчетных зон затенения, применяемых при организации управления системами сбора информации.The aim of the present invention is to reduce the time to build settlement zones of shading used in organizing the management of information collection systems.
В основу изобретения поставлена задача разработки нового способа построения расчетных зон затенения, позволяющего исключить повторные вычисления при прохождении соседних лучей прямой видимости через одни и те же элементарные квадраты цифровой матрицы высот рельефа местности.The basis of the invention is the task of developing a new method of constructing settlement shading zones, which eliminates repeated calculations when passing adjacent lines of direct visibility through the same elementary squares of the digital matrix of the elevation of the terrain.
Решение задачи достигается тем, что построение расчетных зон затенения основывается на цифровом блоке управления, который хранит информацию о наличии дублируемых участков и позволяет расчетные действия заменить управляемым выбором, а также исключить повторный анализ углов α в ходе определения прямой видимости наблюдаемых объектов, находящихся в расположенных рядом квадратах.The solution to the problem is achieved by the fact that the construction of the estimated shading zones is based on a digital control unit that stores information about the presence of duplicated sections and allows the calculation of actions to be replaced by a controlled choice, and also to exclude a reanalysis of the angles α in determining the direct visibility of the observed objects located nearby squares.
Цифровой блок управления создается до начала обработки информации цифровой матрицы высот рельефа местности для формирования зон затенения в соответствии с:A digital control unit is created before processing the information of the digital matrix of the elevation of the terrain for the formation of shading zones in accordance with:
1. размерами площади анализируемой территории;1. the size of the area of the analyzed territory;
2. характеристиками цифровой матрицы высот рельефа;2. characteristics of a digital matrix of elevation of a relief;
3. принятой последовательностью просмотра квадратов цифровой матрицы высот рельефа на наличие прямой видимости объекта, находящегося в пределах периметра выбранного квадрата на поверхности Земли, или на некотором превышении из точки наблюдения, также привязанной к конкретному квадрату цифровой матрицы высот рельефа.3. the adopted sequence of viewing the squares of the digital elevation matrix for the presence of direct visibility of an object located within the perimeter of the selected square on the Earth's surface, or at some excess from the observation point, also tied to a specific square of the digital elevation matrix.
Основу блока управления составляют две матрицы: матрица управления вычислениями и матрица хранения промежуточных результатов вычислений.The basis of the control unit is made up of two matrices: a matrix for managing calculations and a matrix for storing intermediate results of calculations.
Матрица управления вычислениями представляет собой матрицу Y[n×n], где n - сторона наибольшего измерения цифровой матрицы высот, накладываемой на анализируемый участок местности. В ней содержатся индексы квадратов для перехода алгоритма на новый этап вычислений.The computational control matrix is a Y [n × n] matrix, where n is the side of the largest dimension of the digital elevation matrix superimposed on the analyzed area. It contains square indices for the transition of the algorithm to a new stage of calculations.
Матрица хранения промежуточных результатов вычислений представляет собой матрицу P[n×n], где n - сторона наибольшего измерения цифровой матрицы высот, накладываемой на анализируемый участок местности. В ней хранятся углы закрытия, определенные в ходе вычислений, которые могут быть использованы на последующих этапах вычислений для исключения повторений из-за прохождения луча, определяющего "срез" рельефа, по одним и тем же квадратам цифровой карты местности.The matrix for storing intermediate results of calculations is a matrix P [n × n], where n is the side of the largest dimension of the digital matrix of heights superimposed on the analyzed area. It stores the closing angles determined during the calculations, which can be used at subsequent stages of the calculations to exclude repetitions due to the passage of the beam that determines the "slice" of the relief, on the same squares of a digital terrain map.
Сущность функционирования блока управления состоит в том, что перед началом анализа местности на наличие прямой видимости определяется последовательность просмотра квадратов цифровой карты местности относительно квадрата расположения объекта наблюдения. Данное действие необходимо ввиду того, что от последовательности анализируемых квадратов зависят состав квадратов и значения углов закрытия для квадратов, предшествующих очередному анализируемому квадрату, по которым определяется дальнейшая последовательность вычислений.The essence of the operation of the control unit is that before starting the analysis of the area for the presence of direct visibility, the sequence of viewing the squares of a digital map of the area relative to the square of the location of the object of observation is determined. This action is necessary because the composition of the squares and the values of the closing angles for the squares preceding the next analyzed square depend on the sequence of the analyzed squares, which determines the subsequent sequence of calculations.
В ходе вычислений, начиная с первого выбранного для анализа квадрата, каждый очередной угол закрытия δ сохраняется в матрице хранения промежуточных результатов вычислений и запоминаются его индексы для последующего выбора.During the calculations, starting from the first square selected for analysis, each successive closing angle δ is stored in the storage matrix of the intermediate results of the calculations and its indices are stored for subsequent selection.
При анализе "среза" рельефа местности очередного квадрата перед началом сравнения углов α делается проверка информации об углах закрытия ранее рассмотренных квадратов, которая проводится следующим образом.When analyzing the "slice" of the terrain of the next square, before comparing the angles α, information on the closing angles of the previously considered squares is checked, which is carried out as follows.
Пусть в очередном случае (например, для k75) "срез" рельефа проходит по квадратам в следующей последовательности:Let in the next case (for example, for k 75 ) the “slice” of the relief passes along the squares in the following sequence:
Если среди результатов предыдущих вычислений найдется хоть одна последовательность из следующего набора последовательностей:If among the results of previous calculations there is at least one sequence from the following set of sequences:
имеющая для крайнего правого квадрата значение угла закрытия δ, то в зависимости от найденной последовательности вычисление углов в данном очередном случае выполняется только для оставшихся членов последовательности (6), а значение исходного угла α для крайнего правого квадрата принимается равным δ.having a closing angle δ for the far right square, then, depending on the sequence found, the calculation of the angles in this next case is performed only for the remaining members of the sequence (6), and the value of the initial angle α for the far right square is taken equal to δ.
Следует отметить, что при наличии двух и более таких последовательностей выбирается последовательность с большим количеством членов k.It should be noted that in the presence of two or more such sequences, a sequence with a large number of members k is selected.
Информация о наиболее продолжительной последовательности для каждого квадрата определяется заранее и хранится в матрице управления вычислениями, а δ в матрице хранения промежуточных результатов вычислений.Information about the longest sequence for each square is determined in advance and stored in the calculation control matrix, and δ in the storage matrix of intermediate calculation results.
Таким образом, при обработке данных цифровой карты местности с целью построения расчетной зоны затенения циклическая обработка данных перебираемых квадратов "среза" рельефа местности частично заменяется управляемым выбором значений.Thus, when processing data from a digital terrain map in order to construct a calculated shading zone, the cyclic processing of the data of sorted squares of the “slice” of the terrain is partially replaced by a controlled choice of values.
Сопоставительный анализ заявляемого способа построения расчетных зон затенения управляемым перебором квадратов цифровой карты местности с прототипом, выполняющим циклическую обработку данных перебираемых квадратов "среза" рельефа местности, позволяет сделать вывод, что заявляемый способ базируется на новых информационных элементах, входящих в цифровой блок управления: матрица управления вычислениями и матрица хранения промежуточных результатов вычислений, поэтому заявляемый способ соответствует критерию "новизна".A comparative analysis of the proposed method for constructing the estimated shading zones by controlled enumeration of squares of a digital terrain map with a prototype performing cyclic data processing of the selected squares of the "slice" of the terrain allows us to conclude that the claimed method is based on new information elements included in the digital control unit: control matrix calculations and the storage matrix of intermediate results of calculations, therefore, the inventive method meets the criterion of "novelty."
Известные способы построения расчетных зон затенения не позволяют исключить избыточность вычислений, вызванную повторным анализом углов α при вычислении видимости наблюдаемых объектов в соседних квадратах. Введенные элементы за счет своего информационного состояния, сформированного до начала построения расчетных зон затенения, позволяют по-новому осуществлять воздействие на информационное изделие цифровой карты местности, исключают дублирование и в итоге достигается требуемый результат, но с меньшими временными затратами, чем при применении известных способов. Этот эффект, достигаемый применением заявляемого способа, позволяет сделать вывод о его соответствии критерию "существенные отличия".Known methods for constructing design shading zones do not eliminate the redundancy of calculations caused by the repeated analysis of the angles α when calculating the visibility of the observed objects in adjacent squares. The introduced elements, due to their informational state, formed before the construction of the estimated shading zones, allow for a new effect on the information product of the digital terrain map, eliminate duplication and ultimately achieve the desired result, but with less time than using known methods. This effect, achieved by the application of the proposed method, allows us to conclude that it meets the criterion of "significant differences".
Введение указанных отличительных признаков позволило разработать новый способ построения расчетных зон затенения управляемым перебором квадратов цифровой карты местности, свободного от повторных вычислений при прохождении соседних лучей прямой видимости через одни и те же элементарные квадраты цифровой матрицы высот рельефа местности, получить сокращение времени на построение расчетных зон затенения, применяемых при организации управления системами сбора информации, и тем самым достичь цели изобретения.The introduction of these distinguishing features made it possible to develop a new method for constructing shading calculation zones by controlled enumeration of squares of a digital terrain map, free from repeated calculations when passing adjacent line of sight lines through the same elementary squares of a digital matrix of elevations of the terrain, and to obtain a reduction in the time for building settlement shading zones used in organizing the management of information collection systems, and thereby achieve the objective of the invention.
Для примера и последующей оценки эффективности способа построения расчетных зон затенения управляемым перебором квадратов цифровой карты приводится последовательность, когда перебор квадратов осуществляется по секторам 1-8, равным 45 градусам от ближнего квадрата к дальнему по вертикали и от ближней вертикали к дальней по горизонтали.For an example and a subsequent evaluation of the effectiveness of the method of constructing the estimated shading zones by controlled enumeration of the squares of the digital map, a sequence is given when the enumeration of the squares is carried out in sectors 1-8, equal to 45 degrees from the near square to the far vertical and from the near vertical to the far horizontal.
Расположение секторов относительно объекта наблюдения V с координатами (а,b) представлено на Фиг.10. Область оценки зоны затенения в примере ограничена совокупностью квадратов с индексами в пределах: по горизонтали от с до d; по вертикали от е до g.The location of the sectors relative to the object of observation V with coordinates ( a , b) is presented in FIG. 10. The evaluation area of the shading zone in the example is limited to a set of squares with indices in the range: horizontally from c to d; vertically from e to g.
Последовательность перебора в секторах одна и та же. На примере сектора 1 квадраты выбираются следующим образом.The search sequence in sectors is the same. On the example of
1. Создается дополнительная переменная w=0.1. An additional variable w = 0 is created.
2. Задаются индексы начального квадрата i=а+1, j=h+w+2.2. The indices of the initial square i = a +1, j = h + w + 2 are specified.
3. Последующие квадраты выбираются путем увеличения индекса j на единицу до тех пор, пока j≤е.3. Subsequent squares are selected by increasing the index j by one until j≤e.
4. Значение i увеличивается на единицу.4. The value of i increases by one.
5. Если i<d, значение w увеличивается на единицу, а значение индекса j=b+w+2, осуществляется переход к пункту 3, иначе процесс перебора квадратов в секторе прекращается.5. If i <d, the value of w increases by one, and the value of the index j = b + w + 2, go to
Для описанной выше последовательности перебора квадратов перед началом вычислений матрица управления вычислениями будет иметь вид, показанный на Фиг.11. Матрица хранения промежуточных результатов вычислений будет иметь аналогичные размеры, однако перед началом вычислений значения массива заполняются нулями.For the above-described sequence of enumerating the squares before starting the calculations, the calculation control matrix will have the form shown in FIG. 11. The matrix for storing intermediate results of calculations will have similar dimensions, however, before starting the calculations, the array values are filled with zeros.
В матрице управления вычислениями нулями обозначены квадраты, для которых при анализе среза рельефа местности вообще отсутствуют повторяющиеся последовательности. В этом случае управляющие воздействия сводятся к выбору начального квадрата на срезе рельефа местности. Исключение составляют диагонали, вертикали и горизонтали.In the computational control matrix, zeros denote the squares for which when analyzing a slice of the terrain there are no repeating sequences at all. In this case, the control actions are reduced to the choice of the initial square on the slice of the terrain. The exception is diagonals, verticals and horizontals.
Если в квадрате записано число, отличное от нуля, то это значит, что данный квадрат уже имеет последовательность, которая при выбранном способе перебора ранее уже рассматривалась и имеет значение угла закрытия δ. Данное значение ранее, после анализа последовательности, было помещено в матрицу хранения промежуточных результатов вычислений.If a number other than zero is written in the square, this means that this square already has a sequence that, with the selected search method, has already been considered and has a closing angle value δ. This value earlier, after analysis of the sequence, was placed in the storage matrix of intermediate results of calculations.
Индексы квадрата, который является последним в ранее рассмотренной последовательности наибольшей продолжительности, а также значение угла закрытия δ определяется параметрами матрицы управления вычислениями.The indices of the square that is the last in the previously considered sequence of the greatest duration, as well as the value of the closing angle δ, are determined by the parameters of the calculation control matrix.
Порядок определения индексовIndex Definition Procedure
Начальное значение (0,0) совмещается с квадратом (a,b) матрицы высот рельефа, в котором располагается объект наблюдения, Фиг 12.The initial value (0,0) is combined with the square ( a , b) of the matrix of heights of the relief in which the object of observation is located, Fig 12.
Очередному выбранному квадрату (sk,zk) матрицы высот рельефа в этом случае соответствует квадрат КU с координатами (i,j), а также симметричный квадрат К-U с координатами (j,i) матрицы управления вычислениями. При этом для данных условий i=sk-a-1, j=zk-b-1.The next selected square (s k , z k ) of the elevation matrix in this case corresponds to the square K U with coordinates (i, j), as well as the symmetrical square K- U with coordinates (j, i) of the calculation control matrix. Moreover, for these conditions i = s k - a -1, j = z k -b-1.
Индексы (х,у) квадрата Kδ матрицы хранения промежуточных результатов вычислений, в котором хранится значение угла закрытия δ последовательности наибольшей продолжительности для рассматриваемого квадрата, определяются как:The indices (x, y) of the square K δ of the storage matrix of the intermediate calculation results, in which the value of the closing angle δ of the sequence of the longest duration for the square in question is stored, are defined as:
таким образом,in this way,
Индексы (sp,zp) квадрата матрицы высот рельефа, который является последним в последовательности наибольшей продолжительности, получаются с учетом проведенных выше вычислений как:The indices (s p , z p ) of the square of the elevation matrix of the relief, which is the last in the sequence of the greatest duration, are obtained taking into account the above calculations as:
Очевидно, что изобретение не ограничивается вышеописанным примером его осуществления, исходя из него могут быть предусмотрены и другие варианты, не выходящие за рамки предмета изобретения, например выбор другой последовательности перебора квадратов и построение для нее матрицы управления вычислениями.Obviously, the invention is not limited to the above-described example of its implementation, on the basis of it other options may be provided that do not go beyond the scope of the subject invention, for example, selecting another sequence of enumerating squares and constructing a calculation control matrix for it.
Реализация предложенного способа технических трудностей не представляет и может быть выполнена на существующих типах ЭВМ.Implementation of the proposed method does not present technical difficulties and can be performed on existing types of computers.
Технико-экономическая оценка эффективности предлагаемого способа по сравнению с прототипом и другими известными способами проводилась экспериментальным методом. Алгоритм работы был реализован на ЭВМ.Feasibility study of the effectiveness of the proposed method compared with the prototype and other known methods was carried out by an experimental method. The algorithm of work was implemented on a computer.
Исследования показали, что применение предлагаемого способа по сравнению с прототипом позволяет в 2-5 раз сократить время построения расчетных зон затенения.Studies have shown that the application of the proposed method in comparison with the prototype allows 2-5 times to reduce the time it takes to build the estimated shading zones.
Источники информацииInformation sources
1. Шишкин А.В., Боресков А.В. Компьютерная графика. Полигонные модели. М.: ДИАЛОГ-МИФИ, 2000. - 464 с.1. Shishkin A.V., Boreskov A.V. Computer graphics. Polygon models. M .: DIALOGUE-MEPhI, 2000 .-- 464 p.
2. Сети мобильной связи. Частотно-территориальное планирование. Учебное пособие для вузов. / В.Ю.Бабков, М.А.Вознюк, П.А.Михайлов. - 2-е изд., испр. - М.: Горячая линия - Телеком, 2007. - 224 с.: ил.2. Mobile networks. Frequency-territorial planning. Textbook for universities. / V.Yu. Babkov, M.A. Voznyuk, P.A. Mikhailov. - 2nd ed., Rev. - M .: Hot line - Telecom, 2007 .-- 224 p.: Ill.
3. ГОСТ 21667-76. Картография.3. GOST 21667-76. Cartography.
4. ГОСТ 28441-90. Картография цифровая.4. GOST 28441-90. Digital cartography.
5. ГОСТ Р 50828-95. Геоинформационное картографирование. Пространственные данные, цифровые и электронные карты. Общие требования.5. GOST R 50828-95. Geoinformation mapping. Spatial data, digital and electronic maps. General requirements.
6. ГОСТ 51605-2000. Карты цифровые топографические. Общие требования.6. GOST 51605-2000. Digital topographic maps. General requirements.
7. ГОСТ 51606-2000. Карты цифровые топографические. Система классификации и кодирования цифровой картографической информации. Общие требования.7. GOST 51606-2000. Digital topographic maps. Classification and coding system for digital cartographic information. General requirements.
8. ГОСТ 51607-2000. Карты цифровые топографические. Правила цифрового описания картографической информации. Общие требования.8. GOST 51607-2000. Digital topographic maps. Rules for the digital description of cartographic information. General requirements.
9. Положение по информационному картографическому обеспечению топографической службы ВС РФ. - М.: ВТУ ГШ, 2001.9. Regulation on information cartographic support of the topographic service of the Armed Forces of the Russian Federation. - M.: VTU GSH, 2001.
10. ГОСТ Р 52055-2003. Геоинформационное картографирование. Пространственные модели местности. Общие требования.10. GOST R 52055-2003. Geoinformation mapping. Spatial terrain models. General requirements.
11. Справочник по радиорелейной связи. Каменский Н.Н., Модель A.M., Надененко Б.С. и др. Под ред. С.В.Бородича. - Изд-е 2-е, перераб. и доп. - М.: Радио и связь, 1981. - 416 с., ил.11. Handbook of radio relay communication. Kamensky N.N., Model A.M., Nadenenko B.S. et al. Ed. S.V. Borodich. -
12. Сети телевизионного и звукового ОВЧ ЧМ-вещания: Справочник / М.Г.Локшин, А.А.Шур, А.В.Кокорев, Р.А.Краснощеков. - М.: Радио и связь, 1988.12. Networks of television and sound VHF FM broadcasts: Reference book / M.G. Lokshin, A.A. Shur, A.V. Kokorev, R.A. Krasnoshchekov. - M.: Radio and Communications, 1988.
13. Буга Н.Н. и др. Электромагнитная совместимость радиоэлектронных средств: Учебное пособие для вузов. / Н.Н.Буга, В.Я.Конторович, В.И.Носов. Под ред. Н.Н.Буги. - М.: Радио и связь, 1993.13. Bug N.N. et al. Electromagnetic compatibility of electronic equipment: Textbook for universities. / N.N. Buga, V.Ya. Kontorovich, V.I. Nosov. Ed. N.N. Boogie. - M .: Radio and communications, 1993.
14. Рекомендация ITU-R Р.370-7. Графики распространения сигналов диапазона ОВЧ и УВЧ на частотах 30-1000 МГц.14. Recommendation ITU-R P.370-7. VHF and UHF signals distribution graphs at frequencies of 30-1000 MHz.
15. Рекомендация ITU-R Р.529-2. Методы прогнозирования для наземных сухопутных мобильных радиослужб в диапазонах ОВЧ и УВЧ.15. Recommendation ITU-R P.529-2. Prediction methods for terrestrial land mobile radio services in the VHF and UHF bands.
16. Рекомендация ITU-R PN.833-1. Ослабление растительностью.16. Recommendation ITU-R PN.833-1. Attenuation by vegetation.
17. Рекомендация ITU-R P.526-5. Распространение при дифракции.17. Recommendation ITU-R P.526-5. Propagation by diffraction.
Claims (1)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2011126018/08A RU2475852C1 (en) | 2011-06-24 | 2011-06-24 | Method of constructing calculated shadow zones by controlled search of relief height matrix elements |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| RU2011126018/08A RU2475852C1 (en) | 2011-06-24 | 2011-06-24 | Method of constructing calculated shadow zones by controlled search of relief height matrix elements |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| RU2011126018A RU2011126018A (en) | 2012-12-27 |
| RU2475852C1 true RU2475852C1 (en) | 2013-02-20 |
Family
ID=49121138
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| RU2011126018/08A RU2475852C1 (en) | 2011-06-24 | 2011-06-24 | Method of constructing calculated shadow zones by controlled search of relief height matrix elements |
Country Status (1)
| Country | Link |
|---|---|
| RU (1) | RU2475852C1 (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2182713C2 (en) * | 2000-03-28 | 2002-05-20 | Григорьев Владимир Григорьевич | Way of remote determination of coordinates of position of ground object |
| RU2364887C2 (en) * | 2007-09-25 | 2009-08-20 | ОАО "Уральское проектно-конструкторское бюро " Деталь" | Method for navigation of aircraft by radar images of earth surface with application of digital area models |
| TW201000934A (en) * | 2008-06-16 | 2010-01-01 | Chunghwa Telecom Co Ltd | Method of detecting traffic bottleneck points based on space environment data mining technology |
| JP2010068719A (en) * | 2008-09-16 | 2010-04-02 | Hitachi Software Eng Co Ltd | System for monitoring green space, and delivery method for monitoring green space |
-
2011
- 2011-06-24 RU RU2011126018/08A patent/RU2475852C1/en active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2182713C2 (en) * | 2000-03-28 | 2002-05-20 | Григорьев Владимир Григорьевич | Way of remote determination of coordinates of position of ground object |
| RU2364887C2 (en) * | 2007-09-25 | 2009-08-20 | ОАО "Уральское проектно-конструкторское бюро " Деталь" | Method for navigation of aircraft by radar images of earth surface with application of digital area models |
| TW201000934A (en) * | 2008-06-16 | 2010-01-01 | Chunghwa Telecom Co Ltd | Method of detecting traffic bottleneck points based on space environment data mining technology |
| JP2010068719A (en) * | 2008-09-16 | 2010-04-02 | Hitachi Software Eng Co Ltd | System for monitoring green space, and delivery method for monitoring green space |
Also Published As
| Publication number | Publication date |
|---|---|
| RU2011126018A (en) | 2012-12-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| KR102024695B1 (en) | System and method for modeling surveillance camera layout | |
| CN110782708A (en) | Unmanned aerial vehicle flight network modeling method based on low-altitude airspace limiting conditions | |
| Tian et al. | Modeling urban expansion policy scenarios using an agent-based approach for Guangzhou Metropolitan Region of China | |
| CN113012063B (en) | Dynamic point cloud repairing method and device and computer equipment | |
| CN113916130B (en) | Building position measuring method based on least square method | |
| CN110888132B (en) | Bridge deformation analysis method and system based on InSAR monitoring | |
| US7778808B2 (en) | Geospatial modeling system providing data thinning of geospatial data points and related methods | |
| Sahraoui et al. | Integrated GIS software for computing landscape visibility metrics | |
| US20190094196A1 (en) | Method and apparatus for analyzing communication environments and designing networks in consideration of trees | |
| CN112926465A (en) | Coastline property identification method and device based on point cloud type | |
| Cervilla et al. | Total 3D‐Viewshed Map: Quantifying the Visible Volume in Digital Elevation Models | |
| CN113822892B (en) | Evaluation method, device and equipment of simulated radar and computer storage medium | |
| Hao et al. | Integration of multinomial-logistic and Markov-chain models to derive land-use change dynamics | |
| RU2475852C1 (en) | Method of constructing calculated shadow zones by controlled search of relief height matrix elements | |
| CN117470246A (en) | Path planning method and device, storage medium and electronic equipment | |
| KR102401554B1 (en) | Method and device for generating features of spatial data | |
| CN117176918A (en) | Visual field determining method, device, electronic equipment and system of high-point monitoring equipment | |
| CN114543808B (en) | Indoor repositioning method, device, equipment and storage medium | |
| CN115311574A (en) | Building monitoring method, equipment and medium | |
| CN115271564A (en) | Highway slope disaster space danger zoning method and terminal | |
| Xu | [Retracted] Application of Remote Sensing Image Data Scene Generation Method in Smart City | |
| JP5071679B2 (en) | Inclination analyzer, method and program | |
| CN109407087B (en) | Target feature vector diagram matching method and system | |
| Achilleos et al. | Visibility and viewshed algorithms in an information system for environmental management | |
| Kupiec et al. | Visibility analysis as a tool in visual impact analysis procedures for environmental impact assessment in agricultural landscapes |