Project

General

Profile

« Previous | Next » 

Revision 9fd6a7de

Added by Jim Regetz over 13 years ago

  • ID 9fd6a7de83e8f400c9292a8e4c83cae43ef847e7

spell-checked boundary analysis doc (aspell)

View differences:

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