-
Notifications
You must be signed in to change notification settings - Fork 56
LRO LOLA and LROC Example
NASA's Lunar Reconnaissance Orbiter (LRO) carries the Lunar Orbiter Laser Altimeter (LOLA) and the Lunar Reconnaissance Orbiter Camera (LROC). These instruments provide the data products for a lunar digital terrain model.
Data Credit: NASA/GSFC/ASU
1.) Download the following data products:
NAC_ROI_GDANOBRUHIA_E360N1028.IMG
2.) Create a .prj by running gdalinfo on the orthoimage and copying and pasting the projection into a file called ortho.prj. The file should look like this:
PROJCS["EQUIRECTANGULAR MOON",
GEOGCS["GCS_MOON",
DATUM["D_MOON",
SPHEROID["MOON_localRadius",1737400,0]],
PRIMEM["Reference_Meridian",0],
UNIT["degree",0.0174532925199433]],
PROJECTION["Equirectangular"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",102.8],
PARAMETER["standard_parallel_1",36],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0]]
3.) Reproject the DEM to the orthoimage projection and resample it so each pixel is 16 times the size of an orthoimage pixel.
gdalwarp -t_srs ortho.prj -r bilinear -tr 19.2 19.2 ldem_1024_30n_45n_090_120.lbl ldem_1024_30n_45n_090_120.tif
4.) Crop the DEM to cover the same physical area as the orthoimage (use gdalinfo to find the bounds of the orthoimage):
gdal_translate -mo "AREA_OR_POINT=Point" -mo "SCALE=0.5" -strict -projwin -28477.2 1117780.8 28336.8 1063191.6 ldem_1024_30n_45n_090_120.tif ldem_1024_30n_45n_090_120_cropped.tif
5.) Create a landscape directory called GiordanoBrunoCrater.
6.) Create the elevation layer:
layerfactory -type=elevation -file=ldem_1024_30n_45n_090_120_cropped.tif -landscape=GiordanoBrunoCrater -tilesize=128 -globe=moon
7.) Create the orthoimage layer:
layerfactory -type=grayimage -file=NAC_ROI_GDANOBRUHIA_E360N1028.IMG -landscape=GiordanoBrunoCrater -tilesize=2048 -globe=moon
8.) View the landscape with DERT.