[Geoserver-users] WCS: "Unable to get the resolution" during a getMap request

GeoServer GeoServer User
1

Hi list:

Sometimes, when I extend the max y coordinate in BBOX in a GetMap
request, the coverages fail to show, and I get the following console
output:

2826472 [WARNING] org.geotools.renderer.lite.StreamingRenderer -
Transformation doesn't convergence.
2826474 [WARNING] org.geotools.renderer.lite.StreamingRenderer -
Unable to get the resolution
2826658 [INFO] org.vfny.geoserver.servlets.AbstractService - Service handled

Vector layers, however, continue to display fine. I run into this
situation when, for example I want to display the northern most tip of
Canada.

If I diminish the Max y value, the coverage(s) show up on the map again.

Here is the image where the raster shows well:
http://smoke-fire.us/img/normal.png

Here is the image where the raster doesn't show in the desired extend:
http://smoke-fire.us/img/problem.png

The request for normal.png:
http://localhost:8080/geoserver/wms?bbox=-4578324,-2500333,3189008,4311929&styles=raster,polygon&Format=image/png&request=GetMap&layers=test,countries&width=800&height=600&srs=EPSG:41111

The request for problem.png (only the max y value for the BBOX is changed):
http://localhost:8080/geoserver/wms?bbox=-4578324,-2500333,3189008,4311929&styles=raster,polygon&Format=image/png&request=GetMap&layers=test,countries&width=800&height=600&srs=EPSG:41111

EPSG:41111 is an Albers projection I had to define. Here is how it's
listed in srsHelp.do:

  41111 PROJCS["Albers Equal area", 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 latitude", NORTH],
AXIS["Geodetic longitude", EAST], AUTHORITY["EPSG","4326"]],
PROJECTION["Albers_Conic_Equal_Area"], PARAMETER["central_meridian",
-96.0], PARAMETER["latitude_of_origin", 37.5],
PARAMETER["standard_parallel_1", 29.833333333333336],
PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0],
PARAMETER["standard_parallel_2", 45.833333333333336], UNIT["m", 1.0],
AXIS["Easting", EAST], AXIS["Northing", NORTH],
AUTHORITY["EPSG","41111"]]

I would like to show the coverages while the entire vector layer is
visible also.

Thanks for any insight,
Alex

2

Alexander Petkov ha scritto:

Hi list:

Sometimes, when I extend the max y coordinate in BBOX in a GetMap
request, the coverages fail to show, and I get the following console
output:

2826472 [WARNING] org.geotools.renderer.lite.StreamingRenderer -
Transformation doesn't convergence.
2826474 [WARNING] org.geotools.renderer.lite.StreamingRenderer -
Unable to get the resolution
2826658 [INFO] org.vfny.geoserver.servlets.AbstractService - Service handled

Vector layers, however, continue to display fine. I run into this
situation when, for example I want to display the northern most tip of
Canada.

Hmm... this usually happens if anything on the map goes beyond the area of definition for the current
projection... yet it seems the visible part of the coverage is well within the area...
Can you provide me with a sample data set?

Cheers
Andrea

3

Alexander Petkov ha scritto:

Hi list:

EPSG:41111 is an Albers projection I had to define. Here is how it's
listed in srsHelp.do:

  41111 PROJCS["Albers Equal area", 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 latitude", NORTH],
AXIS["Geodetic longitude", EAST], AUTHORITY["EPSG","4326"]],
PROJECTION["Albers_Conic_Equal_Area"], PARAMETER["central_meridian",
-96.0], PARAMETER["latitude_of_origin", 37.5],
PARAMETER["standard_parallel_1", 29.833333333333336],
PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0],
PARAMETER["standard_parallel_2", 45.833333333333336], UNIT["m", 1.0],
AXIS["Easting", EAST], AXIS["Northing", NORTH],
AUTHORITY["EPSG","41111"]]

I opened an issue so that I wont' forget fixing this:
http://jira.codehaus.org/browse/GEOS-765

Now, it occurred to me now that there's a logical problem here that
won't be easy to solve.
Projections defined in WKT do not usually have a validity envelope,
and this may be the source of the problem, because all resolution/scale
computations are performed against the coverge bbox, intersected
with the validity area if available...

So my wild guess is that your coverage envelope goes beyond the
validity area of Albers equal area, and poor gt2 does not know what
to do about it.

Wondering if there's any way to specify the validity area in wkt...
can you check?

Cheers
Andrea Aime

4

Wondering if there's any way to specify the validity area in wkt...
can you check?

Hello Andrea:

I didn't see a way to specify a validity area in WKT, as defined by OGC.

I will define my projection in the HSQL database and see if I still
have the same issue.
Alex

5

On 10/24/06, Alexander Petkov <greenkov@anonymised.com> wrote:

> Wondering if there's any way to specify the validity area in wkt...
> can you check?

Hello Andrea:

I didn't see a way to specify a validity area in WKT, as defined by OGC.

I will define my projection in the HSQL database and see if I still
have the same issue.
Alex

That didn't help :frowning:
The same problem persists...

Alex