[Geoserver-users] Colormap on Float64 raster

Hi,

Thank you for an excellent description and attachments. I can confirm your experience with Geoserver 2.7.0. Output is all white. From the logs I can find this:

2015-03-27 19:47:36,801 DEBUG [org.geotools.util] - InterpolationConverterFactory can be applied from Strings to Interpolation only.
2015-03-27 19:47:36,802 DEBUG [org.geotools.util] - CRSConverterFactory can be applied from Strings to CRS only.
2015-03-27 19:47:36,802 DEBUG [org.geotools.util] - InterpolationConverterFactory can be applied from Strings to Interpolation only.
2015-03-27 19:47:36,802 DEBUG [org.geotools.util] - CRSConverterFactory can be applied from Strings to CRS only.
2015-03-27 19:47:36,835 TRACE [org.geotools.factory] - ENTRY (GridCoverageFactory)
2015-03-27 19:47:36,835 TRACE [org.geotools.factory] - RETURN (GridCoverageFactory): found implementation GridCoverageFactory.
2015-03-27 19:47:36,836 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Rendering using interpolation javax.media.jai.InterpolationNearest@anonymised.com...
2015-03-27 19:47:36,836 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Rendering using interpolation javax.media.jai.InterpolationNearest@anonymised.com...
2015-03-27 19:47:36,837 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Rendering using interpolation javax.media.jai.InterpolationNearest@anonymised.com...
2015-03-27 19:47:36,840 INFO [org.geotools.coverage.processing] - Apply operation "CoverageCrop" on coverage "float_64" with interpolation "Nearest".
2015-03-27 19:47:36,841 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Drawing coverage GridCoverage2D["float_64", GeneralEnvelope[(149.34300000000002, -35.408), (152.157, -32.592)], DefaultGeographicCRS["WGS 84"]]
? Image=RenderedOp["GTCrop"]

2015-03-27 19:47:36,842 DEBUG [org.geotools.renderer.lite.gridcoverage2d] - Unable to create a granuleDescriptor org.geotools.renderer.lite.gridcoverage2d.GridCoverageRenderer@anonymised.com due to jai scale bug
2015-03-27 19:47:36,843 INFO [org.geoserver.wms] -
Request: getMap
  Time =
  Sld = null
  Buffer = 0
  Format = image/jpeg
  Width = 512
  Filters = null
  Filter = null
  StartIndex = null
  FeatureId = null
  RemoteOwsType = null
  Elevation =
  CQLFilter = null
  FormatOptions = {}
  RemoteOwsURL = null
  ViewParams = null
  Height = 511
  Transparent = true
  StyleFormat = sld
  SldBody = null
  StyleUrl = null
  TilesOrigin = null
  ValidateSchema = false
  Tiled = false
  Palette = null
  FeatureVersion = null
  SRS = EPSG:4326
  Styles = [StyleImpl[ name=dem]]
  Layers = [org.geoserver.wms.MapLayerInfo@anonymised.com]
  MaxFeatures = null
  Crs = GEOGCS["WGS 84",
  DATUM["World Geodetic System 1984",
    SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],
    AUTHORITY["EPSG","6326"]],
  PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
  UNIT["degree", 0.017453292519943295],
  AXIS["Geodetic longitude", EAST],
  AXIS["Geodetic latitude", NORTH],
  AUTHORITY["EPSG","4326"]]
  Bbox = SRSEnvelope[149.34375 : 152.15625, -35.40350341796875 : -32.59649658203125]
  Env = {}
  Angle = 0.0
  BgColor = java.awt.Color[r=255,g=255,b=255]
  Exceptions = SE_XML
  StyleBody = null
  SldVersion = null
  StyleVersion = null
  ScaleMethod = null
  Version = 1.1.1
  BaseUrl = http://localhost:8080/geoserver/
  Get = true
  Request = GetMap
  RawKvp = {REQUEST=GetMap, FORMAT=image/jpeg, SRS=EPSG:4326, BBOX=149.34375,-35.40350341796875,152.15625,-32.59649658203125, VERSION=1.1.1, STYLES=, SERVICE=WMS, WIDTH=512, HEIGHT=511, TRANSPARENT=true, LAYERS=topp:float_64}
  RequestCharset = UTF-8
