Titan Data Atlas

Map viewer

Worked example

Georeferencing one MoonKAM frame

The Titan samples only needed a format change. A raw camera frame needs more: it has no coordinates until you give it some. Here is one frame from the GRAIL MoonKAM cameras, tied to a reference map of the Moon by hand.

ApproximateMoon, sphere of radius 1,737.4 km

The original MoonKAM frame, before georeferencing
The original MoonKAM frame, before georeferencing
Side by side: the frame with numbered control points, and the same points on the LRO WAC mosaic
Side by side: the frame with numbered control points, and the same points on the LRO WAC mosaic

On the map

The warped frame on a plain longitude and latitude grid of the Moon.

Layers

Move over the map to read coordinates

Do it yourself, step by step

How this frame was placed on a map of the Moon by hand. You can repeat it in QGIS or ArcGIS Pro with the files on this page. Each step is one action.

What you need

Step 1. Choose a frame

  1. Start from a set of MoonKAM frames sorted by date.
  2. Skip frames that show black sky or a curved horizon. Those views are very tilted.
  3. Skip frames with stripes, colored noise or blocks. Those are transmission errors.
  4. Choose a frame with one large, clear feature, such as a crater with terraces.
  5. Write down the date and time from the file name.

The file name GRAIL_A_CAM1_0WEY0WH100000_2012-089T14 38 20.494.bmp means GRAIL-A, camera 1, day 089 of 2012 (29 March), 14:38:20 UT.

