由地面两点计算是否通视

zq_6389

于 2023-04-25 20:00:20 发布

阅读量208

点赞数

文章标签：笔记经验分享

本文链接：https://blog.csdn.net/zhengquanpu/article/details/130372115

版权

由高度角和发射机和接收机之间的角度进行比较，得出是否能够通视

function [isvis, canShowObstruction, obsLoc, obsEl] = computeLOSUsingTerrainProfile(observerSite, observerCoords, criticalAngle, lats, lons, map)

% Initialize outputs to correspond to no-obstruction case
isvis = true;
canShowObstruction = false;
obsLoc = [];
obsEl = [];

% Return early if not enough data points. Analysis requires at least three
% points: two end points and one middle point.
if numel(lats) <= 2
    return
end

% Compute visibility status by looking for first obstruction along terrain
% between observer and target, which is located where elevation angle from
% observer is greater than angle between observer and target. 通过查找沿地形的第一个障碍物来计算可见性状态在观察点和目标之间，该仰角大于观察点和目标之间的角度。
rxprofile = rxsite('Name', 'internal.terrainsamplesite', ... % Specify to avoid default site naming
    'AntennaHeight', 0, ...
    'Latitude', lats(2:end-1), ...
    'Longitude', lons(2:end-1));
[~, el] = angle(observerSite, rxprofile, 'euclidean', 'Map', map, ...
    'SourceAntennaSiteCoordinates', observerCoords);
obstructionInd = find(el >= criticalAngle, 1);
isvis = isempty(obstructionInd);

% If first sample is obstruction, assume there is complete obstruction
% (i.e. no point to show as visible)
canShowObstruction = obstructionInd ~= 1;

% Compute obstruction point
if ~isvis
    obstructingSite = rxprofile(obstructionInd);
    obstruction = rfprop.internal.AntennaSiteCoordinates.createFromAntennaSites(...
        obstructingSite, map);
    obsLoc = obstruction.DisplayLatitudeLongitude;
    obsEl = obstruction.SurfaceHeightAboveTerrainReference;
end
end