A Two-Step Block Adjustment Method for DSM Accuracy Improvement with Elevation Control of ICESat-2 Data

: Digital surface models (DSMs) have been widely utilized in various applications as fundamental geographic information data. Block adjustment is normally performed on satellite images to enhance the geometric accuracy and DSMs are then generated by stereo mapping. However, new errors may be introduced during the stereo mapping processing and geometric discrepancies between DSMs may still exist. In particular, block adjustment is difﬁcult for multisource satellite images. Therefore, this paper presents a two-step block adjustment approach directly performed on DSMs, with high-accuracy ICESat-2 laser altimetry data used as elevation control. In the method, DSM tie-point matching, elevation control/check point selection from ICESat-2 laser points, and planar and elevation block adjustments are performed in sequence. In the experiments, ZY-3 satellite stereo images and corresponding generated DSMs, as well as SRTM and ALOS DSMs, are used for veriﬁcation. The experimental results show that the absolute elevation accuracy and the relative geometric consistency of the DSMs are both signiﬁcantly improved after two-step DSM block adjustment and it can efﬁciently improve the accuracy, not only for DSMs acquired by the same sensor type, but also for DSMs acquired by different sensor types, which demonstrates the feasibility and advantage of the proposed method.


Introduction
Digital surface models (DSMs) have wide applications in the fields of surveying and civil engineering and are extremely important data for geographic information systems and remote sensing technology [1]. Stereo mapping using remote sensing images is a common way to generate DSMs [2]. However, the process of acquiring remotely sensed images is affected by the observation mode, atmospheric transmission, topographical fluctuations, various sensor factors, etc., all directly impacting the quality of the generated DSM [3,4]. The common approach used to improve an image's geometric accuracy for DSM generation is to perform block adjustment based on a rational function model (RFM) [5][6][7]. For instance, the block adjustment for the ZY-3 images is performed using laser altimetry data for the elevation constraint [8,9]. Auxiliary topographic data such as elevations extracted from the Shuttle Radar Topography Mission (SRTM) DEM can be taken as the initial value of the object elevation of the image tie points, achieving an overall improvement of the elevation positioning accuracy in large-scale image block adjustment [10,11]. Image-based block adjustment is generally suitable for homologous images. However, there are likely to be systematic deviations between DSMs generated by multisource images. Additionally, errors can be propagated throughout the workflows of image stereo mapping processing for DSM generation. Therefore, many studies have been conducted to improve the accuracy of DSMs. mapping processing for DSM generation. Therefore, many studies have been conducted to improve the accuracy of DSMs.
ASTER GDEM data are corrected by using ICESat (Ice, Cloud, and land Elevation Satellite) laser altimetry data [12], where the DEM artefacts and anomalies are removed using a segment-based outlier detection algorithm in order to correct the mosaiced DEM. For the fusion of multisource DEMs, the geometric consistency of DEMs is improved by using the elevation-error-based weights to determine the influence of each DEM in the fusion process [13,14], for which a priori accuracy information is needed. A rapid mosaicking method for InSAR-generated DSMs was proposed [15] that adopted the elevation correction model, which was determined based on error characteristics during interferometric processing [16]. For the TanDEM-X mission, block adjustment was performed using ICESat laser altimetry data as the constraint [17] and, based on analysis of the DSM elevation error generated in the interferometric processing workflow, the error correction model for block adjustment was set to refine the DSMs. However, if the DSMs have large initial location error, block adjustment using only elevation adjustments is insufficient due to the misalignment of the laser points and the DSM.
Compared with the full-waveform lidar system carried by ICESat, the Advanced Topographic Laser Altimeter System (ATLAS) carried by ICESat-2, which was launched in 2018, has higher accuracy and a smaller laser footprint [18,19]. High-accuracy elevation control points were extracted from ICESat-2 ATL08 data, which can be utilized as highprecision elevation reference data [20][21][22]. This paper addresses the fusion and accuracy improvement of multisensory DSMs and presents a two-step block adjustment method with ICESat-2 laser altimetry data as the elevation control to improve accuracy in both planar and elevation measurements.

