Видимость угла обзора между двумя точками в ландшафте
vis = los2(Z,R,lat1,lon1,lat2,lon2)
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1)
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1,alt2)
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1,alt2,alt1opt)
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1,alt2,alt1opt,alt2opt)
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1,alt2,alt1opt,
...
alt2opt,actualradius)
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1,alt2,alt1opt,
...
alt2opt,actualradius,effectiveradius)
[vis,visprofile,dist,H,lattrk,lontrk]
= los2(...)
los2(...)
los2 вычисляет взаимную видимость между двумя точками на отображенной цифровой карте повышения. los2 использует текущий объект, если это - обычная сетка данных или первая обычная сетка данных, найденная на текущей системе координат. zdata сетки используется для профиля. Цветные данные используются в отсутствие данных в z. Две точки выбраны путем нажатия на карту. Результат отображен в новой фигуре. Маркеры указывают на видимые и затененные точки вдоль профиля. Профиль показывают в Декартовой системе координат с источником в местоположении наблюдателя. Отображенная координата z составляет повышение ландшафта и искривление тела.
vis = los2(Z,R,lat1,lon1,lat2,lon2) вычисляет взаимную видимость между парами точек на цифровой карте повышения. Повышения обеспечиваются как обычная сетка данных Z, содержащий повышения в модулях метров. Две точки обеспечиваются как векторы широт и долгот в модулях степеней. Получившийся vis логической переменной равен тому, когда две точки видимы друг другу и нулю, когда угол обзора затенен ландшафтом. Если какой-либо из входных параметров пуст, los2 пытается собрать данные из текущей системы координат. С одним или несколькими выходными аргументами не создаются никакие фигуры, и только данные возвращены.
R может быть географическим растровым объектом ссылки, вектором ссылки или матрицей ссылки. Если R является географическим растровым объектом ссылки, его свойство RasterSize должно быть сопоставимо с size(Z).
Если R является вектором ссылки, это должно быть 1 3 с элементами:
[cells/degree northern_latitude_limit western_longitude_limit]
Если R является матрицей ссылки, это должно быть 3 2 и преобразовать растровые индексы строки и столбца к или от географических координат согласно:
[lon lat] = [row col 1] * R
Если R является матрицей ссылки, он должен задать (невращательный, нескошенный) отношение, в котором каждый столбец сетки данных падает вдоль меридиана и каждой строки падения вдоль параллели. Интерполяция ближайшего соседа используется по умолчанию. NaN возвращен для точек вне пределов сетки или для которого lat или lon содержат NaN. Все углы находятся в модулях степеней.
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1) помещает первую точку на заданной высоте в метрах выше поверхности (на башне, например). Это эквивалентно помещению точки на башне. Если не использовано, точка 1 принята, чтобы быть на поверхности. alt1 может быть или скаляром или вектором с той же длиной как lat1, lon1, lat2 и lon2.
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1,alt2) места обе точки в заданные высоты в метрах выше поверхности. alt2 может быть или скаляром или вектором с той же длиной как lat1, lon1, lat2 и lon2. Если alt2 не использован, точка 2 принята, чтобы быть на поверхности.
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1,alt2,alt1opt) управляет интерпретацией alt1 как любой относительная высота (alt1opt равняется 'AGL', значение по умолчанию) или абсолютная высота (alt1opt равняется 'MSL'). Если высотной опцией является 'AGL', alt1 интерпретирован как высота точки 1 в метрах выше ландшафта (“над уровнем земли”). Если alt1opt является 'MSL', alt1 интерпретирован как высота выше нуля (“средний уровень моря”).
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1,alt2,alt1opt,alt2opt) управляет интерпретацией ALT2.
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1,alt2,alt1opt,
... делает вычисление видимости на сфере с заданным радиусом. Если не использовано, радиус земли в метрах принят. Высоты, повышения и радиус должны быть в тех же модулях. Эта форма вызова является самой полезной для вычислений на телах кроме земли.
alt2opt,actualradius)
vis = los2(Z,R,lat1,lon1,lat2,lon2,alt1,alt2,alt1opt,
... принимает больший радиус для распространения угла обзора. Это может составлять искривление пути прохождения сигнала из-за преломления в атмосфере. Например, радио-распространение в атмосфере обычно обрабатывается как распространение прямой линии на сфере с 4/3 радиус земли. В этом случае последними двумя аргументами был бы
alt2opt,actualradius,effectiveradius)R_e и 4/3*R_e, где R_e является радиусом земли. Используйте Inf в качестве эффективного радиуса для плоских наземных вычислений видимости. Высоты, повышения и радиусы должны быть в тех же модулях.
[vis,visprofile,dist,H,lattrk,lontrk]
= los2(...), для скалярных входных параметров (lat1, lon1, и т.д.), возвращает векторы точек вдоль пути между двумя точками. visprofile является логическим вектором, содержащим верный (logical(1)), где промежуточные точки видимы и ложь (logical(0)) в противном случае. dist является расстоянием вдоль пути (в метрах или модулях радиуса). H содержит профиль ландшафта относительно вертикальной данной величины вдоль пути. lattrk и lontrk являются широтами и долготами точек вдоль пути. Для векторных входных параметров los2 возвращает visprofile, dist, H, lattrk и lontrk как массивы ячеек, с одной ячейкой на элемент lat1,lon1, и т.д.
los2(...), без выходных аргументов, отображает профиль видимости между двумя точками в новой фигуре.
Найдите угол повышения точки 90 степенями от наблюдателя, принимающего, что наблюдатель и цель оба на 1 000 км выше Земли.
lat1 = 0;
lon1 = 0;
alt1 = 1000*1000;
lat2 = 0;
lon2 = 90;
alt2 = 1000*1000;
[az, elev, r] = geodetic2aer(lat2,lon2,alt2,lat1,lon1,alt1,referenceEllipsoid('grs 80'))az = 90
elev = -45
r = 1.0434e+07
Визуально проверяйте результат с помощью функции угла обзора los2. Создайте сетку данных нулей, чтобы представлять поверхность Земли. Функция los2 без выходных аргументов создает фигуру, отображающую геометрию.
Z = zeros(181,361); R = georefpostings([-90 90],[-180 180], size(Z))
R =
GeographicPostingsReference with properties:
LatitudeLimits: [-90 90]
LongitudeLimits: [-180 180]
RasterSize: [181 361]
RasterInterpretation: 'postings'
ColumnsStartFrom: 'south'
RowsStartFrom: 'west'
SampleSpacingInLatitude: 1
SampleSpacingInLongitude: 1
RasterExtentInLatitude: 180
RasterExtentInLongitude: 360
XIntrinsicLimits: [1 361]
YIntrinsicLimits: [1 181]
CoordinateSystemType: 'geographic'
AngleUnit: 'degree'
los2(Z,R,lat1,lon1,lat2,lon2,alt1,alt1);
Z = 500*peaks(100); refvec = [1000 0 0]; [lat1, lon1, lat2, lon2] = deal(-0.027, 0.05, -0.093, 0.042); los2(Z,refvec,lat1,lon1,lat2,lon2,100);
figure;
axesm('globe','geoid',earthRadius('meters'))
meshm(Z, refvec, size(Z), Z); axis tight
camposm(-10,-10,1e6); camupm(0,0)
demcmap('inc', Z, 1000); shading interp; camlight
[vis,visprofile,dist,h,lattrk,lontrk] = ...
los2(Z,refvec,lat1,lon1,lat2,lon2,100);
plot3m(lattrk([1;end]),lontrk([1; end]),...
h([1; end])+[100; 0],'r','linewidth',2)
plotm(lattrk(~visprofile),lontrk(~visprofile),...
h(~visprofile),'r.','markersize',10)
plotm(lattrk(visprofile),lontrk(visprofile),...
h(visprofile),'g.','markersize',10)
