[Geoserver-users] Raster/WMS Display Issue

Hello,

I have a 16 MB Geotiff that I am serving up as a WMS feed on Geoserver 1.7.0. I can only get the raster to display after I zoom in to a 1:6 scale. Well it actually displays prior to that but only values greater than one inch (everything else is blank). It seems to be an intermittent problem (the raster is recreated every 20 minutes).

The details on the raster are:

Driver: GTiff/GeoTIFF

Files: rain.tif

Size is 5000, 5000

Coordinate System is:

GEOGCS[“WGS 84”,

DATUM[“WGS_1984”,

SPHEROID[“WGS 84”,6378137,298.2572235629972,

AUTHORITY[“EPSG”,“7030”]],

AUTHORITY[“EPSG”,“6326”]],

PRIMEM[“Greenwich”,0],

UNIT[“degree”,0.0174532925199433],

AUTHORITY[“EPSG”,“4326”]]

Origin = (-126.000000000000000,22.000000000000000)

Pixel Size = (0.012400000000000,0.005600000000000)

Metadata:

AREA_OR_POINT=Area

Image Structure Metadata:

COMPRESSION=LZW

INTERLEAVE=BAND

Corner Coordinates:

Upper Left (-126.0000000, 22.0000000) (126d 0’0.00"W, 22d 0’0.00"N)

Lower Left (-126.0000000, 50.0000000) (126d 0’0.00"W, 50d 0’0.00"N)

Upper Right ( -64.0000000, 22.0000000) ( 64d 0’0.00"W, 22d 0’0.00"N)

Lower Right ( -64.0000000, 50.0000000) ( 64d 0’0.00"W, 50d 0’0.00"N)

Center ( -95.0000000, 36.0000000) ( 95d 0’0.00"W, 36d 0’0.00"N)

Band 1 Block=256x256 Type=Float32, ColorInterp=Gray

Computed Min/Max=0.000,2.482

Minimum=0.000, Maximum=2.482, Mean=0.009, StdDev=0.045

Metadata:

STATISTICS_MINIMUM=0

STATISTICS_MAXIMUM=2.4823353290558

STATISTICS_MEAN=0.0089933774132436

STATISTICS_STDDEV=0.044540742213422

The raster style is:

<?xml version="1.0" encoding="ISO-8859-1"?>

<StyledLayerDescriptor version=“1.0.0”

xsi:schemaLocation=“http://www.opengis.net/sld StyledLayerDescriptor.xsd”

xmlns=“http://www.opengis.net/sld

xmlns:ogc=“http://www.opengis.net/ogc

xmlns:xlink=“http://www.w3.org/1999/xlink

xmlns:xsi=“http://www.w3.org/2001/XMLSchema-instance”>

precip_raster

Precip Raster Colormap

Colors for precipitation raster image

1.0

Thanks.

bw

Hi Brian,
Can you share your geotiff so I can try to reproduce your problem on my machine?
If you have a link, an ftp space containing that gtiff or something else, let me know.

Regards,
Daniele.

On Wed, Dec 3, 2008 at 9:35 PM, Brian Walawender <Brian.Walawender@anonymised.com> wrote:

Hello,

I have a 16 MB Geotiff that I am serving up as a WMS feed on Geoserver 1.7.0. I can only get the raster to display after I zoom in to a 1:6 scale. Well it actually displays prior to that but only values greater than one inch (everything else is blank). It seems to be an intermittent problem (the raster is recreated every 20 minutes).

The details on the raster are:

Driver: GTiff/GeoTIFF

Files: rain.tif

Size is 5000, 5000

Coordinate System is:

GEOGCS[“WGS 84”,

DATUM[“WGS_1984”,

SPHEROID[“WGS 84”,6378137,298.2572235629972,

AUTHORITY[“EPSG”,“7030”]],

AUTHORITY[“EPSG”,“6326”]],

PRIMEM[“Greenwich”,0],

UNIT[“degree”,0.0174532925199433],

AUTHORITY[“EPSG”,“4326”]]

Origin = (-126.000000000000000,22.000000000000000)

Pixel Size = (0.012400000000000,0.005600000000000)

Metadata:

AREA_OR_POINT=Area

Image Structure Metadata:

COMPRESSION=LZW

INTERLEAVE=BAND

Corner Coordinates:

Upper Left (-126.0000000, 22.0000000) (126d 0’0.00"W, 22d 0’0.00"N)

Lower Left (-126.0000000, 50.0000000) (126d 0’0.00"W, 50d 0’0.00"N)

Upper Right ( -64.0000000, 22.0000000) ( 64d 0’0.00"W, 22d 0’0.00"N)

Lower Right ( -64.0000000, 50.0000000) ( 64d 0’0.00"W, 50d 0’0.00"N)

Center ( -95.0000000, 36.0000000) ( 95d 0’0.00"W, 36d 0’0.00"N)

Band 1 Block=256x256 Type=Float32, ColorInterp=Gray

Computed Min/Max=0.000,2.482

Minimum=0.000, Maximum=2.482, Mean=0.009, StdDev=0.045

Metadata:

STATISTICS_MINIMUM=0

STATISTICS_MAXIMUM=2.4823353290558

STATISTICS_MEAN=0.0089933774132436

STATISTICS_STDDEV=0.044540742213422

The raster style is:

<?xml version="1.0" encoding="ISO-8859-1"?>

<StyledLayerDescriptor version=“1.0.0”