Methodology
The flowchart of the proposed method is shown in Figure 1. For input DSMs, tiepoint matching is conducted first. It is worth noting that the original positioning error of the DSM will result in misalignment between DSM pixels and ICESat-2 data. Therefore, planar block adjustment must be conducted first, followed by elevation control data extraction and elevation block adjustment. The elevation control and check points for the block adjustment are selected from ICESat-2 data according to multiple criteria.

DSM Tie-Point Extraction
The acquired DSMs contain elevation and position errors. The extraction of tie-points needs to consider the influence of these two errors. The elevation contained in the DSM raster represents the terrain information, which can be utilized for tie-point extraction. A moving window is applied to extract the initial tie-points in the overlapping region of DSMs. Then, the elevation's standard deviation, as a relatively robust and efficient descriptor, is adopted to achieve extraction of tie-points [15]. Moreover, bilinear interpolation is used to obtain subpixel-level tie-points. Ultimately, the random sampling consistency (RANSAC) algorithm is performed to eliminate the mismatched tie-points.

Elevation Control/Check Points Selection from ICESat-2 Laser Altimetry Data
The ATL08 data are the Level-3A product of the ICESat-2 altimetry data and contain the land and canopy elevations with elevation accuracies better than 1 m after selection [18]. In the literature [21,22], the impact of various evaluation labels (atmospheric attenuation, surface reflectance, terrain slope, etc.) on the accuracy of ATL08 terrain elevation was analysed, providing reference regarding which label should be considered for elevation control points selection. In our method, the ICESat-2 ATL08 data are first filtered to match the region of the target DSM and then filtered further according to the comprehensive elevation attributes that reflect the elevation precision. The ICESat-2 ATL08 data filtering method in [20] is adopted for first-step selection. It is worth mentioning that the selection criterion in our method is more stringent so as to ensure the high quality of the selected laser points, as shown in Table 1. The points with values for the Cloud_flag_atm attribute less than 1 are preserved, and the height difference between the median terrain elevation and DSM is limited to less than 30 meters. Then, the surface slope, fitted by the photon elevations in the ATL08 segments, is set to less than 2 • so that the retained data should locate at flat terrain. For block adjustment, excessively dense points do not significantly improve the block adjustment results [9]. Therefore, a grid is set to evenly distribute the ATL08 data, and only one ATL08 laser point is ultimately retained in each grid cell. The grid size setting has a large tolerance, and in our method the grid size is set to 1/20 of the smaller of the width or height for input DSMs (For example, we set the grid size 2.85 km for the ZY-3 DSM with a size of 61.8 km × 57 km). Moreover, considering the uncertainty in DSM horizontal accuracy, for the DSM pixel containing the control point location, we count the elevation standard deviation between the target pixel and the surrounding pixels within a window size defined according to the initial horizontal error. The threshold of the standard elevation deviation is set as 0.7 (equal to tan (2 • ) × 20 m) for a DSM with 20 m resolution so that the DSM pixel is in a relatively flat region, as illustrated in the dashed box of Figure 1.

