STC will provide a three-dimensional reconstruction of Mercury surface. The generation of a DTM of the observed features is based on the processing of the acquired images and on the knowledge of the intrinsic and extrinsic parameters of the optical system. The new stereo concept developed for STC needs a pre-flight verification of the actual capabilities to obtain elevation information from stereo couples: for this, a stereo validation setup to get an indoor reproduction of the flight observing condition of the instrument would give a much greater confidence to the developed instrument design. STC is the first stereo satellite camera with two optical channels converging in a unique sensor. Its optical model is based on a brand new concept to minimize mass and volume and to allow push-frame imaging. This model imposed to define a new calibration pipeline to test the reconstruction method in a controlled ambient. An ad-hoc indoor set-up has been realized for validating the instrument designed to operate in deep space, i.e. in-flight STC will have to deal with source/target essentially placed at infinity. This auxiliary indoor setup permits on one side to rescale the stereo reconstruction problem from the operative distance in-flight of 400 km to almost 1 meter in lab; on the other side it allows to replicate different viewing angles for the considered targets. Neglecting for sake of simplicity the Mercury curvature, the STC observing geometry of the same portion of the planet surface at periherm corresponds to a rotation of the spacecraft (SC) around the observed target by twice the 20° separation of each channel with respect to nadir. The indoor simulation of the SC trajectory can therefore be provided by two rotation stages to generate a dual system of the real one with same stereo parameters but different scale. The set of acquired images will be used to get a 3D reconstruction of the target: depth information retrieved from stereo reconstruction and the known features of the target will allow to get an evaluation of the stereo system performance both in terms of horizontal resolution and vertical accuracy. To verify the 3D reconstruction capabilities of STC by means of this stereo validation set-up, the lab target surface should provide a reference, i.e. should be known with an accuracy better than that required on the 3D reconstruction itself. For this reason, the rock samples accurately selected to be used as lab targets have been measured with a suitable accurate 3D laser scanner. The paper will show this method in detail analyzing all the choices adopted to lead back a so complex system to the indoor solution for calibration. |
I.INTRODUCTIONThe Stereo Camera (STC) [2] is the first stereoptic satellite push-frame camera designed for the study of Mercury. STC is onboard of the satellite BepiColombo. Scope of the mission is the study the surface composition and morphology, the geology and the magnetosphere of the planet. Many objectives will be provided by BepiColombo Mission: global mapping with stereo imaging and spectroscopy, color mapping with multi-band filters and hyperspectral imaging of selected regions and high spatial resolution imaging of at least 20% of the planet surface. II.THE STEREO CAMERA (STC)STC has a brand-new optical design (see Fig. 1), it’s a catadrioptic solution which allows to supply two optical wide angle channels with the orientation of +20° and -20° from the nadir direction maintaining limited mass and power. The solution to this minimization is the use of a common detector for both the channels and the adoption of mirrors to have the two different perspectives re-directed almost on the same optical path. [3] On the other side, five different filters in 410-930 nm spectral range permit multi-band acquisitions. Panchromatic ones are used for acquire stereo-pair guaranteeing resulting DTM (Digital Terrain Models) with a stereo vertical accuracy of 80 meters and a grid size of 50-110 m/pixel at periherm over the equator. A.Optical designSix months will be the maximum period in which global stereo mapping, first scientific goal of STC, will be completed. To reach this objective the optical design (proprieties are shown in Table 1) was conceived to reach a coverage of 40×19 km2 for the panchromatic filter and 40×3 km2 for the color ones (periherm configuration on equator). Table 1Optical and on-ground propriety of STC instrument.
Two independent elements compose the STC optical solution [4]: two folding mirrors per each channel, and a common telescope unit. The couple of folding mirrors redirects the ±20° wide beams along direction much closer to the system optical axis (±3.75°). The optical path is in common for both the channel and the rays are focused on a 10 μm pixel size SiPIN CMOS detector, useful in terms of radiation hardness and for the capability of snapshot image acquisition, allowing very short exposure times (< 1 ms). Result of this innovative design is a light instrument which cuts down on energy and permits stereo reconstruction with a narrow focal length of 95 mm to reach the 50 m/pixel scale factor at periherm. B.Stereo Validation Set-up (SVS)Considering the brand-new design adopted a pre-flight verification was needed. Base of the pipeline for the generation of the DTM of Mercury surface is the knowledge of the intrinsic and extrinsic parameters of the optical system. With this aim a stereo validation setup was realized for an indoor reproduction of the flight observing condition [5] [6]. See Fig. 2 for a scheme of the stereo validation set-up concept. Two principal problems were engaged: the need of acquire satellite images in a limited indoor space and the necessity of supply two different perspectives of a surface to use both the optical channels. First hindrance was solved by the use of an auxiliary optical system to observe with the STC instrument not a planetary surface from an altitude greater than 400 km but a selected small stone target at about 1 meter distance: a lens (collimator) that collimates the light rays coming from the target (essentially placing it at infinity) has been foreseen. Optical analysis by means of ray-tracing software led to the choice of an achromatic doublet with focal length of 1000 mm with a diameter of 80 mm. Neglecting for the sake of simplicity the Mercury curvature, multi-perspective acquisition was supplied by the use of two high precision rotational stages (see Fig. 3), which permit the acquisition of images from two different points of view. These stages allow to rotate the instrument and the target without the need of additional space not available in the close clean room environment. The first rotation stage, where the STC instrument (in a thermo-vacuum chamber) is positioned, has the function of aligning the optical axis of the active STC channel with the collimated light coming from the collimator lens; the second rotation stage has the function of rotating a target positioned on the collimator lens focusing plane, reproducing the angle between a pair of homologous rays. With the aim of reproduce different light conditions, a narrow-spot halogen lamp coupled to a collimator has been assembled on a curve rail system. Considering that the target illumination condition has to be maintained in the two acquisitions with STC channels, the rail was mounted on the target rotational stage in order to be moved together with the target. C.Scientific requirements scalingThe ratio between the observable feature size (horizontal and vertical accuracy) in the laboratory images and in the in-flight ones are proportional to the ratio between the collimator focal length and the distance of the SC from the Mercury surface during the acquisition. At periherm the observing distance is about 400 km: considering a collimator focal length of 1 m the scale ratio will be 1:400000. So, considering that the area of the Mercury surface that will be imaged by STC is approximately 40x20 km2, the corresponding target area that should be imaged with this setup in the laboratory is of approximately 9x4 cm2, referring to the panchromatic FoV. Rescaling parameters of Table 1 to obtain the spatial resolution corresponding to that of STC (50 m) with a 1 m collimator focal length, the target surface should be sampled with a step of 105 μm and the elevation should be determined with an accuracy of 190 μm. D.Rock target definitionScaling factor obtained permitted to choose a correct dimension for the target. To reproduce the same albedo condition of the ones attended from the Mercury surface basalt and anorthosite were chosen as composition material of the target stones. These rock samples were acquired by scanning by a CAM2® FaroArm Platinum, a portable Coordinate Measurement Machine (CMM) with a Z axis accuracy of 0.02 mm which supplied a DTM with point spacing of 0.05 mm, i.e. about half of the image resolution [5][6]. This permits to have a reference DTM with an accuracy better than the one required by scientific requirements. DTM obtained by laser scanner is than use to compare stereo reconstruction and to define the effective validity of the instrument approach. III.INTRINSIC AND EXTRINSIC CALIBRATIONAs explained before, if laser scanner data are available directly for stereo reconstruction, intrinsic and extrinsic coordinates are needed on one side to reduce time and error in reconstruction by using epipolar geometry, on the other side to obtain system geometry or triangulate points and to pass from image coordinates to metrical data. A first estimation of these parameters is obtained by an “ad hoc” calibration procedure which has been used for the two channels. Two important limits have to be taken in account in this process: on one side, the channels have a very narrow field of view (about 5.3°x2.3°), and in addition the limited depth of focus (reduced by the fact that any target has to be rotated by the rotational stage) makes difficult the use of strongly convergent images of a calibration target. A standard analytical calibration procedure would therefore lead to large estimation uncertainties as well as to strong correlations between intrinsic and extrinsic parameters. To partially cope with the problem, the principal point position of each channel (defined in this particular race by the collimator axis) has been optically measured using a laser target; moreover, since no significant radial distortion was apparent in the images of the calibration target (see Fig. 4), the corresponding parameters were neglected. Zhang’s algorithm [7] has been used to compute, as the only unknown, the focal length; the estimated value, although with a large uncertainty, turned out to be rather close to the expected nominal value of the focal length of the STC camera – collimator combination. As far as image orientation is concerned, two different approaches have been followed in cascade: Zhang camera calibration and Bundle adjustment refinements. A first set of orientation parameters has been obtained by using the stereo rig module of the well-known Camera Calibration Toolbox [7], introducing part of intrinsic parameters of each camera as known values. A.Zhang camera calibrationA sequence of about 50 image pairs of the chessboard target has been used to compute the relative orientation (extrinsic parameters) of the camera pair for STC FM. Since the relative geometry between the cameras is nominally unchanged, such parameters can be used in the derivation of the DTM, when the stone samples are substituted for the calibration target; however, a not trivial co-registration of the reconstructed DTM with the laser measured one is necessary to compare the two surfaces. Results of the Zhang calibration impose the rotation between the two optical axes of the channels angle are equal to 42.28° (with a relative error around 2%) and stereo base line equal to 1.23 m. These data correspond to a distance pupil-target equal to 1.59 m. In the case of intrinsic parameters while principal point and skew factors are determined, both the channels present a focal length around 164 mm, but the error band is greater and reach the 8%. It has to be considered that these parameters identify just a starting point for bundle adjustment refinement. B.Bundle adjustment refinementsIn the second part of the pipeline, the orientation is computed by means of a bundle block adjustment of a series of automatically identified tie points and some Ground Control Points (GCP) extracted from the reference DTM. A stereo pair of the rock target in the region of the tie-points is shown in Fig. 5. Tie points are automatically extracted and matched by structure from motion algorithms [8], [9]. GCP’s are object points of known coordinates, independently determined with an accuracy better than that obtainable by photogrammetric techniques; using them in a bundle adjustment allows to fix the reference system to the object and to control residual errors in intrinsic parameters. Tie points have two important scope: on one side they are well-recognizable both on images, DTM results and on points cloud results of the laser scanner acquisition; on the other side they permit a fast and precise comparison between stereo results and laser scanner acquisition without resorting to closest point algorithms. IV.RECONSTRUCTION RESULTSAs intrinsic and extrinsic parameters are recomputed by bundle adjustment, the DTM generation of the rock samples has been performed with the program Dense Matcher being developed at the University of Parma since 2006 [10]; this software implements the NCC (Normalized Cross Correlation) method, the Least Squares Matching (LSM) method [11] and the Multi-photo Geometrically Constrained Matching (MGCM) method [12]. LSM is a local area-based method, where a linear radiometric transformation and the affine mapping are introduced to minimize the squared sum of the grey values differences between the patch image and the template image. As input data the program requires interior and exterior orientation parameters of the images. The matching stage is embedded in a multi-resolution approach where three levels of an image pyramid are processed. A parallel dense matching procedure is implemented where the initial disparity map is computed by performing a feature based matching followed by interpolation on a grid. As described in Fig. 6 different tests were followed out varying the shape function in the LSM (Affine, Projective and Polynomial) and the template size (T17: 17×17 pix, T21: 21×21 pix, T25:25×25 pix). An example of the results of the different choices undertaken is shown in Fig. 7. The image shows that high spatial frequency are better reconstructed by the use of Polynomial shape function (more details are visible in Fig. 6b). Metrical analysis otherwise demonstrates how convergence to the laser scanner acquisitions is performed by all the models in the same way making useless the implementation of a higher degree functional mode. Several target stereo couples of images obtained with the STC functional breadboard and the SVS are currently under analysis. The point cloud produced by Dense Matcher has been imported, together with the reference DTM, in a 3D modeling software for the comparison. V.CONCLUSIONSAn innovative setup for the validation of stereo reconstruction capabilities of the SIMBIO-SYS STC has been developed and tested. The preliminary results obtained presented an overall reconstruction error that is compatible with the specifications. Aim of the work was in fact the achievement of a stereo accuracy less than 80 m on the Mercury Surface derived from the instrument scientific requirements. This upper bound was scaled in the controlled ambient of a clean room to 190 μm. A validation set-up to reproduce geometrical and stereo proprieties of the observatory strategy of BepiColombo was produced respecting light proprieties and material albedo characteristics. A step by step pipe-line was implemented with the aim of calibrate the system and obtain DTM reconstruction of the target stone reaching a RMS (compared to more fitted and precise laser reconstruction) less than 70 μm. The fact that RMS obtained by this method is less than half the required stereo accuracy demonstrated both the capabilities of the testing method and the validity of STC project. Work is still in progress in order both to better refine the stereo reconstruction techniques and also to use a novel approach to stereo reconstruction [13], [14]; moreover, additional sets of data will be acquired, to evaluate the impact of different sample morphologies and lighting conditions on stereo performance. The set up validated the reconstruction method for STC and supplied a bench of stereo images with different light condition that will permits to test different sun direction condition on Mercury and the appropriate parameter for stereo algorithm. These bench will be used even to improve the semi-global SIM algorithm [15] projected for this mission. ACKNOWLEDGMENTThis activity has been realized under the BepiColombo Agenzia Spaziale Italiana (ASI) contract to the Istituto Nazionale di Astrofisica (INAF) no I/022/10/0 in collaboration with Selex Electronic Systems S.p.A, Campi Bisenzio (FI). REFERENCESE. Flamini,
“SIMBIO-SYS: the spectrometer and imagers integrated observatory system for the BepiColombo planetary orbiter,”
Planet. Space Sci., 58 125
–143
(2010). Google Scholar
G. Cremonese,
“The stereo camera on the BepiColombo ESA/JAXA mission: a novel approach,”
Advances in Geosciences, 15 305
(2009). Google Scholar
V. Da Deppo, G.Naletto, G.Cremonese, L.Calamai,
“Optical design of the single-detector planetary stereo camera for the BepiColombo European Space Agency mission to Mercury,”
Appl. Opt., 49 2910
–2919
(2010). Google Scholar
V. Da Deppo,
“Optical Design performance of the stereo channel for SIMBIOSYS onboard the BepiColombo ESA mission,”
in Proceedings of the International Conference on Space Optics (ICSO2010),
(2010). Google Scholar
G. Naletto,
“Innovative optical setup for testing a stereo camera for space applications,”
in Proc. SPIE 8442, Space Telescope and Instrumentation 2012: Optical, Infrared and Millimeter Wave,
(2012). Google Scholar
E.Simioni, V.Da Deppo, G.Cremonese, C.Re, G.Naletto,
“Stereo Camera for satellite application: A new testing method,”
IEEE Workshop on Metrology for AeroSpace, in press, Google Scholar
Z. Zhang,
“A flexible new technique for camera calibration,”
IEEE Transactions on Pattern Analysis and Machine Intelligence, 22
(11), 1330
–1334
(2000). Google Scholar
Bouguet J.-Y.,
“Camera Calibration Toolbox for Matlab,”
(2010) http://www.vision.caltech.edu/bouguetj/calib Google Scholar
Bay H., T. Tuytelaars, and L. Van Gool,
“SURF: Speeded Up Robust Features,”
CVIU, 110
(3), 346
–359
(2008). Google Scholar
Nister D.,
“An efficient solution to the five-point relative pose problem,”
IEEE Transactions on Pattern Analysis and Machine Intelligence, 26
(6), 756
–7701
(2004). Google Scholar
Re C., R. Roncella, G. Forlani, G. Cremonese, G. Naletto,
“Evaluation of area-based image matching applied to DTM generation with Hirise images,”
ISPRS Annals of the Photogrammetry, Remote Sensing and Spatial Information Sciences, I-4
(2012). Google Scholar
Grün A.,
“Adaptive least squares correlations: a powerful matching techniques,”
South African J. Photogramm. Remote Sensing and Cartography, 14
(3), 175
–187
(1985). Google Scholar
A. Grün, E.P. Baltsavias,
“Geometrically constrained multiphoto matching,”
PE&RS, 54
(5), 633
–641
(1988). Google Scholar
M. Massironi,
“Simulations using terrestrial geological analogues to assess interpretability of potential geological features of the Hermean surface restituted by the STereo imaging Camera of the SIMBIOSYS package (BepiColombo mission),”
Planetary and Space Science, 56
(8), 1079
–1092
(2008). Google Scholar
E. Simioni,
“A new stereo algorithm based on snakes,”
Photogrammetric Engineering & Remote Sensing, 77
(5), 495
–507
(2011). Google Scholar
|