Galileo Controlled DTMs Created by the USGS
These data can be accessed in three ways:Knowledge of a planetary surface’s topography is necessary to understand its geology and enable landed mission operations. The Solid State Imager (SSI) on board NASA’s Galileo spacecraft ( Citation: Belton, Head & al., 1996 Belton, M., Head, J., Ingersoll, A., Greeley, R., McEwen, A., Klaasen, K., Senske, D., Pappalardo, R., Collins, G., Vasavada, A., Sullivan, R., Simonelli, D., Geissler, P., Carr, M., Davies, M., Veverka, J., Gierasch, P., Banfield, D., Bell, M., Chapman, C., Anger, C., Greenberg, R., Neukum, G., Pilcher, C., Beebe, R., Burns, J., Fanale, F., Ip, W., Johnson, T., Morrison, D., Moore, J., Orton, G., Thomas, P. & , R. (1996). Galileo’s First Images of Jupiter and the Galilean Satellites. Science, 274(5286). 377–385. https://doi.org/10.1126/science.274.5286.377 ) acquired more than 700 images of Jupiter’s moon Europa. Although moderate- and high-resolution coverage is extremely limited, repeat coverage of a small number of sites enables the creation of digital terrain models (DTMs) via stereophotogrammetry. Here we provide stereo-derived DTMs of five sites on Europa. The sites are the bright band Agenor Linea, the crater Cilix, the crater Pwyll, pits and chaos adjacent to Rhadamanthys Linea, and ridged plains near Yelland Linea. We generated the DTMs using BAE’s SOCET SET® software and each was manually edited to correct identifiable errors from the automated stereo matching process ( Citation: Bland, Kirk & al., 2021 Bland, M., Kirk, R., Galuszka, D., Mayer, D., Beyer, R. & Fergason, R. (2021). How Well Do We Know Europa’s Topography? An Evaluation of the Variability in Digital Terrain Models of Europa. Remote Sensing, 13(24). ) . Additionally, we used the recently updated image pointing information provided by the U.S. Geological Survey ( Citation: Bland, Weller & al., 2021 Bland, M., Weller, L., Archinal, B., Smith, E. & Wheeler, B. (2021). Improving the Usability of Galileo and Voyager Images of Jupiter’s Moon Europa. Earth and Space Science, 8(12). e2021EA001935. https://doi.org/https://doi.org/10.1029/2021EA001935 ) , which ties the DTMs to an existing horizontal datum and enables the DTMs to easily be used in coordination with that globally controlled image set. The DTMs are of the highest quality achievable with Galileo data and are therefore suitable for most scientific analysis. However, there are inherent uncertainties in the DTMs including horizontal resolutions that are typically 1–2 km (~10x the image pixel scale) and expected vertical precision (the root mean square (RMS) uncertainty in a point elevation) of 10s–100s of m. The DTMs and their uncertainties are discussed in detail in ( Citation: Bland, Kirk & al., 2021 Bland, M., Kirk, R., Galuszka, D., Mayer, D., Beyer, R. & Fergason, R. (2021). How Well Do We Know Europa’s Topography? An Evaluation of the Variability in Digital Terrain Models of Europa. Remote Sensing, 13(24). ) .
Raw images were downloaded from NASA’s Planetary Data System (PDS) in .img format and converted to ISIS (Integrated Software for Imagers and Spectrometers) image cubes using the ISIS command gllssi2isis. Position and pointing information (i.e., Navigation and Ancillary Information Facility Spacecraft, Planet, Instrument, Camera-matrix, Events (NAIF SPICE)) was applied, including the updated pointing kernels (CK) and associated planetary constants kernels (PCK) provided by the USGS for Europa Galileo images ( Citation: Bland, Weller & al., 2021 Bland, M., Weller, L., Archinal, B., Smith, E. & Wheeler, B. (2021). Improving the Usability of Galileo and Voyager Images of Jupiter’s Moon Europa. Earth and Space Science, 8(12). e2021EA001935. https://doi.org/https://doi.org/10.1029/2021EA001935 ) . Radiometric calibration was applied (ISIS’ gllssical), and then the images were transferred to BAE’s SOCET SET® software where an additional bundle adjustment was performed using SOCET SET®’s Multi-Sensor Triangulation (MST). In MST, tie-points, and in some cases “z” (or height) constrained ground control points are measured, followed by iterative bundle adjustment to achieve a solution with sub-pixel accuracy. The DTMs were then extracted using the known camera orientation and automated high-density area and feature-based matching using SOCET SET®’s Next Generation Automatic Terrain Extraction (NGATE) ( Citation: Zhang, 2006 Zhang, B. (2006). Toward a higher level of automation in softcopy photogrammetry: NGATE and LIDAR processing in SOCET SET®.. ) algorithm with a low smooth option. In cases that used multiple images, the DTM from each stereo model was extracted separately and then merged. We assessed each DTM to ensure that proper alignment with the globally controlled image set of ( Citation: Bland, Weller & al., 2021 Bland, M., Weller, L., Archinal, B., Smith, E. & Wheeler, B. (2021). Improving the Usability of Galileo and Voyager Images of Jupiter’s Moon Europa. Earth and Space Science, 8(12). e2021EA001935. https://doi.org/https://doi.org/10.1029/2021EA001935 ) was maintained. Each DTM was then manually edited to remove artifacts and other errors. In the manual editing process, a technician uses a stereo viewer to manually identify errors due to incorrect automated stereo-matching and adjust the height of those posts closer to the actual three-dimensional surface. This significantly improves the DTM accuracy where errors due to incorrect matching have occurred.
The only significant procedural difference between the five DTMs was in the details of how MST was run. For the pits and chaos near Rhadamanthys Linea and for Cilix crater (single stereo pairs), the nadir-most image was “held” (prevented from adjusting its position in space), resulting in a so-called “dependent-relative solution” using only tie-points. Pwyll crater (three images, two stereo pairs) was similarly controlled using only tie points while holding the single overlapping image and allowing the remaining two images to adjust. The DTM of ridged terrain near Yelland Linea (three images, two stereo pairs) required establishing three XYZ-control points in a wide triangle pattern across the stereo model mosaic in addition to a set of tie-points. The coordinates for the control points (longitude, latitude, height) were obtained from the a priori position and pointing of the updated images directly upon import into SOCET SET®. For Agenor Linea (nine images, seven stereo pairs), we applied four XYZ-control points using a priori image coordinates from the updated images. We established these four control points near the edges of the stereo model mosaic in addition to a set of tie-points. For both the ridged terrain near Yelland Linea and Agenor Linea all images were allowed to adjust, resulting in a so-called “independent relative solution.”
Five DTMs are included in this release:
- Cilix crater
- Chaos and pits near Rhadamanthys Linea
- Pwyll crater
- Ridged plains near Yelland Linea
- A portion of Agenor Linear
The image and stereo criteria for each DTM is shown in Table 1. Each DTM is accompanied by a figure of merit (FOM) and colorized confidence map (ClrCONF), which are described below under Accuracy, Issues, and Errors.
- DTM: The USGS created DTM
- Hillshade: A shaded-relief image derived from the DTM
- FOM: Figure of merit image showing pixel provenance for each pixel in the DTM.
- Confidence: Confidence images that group FOM classes into qualitative classes.
The following description is based on the analysis described in ( Citation: Bland, Kirk & al., 2021 Bland, M., Kirk, R., Galuszka, D., Mayer, D., Beyer, R. & Fergason, R. (2021). How Well Do We Know Europa’s Topography? An Evaluation of the Variability in Digital Terrain Models of Europa. Remote Sensing, 13(24). ) ; however, the five DTMs provided here are not identical to those described in that work. Subsequent to publication of that work, updated image locations were made available for Galileo images of Europa that substantially improve the horizontal registration of that dataset ( Citation: Bland, Weller & al., 2021 Bland, M., Weller, L., Archinal, B., Smith, E. & Wheeler, B. (2021). Improving the Usability of Galileo and Voyager Images of Jupiter’s Moon Europa. Earth and Space Science, 8(12). e2021EA001935. https://doi.org/https://doi.org/10.1029/2021EA001935 ) . We therefore re-generated the DTMs provided here using the updated image locations so that the DTMs are co-aligned with the image data. These DTMs incorporate small improvements in our processing approach and are thus effectively version 2 (whereas those released alongside ( Citation: Bland, Kirk & al., 2021 Bland, M., Kirk, R., Galuszka, D., Mayer, D., Beyer, R. & Fergason, R. (2021). How Well Do We Know Europa’s Topography? An Evaluation of the Variability in Digital Terrain Models of Europa. Remote Sensing, 13(24). ) are version 1). The tweaks to our technical approach and the manual editing processes result in modest differences in pixel values between version 1 and 2. However, the general findings of ( Citation: Bland, Kirk & al., 2021 Bland, M., Kirk, R., Galuszka, D., Mayer, D., Beyer, R. & Fergason, R. (2021). How Well Do We Know Europa’s Topography? An Evaluation of the Variability in Digital Terrain Models of Europa. Remote Sensing, 13(24). ) regarding DTM accuracy and errors, which we summarized next, still hold. The version 2 DTMs provided here should be comparable or superior to version 1. The descriptions of DTM quality below are therefore conservative.
DTM accuracy is summarized in Table 1. For each DTM, we provide the minimum effective resolution based on the smallest observable features in the DTM and hillshade. However, we emphasize that this measure of resolution does not indicate that all surface features at that scale can be observed. For example, in several cases (especially Pwyll crater) the “smallest” observable features are double ridges that are just two DTM pixels wide, but extend across large sections of the DTM. We also show two different values for the expected vertical precision (EP). One, which we simply designate EP, was calculated from Root Mean Square (RMS) ground sample distance (GSD), the image pair’s base-height ratio, and assuming a pixel matching accuracy (ρ) of 0.3 (i.e., this is the “formal” EP). The other, designated EPderived, is based on a detailed evaluation of the DTMs themselves and does not assume ρ is 0.3 pixels. ( Citation: Bland, Kirk & al., 2021 Bland, M., Kirk, R., Galuszka, D., Mayer, D., Beyer, R. & Fergason, R. (2021). How Well Do We Know Europa’s Topography? An Evaluation of the Variability in Digital Terrain Models of Europa. Remote Sensing, 13(24). ) showed that ρ~1 for Galileo images of Europa, so we consider EPderived a better indication of the vertical precision. We emphasize that in either case EP (or EPderived) is not an “uncertainty.” Rather, it is the RMS uncertainty in a point elevation sampled from the DTM. That is, it is a statistical measure of the quasi-normally distributed error. Thus errors several times larger than EP are expected for some small fraction of points. In no case should EP be thought of as an upper bound on the error within a given DTM. Although errors are often close to normally distributed such that EPRMS is similar to the standard deviation, this is not always the case. For example, the error distribution for the Pwyll crater DTM has extremely long tails and the standard deviation is much larger than EPRMS. While no obvious DTM blunders or gaps occur (they are removed by the manual editing process), it is likely that textures near the resolution limit are artifacts rather than a reflection of Europa’s topography.
Imagea | GSD (m/pixel) |
GSDRMSb (m/pixel) |
ic | ed | DSTe | (p/h)f | |
Cilix crater | s0449965000 | 110 | 90 | 31.8° | 25.9° | 0.044 | 1.18 |
s0449967535 | 63 | 36.8° | 29.9° | ||||
Chaos and pits | s0449961826 | 232.1 | 170 | 79.8° | 38.8° | 0.14 | 1.32 |
s0484888726 | 63.3 | 79.9° | 47.6° | ||||
Pwyll crater | s0383715500 | 243.8 | 194 | 76.2° | 40.3° | 4.1 | 1.8 |
s0383715504 | 244.7 | 80.0° | 43.9° | ||||
s0426268700 | 125.8 | 57.9° | 48.3° | ||||
Ridged Plains | s0426272642 | 14 | 21 | 29.8° | 37.8° | 0.004 | 0.55 |
s0426272646 | 14 | 29.8° | 37.8° | ||||
s0426272821 | 25.6 | 29.9° | 47.7° | ||||
Agenor Linea | s0466669500 | 50 | 156 | 75-80° | 37-39° | 2.2 | 0.88 |
s0466669513 | |||||||
s0466669526 | |||||||
s0466669539 | |||||||
s0466669552 | |||||||
s0466669565 | |||||||
s0466669578 | |||||||
s0466664665 | 210-223 | 71-80° | 31-40° | ||||
s0466664326 |
bGSDRMS is the root mean square (RMS) ground sample distance (GSD) or pixel scale.
ci is the incidence angle
de is the emission angle
eDST is shadow tip difference. Low values are preferable and indicate similar illumination. See ( Citation: Becker, Archinal & al., 2015 Becker, K., Archinal, B., Hare, T., Kirk, R., Howington-Kraus, E., Robinson, M. & Rosiek, M. (2015). Criteria for Automated Identification of Stereo Image Pairs. ) .
f(p/h) is the base-height ratio. Ideal values are from 0.4 to 0.6, which correspond to convergence angles of 20– 30o.
For each DTM, we also provide a “Figure of Merit” and a colorized confidence map. The FOM is a map of numerical values generated by SOCET SET® that indicates one of three things for a given post measurement:
- It may be an error flag value (e.g., for poor automatic measurement).
- It may indicate a successful measurement.
- It may be an edit flag value, indicating the type of editing used. FOM values ranging from 0 to 39 indicate error or edit flags. FOM codes greater than 40 indicate a successful automatic correlation, with larger values indicating a better measurement. Table 3 lists the meaning of each FOM value.
Because the FOM codes can be difficult to parse, we also provide a confidence map for each DTM. The confidence map is derived from the FOM and simply groups different types of FOM values together. For example, all nine types of edit flags are grouped together as simply “manually interpolated.” The map is especially useful for understanding which posts were edited, as well as the quality of the correlation. The confidence map legend is provided below.
The DTMs were generated from image data with improved locations relative to the global imaging data set ( Citation: Bland, Weller & al., 2021 Bland, M., Weller, L., Archinal, B., Smith, E. & Wheeler, B. (2021). Improving the Usability of Galileo and Voyager Images of Jupiter’s Moon Europa. Earth and Space Science, 8(12). e2021EA001935. https://doi.org/https://doi.org/10.1029/2021EA001935 ) . The relative and absolute horizontal accuracy of the DTMs is therefore identical to that of the image data set. That data set have an image-to-image (i.e., relatively alignment) RMS uncertainty of 246.6 m in latitude, 307.0 m in longitude, and 70.5 m in radius. The image data are also tied to Europa’s geodetic reference system. The center of the crater Cilix, which defines the reference system at 182º W (178º E), is located at 181.9991414ºW (178.0008586ºE). The difference is approximately equivalent to 23 m, which is less than the pixel scale of the highest resolution Galileo images of Cilix. The location is therefore “known” to a fraction of a pixel and actually exceeds the precision with which the center of Cilix (a 19-km-diameter crater) is known. However, given that the data set is tied to the reference frame only at a single point, the certainty with which absolute coordinates are known degrades with distance from Cilix (generally as the √N where N is the number of images away from Cilix).
The manual editing process ensures that these DTMs are of extremely high quality with the vast majority of blunders corrected. They are generated from, and therefore well-aligned with, the improved Galileo and Voyager image data set of ( Citation: Bland, Weller & al., 2021 Bland, M., Weller, L., Archinal, B., Smith, E. & Wheeler, B. (2021). Improving the Usability of Galileo and Voyager Images of Jupiter’s Moon Europa. Earth and Space Science, 8(12). e2021EA001935. https://doi.org/https://doi.org/10.1029/2021EA001935 ) . We therefore consider them to be the highest quality possible from stereophotogrammetry of Galileo images, and they are therefore fit for most scientific uses, including geologic and morphometric analysis. However, the user should be aware of their inherent limitations. ( Citation: Bland, Kirk & al., 2021 Bland, M., Kirk, R., Galuszka, D., Mayer, D., Beyer, R. & Fergason, R. (2021). How Well Do We Know Europa’s Topography? An Evaluation of the Variability in Digital Terrain Models of Europa. Remote Sensing, 13(24). ) found that the horizontal resolution is typically ~10x the DTM pixel scale (m/post), which is consistent with recent studies of DTMs generated for Mars ( ( Citation: Kirk, Mayer & al., 2021 Kirk, R., Mayer, D., Fergason, R., Redding, B., Galuszka, D., Hare, T. & Gwinner, K. (2021). Evaluating Stereo Digital Terrain Model Quality at Mars Rover Landing Sites with HRSC, CTX, and HiRISE Images. Remote Sensing, 13(17). https://doi.org/10.3390/rs13173511 ) found 10–20 km). Features below the resolution given in Table 2 (typically 1–2 km) are almost certainly present on Europa. Furthermore, evaluation of slopes are also limited by the DTM pixel scale (100s of m), and steeper slopes may be present at small scaler scales. The expected vertical precision given in Table 2 is a statistical measure of the DTMs error distribution, and the uncertainty at a single point may be several times larger. These issues are discussed at length in ( Citation: Bland, Kirk & al., 2021 Bland, M., Kirk, R., Galuszka, D., Mayer, D., Beyer, R. & Fergason, R. (2021). How Well Do We Know Europa’s Topography? An Evaluation of the Variability in Digital Terrain Models of Europa. Remote Sensing, 13(24). ) , and we encourage the user of these DTMs to read that manuscript to obtain a better understanding of uncertainties inherent in these data products.
DTM | Post Spacing | Resolution | EP | EPderived |
---|---|---|---|---|
Cilix crater | 250 m | 1.1 km | 23 m | 36 m |
Pwyll crater | 750 m | 1.5 km | 32 m | 112 m |
Chaos and Pits | 200 m | 1.8 km | 39 m | 66 m |
Ridged Plains | 60 m | 180 m | 11 m | 14–21 m |
Agenor Linea | 450 m | 3.0 km | 52 m | ~100 m |
The DTMs were generated from Galileo images with improved geometry ( Citation: Bland, Weller & al., 2021 Bland, M., Weller, L., Archinal, B., Smith, E. & Wheeler, B. (2021). Improving the Usability of Galileo and Voyager Images of Jupiter’s Moon Europa. Earth and Space Science, 8(12). e2021EA001935. https://doi.org/https://doi.org/10.1029/2021EA001935 ) . Those images and accompanying documentation are available via a SpatioTemporal Asset Catalog at https://psdi.astrogeology.usgs.gov/europa/data/data_products/ .
These data are released under the CC0-1.0 license , meaning you can copy, modify, and distribution these data without permissions. We ask that you cite these data if you make use of them. The citation to be used is:
Bland, M.T.; Kirk, R.L.; Galuszka, D.M.; Mayer, D.P.; Beyer, R.A.; Fergason, R.L. How Well Do We Know Europa’s Topography? An Evaluation of the Variability in Digital Terrain Models of Europa. Remote Sens. 2021, 13, 5097. https://doi.org/10.3390/rs13245097
- Becker, K., Archinal, B., Hare, T., Kirk, R., Howington-Kraus, E., Robinson, M. & Rosiek, M. (2015). Criteria for Automated Identification of Stereo Image Pairs.
- Belton, M., Head, J., Ingersoll, A., Greeley, R., McEwen, A., Klaasen, K., Senske, D., Pappalardo, R., Collins, G., Vasavada, A., Sullivan, R., Simonelli, D., Geissler, P., Carr, M., Davies, M., Veverka, J., Gierasch, P., Banfield, D., Bell, M., Chapman, C., Anger, C., Greenberg, R., Neukum, G., Pilcher, C., Beebe, R., Burns, J., Fanale, F., Ip, W., Johnson, T., Morrison, D., Moore, J., Orton, G., Thomas, P. & , R. (1996). Galileo’s First Images of Jupiter and the Galilean Satellites. Science, 274(5286). 377–385. https://doi.org/10.1126/science.274.5286.377
- Bland, M., Weller, L., Archinal, B., Smith, E. & Wheeler, B. (2021). Improving the Usability of Galileo and Voyager Images of Jupiter’s Moon Europa. Earth and Space Science, 8(12). e2021EA001935. https://doi.org/https://doi.org/10.1029/2021EA001935
- Bland, M., Kirk, R., Galuszka, D., Mayer, D., Beyer, R. & Fergason, R. (2021). How Well Do We Know Europa’s Topography? An Evaluation of the Variability in Digital Terrain Models of Europa. Remote Sensing, 13(24).
- Kirk, R., Mayer, D., Fergason, R., Redding, B., Galuszka, D., Hare, T. & Gwinner, K. (2021). Evaluating Stereo Digital Terrain Model Quality at Mars Rover Landing Sites with HRSC, CTX, and HiRISE Images. Remote Sensing, 13(17). https://doi.org/10.3390/rs13173511
- Zhang, B. (2006). Toward a higher level of automation in softcopy photogrammetry: NGATE and LIDAR processing in SOCET SET®..