Two-Step DSM Block Adjustment
We propose a two-step DSM block adjustment approach to enhance the position accuracy of DSMs. First, planar block adjustment is performed by utilising the constraints of the consistency of the geographic coordinates of the tie-points. For DSMs, the geographic coordinates of each grid cell are determined by the geodetic coordinates of the corner point and the resolution in the lateral and longitudinal directions. The six-parameter affine transformation model is adopted to optimize the positioning parameters of each DSM, as shown in Equation (1).
where X and Y correspond to the geographic coordinates of a certain grid cell of the DSM; x and y are the corresponding row and column numbers; dx and dy are the latitude and longitude of the corner points of the DSM. That is, the planar block adjustment will optimize the location of the DSM, including the corner point coordinates and the A matrix (representing the rotation parameter and the scale parameter). Second, DSM elevation block adjustment is conducted after planar block adjustment to correct the relative elevation inaccuracy between adjacent DSMs. The parameters of the elevation error model for each DSM are optimized by the constraint of the elevation control points. Simultaneously, tie-points of the overlapping region are utilized to establish the consistency restrictions of DSM pixels. For satellite images, there are biases in the RFM coefficients owing to observation or sensor orientation errors [2]. The discrepancies between the measured and nominal line and sample coordinates are generally characterized as polynomial models of the image line and sample coordinates [5]. For input DSMs, the initial elevation error is first evaluated by the ATL08 check points. Then, the order of the polynomial model is determined by the elevation error distribution calculated by the extracted ATL08 data, as shown in Equation (2).
where g(x i , y i ) represents the elevation correction value and a i and b i are the model parameters corresponding to the DSM row and column numbers. For a certain tie-point pair x J , y J and (x K , y K ) on DSM J and K, respectively, the block adjustment constraint is that their final elevations are equivalent after corresponding elevation correction, as shown in Equation (3).
where g J x J , y J is the elevation correction value of DSM J, g K (x K , y K ) is the corresponding value of DSM K. h DEM_J x J , y J is the initial elevation of DSM J at location x J , y J , and h DEM_K (x K , y K ) is the elevation of DSM K at location (x K , y K ). Combined with the constraint on the control points, the observation equation can be written as Equation (4).
The first two lines (H GCP_J and H GCP_K ) of matrix A indicate that ICESat-2 control data are utilized to improve DSM J and K, respectively. The third line (H TP_J and H TP_K ) indicates the elevation constraints between the tie-points of adjacent DSMs. These matrices have the same form consisting of different orders of row and column numbers as follows: Then, the coefficients of the error correction model can be obtained by using the least square principle to eliminate the elevation system error of each DSM.

Accuracy Assessment
First, to intuitively compare the results between the block adjustment performed on the image space and the method directly conducted on DSMs proposed in this paper, corresponding DSMs are generated from the refined images by stereo mapping after the RFM-based image block adjustment [8,9,23,24]. The elevation check points are utilized to evaluate the elevation accuracy of DSMs generated by the two methods.
Then, the accuracy evaluation of the experimental results was carried out for two aspects for the ICESat-2 data-aided DSM block adjustment. On the one hand, plane block adjustment optimizes the coordinates of the corner points and the rotation parameters of DSMs. Therefore, the difference between the geographic latitude and longitude coordinates of the tie-points will be estimated before and after the block adjustment, corresponding to the relative planar error, and the height difference of the tie points correspond to the relative elevation error. On the other hand, the overall elevation accuracy of the DSM is enhanced through the constraint of the elevation consistency of the tie-points (improve the relative elevation accuracy) and the constraint of the ATL08 control points. Therefore, to assess the height accuracy evaluation of the DSMs before and after the block adjustment, the RMSE of the elevation difference between the check points (randomly chosen from the retained ICESat-2 data) and the corresponding DSM grid is calculated.

Study Area and Experimental Data
Based on the collected ZY-3 images acquired in April 2016 over the Taihu Lake Basin, which includes urban and suburban areas and where the terrain is mainly plain with some hilly areas, the corresponding DSMs with a resolution of 10 m are generated using the ZY-3 backward and forward images and the coverage is illustrated in Figure 2. To demonstrate the effectiveness of the proposed block adjustment method for multisensory DSM data, elevation data are 1extracted from SRTM data, with a 30 m resolution, and ALOS PALSAR DSM data, with a resolution of 12.5 m, both of which have absolute elevation accuracy within 10 m [25,26]. The coordinate system of all the experimental DSMs is based on the WGS84 datum. Moreover, the ICESat-2 ATL08 data that overlaps with the area of the DSMs acquired in a period from October 2018 are acquired and filtered to serve as elevation control and check points, which is the available earliest data since its launch, so that the altimetry data time is as close as possible to the DSM acquisition time to avoid, as much as possible, terrain change due to human activities. The height accuracy of the ATL08 data is 0.2 m in the plain area [20,27]. mote Sens. 2022, 14, x FOR PEER REVIEW 6 of 14 much as possible, terrain change due to human activities. The height accuracy of the ATL08 data is 0.2 m in the plain area [20,27].