2015-03-27 19:47:36,849 DEBUG [org.geotools.image] - Encoding input image to write out as JPEG.
2015-03-27 19:47:36,852 DEBUG [org.geotools.image] - Getting a JPEG writer and configuring it.
2015-03-27 19:47:36,852 DEBUG [org.geotools.image] - Writing out...
2015-03-27 19:47:36,889 DEBUG [org.geotools.image] - Writing out... Done!

-Jukka Rahkonen-

-----Alkuperäinen viesti-----
Lähettäjä: CGA [mailto:cgallay.artal@anonymised.com]
Lähetetty: 27. maaliskuuta 2015 17:49
Vastaanottaja: geoserver-users@lists.sourceforge.net
Aihe: [Geoserver-users] Colormap on Float64 raster

Hello,

I'm trying to user GeoServer to publish rasters stored in Float64 type.
For the data to be visually understandable, I put in place an SLD file (using uDig).

Rendering these rasters with uDig and even QGis works well, but when it comes to rendering with geoserver, all i get is an empty image.

Here are the screenshots :
Rendering in QGIS directly from file and created SLD file <http://osgeo-org.1560.x6.nabble.com/file/n5195778/from-file-with-sld.png&gt;
Result of getmap request : all white
<http://osgeo-org.1560.x6.nabble.com/file/n5195778/with-wms.png&gt;

Just to be sure, I tried using WCS to get the raw data and everything seems fine <http://osgeo-org.1560.x6.nabble.com/file/n5195778/with-wcs.png&gt;

Here is the SLD used :

<?xml version="1.0" encoding="UTF-8"?><sld:StyledLayerDescriptor
xmlns="http://www.opengis.net/sld&quot; xmlns:sld="http://www.opengis.net/sld&quot;
xmlns:ogc="http://www.opengis.net/ogc&quot;
xmlns:gml="http://www.opengis.net/gml&quot; version="1.0.0">
    <sld:UserLayer>
        <sld:LayerFeatureConstraints>
            <sld:FeatureTypeConstraint/>
        </sld:LayerFeatureConstraints>
        <sld:UserStyle>
            <sld:Name>1</sld:Name>
            <sld:Title>Default Raster</sld:Title>
            <sld:Abstract>A sample style that draws a raster, good for displaying imagery</sld:Abstract>
            <sld:FeatureTypeStyle>
                <sld:Name>name</sld:Name>
                <sld:FeatureTypeName>Feature</sld:FeatureTypeName>
                <sld:Rule>
                    <sld:Name>rule1</sld:Name>
                    <sld:Title>Opaque Raster</sld:Title>
                    <sld:Abstract>A raster with 100% opacity</sld:Abstract>
                    <sld:RasterSymbolizer>
                        <sld:Geometry>
                            <ogc:PropertyName>geom</ogc:PropertyName>
                        </sld:Geometry>
                        <sld:ColorMap>
                            <sld:ColorMapEntry color="#D7191C" opacity="1"
quantity="0" label="1"/>
                            <sld:ColorMapEntry color="#FDAE61" opacity="1"
quantity="0.002383" label="2"/>
                            <sld:ColorMapEntry color="#FFFFBF" opacity="1"
quantity="0.004766" label="3"/>
                            <sld:ColorMapEntry color="#ABDDA4" opacity="1"
quantity="0.007149" label="4"/>
                            <sld:ColorMapEntry color="#2B83BA" opacity="1"
quantity="0.009532" label="5"/>
                        </sld:ColorMap>
                    </sld:RasterSymbolizer>
                </sld:Rule>
            </sld:FeatureTypeStyle>
        </sld:UserStyle>
    </sld:UserLayer>
</sld:StyledLayerDescriptor>

Here is the raw raster used :
1.tif <http://osgeo-org.1560.x6.nabble.com/file/n5195778/1.tif&gt;

Cyril

--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Colormap-on-Float64-raster-tp5195778.html
Sent from the GeoServer - User mailing list archive at Nabble.com.

------------------------------------------------------------------------------
Dive into the World of Parallel Programming The Go Parallel Website, sponsored by Intel and developed in partnership with Slashdot Media, is your hub for all things parallel software development, from weekly thought leadership blogs to news, videos, case studies, tutorials and more. Take a look and join the conversation now. http://goparallel.sourceforge.net/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users