[Geoserver-users] Reprojected point disappears past -67 south [Sec=Unclassified]

-----Original Message-----
From: Andrea Aime [mailto:aaime@anonymised.com] >
Miles Jordan ha scritto:
> Hi all,
>
> I've come across something strange with geoserver 1.5.1 Latest war.
>
> I am requesting a map using the EPSG:3031 antarctic polar
> stereographic projection (0,0 is the south pole), showing
the whole of
> antarctica. I am using a shapefile of the antarctic
continent, which
> is also in EPSG:3031, to view the entire antarctic
continent on a map.
> This all works fine.
>
> I am using an InlineFeature to display a yellow dot on the map. The
> location of the dot is specified using EPSG:4326 lat/lon decimal
> coordinates, and reprojected on the fly (which is a great
feature, by
> the way). The dot displays on any longitude, but disappears
when it is
> placed at a latitude of -68 or higher.

Ah ha! Yeah, I already faced this issue one month ago,
shortly after the 1.5.1 release. The renderer is doing a bbox
query against the data to reduce the amount of data that
needs to be loaded, and the bbox reprojection fails to take
into account the middle of the box, but
considers only the perimeter, which lies more or less at the
latitude you mentioned for Antarctica.

We fixed it, so you may want to try out a nightly of the 1.5.2 series
here:
http://geo.openplans.org/nightly/1.5.x/

If the problem still persists, it means not all code paths
have been taken care of, and I'll need some sample data and
the full request.

Cheers
Andrea

Hi Andrea,

I'm afraid the problem does still exist. I tried the 071007 build.
Here's the request:
http://YOUR-GEOSERVER-HOST/geoserver/wms?bgcolor=0xf2AAAFF&bbox=-2661827
8078423585,-2489793.885787334,2750918.460493497,2321607.9999959785&SLD=
http://milesjordan.com/sld/geos-test.sld&Format=image/png&request=GetMap
&width=600&height=600&SRS=EPSG:3031

The above should show a blue map with yellow dots in opposite corners,
one at -45,-67 and the other at 135,-80, however only the one at -45,-67
appears. I took out the reference to the antarctic coastline layer
because I can't redistribute the data (although it it available from the
Antarctic Digital Database) - it's not needed anyway, as long as the SRS
in the request is set to EPSG:3031.

I noticed 3031 snuck it's way out of the EPSG codes that come with
geoserver releases again, but assume you already have it defined as you
seem to have been doing some work with it!

I see no difference at all with the nightly build - is it possible that
the changes aren't in it yet? Thanks for your help.

Regards,

Miles

___________________________________________________________________________

    Australian Antarctic Division - Commonwealth of Australia
IMPORTANT: This transmission is intended for the addressee only. If you are not the
intended recipient, you are notified that use or dissemination of this communication is
strictly prohibited by Commonwealth law. If you have received this transmission in error,
please notify the sender immediately by e-mail or by telephoning +61 3 6232 3209 and
DELETE the message.
        Visit our web site at http://www.antarctica.gov.au/
___________________________________________________________________________

Miles Jordan ha scritto:

I noticed 3031 snuck it's way out of the EPSG codes that come with
geoserver releases again, but assume you already have it defined as you
seem to have been doing some work with it!

Yes, Geoserver 1.5.x cannot support it because of the peculiar definition of an axis direction, only 1.6.x can work with it correctly.
In fact I'm curious about how you're making your tests... on my box
I cannot have 3031 parsed at all (even if I put the wkt in the user
projections file).

Cheers
Andrea

Miles Jordan ha scritto:

I'm afraid the problem does still exist. I tried the 071007 build.
Here's the request:
http://YOUR-GEOSERVER-HOST/geoserver/wms?bgcolor=0xf2AAAFF&bbox=-2661827
8078423585,-2489793.885787334,2750918.460493497,2321607.9999959785&SLD=
http://milesjordan.com/sld/geos-test.sld&Format=image/png&request=GetMap
&width=600&height=600&SRS=EPSG:3031

Soo... I've found quite a few issues, and I think I fixed them all,
but only on trunk, on 1.5.x I could to make 3031 work (if you have any
magic recipe of your own, please share, it's probably just an issue
of setting up the proper wkt definition :)).

One thing, the srsName you provided in the SLD is invalid, the
correct form is srsName="srsName="http://www.opengis.net/gml/srs/epsg.xml#4326"

If you try out tomorrow's trunk nightly build you should get the
same image as the one attached with the following request (fix
the sld url to point to your own):

http://localhost:8080/geoserver/wms?bgcolor=0xf2AAAFF&bbox=-2761827.8078423585,-2589793.885787334,2850918.460493497,2421607.9999959785&SLD=http://localhost:8080/geoserver/www/test.sld&Format=image/png&request=GetMap&width=600&height=600&SRS=EPSG:3031

Cheers
Andrea

wms.png