Visualize Coverage Maps over Lunar Terrain Using Ray Tracing
Display viewsheds and coverage maps for a base station on the surface of the Moon using terrain data and ray tracing analysis. A viewshed includes areas that are within line of sight (LOS) of a set of points, and a coverage map includes received power for the areas within range of a set of wireless transmitters. You can use a viewshed as a first-order approximation of a coverage map.
The goal of this example is to analyze radio coverage for a lunar module during several extra-vehicular activities (EVAs). During these EVAs, a moon rover travels from the lunar module to nearby hill clusters. The trajectories of the moon rover are based on the 1971 Apollo 15 mission, which was the first Apollo mission to use a moon rover. Because the moon rover enabled EVAs of longer distances than previous Apollo missions, maintaining communications between the rover, the module, and Earth was critical to the mission. During the mission, the rover frequently drove through areas with no radio coverage, inhibiting communication with the module.
When using geographic coordinates, the ray tracing capabilities in Antenna Toolbox™ and Communications Toolbox™ use an Earth-based region of interest (ROI). By combining the ray tracing capabilities with the mapping capabilities in Mapping Toolbox™, you can visualize coverage over generic planetary terrain.
This example shows you how to:
Read lunar terrain data and a reference spheroid from a GeoTIFF file.
Get a reference ellipsoid for the Moon.
Generate the trajectory of a moon rover.
View terrain, coverage maps, and viewsheds in 2-D.
View terrain and propagation paths in 3-D.
Read Terrain and Reference Spheroid
The Apollo 15 mission took place in a region surrounding the Hadley Plains. Read lunar terrain elevation data [1] for the Hadley Plains from a GeoTIFF file. The file reports elevations in meters with respect to a reference spheroid that represents the shape of the Moon. You can download the data in JPEG 2000 format from the Goddard Space Flight Center [2]. The file used in this example has been cropped from the original data set and exported to a GeoTIFF file.
filename = "SLDEM2015_Hadley_Plains_Apollo15.tif"; [Z,R] = readgeoraster(filename,OutputType="double");
The readgeoraster function imports the data into the workspace as an array and a raster reference object in geographic coordinates. The reference object stores information such as the latitude limits, the longitude limits, and the geographic coordinate reference system (CRS). Display the properties of the reference object.
disp(R)
GeographicCellsReference with properties:
LatitudeLimits: [24.6328125 27.630859375]
LongitudeLimits: [2.134765625 5.1328125]
RasterSize: [307 307]
RasterInterpretation: 'cells'
ColumnsStartFrom: 'north'
RowsStartFrom: 'west'
CellExtentInLatitude: 0.009765625
CellExtentInLongitude: 0.009765625
RasterExtentInLatitude: 2.998046875
RasterExtentInLongitude: 2.998046875
XIntrinsicLimits: [0.5 307.5]
YIntrinsicLimits: [0.5 307.5]
CoordinateSystemType: 'geographic'
GeographicCRS: [1×1 geocrs]
AngleUnit: 'degree'
The geographic CRS stores the reference spheroid for the Moon. You must use the reference spheroid to generate the trajectory of the moon rover, to perform the ray tracing analysis, and to create maps of the region. Extract the reference spheroid from the CRS.
moonSpheroid = R.GeographicCRS.Spheroid;
Define Lunar Module, Moon Rover, and Trajectory of Rover
Specify the latitude and longitude coordinates [3] of the lunar module. Approximate the height of the lunar module in meters.
moduleLat = 26.132; moduleLon = 3.634; moduleHeight = 7;
The lunar module uses the Unified-S band radio system for tracking and communication. Specify the frequency in hertz and the power in watts [4].
moduleFreq = 2.82e9; modulePower = 20;
Approximate the height of the moon rover in meters.
roverHeight = 1.5;
Approximate the trajectories of three EVAs: one to the Hadley Rille, another to the Mons Hadley Delta hill complex, and the last to the South Cluster hill complex. For each trajectory:
Specify the approximate coordinates of the rover destination.
Calculate the trajectory by finding the coordinates of points between the lunar module and the rover destination. Get points along a geodesic on the Moon by using the
track2function and specifying the Moon reference spheroid.
eva1Lat = 26.15; % 1st EVA eva1Lon = 3.56; [traj1Lat,traj1Lon] = track2(moduleLat,moduleLon,eva1Lat,eva1Lon,moonSpheroid); eva2Lat = 26.01; % 2nd EVA eva2Lon = 3.59; [traj2Lat,traj2Lon] = track2(moduleLat,moduleLon,eva2Lat,eva2Lon,moonSpheroid); eva3Lat = 25.98; % 3rd EVA eva3Lon = 3.697; [traj3Lat,traj3Lon] = track2(moduleLat,moduleLon,eva3Lat,eva3Lon,moonSpheroid);
To facilitate plotting, collect the rover destination and trajectory variables for each EVA into combined variables. Use one variable each for these attributes: destination latitudes, destination longitudes, destination heights, trajectory latitudes, trajectory longitudes, and trajectory heights.
evaLat = [eva1Lat; eva2Lat; eva3Lat]; evaLon = [eva1Lon; eva2Lon; eva3Lon]; evaHeight = ones(length(evaLat),1) * roverHeight; trajLat = [traj1Lat; traj2Lat; traj3Lat]; trajLon = [traj1Lon; traj2Lon; traj3Lon]; trajHeight = ones(length(trajLat),1) * roverHeight;
Specify the colors to use when plotting the trajectories, the lunar module, and the rover destinations. Use white for the trajectories, red for the lunar module, yellow for the rover destinations, and black for the marker outlines.
trajColor = "w"; moduleColor = "r"; roverColor = "y"; outlineColor = "k";
Visualize Terrain and Viewshed in 2-D
Visualize the terrain for the region and the viewshed from the lunar module by using axesm-based maps. axesm-based maps enable you to display geographic data using a map projection.
Visualize Terrain
Visualize the terrain for the region.
Create an axesm-based map by using an equidistant cylindrical projection and the Moon reference spheroid. The equidistant cylindrical projection is a reasonable choice for the region because it minimizes distortion near the equator. Set the latitude and longitude limits of the map using the latitude and longitude limits of the terrain data. Customize the appearance of the map by setting the positions of the meridian and parallel labels, specifying the number of significant digits to include in the labels, specifying the color and line style of the grid, and removing the frame. Reduce the white space around the map by using the tightmap function.
figure axesm("eqdcylin",Geoid=moonSpheroid, ... % projection and spheroid MapLatLimit=R.LatitudeLimits,MapLonLimit=R.LongitudeLimits, ... % latitude and longitude limits MeridianLabel="on",MLineLocation=0.5,MLabelParallel="south", ... % meridian label position ParallelLabel="on",PLineLocation=0.5, ... % parallel label position MlabelRound=-1,PlabelRound=-1, ... % significant digits for labels Grid="on",GColor=[0.8 0.8 0.8],GLineStyle="-.",GAltitude=Inf, ... % grid appearance Frame="off") % frame tightmap
Display the terrain data on a map using a surface. Add a labeled color bar.
geoshow(Z,R,DisplayType="texturemap") cb = colorbar; cb.Label.String = "Elevation in meters (with respect to reference spheroid)";
Display the trajectories on the map using white circle markers.
geoshow(trajLat,trajLon,DisplayType="point",ZData=trajHeight, ... Marker="o",MarkerSize=1,MarkerFaceColor=trajColor,MarkerEdgeColor=trajColor)
Display the lunar module on the map using a red circle marker. Prepare to include the module in a legend by specifying a display name and setting the plot object to a variable.
modulePt = geoshow(moduleLat,moduleLon,DisplayType="point",ZData=moduleHeight, ... Marker="o",MarkerSize=7.5,MarkerFaceColor=moduleColor, ... LineWidth=1.2,MarkerEdgeColor=outlineColor,DisplayName="Lunar Module");
Display the rover destinations on the map using yellow triangle markers. Prepare to include the destinations in the legend by specifying a display name and setting the plot object to a variable.
roverPt = geoshow(evaLat,evaLon,DisplayType="point",ZData=evaHeight, ... Marker="^",MarkerSize=7.5,MarkerFaceColor=roverColor, ... LineWidth=1.2,MarkerEdgeColor=outlineColor,DisplayName="Rover");
Add a title. Include the lunar module and the rover destinations in a legend.
title("Lunar Terrain for Hadley Plains")
legend([modulePt roverPt])
Visualize Terrain for ROI
Visualize the terrain for only the area surrounding the lunar module by cropping the terrain data to a smaller region.
Specify an ROI by finding the coordinates of points that are 5 kilometers north, east, south, and west of the lunar module. Calculate the latitude and longitude limits of the ROI.
[boundinglat,boundinglon] = reckon(moduleLat,moduleLon,5000,[0; 90; 180; 270],moonSpheroid); [latlim,lonlim] = geoquadpt(boundinglat,boundinglon);
Crop the terrain elevation data to the ROI. Resize the cropped data so that the array has 200 rows and 200 columns.
[ZROI,RROI] = geocrop(Z,R,latlim,lonlim); gridSize = [200 200]; latscale = gridSize(1)/RROI.RasterSize(1); lonscale = gridSize(2)/RROI.RasterSize(2); [ZROI,RROI] = georesize(ZROI,RROI,latscale,lonscale);
Create a new axesm-based map. Then, display the terrain data on the map.
figure axesm("eqdcylin",Geoid=moonSpheroid, ... % projection and spheroid MapLatLimit=RROI.LatitudeLimits,MapLonLimit=RROI.LongitudeLimits, ... % latitude and longitude limits MeridianLabel="on",MLineLocation=0.1,MLabelParallel="south", ... % meridian label position ParallelLabel="on",PLineLocation=0.1, ... % parallel label position MlabelRound=-1,PlabelRound=-1, ... % significant digits for labels Grid="on",GColor=[0.8 0.8 0.8],GLineStyle="-.",GAltitude=Inf, ... % grid appearance Frame="off") % frame tightmap geoshow(ZROI,RROI,DisplayType="texturemap")
Specify limits for the colormap using the minimum and maximum elevation values. Add a labeled color bar.
[cmin,cmax] = bounds(ZROI,"all"); clim([cmin cmax]) cb = colorbar; cb.Label.String = "Elevation in meters (with respect to reference spheroid)";
Display the trajectories, the lunar module, and the rover destinations on the map.
geoshow(trajLat,trajLon,DisplayType="point",ZData=trajHeight, ... Marker="o",MarkerSize=2,MarkerFaceColor=trajColor,MarkerEdgeColor=trajColor) modulePt = geoshow(moduleLat,moduleLon,DisplayType="point",ZData=moduleHeight, ... Marker="o",MarkerSize=12,MarkerFaceColor=moduleColor, ... LineWidth=1.2,MarkerEdgeColor=outlineColor,DisplayName="Lunar Module"); roverPt = geoshow(evaLat,evaLon,DisplayType="point",ZData=evaHeight, ... Marker="^",MarkerSize=12,MarkerFaceColor=roverColor, ... LineWidth=1.2,MarkerEdgeColor=outlineColor,DisplayName="Rover");
Add a title. Include the lunar module and the rover destinations in a legend.
legend([modulePt roverPt])
title("Lunar Terrain for ROI")
Visualize Viewshed for ROI
Calculate the viewshed for the lunar module. Indicate that the heights of the lunar module and the grid coordinates are referenced to the terrain by specifying the height references as "AGL". Specify the radius of the reference sphere using the Moon spheroid.
visBase = viewshed(ZROI,RROI,moduleLat,moduleLon, ... moduleHeight,roverHeight,"AGL","AGL",moonSpheroid.Radius);
The visBase array returned by the viewshed function indicates visible areas using values of 1 and obscured areas using values of 0. To avoid plotting the obscured areas, replace the 0 values with NaN values. Define a colormap for the viewshed that uses turquoise for visible areas.
visBase(visBase == 0) = NaN; covcolor = [0.1034 0.8960 0.7150];
Create a new axesm-based map. Then, display the terrain data on the map.
figure axesm("eqdcylin",Geoid=moonSpheroid, ... % projection and spheroid MapLatLimit=RROI.LatitudeLimits,MapLonLimit=RROI.LongitudeLimits, ... % latitude and longitude limits MeridianLabel="on",MLineLocation=0.1,MLabelParallel="south", ... % meridian label position ParallelLabel="on",PLineLocation=0.1, ... % parallel label position MlabelRound=-1,PlabelRound=-1, ... % significant digits for labels Grid="on",GColor=[0.8 0.8 0.8],GLineStyle="-.",GAltitude=Inf, ... % grid appearance Frame="off") % frame tightmap geoshow(ZROI,RROI,DisplayType="texturemap")
Apply a gray colormap to the plotted terrain data. Specify the limits for the colormap using the minimum and maximum elevation values. Add a labeled color bar.
colormap("gray") [cmin,cmax] = bounds(ZROI,"all"); clim([cmin cmax]) cb = colorbar; cb.Label.String = "Elevation in meters (with respect to reference spheroid)";
Display the viewshed for the lunar module over the terrain data. Prepare to include the viewshed in a legend by specifying a display name and setting the plot object to a variable.
viewshedSurf = geoshow(visBase,RROI,EdgeColor="none", ... DisplayType="surface",FaceColor=covcolor,DisplayName="Viewshed");
Display the trajectories, the lunar module, and the rover destinations on the map.
geoshow(trajLat,trajLon,DisplayType="point",ZData=trajHeight, ... Marker="o",MarkerSize=2,MarkerFaceColor=trajColor,MarkerEdgeColor=trajColor) modulePt = geoshow(moduleLat,moduleLon,DisplayType="point",ZData=moduleHeight, ... Marker="o",MarkerSize=12,MarkerFaceColor=moduleColor, ... LineWidth=1.2,MarkerEdgeColor=outlineColor,DisplayName="Lunar Module"); roverPt = geoshow(evaLat,evaLon,DisplayType="point",ZData=evaHeight, ... Marker="^",MarkerSize=12,MarkerFaceColor=roverColor, ... LineWidth=1.2,MarkerEdgeColor=outlineColor,DisplayName="Rover");
Add a title. Include the viewshed, the lunar module, and the rover destinations in a legend.
legend([viewshedSurf modulePt roverPt])
title("Viewshed of Apollo 15 Lunar Module")
Visualize Coverage in 2-D
A viewshed is comparable to a coverage map that includes only line-of-sight (LOS) propagation. You can create more detailed coverage maps and calculate received power by using antenna sites and ray tracing propagation models.
Create Antenna Sites
Create a base station on top of the lunar module. Specify the location of the base station using geographic coordinates, and transform the coordinates to local ENU coordinates by using the helperLatLonToENU helper function. Represent the base station using a transmitter site. By default, transmitter sites use isotropic antennas.
txlat = moduleLat; txlon = moduleLon; originlat = moduleLat; originlon = moduleLon; originht = 0; [txx,txy,txz] = helperLatLonToENU(Z,R,txlat,txlon,originlat,originlon,originht,moonSpheroid); txz = txz + moduleHeight; % place antenna on top of lunar module tx = txsite("cartesian", ... AntennaPosition=[txx; txy; txz], ... TransmitterFrequency=moduleFreq, ... TransmitterPower=modulePower);
Extract a grid of geographic coordinates from the reference object for the ROI. Then, prepare to create a coverage map for the ROI by creating a receiver site for each point in the grid of geographic coordinates. Transform the coordinates to local ENU coordinates by using the helperLatLonToENU helper function. Specify the height values so that each receiver site sits on top of the moon rover.
[latROI,lonROI] = geographicGrid(RROI); [rxx,rxy,rxz] = helperLatLonToENU(Z,R,latROI(:),lonROI(:),originlat,originlon,originht,moonSpheroid); rxz = rxz + roverHeight; % place antennas on top of moon rover rxs = rxsite("cartesian",AntennaPosition=[rxx rxy rxz]');
Create Triangulation of Lunar Surface
Convert the terrain elevation data to a triangulation object. A triangulation object enables you to create coverage maps using ray tracing propagation model objects.
To use the ray tracing propagation model with a triangulation object, you must use a triangulation object referenced to Cartesian coordinates. Extract a grid of geographic coordinates from the reference object for the terrain data.
[gridlat,gridlon] = geographicGrid(R);
Convert the coordinates to Cartesian coordinates in an east-north-up (ENU) local coordinate system. Specify the origin of the ENU system using the coordinates of the lunar module.
[x,y,z] = geodetic2enu(gridlat,gridlon,Z,originlat,originlon,originht,moonSpheroid);
Convert the Cartesian coordinates to a patch with triangular faces. The surf2patch function represents the patch using two variables: V contains the patch vertices, and F determines which vertices to connect to form each triangular face. Create a triangulation object from the patch.
[F,V] = surf2patch(x,y,z,"triangles");
tri = triangulation(F,V);Create Ray Tracing Propagation Model
Create a ray tracing propagation model, which MATLAB® represents using a RayTracing object. Specify the surface material using the real relative permittivity and the conductivity of the Moon surface [5]. By default, the model uses the shooting-and-bouncing rays (SBR) method.
epsR = 3.7; % real relative permittivity of Moon surface lossTan = 0.1; % dielectric loss tangent epsilon = 8.8541878128e-12; % permittivity of free space (vacuum) cond = lossTan*(2*pi*moduleFreq)*epsR*epsilon; % conductivity of Moon surface pm = propagationModel("raytracing",CoordinateSystem="cartesian", ... SurfaceMaterial="custom",SurfaceMaterialPermittivity=epsR,SurfaceMaterialConductivity=cond);
Visualize Line-of-Sight Coverage
Configure the propagation model to find LOS paths. Then, calculate the received power, in dBm, at each receiver site. Specify the terrain using the triangulation object.
pm.MaxNumReflections = 0; sslos = sigstrength(rxs,tx,pm,Map=tri);
Prepare to display the coverage map.
Specify power levels of interest between –120 dBm and –50 dBm.
Define a colormap that uses 5 colors per 10 dBm.
Specify colormap limits.
Reshape the power values into a grid, and replace values below the minimum power threshold with
NaNvalues.
levels = -120:10:-50; cmap = turbo(5*length(levels)); climLevels = [min(levels)+20 max(levels)+30]; ss = sslos; ss = reshape(ss,gridSize); ss (ss < min(levels)) = NaN;
To display the terrain data and the received power values on the same map, use a tiled chart layout with one tile.
figure t = tiledlayout(1,1);
Create an axesm-based map in the tile. Display the terrain data, and apply a gray colormap.
ax1 = axes(t); axesm("eqdcylin",Geoid=moonSpheroid, ... % projection and spheroid MapLatLimit=RROI.LatitudeLimits,MapLonLimit=RROI.LongitudeLimits, ... % latitude and longitude limits MeridianLabel="on",MLineLocation=0.1,MLabelParallel="south", ... % meridian label position ParallelLabel="on",PLineLocation=0.1, ... % parallel label position MlabelRound=-1,PlabelRound=-1, ... % significant digits for labels Grid="on",GColor=[0.8 0.8 0.8],GLineStyle="-.",GAltitude=Inf, ... % grid appearance Frame="off") % frame tightmap geoshow(ax1,ZROI,RROI,DisplayType="texturemap") colormap(ax1,"gray")
Create a second axesm-based map in the same tile. Use the same projection, reference spheroid, and limits as the first map.
ax2 = axes(t); abm2 = axesm("eqdcylin",Frame="off",Geoid=moonSpheroid, ... MapLatLimit=RROI.LatitudeLimits,MapLonLimit=RROI.LongitudeLimits); abm2.Visible = "off"; tightmap
Display the received power values as a surface. Apply the power level colormap and set the color limits. Add a labeled color bar.
geoshow(ax2,ss,RROI,DisplayType="surface") colormap(ax2,cmap) clim(ax2,climLevels) cb = colorbar; cb.Label.String = "Power (dBm)";
Display the trajectories, the lunar module, and the rover destinations in the second map.
geoshow(trajLat,trajLon,DisplayType="point",ZData=trajHeight, ... Marker="o",MarkerSize=2,MarkerFaceColor=trajColor,MarkerEdgeColor=trajColor) modulePt = geoshow(moduleLat,moduleLon,DisplayType="point",ZData=moduleHeight, ... Marker="o",MarkerSize=12,MarkerFaceColor=moduleColor, ... LineWidth=1.2,MarkerEdgeColor=outlineColor,DisplayName="Lunar Module"); roverPt = geoshow(evaLat,evaLon,DisplayType="point",ZData=evaHeight, ... Marker="^",MarkerSize=12,MarkerFaceColor=roverColor, ... LineWidth=1.2,MarkerEdgeColor=outlineColor,DisplayName="Rover");
Link the axes objects that are stored in the tiles so that the axes objects share the same axis limits.
linkaxes([ax1 ax2],"xy")Add a title and subtitle. Include the lunar module and the rover destinations in a legend.
title("Coverage Map for Apollo 15 Lunar Module",Visible="on") subtitle("LOS Paths",Visible="on") legend([modulePt roverPt])
Visualize Coverage with One Reflection
Configure the propagation model to find paths with up to one surface reflection. Then, recalculate the received power, in dBm, at each receiver site.
pm.MaxNumReflections = 1; ssref = sigstrength(rxs,tx,pm,Map=tri);
Prepare to display the coverage map by reshaping the power values into a grid and by replacing values below the minimum threshold with NaN values.
ss = ssref; ss = reshape(ss,gridSize); ss (ss < min(levels)) = NaN;
Create a new tiled chart layout with one tile.
figure t = tiledlayout(1,1);
Create an axesm-based map in the tile. Display the terrain data and apply a gray colormap.
ax1 = axes(t); axesm("eqdcylin",Geoid=moonSpheroid, ... % projection and spheroid MapLatLimit=RROI.LatitudeLimits,MapLonLimit=RROI.LongitudeLimits, ... % latitude and longitude limits MeridianLabel="on",MLineLocation=0.1,MLabelParallel="south", ... % meridian label position ParallelLabel="on",PLineLocation=0.1, ... % parallel label position MlabelRound=-1,PlabelRound=-1, ... % significant digits for labels Grid="on",GColor=[0.8 0.8 0.8],GLineStyle="-.",GAltitude=Inf, ... % grid appearance Frame="off") % frame tightmap geoshow(ax1,ZROI,RROI,DisplayType="texturemap") colormap(ax1,"gray")
Create a second axesm-based map in the same tile. Use the same projection, reference spheroid, and limits as the first map.
ax2 = axes(t); abm2 = axesm("eqdcylin",Frame="off",Geoid=moonSpheroid, ... MapLatLimit=RROI.LatitudeLimits,MapLonLimit=RROI.LongitudeLimits); abm2.Visible = "off"; tightmap
Display the received power values as a surface. Apply the power level colormap and set the color limits. Add a labeled color bar.
geoshow(ax2,ss,RROI,DisplayType="surface") colormap(ax2,cmap) clim(ax2,climLevels) cb = colorbar; cb.Label.String = "Power (dBm)";
Display the trajectories, the lunar module, and the rover destinations in the second tile.
geoshow(trajLat,trajLon,DisplayType="point",ZData=trajHeight, ... Marker="o",MarkerSize=2,MarkerFaceColor=trajColor,MarkerEdgeColor=trajColor) modulePt = geoshow(moduleLat,moduleLon,DisplayType="point",ZData=moduleHeight, ... Marker="o",MarkerSize=12,MarkerFaceColor=moduleColor, ... LineWidth=1.2,MarkerEdgeColor=outlineColor,DisplayName="Lunar Module"); roverPt = geoshow(evaLat,evaLon,DisplayType="point",ZData=evaHeight, ... Marker="^",MarkerSize=12,MarkerFaceColor=roverColor, ... LineWidth=1.2,MarkerEdgeColor=outlineColor,DisplayName="Rover");
Link the axes objects that are stored in the tiles so that the axes objects share the same axis limits.
linkaxes([ax1 ax2],"xy")Add a title and subtitle. Include the lunar module and the rover destinations in a legend.
title("Coverage Map for Apollo 15 Lunar Module",Visible="on") subtitle("LOS Paths and Paths with 1 Reflection",Visible="on") legend([modulePt roverPt])
While the addition of paths with one reflection does not change the general shape of the coverage area, the interference between the line-of-sight paths and the reflected paths cause speckles in the power profile.
Visualize Terrain and Propagation Paths in 3-D
Visualize the terrain and propagation paths in 3-D using Site Viewer.
Import the triangulation object that represents the terrain model into Site Viewer. View the terrain for the region by zooming out and rotating the model.
sv = siteviewer(SceneModel=tri,ShowOrigin=false);
Create a receiver site that represents the destination for the EVA to the Mons Hadley Delta hill complex. Convert the geographic coordinates to local ENU coordinates by using the helperLatLonToENU function.
[rxx,rxy,rxz] = helperLatLonToENU(Z,R,eva2Lat,eva2Lon,originlat,originlon,originht,moonSpheroid);
rxz = rxz + roverHeight;
rx = rxsite("cartesian",AntennaPosition=[rxx rxy rxz]');Display the transmitter site, the receiver site, and the LOS visibility status from the lunar module to the rover destination. The green line indicates that the rover is within LOS of the lunar module.
show(tx) show(rx) los(tx,rx)
Calculate and display the LOS and single-reflection propagation paths between transmitter and the receiver.
raytrace(tx,rx,pm,ColorLimits=[-100 -5])
References
[1] Barker, M.K., E. Mazarico, G.A. Neumann, M.T. Zuber, J. Haruyama, and D.E. Smith. "A New Lunar Digital Elevation Model from the Lunar Orbiter Laser Altimeter and SELENE Terrain Camera." Icarus 273 (July 2016): 346–55. https://doi.org/10.1016/j.icarus.2015.07.039.
[2] NASA's Planetary Geology, Geophysics and Geochemistry Laboratory. "PGDA - High-resolution Lunar Topography (SLDEM2015)." NASA Goddard Space Flight Center. Accessed August 28, 2023. https://pgda.gsfc.nasa.gov/products/54.
[3] NASA National Space Science Data Center. "Apollo 15 Lunar Module /ALSEP." Accessed September 1, 2023. https://nssdc.gsfc.nasa.gov/nmc/spacecraft/display.action?id=1971-063C.
[4] Peltzer, K.E. “Apollo Unified S-Band System.” Technical Memorandum (TM). Greenbelt, Maryland: NASA Goddard Space Flight Center, April 1, 1966. https://ntrs.nasa.gov/citations/19660018739.
[5] Heiken, Grant, David Vaniman, and Bevan M. French, eds. Lunar Sourcebook: A User’s Guide to the Moon. Cambridge [England]; New York: Cambridge University Press, 1991.
Helper Functions
The helperLatLonToENU helper function converts the geographic coordinates lat and lon along the terrain specified by Zgeo and Rgeo to local ENU coordinates. The input arguments lat0, lon0, and h0 specify the origin of the ENU system. The input argument spheroid specifies the reference spheroid for the geographic coordinates.
function [x,y,z] = helperLatLonToENU(Zgeo,Rgeo,lat,lon,lat0,lon0,h0,spheroid) % Convert geographic coordinates to local ENU coordinates h = geointerp(Zgeo,Rgeo,lat,lon); [x,y,z] = geodetic2enu(lat,lon,h,lat0,lon0,h0,spheroid); end
See Also
Functions
readgeoraster(Mapping Toolbox) |geodetic2enu(Mapping Toolbox) |sigstrength|raytrace
Objects
siteviewer(Antenna Toolbox) |RayTracing(Antenna Toolbox)
Related Topics
- Visualize Viewsheds and Coverage Maps Using Terrain (Mapping Toolbox)
You can also select a web site from the following list:
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)