Elevation Control and Check Point Selection and Elevation Correction Model Detection
The ATL08 data are first filtered to generate elevation control and check points. Considering the preliminary positioning error of the collected DSMs, we set a window with a size of 7 × 7 pixels to restrict the topographic relief around the target DSM grid cell that corresponds to the control point (this parameter has a large tolerance because we no longer need many elevation control points). The distribution of the ATL08 control and check points after the selection are illustrated in Figure 2. A total of 55,119 points are preliminarily acquired in the DSM region. After setting the threshold and a 280 × 280 grid for evenly distributed control points, 1314 ATL08 control points are reserved, one-third of which are randomly selected as the check points. The first-order and second-order polynomial models are adopted to rectify the elevation error for the near-linear and systematic elevation error distribution, which is in contrast with the ICESat-2 data. The accuracy of these two models has been compared in the block adjustment experiments for DSMs with control points, affirming in advance that the experimental results prove that the accuracy of the second-order polynomial is relatively higher. Therefore, if there is no explanation, we uniformly adopt the second-order polynomial model to suit the elevation error of DSMs.

Comparison of GCP-Free Block Adjustment Conducted on Image and Object Space
To verify the effectiveness of the proposed two-step block adjustment approach carried out on the object space, GCP-free DSM block adjustment is performed for equal comparison with the RFM-based image block adjustment method [7,24]. Through the

Elevation Control and Check Point Selection and Elevation Correction Model Detection
The ATL08 data are first filtered to generate elevation control and check points. Considering the preliminary positioning error of the collected DSMs, we set a window with a size of 7 × 7 pixels to restrict the topographic relief around the target DSM grid cell that corresponds to the control point (this parameter has a large tolerance because we no longer need many elevation control points). The distribution of the ATL08 control and check points after the selection are illustrated in Figure 2. A total of 55,119 points are preliminarily acquired in the DSM region. After setting the threshold and a 280 × 280 grid for evenly distributed control points, 1314 ATL08 control points are reserved, one-third of which are randomly selected as the check points. The first-order and second-order polynomial models are adopted to rectify the elevation error for the near-linear and systematic elevation error distribution, which is in contrast with the ICESat-2 data. The accuracy of these two models has been compared in the block adjustment experiments for DSMs with control points, affirming in advance that the experimental results prove that the accuracy of the second-order polynomial is relatively higher. Therefore, if there is no explanation, we uniformly adopt the second-order polynomial model to suit the elevation error of DSMs.