Step 2. Find where the frame is

  1. Look up where the spacecraft was at that time in JPL Horizons (https://ssd.jpl.nasa.gov/horizons/), with GRAIL-A (code -177) as the observer and the Moon as the target.
  2. Note the point below the spacecraft. At 14:38:20 it was about 26.9 E, 8.9 N, about 45 km up, flying south.
  3. Open the basemap around that point, about 5 degrees in each direction.
  4. Look for a crater with the same size and shape as the crater in the frame.
  5. Check at least 3 other features around it. Use their positions and sizes, not their shadows.
  6. Find the crater name in the Gazetteer of Planetary Nomenclature (https://planetarynames.wr.usgs.gov/).

The spacecraft position is only a hint. The camera does not always point straight down. In this frame, Carrel is about 50 km north of the point below the spacecraft. The match comes from the shapes in the image.

What matched here: Carrel and its terraced inside wall, the ridge just south of it, the straight rille running south-east, a long mound with a dark pit at its north-west end, and a cluster of small craters east of Carrel.

Step 3. Work out which way is north

  1. Find the shadow inside the big crater.
  2. Find where the Sun was. On 29 March 2012 the Moon was near first quarter, so at 27 E the Sun was in the east.
  3. Remember that the inside wall nearest the Sun is in shadow.
  4. Turn the frame until the shadow is on the east side.

In this frame the shadow is at the bottom, so east is down and north is to the right. Rotate the frame 90 degrees counterclockwise to see it with north up.

The basemap is lit from the west, so its shadows fall on the opposite side. Do not match shadow to shadow.

Step 4. Place control points

  1. Open the frame in the Georeferencer (QGIS), or add it to the map and click Georeference (ArcGIS Pro).
  2. Set the target coordinate system to IAU_2015:30100.
  3. Click the center of a feature in the frame.
  4. Click the center of the same feature on the basemap.
  5. Use crater centers, not crater edges. The center does not move when the light changes.
  6. Spread the points over the whole frame, including the corners.
  7. Place at least 6 points. We placed 9.

A control point joins one pixel in the frame to one longitude and latitude on the basemap. Our 9 points are in the table below.

You can skip the clicking and load our points. In QGIS, load the QGIS Georeferencer points file from the downloads. In ArcGIS Pro, import the control points table (CSV). Pixel x counts from the left edge, pixel y counts down from the top edge (QGIS writes y as a negative number).

Step 5. Fit and check

  1. Choose the transformation "Polynomial 1" (QGIS) or "1st Order Polynomial (Affine)" (ArcGIS Pro).
  2. Read the residual for each point in the table.
  3. Find any point with a residual much larger than the others.
  4. Check that point again on both images. Move it or delete it.
  5. Stop when no single point stands out.

A first-order fit can move, turn, stretch and slant the frame, but it cannot bend it. Residuals of about 1 km here come from the tilted view, not from wrong points.

Step 6. Save the result

  1. Set the resampling method to bilinear.
  2. Set the output to GeoTIFF.
  3. Run the georeferencing.
  4. Add the new GeoTIFF over the basemap.
  5. Set its transparency to 50 percent.
  6. Check that Carrel and the rille line up with the basemap.

Our result covers about 26.22 to 27.59 E and 9.38 to 11.02 N. Compare yours with the georeferenced GeoTIFF in the downloads.

Control points and residuals

PointFeaturePixelLongitude, latitudeResidual (frame pixels)
1Center of crater Carrel567.5, 170.526.7071, 10.66481.8
2Small sharp crater north-west of Carrel661.5, 81.526.4655, 10.89516.2
3Small crater on Carrel's south-east flank484.5, 280.526.99, 10.4779
4Crater B in the eastern crater cluster613.5, 388.527.3321, 10.722815
5Crater A in the eastern crater cluster658.5, 402.527.3735, 10.86256.9
6Dark pit at the north-west end of the long mound153.5, 345.527.1591, 9.745911.4
7Middle of the long bright mound190.5, 383.527.2687, 9.846117
8Isolated fresh crater south of Carrel113.5, 177.526.7461, 9.633713.6
9Small dark crater beside the rille93.5, 387.527.3303, 9.583311.8

RMS error: 11.2 frame pixels, about 847 m. One frame pixel is about 75 m on the ground.

Limits

Treat the result as approximate: good for "where is this", not for measuring.

Technical notes
  • RMS error 11.2 frame pixels, about 847 m, over 9 control points (largest 17.0 px, 1236 m).
  • Approximate. MoonKAM frames are not taken straight down. This one looks back at Carrel from the south at an angle, so the scale changes across the frame. A first-order fit cannot remove that, so features near the west and south-west edges can be off by 1 to 2 km.
  • Control points were picked by eye, and the Sun lit the frame from the east while the mosaic is lit from the west, so shadows fall on opposite sides. Crater centers were used instead of shadow edges.
  • Location confirmed visually: crater Carrel with its terraced interior, the ridge to its south, the long straight rille running south-east, the elongated mound with a dark pit at its end, and the crater cluster east of Carrel all match the LRO WAC mosaic in position and size.
  • The spacecraft position (JPL Horizons, GRAIL-A about 45 km above 26.9 E, 8.9 N at this time) was used only as a hint for where to look.
  • For comparison, a second-order polynomial on the same points gives an RMS of about 534 m, but with 9 points and 6 terms per axis it has little check on itself, so the product uses the first-order fit.
  • The GeoTIFF is a Cloud-Optimized GeoTIFF with nodata 0 outside the frame; the PNG has no alpha so JMARS can load it.

A residual is how far a control point lands from where it should after the warp. Small residuals mean the points agree with each other; they do not prove the points were matched to the right craters. That is why this layer is labeled approximate.

One MoonKAM frame over crater Carrel, georeferenced by hand

rasterApproximate
Coordinate system
Moon 2015 sphere, geographic lon/lat (equirectangular grid) IAU_2015:30100
Bounds
26.224° E to 27.592° E, 9.382° N to 11.018° N
Resolution
0.002 degrees per pixel (about 61 m north-south); the frame itself is about 75 m per pixel
Full coordinate system definition (WKT)
WKT
GEOGCS["Moon (2015) - Sphere / Ocentric",DATUM["Moon (2015) - Sphere",SPHEROID["Moon (2015) - Sphere",1737400,0,AUTHORITY["IAU","30100"]],AUTHORITY["IAU","30100"]],PRIMEM["Reference Meridian",0,AUTHORITY["IAU","30100"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["IAU","30100"]]

Georeferenced from the original 720 x 486 BMP with a first-order polynomial (affine) fit to 9 hand-picked control points, warped with rasterio 1.5.2 / GDAL 3.12.2 on 2026-10-07

MoonKAM imagery credit: NASA / Caltech-JPL / MIT / Sally Ride Science.