<div dir="ltr"><div>That's a great solution. Thanks a lot.</div><div><br></div><div>What about "No data" value ? Is there any way to display this value on the getCapabilities so that a client can have this information ?</div><div><br></div><div>Jean-Christophe</div><div><br><div><br></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Mar 29, 2018 at 7:13 PM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div class="HOEnZb"><div class="h5">On jeudi 29 mars 2018 16:36:34 CEST Jean-Christophe Malapert wrote:<br>
> Hi,<br>
><br>
> I have a DEM with the following information:<br>
> PROJCS["SimpleCylindrical Moon",<br>
>     GEOGCS["GCS_Moon",<br>
>         DATUM["D_Moon",<br>
>             SPHEROID["Moon",1737400,0]],<br>
>         PRIMEM["Reference_Meridian",0]<wbr>,<br>
>         UNIT["degree",0.<wbr>0174532925199433]],<br>
>     PROJECTION["Equirectangular"],<br>
>     PARAMETER["latitude_of_origin"<wbr>,0],<br>
>     PARAMETER["central_meridian",<wbr>0],<br>
>     PARAMETER["standard_parallel_<wbr>1",0],<br>
>     PARAMETER["false_easting",0],<br>
>     PARAMETER["false_northing",0],<br>
>     UNIT["metre",1,<br>
>         AUTHORITY["EPSG","9001"]]]<br>
> Origin = (-5458203.076000000350177,<wbr>2729101.538000000175089)<br>
> Pixel Size = (118.450587586805568,-118.<wbr>450587586805568)<br>
> Metadata:<br>
>   AREA_OR_POINT=Area<br>
> Image Structure Metadata:<br>
>   INTERLEAVE=BAND<br>
> Corner Coordinates:<br>
> Upper Left  (-5458203.076, 2729101.538) (180d 0' 0.00"W, 90d 0' 0.00"N)<br>
> Lower Left  (-5458203.076,-2729101.538) (180d 0' 0.00"W, 90d 0' 0.00"S)<br>
> Upper Right ( 5458203.076, 2729101.538) (180d 0' 0.00"E, 90d 0' 0.00"N)<br>
> Lower Right ( 5458203.076,-2729101.538) (180d 0' 0.00"E, 90d 0' 0.00"S)<br>
> Center      (   0.0000000,   0.0000000) (  0d 0' 0.01"E,  0d 0' 0.01"N)<br>
> Band 1 Block=92160x1 Type=Int16, ColorInterp=Gray<br>
>   NoData Value=-32768<br>
>   Offset: 0,   Scale:0.5<br>
><br>
> From my understanding of the DEM information, each pixel value must be<br>
> multiplied by 0.5. If I want to access to my DEM through the WCS protocol,<br>
> how can I say to mapserver to multiply each value by 0.5 so that my client<br>
> displays the right values ? If it is not possible, is there a way to have<br>
> this information in the getCapabilities ?<br>
<br>
</div></div>Jean-Christophe,<br>
<br>
You can create a VRT that will apply the scaling to raw values (ie multiply<br>
them by 0.5)<br>
<br>
gdal_translate your.dem unscaled.vrt -of VRT -unscale<br>
<br>
and serve unscaled.vrt with MapServer<br>
<br>
If yo can spare the disk space for a full unscale copy, you can also<br>
gdal_translate to GeoTIFF.<br>
<span class="HOEnZb"><font color="#888888"><br>
Even<br>
<br>
<br>
--<br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</font></span></blockquote></div><br></div>