Comparison of GCP-Free Block Adjustment Conducted on Image and Object Space
To verify the effectiveness of the proposed two-step block adjustment approach carried out on the object space, GCP-free DSM block adjustment is performed for equal comparison with the RFM-based image block adjustment method [7,24]. Through the uncontrolled block adjustment of the ZY-3 original images, corresponding DSMs are generated by performing stereo mapping for subsequent accuracy comparison. Simultaneously, uncontrolled block adjustments are performed by the proposed approach on DSMs generated by the original ZY-3 images. The accuracy of the DSMs generated before and after the block adjustment Remote Sens. 2022, 14, 4455 7 of 13 by these two strategies are in contrast, as illustrated in Figure 3. The histograms of the DSM elevation errors from the check data are acquired before and after the GCP-free block adjustments with the two methods, as shown in Figure 4a-c. Evaluated through the checkpoints, the accuracies of these two strategies conducted in image or object space are equivalent, as shown in Table 2, with RMSEs of 7.00 m and 5.85 m. Therefore, it is difficult to confirm that uncontrolled block adjustment clearly improves absolute elevation accuracy. Both block adjustment results indicate the capability to correct the elevation disagreements in the overlapping DSMs, because their relative elevation error was reduced from 8.89 m to less than 2 m after block adjustment. For the original DSMs and corresponding results after block adjustment, the locations of tie-points are calculated, with the RMSE of the planar location error decreasing from 13.45 m initially to 3.26 m (DSMs after uncontrolled image block adjustment) and 3.11 m (after uncontrolled block adjustment directly in object space). For each DSM, since the ZY-3 images are obtained by stitching the data acquired by the block CCD, there will be elevation inconsistencies inside the DSMs, and this internal error is temporarily ignored. In addition, since it is difficult to precisely match the ICESat-2 data and the optical images, and the focus of this paper is to prove the effectiveness of the block adjustment in object space, the block adjustment experiment with control points conducted on images is not performed for comparison with the adjustment in object space.
uncontrolled block adjustment of the ZY-3 original images, corresponding DSMs are generated by performing stereo mapping for subsequent accuracy comparison. Simultaneously, uncontrolled block adjustments are performed by the proposed approach on DSMs generated by the original ZY-3 images. The accuracy of the DSMs generated before and after the block adjustment by these two strategies are in contrast, as illustrated in Figure  3. The histograms of the DSM elevation errors from the check data are acquired before and after the GCP-free block adjustments with the two methods, as shown in Figure 4ac. Evaluated through the checkpoints, the accuracies of these two strategies conducted in image or object space are equivalent, as shown in Table 2, with RMSEs of 7.00 m and 5.85 m. Therefore, it is difficult to confirm that uncontrolled block adjustment clearly improves absolute elevation accuracy. Both block adjustment results indicate the capability to correct the elevation disagreements in the overlapping DSMs, because their relative elevation error was reduced from 8.89 m to less than 2 m after block adjustment. For the original DSMs and corresponding results after block adjustment, the locations of tie-points are calculated, with the RMSE of the planar location error decreasing from 13.45 m initially to 3.26 m (DSMs after uncontrolled image block adjustment) and 3.11 m (after uncontrolled block adjustment directly in object space). For each DSM, since the ZY-3 images are obtained by stitching the data acquired by the block CCD, there will be elevation inconsistencies inside the DSMs, and this internal error is temporarily ignored. In addition, since it is difficult to precisely match the ICESat-2 data and the optical images, and the focus of this paper is to prove the effectiveness of the block adjustment in object space, the block adjustment experiment with control points conducted on images is not performed for comparison with the adjustment in object space.

Block Adjustment Results Conducted on DSMs with Elevation Control
We have demonstrated through uncontrolled block adjustment that the proposed two-step block adjustment conducted on object space can also enhance the elevation uniformity of DSMs. To further make use of the ICESat-2 laser altimetry data, block adjustments with elevation control points are subsequently carried out. We first compared the results of block adjustment with control points by adopting the first-order and secondorder polynomial elevation error models. As shown in Table 2, the block adjustment results obtained from the second-order polynomial are significantly better. Therefore, a second-order polynomial is adopted as the regression model to fit the systemically distributed elevation errors of DSMs. Because the usage of control points will bring extra influence on the parameters of elevation correction model, the deviation of elevation between tie-points after uncontrolled DSM block adjustment is smaller than that with control points. Figure 3a,d contrasts the ZY-3 DSMs before and after block adjustment and the introduction of the elevation control points. The distribution histogram of DSM elevation errors derived from check points is obtained before and after the block adjustments (with and without control points), as shown in Figure 4b,d. After block adjustment with elevation control points using the second-order elevation error model, the mean elevation error of the check points decreased from 3.18 m to 0.55 m, and the RMSE of the elevation error decreased from 7.62 m to 2.60 m. Figure 5a-d illustrates the error distribution of check points in the DSMs corresponding to the results in Figure 3. It can be seen that systematic errors remained after block adjustment from the image and object space without control

