<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";}
@page WordSection1
{size:8.5in 11.0in;
margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-US" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal">How do I construct a map file that can handle different resolutions of terrain, with fallback to lower resolution when a no-value is found? In Google Earth Fusion, I can specify multiple terrain sources and tell it which value(s) represent
missing data, and it will then process the terrain, checking for availability of each type of terrain, falling back if there isn’t any in a geographic area and falling back if there is missing data.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">For simple fallback based on area coverage, I can use a map file like the one below (using 30 Arcsecond SRTM and 1 Arcsecond DTED). Now, what would I need to do to have it use/interpolate data from the 30 Arcsecond values when there is
a missing data/NODATA/hole in the DTED? In the file below, gdaltindex has been used to generate the coverage shapefiles.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">The example below is a simple one, and it would be really helpful to be able to have multiple terrain sources and a way to automatically fill voids the same way that transparency works with raster image layers. Specifying that terrain
values are posts and not pixels would be good too, but lower priority.<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Xander<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">MAP<o:p></o:p></p>
<p class="MsoNormal"> NAME "srtm"<o:p></o:p></p>
<p class="MsoNormal"> EXTENT -180 -90 180 90<o:p></o:p></p>
<p class="MsoNormal"> SIZE 512 512<o:p></o:p></p>
<p class="MsoNormal"> UNITS DD<o:p></o:p></p>
<p class="MsoNormal"> IMAGETYPE bil<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> PROJECTION<o:p></o:p></p>
<p class="MsoNormal"> "init=epsg:4326"<o:p></o:p></p>
<p class="MsoNormal"> END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> OUTPUTFORMAT<o:p></o:p></p>
<p class="MsoNormal"> NAME "bil"<o:p></o:p></p>
<p class="MsoNormal"> DRIVER "GDAL/EHdr"<o:p></o:p></p>
<p class="MsoNormal"> MIMETYPE "image/bil"<o:p></o:p></p>
<p class="MsoNormal"> EXTENSION "bil"<o:p></o:p></p>
<p class="MsoNormal"> IMAGEMODE INT16<o:p></o:p></p>
<p class="MsoNormal"> END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> OUTPUTFORMAT<o:p></o:p></p>
<p class="MsoNormal"> NAME "bil16"<o:p></o:p></p>
<p class="MsoNormal"> DRIVER "GDAL/EHdr"<o:p></o:p></p>
<p class="MsoNormal"> MIMETYPE "application/bil16"<o:p></o:p></p>
<p class="MsoNormal"> EXTENSION "bil"<o:p></o:p></p>
<p class="MsoNormal"> IMAGEMODE INT16<o:p></o:p></p>
<p class="MsoNormal"> END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> OUTPUTFORMAT<o:p></o:p></p>
<p class="MsoNormal"> NAME "gtiff_16"<o:p></o:p></p>
<p class="MsoNormal"> DRIVER "GDAL/GTiff"<o:p></o:p></p>
<p class="MsoNormal"> MIMETYPE "image/tiff"<o:p></o:p></p>
<p class="MsoNormal"> EXTENSION "tif"<o:p></o:p></p>
<p class="MsoNormal"> IMAGEMODE INT16<o:p></o:p></p>
<p class="MsoNormal"> END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> WEB<o:p></o:p></p>
<p class="MsoNormal"> IMAGEPATH "/tmp/"<o:p></o:p></p>
<p class="MsoNormal"> IMAGEURL "/tmp/"<o:p></o:p></p>
<p class="MsoNormal"> LOG "/tmp/srtm.err"<o:p></o:p></p>
<p class="MsoNormal"> METADATA<o:p></o:p></p>
<p class="MsoNormal"> wms_title "SRTM Terrain"<o:p></o:p></p>
<p class="MsoNormal"> END<o:p></o:p></p>
<p class="MsoNormal"> END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> #<o:p></o:p></p>
<p class="MsoNormal"> # Coarse (30 ArcSecond) terrain. There is worldwide coverage, so this<o:p></o:p></p>
<p class="MsoNormal"> # is good for areas uncovered by SRTM data, as well as for high<o:p></o:p></p>
<p class="MsoNormal"> # scale values.<o:p></o:p></p>
<p class="MsoNormal"> #<o:p></o:p></p>
<p class="MsoNormal"> LAYER<o:p></o:p></p>
<p class="MsoNormal"> NAME "srtm"<o:p></o:p></p>
<p class="MsoNormal"> TILEINDEX "SRTM30/bathmetry.shp"<o:p></o:p></p>
<p class="MsoNormal"> TYPE RASTER<o:p></o:p></p>
<p class="MsoNormal"> PROCESSING "RESAMPLE=BICUBIC"<o:p></o:p></p>
<p class="MsoNormal"> STATUS ON<o:p></o:p></p>
<p class="MsoNormal"> <o:p></o:p></p>
<p class="MsoNormal"> METADATA<o:p></o:p></p>
<p class="MsoNormal"> "wms_title" "SRTM 30 ArcSecond Bathmetry"<o:p></o:p></p>
<p class="MsoNormal"> "wms_group_title" "srtm30"<o:p></o:p></p>
<p class="MsoNormal"> "wms_srs" "epsg:4326"<o:p></o:p></p>
<p class="MsoNormal"> "wms_style" "Default"<o:p></o:p></p>
<p class="MsoNormal"> "wms_extent" "-180 -90 180 90"<o:p></o:p></p>
<p class="MsoNormal"> END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> MAXSCALEDENOM 1000000000<o:p></o:p></p>
<p class="MsoNormal"> MINSCALEDENOM 10<o:p></o:p></p>
<p class="MsoNormal"> END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> #<o:p></o:p></p>
<p class="MsoNormal"> # Moderate resolution (1 ArcSecond) terrain. Coverage is available worldwide<o:p></o:p></p>
<p class="MsoNormal"> # for most land areas between latitudes -60 and +60.<o:p></o:p></p>
<p class="MsoNormal"> #<o:p></o:p></p>
<p class="MsoNormal"> LAYER<o:p></o:p></p>
<p class="MsoNormal"> NAME "srtm"<o:p></o:p></p>
<p class="MsoNormal"> TILEINDEX "srt2f_1/srt2f_1.shp"<o:p></o:p></p>
<p class="MsoNormal"> TYPE RASTER<o:p></o:p></p>
<p class="MsoNormal"> PROCESSING "RESAMPLE=BICUBIC"<o:p></o:p></p>
<p class="MsoNormal"> STATUS ON<o:p></o:p></p>
<p class="MsoNormal"> <o:p></o:p></p>
<p class="MsoNormal"> METADATA<o:p></o:p></p>
<p class="MsoNormal"> "wms_title" "SRTM 1 ArcSecond Bathmetry"<o:p></o:p></p>
<p class="MsoNormal"> "wms_group_title" "srtm30"<o:p></o:p></p>
<p class="MsoNormal"> "wms_srs" "epsg:4326"<o:p></o:p></p>
<p class="MsoNormal"> "wms_style" "Default"<o:p></o:p></p>
<p class="MsoNormal"> "wms_extent" "-180 -60 180 60"<o:p></o:p></p>
<p class="MsoNormal"> END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"> MAXSCALEDENOM 1000000<o:p></o:p></p>
<p class="MsoNormal"> END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">END<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</body>
</html>