Revision 9fd6a7de
Added by Jim Regetz over 13 years ago
- ID 9fd6a7de83e8f400c9292a8e4c83cae43ef847e7
terrain/boundary/boundary.Rnw | ||
---|---|---|
154 | 154 |
additional approaches to fusing the layers, but do not include further |
155 | 155 |
assessment of these here. |
156 | 156 |
|
157 |
\subparagraph{Simple fusion} Naive concatentation of SRTM below 60\N with
|
|
157 |
\subparagraph{Simple fusion} Naive concatenation of SRTM below 60\N with |
|
158 | 158 |
ASTER above 60\N, without applying any modifications to deal with |
159 | 159 |
boundary artifacts (Figure \ref{blend-simple}). |
160 | 160 |
|
... | ... | |
199 | 199 |
|
200 | 200 |
\emph{Fused with exponential ramp north of 60\N.} |
201 | 201 |
The first step was to take the pixel-wise difference between SRTM and |
202 |
ASTER in the row immediately below 60\N (i.e., the northmost extent of |
|
202 |
ASTER in the row immediately below 60\N (i.e., the northernmost extent of
|
|
203 | 203 |
SRTM). An exponentially declining fraction of this difference was then |
204 | 204 |
then added back into the ASTER values north of 60\N. This does a fine |
205 | 205 |
job of eliminating the artificial shelf and thus the appearance of a |
... | ... | |
300 | 300 |
than SRTM for \Sexpr{round(100*mean(d.delta.vals>0))}\% of pixels (see |
301 | 301 |
Figure \ref{aster-srtm-scatter}). Thus, although adding a constant |
302 | 302 |
offset of \Sexpr{-delta.median} meters to the ASTER DEM would clearly |
303 |
center it with respect to the STRM (at least in the Canada focal
|
|
303 |
center it with respect to the SRTM (at least in the Canada focal
|
|
304 | 304 |
region), appreciable differences would remain. Figure |
305 | 305 |
\ref{aster-srtm-scatter} also highlights the existence of several |
306 | 306 |
obviously spurious ASTER spikes of >1000m; although not shown here, |
... | ... | |
310 | 310 |
Not surprisingly, simple fusion produces an artificial $\sim$12m cliff |
311 | 311 |
in the mean elevation profile (Figure \ref{mean-elevation}). At least in |
312 | 312 |
terms of mean elevation, this artifact is completely removed by both the |
313 |
multiresolution spline and gaussian weighted average methods. The
|
|
313 |
multiresolution spline and Gaussian weighted average methods. The
|
|
314 | 314 |
transition is, to the eye, slightly smoother in the former case, |
315 | 315 |
although ultimately this would depend on the chosen zone of overlap and |
316 | 316 |
on the exact parameterization of the weighting function. |
... | ... | |
320 | 320 |
degree (Figure \ref{mean-slope}). However, the shape of the profile |
321 | 321 |
itself is nearly identical between the two. Although this may partly |
322 | 322 |
reflect inherent SRTM vs ASTER differences, my guess is that CGIAR |
323 |
postprocessing of the particular SRTM product we're using has removed |
|
323 |
post-processing of the particular SRTM product we're using has removed
|
|
324 | 324 |
some of the high frequency ``noise'' that remains in ASTER? |
325 | 325 |
\textbf{\color{red}[todo: check!]} |
326 | 326 |
|
... | ... | |
386 | 386 |
expect in the presence of a cliff artifact at the seam. |
387 | 387 |
|
388 | 388 |
Interestingly, the aspect layers derived from the two blended DEMs |
389 |
(muliresolution spline and Gaussian weighted average) exhibit a |
|
389 |
(multiresolution spline and Gaussian weighted average) exhibit a
|
|
390 | 390 |
consistent mean northward inclination at all latitudes in their |
391 | 391 |
respective fusion zones. This pattern is visually obvious at latitudes |
392 | 392 |
between 59.95\N and 60\N in the bottom two panels of Figure |
... | ... | |
439 | 439 |
Gaussian weighted averaging produces a layer that exhibits a gradual |
440 | 440 |
transition from ASTER to SRTM, whereas the multiresolution spline yields |
441 | 441 |
an abrupt transition. Not surprisingly, the simple fused layer is even |
442 |
worse, producing not only a sudden transiton but also abberant values
|
|
442 |
worse, producing not only a sudden transition but also aberrant values
|
|
443 | 443 |
at the fusion seam itself; note downward (upward) spikes in correlation |
444 | 444 |
(RMSE) at 60\N in the first column of plots in Figures \ref{corr-slope} |
445 | 445 |
and \ref{rmse-slope}. |
... | ... | |
511 | 511 |
unexpected (and as-yet unexplained) negative spike at $\sim$59.95\N, the |
512 | 512 |
correlation profiles are fairly well-behaved for both blended layers, |
513 | 513 |
and don't show the same odd behavior as was the case for aspect. Again |
514 |
it is clear that the multiresultion spline results in a much more abrupt
|
|
514 |
it is clear that the multiresolution spline results in a much more abrupt
|
|
515 | 515 |
transition than does the Gaussian weighted average. The RMSE flow |
516 | 516 |
direction profiles echo this pattern (Figure \ref{rmse-flowdir}), and |
517 | 517 |
indeed look almost the same as those computed using aspect (Figure |
Also available in: Unified diff
spell-checked boundary analysis doc (aspell)