Block Adjustment Results Conducted on DSMs with Elevation Control
We have demonstrated through uncontrolled block adjustment that the proposed twostep block adjustment conducted on object space can also enhance the elevation uniformity of DSMs. To further make use of the ICESat-2 laser altimetry data, block adjustments with elevation control points are subsequently carried out. We first compared the results of block adjustment with control points by adopting the first-order and second-order polynomial elevation error models. As shown in Table 2, the block adjustment results obtained from the second-order polynomial are significantly better. Therefore, a second-order polynomial is adopted as the regression model to fit the systemically distributed elevation errors of DSMs. Because the usage of control points will bring extra influence on the parameters of elevation correction model, the deviation of elevation between tie-points after uncontrolled DSM block adjustment is smaller than that with control points. Figure 3a,d contrasts the ZY-3 DSMs before and after block adjustment and the introduction of the elevation control points. The distribution histogram of DSM elevation errors derived from check points is obtained before and after the block adjustments (with and without control points), as shown in Figure 4b,d. After block adjustment with elevation control points using the second-order elevation error model, the mean elevation error of the check points decreased from 3.18 m to 0.55 m, and the RMSE of the elevation error decreased from 7.62 m to 2.60 m. Figure 5a-d illustrates the error distribution of check points in the DSMs corresponding to the results in Figure 3. It can be seen that systematic errors remained after block adjustment from the image and object space without control points, whilst Figure 5d indicates that the overall elevation accuracy is considerably improved after the additional elevation constraint.

Block Adjustment Results Are Conducted on Multisource DSM with Elevation Control
To further confirm the suitability of our proposed method for multisensory DSMs, block adjustment experiments among SRTM, ALOS and ZY-3 DSMs are conducted. For the multisensory DSM block adjustment processing, we resampled SRTM DSMs and ALOS DSMs to make the resolution consistent with ZY-3 DSMs. Figures 6 and 7 illustrate the block adjustment results for the SRTM DSM from the aspect of the overall elevation consistency of the DSMs and the elevation error histogram of the check points, respectively. Figure 8 illustrates the error distribution of check points in the DSMs corresponding to the results in Figure 6. To avoid confusion with the check points of ZY-3, only the check points on the SRTM are exhibited. A systematic elevation deviation of more than 5 meters is shown, which is well corrected after block adjustment.
Remote Sens. 2022, 14, x FOR PEER REVIEW 9 of 14 points, whilst Figure 5d indicates that the overall elevation accuracy is considerably improved after the additional elevation constraint.4.4. Block adjustment results are conducted on multisource DSM with elevation control. To further confirm the suitability of our proposed method for multisensory DSMs, block adjustment experiments among SRTM, ALOS and ZY-3 DSMs are conducted. For the multisensory DSM block adjustment processing, we resampled SRTM DSMs and ALOS DSMs to make the resolution consistent with ZY-3 DSMs. Figures 6 and 7 illustrate the block adjustment results for the SRTM DSM from the aspect of the overall elevation consistency of the DSMs and the elevation error histogram of the check points, respectively. Figure 8 illustrates the error distribution of check points in the DSMs corresponding to the results in Figure 6. To avoid confusion with the check points of ZY-3, only the check points on the SRTM are exhibited. A systematic elevation deviation of more than 5 meters is shown, which is well corrected after block adjustment.   To further confirm the suitability of our proposed method for multisensory DSMs, block adjustment experiments among SRTM, ALOS and ZY-3 DSMs are conducted. For the multisensory DSM block adjustment processing, we resampled SRTM DSMs and ALOS DSMs to make the resolution consistent with ZY-3 DSMs. Figures 6 and 7 illustrate the block adjustment results for the SRTM DSM from the aspect of the overall elevation consistency of the DSMs and the elevation error histogram of the check points, respectively. Figure 8 illustrates the error distribution of check points in the DSMs corresponding to the results in Figure 6. To avoid confusion with the check points of ZY-3, only the check points on the SRTM are exhibited. A systematic elevation deviation of more than 5 meters is shown, which is well corrected after block adjustment.      Table 3, the ALOS DSMs with higher resolution have a higher initial elevation accuracy of 2.06 m, and the overall improvement of 41% is less noticeable than that of SRTM, 61%. Both the relative elevation errors calculated by the tie-points are reduced to less than 1.5 m after block adjustment. The accuracy of planar location also has a significant improvement after block adjustment, which is within a pixel. The error of check points for ALOS is presented in Figure 11, which also proves the improvement of elevation accuracy after block adjustment, and it can be seen that the initial elevation accuracy of ALOS is relatively higher than SRTM. However, due to the difference in sensors, the elevation of the DSM overlapping region is not as consistent as that of the homologous-source DSMs. Additionally, due to the lack of texture and various imaging  Table 3, the ALOS DSMs with higher resolution have a higher initial elevation accuracy of 2.06 m, and the overall improvement of 41% is less noticeable than that of SRTM, 61%. Both the relative elevation errors calculated by the tie-points are reduced to less than 1.5 m after block adjustment. The accuracy of planar location also has a significant improvement after block adjustment, which is within a pixel. The error of check points for ALOS is presented in Figure 11, which also proves the improvement of elevation accuracy after block adjustment, and it can be seen that the initial elevation accuracy of ALOS is relatively higher than SRTM. However, due to the difference in sensors, the elevation of the DSM overlapping region is not as consistent as that of the homologoussource DSMs. Additionally, due to the lack of texture and various imaging mechanisms, the elevation consistency in the lake area between the original ZY-3 and SAR DSMs is relatively obvious, as illustrated in the lower left corners of Figures 6 and 9.