xsi:schemaLocation=“http://www.opengis.net/sld StyledLayerDescriptor.xsd”

xmlns=“http://www.opengis.net/sld

xmlns:ogc=“http://www.opengis.net/ogc

xmlns:xlink=“http://www.w3.org/1999/xlink

xmlns:xsi=“http://www.w3.org/2001/XMLSchema-instance”>

precip_raster

Precip Raster Colormap

Colors for precipitation raster image

1.0

Thanks.

bw


This SF.Net email is sponsored by the Moblin Your Move Developer’s challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/


Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

Eng. Daniele Romagnoli
Software Engineer

GeoSolutions S.A.S.
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 328 0559267

http://www.geo-solutions.it


Ciao Brian,
I am pretty sure you have run into a JAI-ImageIO tiff bug (actually a
combination of two bugs) as I reported a while back. As Daniele
mentioned, it would be great to get some sample data to make sure it
is not a bug somewhere in geotools. I guess it is time to push out an
imageio-ext release and to include the improved imageio-tiff
reader/writer, what do you think daniele?

Simone.
-------------------------------------------------------
Ing. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Wed, Dec 3, 2008 at 9:35 PM, Brian Walawender
<Brian.Walawender@anonymised.com> wrote:

Hello,

I have a 16 MB Geotiff that I am serving up as a WMS feed on Geoserver
1.7.0. I can only get the raster to display after I zoom in to a 1:6
scale. Well it actually displays prior to that but only values greater
than one inch (everything else is blank). It seems to be an intermittent
problem (the raster is recreated every 20 minutes).

The details on the raster are:

Driver: GTiff/GeoTIFF

Files: rain.tif

Size is 5000, 5000

Coordinate System is:

GEOGCS["WGS 84",

    DATUM["WGS_1984",

        SPHEROID["WGS 84",6378137,298.2572235629972,

            AUTHORITY["EPSG","7030"]],

        AUTHORITY["EPSG","6326"]],

    PRIMEM["Greenwich",0],

    UNIT["degree",0.0174532925199433],

    AUTHORITY["EPSG","4326"]]

Origin = (-126.000000000000000,22.000000000000000)

Pixel Size = (0.012400000000000,0.005600000000000)

Metadata:

  AREA_OR_POINT=Area

Image Structure Metadata:

  COMPRESSION=LZW

  INTERLEAVE=BAND

Corner Coordinates:

