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

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.

If I request the map using EPSG:4326, the dot will display at any
coordinate within the bbox, and the shapefile is reprojected correctly
so it all looks fine. Antarctica doesn't look like that, though :slight_smile:

Here's what it looks like with the dot at -110,-66
(http://milesjordan.com/images/work/getMap3031_working.png)
And at -110,-80
(http://milesjordan.com/images/work/getMap3031_broken.png)
And at -110,-80 using EPSG:4326
(http://milesjordan.com/images/work/getMap4326.png)

There doesn't seem to be a problem with the SLD it is defined in because
it shouldn't work for one but not the other, but it's accessible at
http://milesjordan.com/sld/inline.sld.

There is a log extract at http://milesjordan.com/sld/geoserver.log
although everything seems fine, except for the odd "WARNING: No CRS set
on inline features default geometry. Assuming the requestor has their
inlinefeatures in the boundingbox CRS." (see line 16). This is strange
because a) the CRS is set for the inline Point and b) if we were to
assume that it was using the same CRS as specified by the bounding box,
it wouldn't even display properly above 67 degrees south.

I've spent most of my day trying to figure this out. Does anyone have
any ideas?

Cheers,

Miles Jordan

Applications Developer
Data Centre
The Australian Antarctic Division
203 Channel Hwy, Kingston TAS 7053

[p] +61 (0) 3 6232 3486
[m] +61 (0) 424 879 668
[e] miles.jordan@anonymised.com
[w] http://www.aad.gov.au
___________________________________________________________________________

    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:

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