Discussion
The method in this paper clearly improves the consistency and position accuracy of multisource DSMs by two-step block adjustment from the object space. Through the experimental results of block adjustment without the constraint of control points, it has been demonstrated that the proposed approach considerably improves the results as compared with the traditional correction method based on block adjustment conducted in image space. In the absence of control points, the two methods do not significantly improve the

Discussion
The method in this paper clearly improves the consistency and position accuracy of multisource DSMs by two-step block adjustment from the object space. Through the experimental results of block adjustment without the constraint of control points, it has been demonstrated that the proposed approach considerably improves the results as compared with the traditional correction method based on block adjustment conducted in image space. In the absence of control points, the two methods do not significantly improve the absolute accuracy of the DSM elevation, but the geometric consistency of the overlapping areas is significantly improved. However, block adjustment from the object space has the following advantages. On the one hand, block adjustment from the object space can directly register the multisource DSMs with the terrain information, avoiding complex matching processing between multisource images. On the other hand, for the ICESat-2 control data, it is easier to effectively match the 3D terrain contained in the DSM data, while the traditional block adjustment in image space requires the referenced footprint images of control data for registration. Moreover, it can be seen from the experiment that our method can improve the absolute accuracy of the multisource DSMs by more than 40% using the ICESat-2 control data as a constraint, and the consistency of the overlapping area is also improved significantly.

Conclusions
This paper proposes a two-step DSM block adjustment method using ICESat-2 ATL08 data as the elevation control. By calculating the elevation difference between the elevation check points and DSMs and the location deviation of the tie points, it has been demonstrated that the elevation accuracy and geometric consistency of the input DSMs have been greatly improved after block adjustment. The SRTM and ALOS DSMs generated from SAR have been utiliz2ed for experiments as multisource DSMs, and the experimental results have verified the feasibility and advantage of the proposed block adjustment method for improving the geometric consistency and quality of multisource DSMs.  Data Availability Statement: Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.