Upper Left (-126.0000000, 22.0000000) (126d 0'0.00"W, 22d 0'0.00"N)

Lower Left (-126.0000000, 50.0000000) (126d 0'0.00"W, 50d 0'0.00"N)

Upper Right ( -64.0000000, 22.0000000) ( 64d 0'0.00"W, 22d 0'0.00"N)

Lower Right ( -64.0000000, 50.0000000) ( 64d 0'0.00"W, 50d 0'0.00"N)

Center ( -95.0000000, 36.0000000) ( 95d 0'0.00"W, 36d 0'0.00"N)

Band 1 Block=256x256 Type=Float32, ColorInterp=Gray

    Computed Min/Max=0.000,2.482

  Minimum=0.000, Maximum=2.482, Mean=0.009, StdDev=0.045

  Metadata:

    STATISTICS_MINIMUM=0

    STATISTICS_MAXIMUM=2.4823353290558

    STATISTICS_MEAN=0.0089933774132436

    STATISTICS_STDDEV=0.044540742213422

The raster style is:

<?xml version="1.0" encoding="ISO-8859-1"?>

<StyledLayerDescriptor version="1.0.0"

    xsi:schemaLocation="http://www.opengis.net/sld
StyledLayerDescriptor.xsd"

    xmlns="http://www.opengis.net/sld&quot;

    xmlns:ogc="http://www.opengis.net/ogc&quot;

    xmlns:xlink="http://www.w3.org/1999/xlink&quot;

    xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance&quot;&gt;

    <!-- a named layer is the basic building block of an sld document -->

    <NamedLayer>

        <Name>precip_raster</Name>

        <UserStyle>

        <!-- they have names, titles and abstracts -->

            <Title>Precip Raster Colormap</Title>

            <Abstract>Colors for precipitation raster image</Abstract>

            <FeatureTypeStyle>

                    <Rule>

                                <RasterSymbolizer>

                                    <Opacity>1.0</Opacity>

                                    <ColorMap type="Intervals">

                                       <ColorMapEntry color="#FFFFFF"
quantity="0.00" opacity="0.0" />

                                       <ColorMapEntry color="#FFFFFF"
quantity="0.0099" opacity="0.0" />

                                       <ColorMapEntry color="#7FFF00"
quantity="0.099"/>

                                       <ColorMapEntry color="#00CD00"
quantity="0.249"/>

                                       <ColorMapEntry color="#008B00"
quantity="0.499"/>

                                       <ColorMapEntry color="#104E8B"
quantity="0.749"/>

                                       <ColorMapEntry color="#1E90FF"
quantity="0.999"/>

                                       <ColorMapEntry color="#00B2EE"
quantity="1.249"/>

                                       <ColorMapEntry color="#00EEEE"
quantity="1.499"/>

                                       <ColorMapEntry color="#8968CD"
quantity="1.749"/>

                                       <ColorMapEntry color="#912CEE"
quantity="1.999"/>

                                       <ColorMapEntry color="#8B008B"
quantity="2.499"/>

                                       <ColorMapEntry color="#8B0000"
quantity="2.999"/>

                                       <ColorMapEntry color="#CD0000"
quantity="3.999"/>

                                       <ColorMapEntry color="#EE4000"
quantity="4.999"/>

                                       <ColorMapEntry color="#FF7F00"
quantity="5.999"/>

                                       <ColorMapEntry color="#CD8500"
quantity="6.999"/>

                                       <ColorMapEntry color="#EEEE00"
quantity="7.999"/>

                                       <ColorMapEntry color="#FFFF00"
quantity="8.999"/>

                                    </ColorMap>

                                </RasterSymbolizer>

                </Rule>

            </FeatureTypeStyle>

        </UserStyle>

    </NamedLayer>

</StyledLayerDescriptor>

Thanks.

bw

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great
prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

Daniele and Simone,

Here is a link to the raster I had the problems with:

http://www.crh.noaa.gov/sharedoc/rain.tif

Let me know if you need more information on how it was created.

Thanks for your quick response.

bw

----- Original Message -----
From: Daniele Romagnoli <dany.geotools@anonymised.com>
Date: Wednesday, December 3, 2008 5:28 pm
Subject: Re: [Geoserver-users] Raster/WMS Display Issue
To: Brian.Walawender@anonymised.com
Cc: geoserver-users@lists.sourceforge.net

Hi Brian,
Can you share your geotiff so I can try to reproduce your problem on my
machine?
If you have a link, an ftp space containing that gtiff or something else,
let me know.

Regards,
Daniele.
On Wed, Dec 3, 2008 at 9:35 PM, Brian Walawender
<Brian.Walawender@anonymised.com>wrote:

> Hello,
>
>
>
> I have a 16 MB Geotiff that I am serving up as a WMS feed on Geoserver
> 1.7.0. I can only get the raster to display after I zoom in to a
1:6
> scale. Well it actually displays prior to that but only values greater
> than one inch (everything else is blank). It seems to be an intermittent
> problem (the raster is recreated every 20 minutes).
>
>
>
> The details on the raster are:
>
>
>
> Driver: GTiff/GeoTIFF
>
> Files: rain.tif
>
> Size is 5000, 5000
>
> Coordinate System is:
>
> GEOGCS["WGS 84",
>
> DATUM["WGS_1984",
>
> SPHEROID["WGS 84",6378137,298.2572235629972,
>
> AUTHORITY["EPSG","7030"]],
>
> AUTHORITY["EPSG","6326"]],
>
> PRIMEM["Greenwich",0],
>
> UNIT["degree",0.0174532925199433],
>
> AUTHORITY["EPSG","4326"]]
>
> Origin = (-126.000000000000000,22.000000000000000)
>
> Pixel Size = (0.012400000000000,0.005600000000000)
>
> Metadata:
>
> AREA_OR_POINT=Area
>
> Image Structure Metadata:
>
> COMPRESSION=LZW
>
> INTERLEAVE=BAND
>
> Corner Coordinates:
>
> Upper Left (-126.0000000, 22.0000000) (126d 0'0.00"W, 22d 0'0.00"N)
>
> Lower Left (-126.0000000, 50.0000000) (126d 0'0.00"W, 50d 0'0.00"N)
>
> Upper Right ( -64.0000000, 22.0000000) ( 64d 0'0.00"W, 22d 0'0.00"N)
>
> Lower Right ( -64.0000000, 50.0000000) ( 64d 0'0.00"W, 50d 0'0.00"N)
>
> Center ( -95.0000000, 36.0000000) ( 95d 0'0.00"W, 36d 0'0.00"N)
>
> Band 1 Block=256x256 Type=Float32, ColorInterp=Gray
>
> Computed Min/Max=0.000,2.482
>
> Minimum=0.000, Maximum=2.482, Mean=0.009, StdDev=0.045
>
> Metadata:
>
> STATISTICS_MINIMUM=0
>
> STATISTICS_MAXIMUM=2.4823353290558
>
> STATISTICS_MEAN=0.0089933774132436
>
> STATISTICS_STDDEV=0.044540742213422
>
>
>
> The raster style is:
>
>
>
> <?xml version="1.0" encoding="ISO-8859-1"?>
>
> <StyledLayerDescriptor version="1.0.0"
>
> xsi:schemaLocation=""
>
> xmlns=""
>
> xmlns:ogc=""
>
> xmlns:xlink=""
>
> xmlns:xsi="">
>
> <!-- a named layer is the basic building block of an sld
document -->
>
>
>
> <NamedLayer>
>
> <Name>precip_raster</Name>
>
> <UserStyle>
>
> <!-- they have names, titles and abstracts -->
>
> <Title>Precip Raster Colormap</Title>
>
> <Abstract>Colors for precipitation raster image</Abstract>
>
> <FeatureTypeStyle>
>
> <Rule>
>
> <RasterSymbolizer>
>
> <Opacity>1.0</Opacity>
>
> <ColorMap type="Intervals">
>
> <ColorMapEntry color="#FFFFFF"
> quantity="0.00" opacity="0.0" />
>
> <ColorMapEntry color="#FFFFFF"
> quantity="0.0099" opacity="0.0" />
>
> <ColorMapEntry color="#7FFF00"
> quantity="0.099"/>
>
> <ColorMapEntry color="#00CD00"
> quantity="0.249"/>
>
> <ColorMapEntry color="#008B00"
> quantity="0.499"/>
>
> <ColorMapEntry color="#104E8B"
> quantity="0.749"/>
>
> <ColorMapEntry color="#1E90FF"
> quantity="0.999"/>
>
> <ColorMapEntry color="#00B2EE"
> quantity="1.249"/>
>
> <ColorMapEntry color="#00EEEE"
> quantity="1.499"/>
>
> <ColorMapEntry color="#8968CD"
> quantity="1.749"/>
>
> <ColorMapEntry color="#912CEE"
> quantity="1.999"/>
>
> <ColorMapEntry color="#8B008B"
> quantity="2.499"/>
>
> <ColorMapEntry color="#8B0000"
> quantity="2.999"/>
>
> <ColorMapEntry color="#CD0000"
> quantity="3.999"/>
>
> <ColorMapEntry color="#EE4000"
> quantity="4.999"/>
>
> <ColorMapEntry color="#FF7F00"
> quantity="5.999"/>
>
> <ColorMapEntry color="#CD8500"
> quantity="6.999"/>
>
> <ColorMapEntry color="#EEEE00"
> quantity="7.999"/>
>
> <ColorMapEntry color="#FFFF00"
> quantity="8.999"/>
>
> </ColorMap>
>
> </RasterSymbolizer>
>
> </Rule>
>
> </FeatureTypeStyle>
>
> </UserStyle>
>
> </NamedLayer>
>
> </StyledLayerDescriptor>
>
>
>
> Thanks.
>
>
>
> bw
>
>
>
>
>
> -------------------------------------------------------------------------
> This SF.Net email is sponsored by the Moblin Your Move Developer's
> challenge
> Build the coolest Linux based applications with Moblin SDK & win great
> prizes
> Grand prize is a trip for two to an Open Source event anywhere in
the world
>
> _______________________________________________
> Geoserver-users mailing list
> Geoserver-users@lists.sourceforge.net
>
>
>

--
-------------------------------------------------------
Eng. Daniele Romagnoli
Software Engineer

GeoSolutions S.A.S.
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 328 0559267

-------------------------------------------------------
-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great
prizes
Grand prize is a trip for two to an Open Source event anywhere in the
world
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net

Simone and Daniele,

One more comment, it looks like the geotiff values are being truncated to
integers at the initial zooms and then being displayed as a float once you
zoom in.

bw

-----Original Message-----
From: Simone Giannecchini [mailto:simone.giannecchini@anonymised.com]
Sent: Wednesday, December 03, 2008 6:01 PM
To: Brian.Walawender@anonymised.com
Cc: geoserver-users@lists.sourceforge.net
Subject: Re: [Geoserver-users] Raster/WMS Display Issue

Ciao Brian,
I am pretty sure you have run into a JAI-ImageIO tiff bug (actually a
combination of two bugs) as I reported a while back. As Daniele
mentioned, it would be great to get some sample data to make sure it
is not a bug somewhere in geotools. I guess it is time to push out an
imageio-ext release and to include the improved imageio-tiff
reader/writer, what do you think daniele?

Simone.
-------------------------------------------------------
Ing. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Wed, Dec 3, 2008 at 9:35 PM, Brian Walawender
<Brian.Walawender@anonymised.com> wrote:

Hello,

I have a 16 MB Geotiff that I am serving up as a WMS feed on Geoserver
1.7.0. I can only get the raster to display after I zoom in to a 1:6
scale. Well it actually displays prior to that but only values greater
than one inch (everything else is blank). It seems to be an intermittent
problem (the raster is recreated every 20 minutes).

The details on the raster are:

Driver: GTiff/GeoTIFF

Files: rain.tif

Size is 5000, 5000

Coordinate System is:

GEOGCS["WGS 84",

    DATUM["WGS_1984",

        SPHEROID["WGS 84",6378137,298.2572235629972,

            AUTHORITY["EPSG","7030"]],

        AUTHORITY["EPSG","6326"]],

    PRIMEM["Greenwich",0],

    UNIT["degree",0.0174532925199433],

    AUTHORITY["EPSG","4326"]]

Origin = (-126.000000000000000,22.000000000000000)

Pixel Size = (0.012400000000000,0.005600000000000)

Metadata:

  AREA_OR_POINT=Area

Image Structure Metadata:

  COMPRESSION=LZW

  INTERLEAVE=BAND

Corner Coordinates:

Upper Left (-126.0000000, 22.0000000) (126d 0'0.00"W, 22d 0'0.00"N)

Lower Left (-126.0000000, 50.0000000) (126d 0'0.00"W, 50d 0'0.00"N)

Upper Right ( -64.0000000, 22.0000000) ( 64d 0'0.00"W, 22d 0'0.00"N)

Lower Right ( -64.0000000, 50.0000000) ( 64d 0'0.00"W, 50d 0'0.00"N)

Center ( -95.0000000, 36.0000000) ( 95d 0'0.00"W, 36d 0'0.00"N)

Band 1 Block=256x256 Type=Float32, ColorInterp=Gray

    Computed Min/Max=0.000,2.482

  Minimum=0.000, Maximum=2.482, Mean=0.009, StdDev=0.045

  Metadata:

    STATISTICS_MINIMUM=0

    STATISTICS_MAXIMUM=2.4823353290558

    STATISTICS_MEAN=0.0089933774132436

    STATISTICS_STDDEV=0.044540742213422

The raster style is:

<?xml version="1.0" encoding="ISO-8859-1"?>

<StyledLayerDescriptor version="1.0.0"

    xsi:schemaLocation="http://www.opengis.net/sld
StyledLayerDescriptor.xsd"

    xmlns="http://www.opengis.net/sld&quot;

    xmlns:ogc="http://www.opengis.net/ogc&quot;

    xmlns:xlink="http://www.w3.org/1999/xlink&quot;

    xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance&quot;&gt;

    <!-- a named layer is the basic building block of an sld document -->

    <NamedLayer>

        <Name>precip_raster</Name>

        <UserStyle>

        <!-- they have names, titles and abstracts -->

            <Title>Precip Raster Colormap</Title>

            <Abstract>Colors for precipitation raster image</Abstract>

            <FeatureTypeStyle>

                    <Rule>

                                <RasterSymbolizer>

                                    <Opacity>1.0</Opacity>

                                    <ColorMap type="Intervals">

                                       <ColorMapEntry color="#FFFFFF"
quantity="0.00" opacity="0.0" />

                                       <ColorMapEntry color="#FFFFFF"
quantity="0.0099" opacity="0.0" />

                                       <ColorMapEntry color="#7FFF00"
quantity="0.099"/>

                                       <ColorMapEntry color="#00CD00"
quantity="0.249"/>

                                       <ColorMapEntry color="#008B00"
quantity="0.499"/>

                                       <ColorMapEntry color="#104E8B"
quantity="0.749"/>

                                       <ColorMapEntry color="#1E90FF"
quantity="0.999"/>

                                       <ColorMapEntry color="#00B2EE"
quantity="1.249"/>

                                       <ColorMapEntry color="#00EEEE"
quantity="1.499"/>

                                       <ColorMapEntry color="#8968CD"
quantity="1.749"/>

                                       <ColorMapEntry color="#912CEE"
quantity="1.999"/>

                                       <ColorMapEntry color="#8B008B"
quantity="2.499"/>

                                       <ColorMapEntry color="#8B0000"
quantity="2.999"/>

                                       <ColorMapEntry color="#CD0000"
quantity="3.999"/>

                                       <ColorMapEntry color="#EE4000"
quantity="4.999"/>

                                       <ColorMapEntry color="#FF7F00"
quantity="5.999"/>

                                       <ColorMapEntry color="#CD8500"
quantity="6.999"/>

                                       <ColorMapEntry color="#EEEE00"
quantity="7.999"/>

                                       <ColorMapEntry color="#FFFF00"
quantity="8.999"/>

                                    </ColorMap>

                                </RasterSymbolizer>

                </Rule>

            </FeatureTypeStyle>

        </UserStyle>

    </NamedLayer>

</StyledLayerDescriptor>

Thanks.

bw

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's

challenge

Build the coolest Linux based applications with Moblin SDK & win great
prizes
Grand prize is a trip for two to an Open Source event anywhere in the

world

http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great
prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

That's exactly the erroneous behavior I have fixed in the SUN tiff reader.

Simone.
-------------------------------------------------------
Ing. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Thu, Dec 4, 2008 at 3:28 PM, Brian Walawender
<Brian.Walawender@anonymised.com> wrote:

Simone and Daniele,

One more comment, it looks like the geotiff values are being truncated to
integers at the initial zooms and then being displayed as a float once you
zoom in.

bw

-----Original Message-----
From: Simone Giannecchini [mailto:simone.giannecchini@anonymised.com]
Sent: Wednesday, December 03, 2008 6:01 PM
To: Brian.Walawender@anonymised.com
Cc: geoserver-users@lists.sourceforge.net
Subject: Re: [Geoserver-users] Raster/WMS Display Issue

Ciao Brian,
I am pretty sure you have run into a JAI-ImageIO tiff bug (actually a
combination of two bugs) as I reported a while back. As Daniele
mentioned, it would be great to get some sample data to make sure it
is not a bug somewhere in geotools. I guess it is time to push out an
imageio-ext release and to include the improved imageio-tiff
reader/writer, what do you think daniele?

Simone.
-------------------------------------------------------
Ing. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Wed, Dec 3, 2008 at 9:35 PM, Brian Walawender
<Brian.Walawender@anonymised.com> wrote:

Hello,

I have a 16 MB Geotiff that I am serving up as a WMS feed on Geoserver
1.7.0. I can only get the raster to display after I zoom in to a 1:6
scale. Well it actually displays prior to that but only values greater
than one inch (everything else is blank). It seems to be an intermittent
problem (the raster is recreated every 20 minutes).

The details on the raster are:

Driver: GTiff/GeoTIFF

Files: rain.tif

Size is 5000, 5000

Coordinate System is:

GEOGCS["WGS 84",

    DATUM["WGS_1984",

        SPHEROID["WGS 84",6378137,298.2572235629972,

            AUTHORITY["EPSG","7030"]],

        AUTHORITY["EPSG","6326"]],

    PRIMEM["Greenwich",0],

    UNIT["degree",0.0174532925199433],

    AUTHORITY["EPSG","4326"]]

Origin = (-126.000000000000000,22.000000000000000)

Pixel Size = (0.012400000000000,0.005600000000000)

Metadata:

  AREA_OR_POINT=Area

Image Structure Metadata:

  COMPRESSION=LZW

  INTERLEAVE=BAND

Corner Coordinates:

Upper Left (-126.0000000, 22.0000000) (126d 0'0.00"W, 22d 0'0.00"N)

Lower Left (-126.0000000, 50.0000000) (126d 0'0.00"W, 50d 0'0.00"N)

Upper Right ( -64.0000000, 22.0000000) ( 64d 0'0.00"W, 22d 0'0.00"N)

Lower Right ( -64.0000000, 50.0000000) ( 64d 0'0.00"W, 50d 0'0.00"N)

Center ( -95.0000000, 36.0000000) ( 95d 0'0.00"W, 36d 0'0.00"N)

Band 1 Block=256x256 Type=Float32, ColorInterp=Gray

    Computed Min/Max=0.000,2.482

  Minimum=0.000, Maximum=2.482, Mean=0.009, StdDev=0.045

  Metadata:

    STATISTICS_MINIMUM=0

    STATISTICS_MAXIMUM=2.4823353290558

    STATISTICS_MEAN=0.0089933774132436

    STATISTICS_STDDEV=0.044540742213422

The raster style is:

<?xml version="1.0" encoding="ISO-8859-1"?>

<StyledLayerDescriptor version="1.0.0"

    xsi:schemaLocation="http://www.opengis.net/sld
StyledLayerDescriptor.xsd"

    xmlns="http://www.opengis.net/sld&quot;

    xmlns:ogc="http://www.opengis.net/ogc&quot;

    xmlns:xlink="http://www.w3.org/1999/xlink&quot;

    xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance&quot;&gt;

    <!-- a named layer is the basic building block of an sld document -->

    <NamedLayer>

        <Name>precip_raster</Name>

        <UserStyle>

        <!-- they have names, titles and abstracts -->

            <Title>Precip Raster Colormap</Title>

            <Abstract>Colors for precipitation raster image</Abstract>

            <FeatureTypeStyle>

                    <Rule>

                                <RasterSymbolizer>

                                    <Opacity>1.0</Opacity>

                                    <ColorMap type="Intervals">

                                       <ColorMapEntry color="#FFFFFF"
quantity="0.00" opacity="0.0" />

                                       <ColorMapEntry color="#FFFFFF"
quantity="0.0099" opacity="0.0" />

                                       <ColorMapEntry color="#7FFF00"
quantity="0.099"/>

                                       <ColorMapEntry color="#00CD00"
quantity="0.249"/>

                                       <ColorMapEntry color="#008B00"
quantity="0.499"/>

                                       <ColorMapEntry color="#104E8B"
quantity="0.749"/>

                                       <ColorMapEntry color="#1E90FF"
quantity="0.999"/>

                                       <ColorMapEntry color="#00B2EE"
quantity="1.249"/>

                                       <ColorMapEntry color="#00EEEE"
quantity="1.499"/>

                                       <ColorMapEntry color="#8968CD"
quantity="1.749"/>

                                       <ColorMapEntry color="#912CEE"
quantity="1.999"/>

                                       <ColorMapEntry color="#8B008B"
quantity="2.499"/>

                                       <ColorMapEntry color="#8B0000"
quantity="2.999"/>

                                       <ColorMapEntry color="#CD0000"
quantity="3.999"/>

                                       <ColorMapEntry color="#EE4000"
quantity="4.999"/>

                                       <ColorMapEntry color="#FF7F00"
quantity="5.999"/>

                                       <ColorMapEntry color="#CD8500"
quantity="6.999"/>

                                       <ColorMapEntry color="#EEEE00"
quantity="7.999"/>

                                       <ColorMapEntry color="#FFFF00"
quantity="8.999"/>

                                    </ColorMap>

                                </RasterSymbolizer>

                </Rule>

            </FeatureTypeStyle>

        </UserStyle>

    </NamedLayer>

</StyledLayerDescriptor>

Thanks.

bw

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's

challenge

Build the coolest Linux based applications with Moblin SDK & win great
prizes
Grand prize is a trip for two to an Open Source event anywhere in the

world

http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great
prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

This is also similar to the behavior I have posted on the list before:

Hi List,

I’m fairly new to Geoserver, and I’m very impressed by how easy it is to
get up and running quickly.

I’ve been playing with serving data from Geoserver’s WMS. I’m using the
stand alone Geoserver 1.7.0-RC4 on linux with jdk 1.6.0_06 and native jai.

I’ve been noticing some strange image rendering behavior when generating
WMS images of radar data. This data contains many polygons (~80,000) and
is in a custom orthographic projection. The polygons near the center of
the data are much smaller than the ones near the perimeter (since it is
radial data).

Many of the small polygons are not rendered or are rendered in a way
that leaves a strange pattern. You can see this in
http://venus1.wrh.noaa.gov/scratch/sutula/wms1.png. That image was
generated with the following WMS request:

http://localhost:8080/geoserver/wms?WIDTH=800&HEIGHT=740&LAYERS=topp:states,topp%3ABRO_20080723_1807&STYLES=&SRS=EPSG%3A500001&FORMAT=image%2Fpng&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&EXCEPTIONS=application%2Fvnd.ogc.se_inimage&FORMAT_OPTIONS=antialias%3Anone&BBOX=-268876.70491178706,-160532.58400467926,134399.955480083,212498.32685780057

Strangely, if I make the output image half the size, the data does not
appear at all as seen in
http://venus1.wrh.noaa.gov/scratch/sutula/wms2.png. That was generated with:

http://localhost:8080/geoserver/wms?WIDTH=400&HEIGHT=370&LAYERS=topp:states,topp%3ABRO_20080723_1807&STYLES=&SRS=EPSG%3A500001&FORMAT=image%2Fpng&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&EXCEPTIONS=application%2Fvnd.ogc.se_inimage&FORMAT_OPTIONS=antialias%3Anone&BBOX=-268876.70491178706,-160532.58400467926,134399.955480083,212498.32685780057

My sld has no scale dependence and only specifies a fill color and
opacity for different values. The behavior happens with or without anti
aliasing enabled.

Any ideas what is going on here?

Thank you,
Aaron

In summary, the data does not render at all when you zoom out a bit, and when it does render, some pixels are not filled in where the features being rendered are relatively small. I see the same behavior in Udig, but the data renders correctly in other software like Mapserver. It would be wonderful to get this fixed.

Aaron

Simone Giannecchini wrote:

Ciao Aaron,
I have fixed this like two months ago, but the problem is that the fix
is in the ImageIO code not in imageio-ext or geotools or geoserver,
hence we can control it, which means that we cannot apply a patch
ourselves to one of our projects and make a new release but we need to
wait for sun to do that. I can send you an update version of the
ImageIO jar so that you could test it while we take a decision on how
to proceed. However I think that we will just ship the update imageio
tiff l ibrary in imageio-ext.

Simone.
-------------------------------------------------------
Ing. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Thu, Dec 4, 2008 at 4:31 PM, Aaron Sutula <Aaron.Sutula@anonymised.com> wrote:

This is also similar to the behavior I have posted on the list before:

Hi List,

I'm fairly new to Geoserver, and I'm very impressed by how easy it is to
get up and running quickly.

I've been playing with serving data from Geoserver's WMS. I'm using the
stand alone Geoserver 1.7.0-RC4 on linux with jdk 1.6.0_06 and native jai.

I've been noticing some strange image rendering behavior when generating
WMS images of radar data. This data contains many polygons (~80,000) and
is in a custom orthographic projection. The polygons near the center of
the data are much smaller than the ones near the perimeter (since it is
radial data).

Many of the small polygons are not rendered or are rendered in a way
that leaves a strange pattern. You can see this in
http://venus1.wrh.noaa.gov/scratch/sutula/wms1.png. That image was
generated with the following WMS request:

http://localhost:8080/geoserver/wms?WIDTH=800&HEIGHT=740&LAYERS=topp:states,topp%3ABRO_20080723_1807&STYLES=&SRS=EPSG%3A500001&FORMAT=image%2Fpng&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&EXCEPTIONS=application%2Fvnd.ogc.se_inimage&FORMAT_OPTIONS=antialias%3Anone&BBOX=-268876.70491178706,-160532.58400467926,134399.955480083,212498.32685780057

Strangely, if I make the output image half the size, the data does not
appear at all as seen in
http://venus1.wrh.noaa.gov/scratch/sutula/wms2.png. That was generated with:

http://localhost:8080/geoserver/wms?WIDTH=400&HEIGHT=370&LAYERS=topp:states,topp%3ABRO_20080723_1807&STYLES=&SRS=EPSG%3A500001&FORMAT=image%2Fpng&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&EXCEPTIONS=application%2Fvnd.ogc.se_inimage&FORMAT_OPTIONS=antialias%3Anone&BBOX=-268876.70491178706,-160532.58400467926,134399.955480083,212498.32685780057

My sld has no scale dependence and only specifies a fill color and
opacity for different values. The behavior happens with or without anti
aliasing enabled.

Any ideas what is going on here?

Thank you,
Aaron

In summary, the data does not render at all when you zoom out a bit, and
when it does render, some pixels are not filled in where the features being
rendered are relatively small. I see the same behavior in Udig, but the data
renders correctly in other software like Mapserver. It would be wonderful to
get this fixed.

Aaron

Simone Giannecchini wrote:

That's exactly the erroneous behavior I have fixed in the SUN tiff reader.

Simone.
-------------------------------------------------------
Ing. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Thu, Dec 4, 2008 at 3:28 PM, Brian Walawender
<Brian.Walawender@anonymised.com> wrote:

Simone and Daniele,

One more comment, it looks like the geotiff values are being truncated to
integers at the initial zooms and then being displayed as a float once you
zoom in.

bw

-----Original Message-----
From: Simone Giannecchini [mailto:simone.giannecchini@anonymised.com]
Sent: Wednesday, December 03, 2008 6:01 PM
To: Brian.Walawender@anonymised.com
Cc: geoserver-users@lists.sourceforge.net
Subject: Re: [Geoserver-users] Raster/WMS Display Issue

Ciao Brian,
I am pretty sure you have run into a JAI-ImageIO tiff bug (actually a
combination of two bugs) as I reported a while back. As Daniele
mentioned, it would be great to get some sample data to make sure it
is not a bug somewhere in geotools. I guess it is time to push out an
imageio-ext release and to include the improved imageio-tiff
reader/writer, what do you think daniele?

Simone.
-------------------------------------------------------
Ing. Simone Giannecchini
GeoSolutions S.A.S.
Owner - Software Engineer
Via Carignoni 51
55041 Camaiore (LU)
Italy

phone: +39 0584983027
fax: +39 0584983027
mob: +39 333 8128928

http://www.geo-solutions.it
http://simboss.blogspot.com/
http://www.linkedin.com/in/simonegiannecchini

-------------------------------------------------------

On Wed, Dec 3, 2008 at 9:35 PM, Brian Walawender
<Brian.Walawender@anonymised.com> wrote:

Hello,

I have a 16 MB Geotiff that I am serving up as a WMS feed on Geoserver
1.7.0. I can only get the raster to display after I zoom in to a 1:6
scale. Well it actually displays prior to that but only values greater
than one inch (everything else is blank). It seems to be an intermittent
problem (the raster is recreated every 20 minutes).

The details on the raster are:

Driver: GTiff/GeoTIFF

Files: rain.tif

Size is 5000, 5000

Coordinate System is:

GEOGCS["WGS 84",

    DATUM["WGS_1984",

        SPHEROID["WGS 84",6378137,298.2572235629972,

            AUTHORITY["EPSG","7030"]],

        AUTHORITY["EPSG","6326"]],

    PRIMEM["Greenwich",0],

    UNIT["degree",0.0174532925199433],

    AUTHORITY["EPSG","4326"]]

Origin = (-126.000000000000000,22.000000000000000)

Pixel Size = (0.012400000000000,0.005600000000000)

Metadata:

  AREA_OR_POINT=Area

Image Structure Metadata:

  COMPRESSION=LZW

  INTERLEAVE=BAND

Corner Coordinates:

Upper Left (-126.0000000, 22.0000000) (126d 0'0.00"W, 22d 0'0.00"N)

Lower Left (-126.0000000, 50.0000000) (126d 0'0.00"W, 50d 0'0.00"N)

Upper Right ( -64.0000000, 22.0000000) ( 64d 0'0.00"W, 22d 0'0.00"N)

Lower Right ( -64.0000000, 50.0000000) ( 64d 0'0.00"W, 50d 0'0.00"N)

Center ( -95.0000000, 36.0000000) ( 95d 0'0.00"W, 36d 0'0.00"N)

Band 1 Block=256x256 Type=Float32, ColorInterp=Gray

    Computed Min/Max=0.000,2.482

  Minimum=0.000, Maximum=2.482, Mean=0.009, StdDev=0.045

  Metadata:

    STATISTICS_MINIMUM=0

    STATISTICS_MAXIMUM=2.4823353290558

    STATISTICS_MEAN=0.0089933774132436

    STATISTICS_STDDEV=0.044540742213422

The raster style is:

<?xml version="1.0" encoding="ISO-8859-1"?>

<StyledLayerDescriptor version="1.0.0"

    xsi:schemaLocation="http://www.opengis.net/sld
StyledLayerDescriptor.xsd"

    xmlns="http://www.opengis.net/sld&quot;

    xmlns:ogc="http://www.opengis.net/ogc&quot;

    xmlns:xlink="http://www.w3.org/1999/xlink&quot;

    xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance&quot;&gt;

    <!-- a named layer is the basic building block of an sld document -->

    <NamedLayer>

        <Name>precip_raster</Name>

        <UserStyle>

        <!-- they have names, titles and abstracts -->

            <Title>Precip Raster Colormap</Title>

            <Abstract>Colors for precipitation raster image</Abstract>

            <FeatureTypeStyle>

                    <Rule>

                                <RasterSymbolizer>

                                    <Opacity>1.0</Opacity>

                                    <ColorMap type="Intervals">

                                       <ColorMapEntry color="#FFFFFF"
quantity="0.00" opacity="0.0" />

                                       <ColorMapEntry color="#FFFFFF"
quantity="0.0099" opacity="0.0" />

                                       <ColorMapEntry color="#7FFF00"
quantity="0.099"/>

                                       <ColorMapEntry color="#00CD00"
quantity="0.249"/>

                                       <ColorMapEntry color="#008B00"
quantity="0.499"/>

                                       <ColorMapEntry color="#104E8B"
quantity="0.749"/>

                                       <ColorMapEntry color="#1E90FF"
quantity="0.999"/>

                                       <ColorMapEntry color="#00B2EE"
quantity="1.249"/>

                                       <ColorMapEntry color="#00EEEE"
quantity="1.499"/>

                                       <ColorMapEntry color="#8968CD"
quantity="1.749"/>

                                       <ColorMapEntry color="#912CEE"
quantity="1.999"/>

                                       <ColorMapEntry color="#8B008B"
quantity="2.499"/>

                                       <ColorMapEntry color="#8B0000"
quantity="2.999"/>

                                       <ColorMapEntry color="#CD0000"
quantity="3.999"/>

                                       <ColorMapEntry color="#EE4000"
quantity="4.999"/>

                                       <ColorMapEntry color="#FF7F00"
quantity="5.999"/>

                                       <ColorMapEntry color="#CD8500"
quantity="6.999"/>

                                       <ColorMapEntry color="#EEEE00"
quantity="7.999"/>

                                       <ColorMapEntry color="#FFFF00"
quantity="8.999"/>

                                    </ColorMap>

                                </RasterSymbolizer>

                </Rule>

            </FeatureTypeStyle>

        </UserStyle>

    </NamedLayer>

</StyledLayerDescriptor>

Thanks.

bw

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's

challenge

Build the coolest Linux based applications with Moblin SDK & win great
prizes
Grand prize is a trip for two to an Open Source event anywhere in the

world

http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great
prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great
prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users

-------------------------------------------------------------------------
This SF.Net email is sponsored by the Moblin Your Move Developer's challenge
Build the coolest Linux based applications with Moblin SDK & win great
prizes
Grand prize is a trip for two to an Open Source event anywhere in the world
http://moblin-contest.org/redirect.php?banner_id=100&url=/
_______________________________________________
Geoserver-users mailing list
Geoserver-users@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/geoserver-users