Best way to calculate elliptical trajectory
2 次查看(过去 30 天)
显示 更早的评论
How do I calculate the ellptical trajectory on a spheroid at constant altitude if I know the center of ellipse and boundary points of imaginary rectangle?
% earth ref spheroid
wgs84 = wgs84Ellipsoid;
% Aircraft location
lat_aircraft = 45; %deg
lon_aircraft = 8; %deg
alt_aircraft = 5000; %m
%% Receiver extent
lon_extent = 0.5;
lat_extent = 1;
%% Find boundry points of rectangle, and add center of the location
lat_range = [lat_aircraft - (lat_extent/2), lat_aircraft, lat_aircraft + (lat_extent/2)];
lon_range = [lon_aircraft - (lon_extent/2), lon_aircraft, lon_aircraft + (lon_extent/2)];
latlonPts = [combvec(lat_range, lon_range)]';
geoplot(latlonPts(:,1), latlonPts(:,2),'*')
%% Create elliptical traj at constant altitude
0 个评论
回答(1 个)
Adam Danz
2021-6-1
编辑:Adam Danz
2021-6-2
Existing code
% earth ref spheroid
wgs84 = wgs84Ellipsoid;
% Aircraft location
lat_aircraft = 45; %deg
lon_aircraft = 8; %deg
alt_aircraft = 5000; %m
%% Receiver extent
lon_extent = 0.5;
lat_extent = 1;
%% Find boundry points of rectangle, and add center of the location
lat_range = [lat_aircraft - (lat_extent/2), lat_aircraft, lat_aircraft + (lat_extent/2)];
lon_range = [lon_aircraft - (lon_extent/2), lon_aircraft, lon_aircraft + (lon_extent/2)];
latlonPts = [combvec(lat_range, lon_range)]';
geoplot(latlonPts(:,1), latlonPts(:,2),'*')
Create ellipse within bounding box
%% Create elliptical traj at constant altitude
theta = linspace(0,2*pi) ;
x = lon_extent/2 * cos(theta) + lon_aircraft;
y = lat_extent/2 * sin(theta) + lat_aircraft;
hold on
geoplot(y,x,'r-')
1 个评论
Adam Danz
2021-6-3
Note, the simple solution above works when the rectangle is parallel with the axes but does not work with rotated rectangles.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Coordinate Systems 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!