Re: [Geoserver-users] Simple WMS query driving me mad!

SourceForge Headquarters 1320 Columbia Street Suite 310 San Diego, CA 92101 +1 (858) 422-6466

Hi guys,

Sorry for the delay, I just got back to town last night.

You can get bursa wolf parameters from the epsg database. The databas=
e
contains a list of transformation methods from one CRS to another CRS
and the associated parameters. There can be more than one transformat=
ion
method for a datum shift between CRS's and more than one set of
parameters for a given transformation method. For example, there are
many sets of bursa-wolf parameters when converting from NAD 27 to WGS=
 84
(a north american example). The best parameters usually depends on yo=
ur
location.

Another source of parameters is appendix B of "Department of Defense
World Geodetic System 1984" (NIMA TR8350.2). I think this is the sour=
ce
of some bursa-wolf parameter information in the epsg database.

The other source of these parameters is your government's mapping
agency. They should have parameters calculated for many different
methods and may even have something more accurate than a simple
bursa-wolf. These parameters are also fed into the epsg database, but
you should always check that they are correct (if accuracy in importa=
nt
for your application).

Note that using parameters of [0,0,0] will not perform any datum shif=
t.
In my above north american example, it will just tell geotools than
NAD27 =3D=3D WGS84. This might be ok if you are mapping at a really s=
mall
scale (ie. a map of canada). The error is about 200 m, so this is not
suitable at larger scales.

If you have geotools jars, you can have geotools use the epsg databas=
e
to create WKT with bursa-wolf parameters. The following command uses
unreleased jars built with maven. It should be the same for the relea=
sed
2.2 jars, but the jar names might be a bit different.
epsg-hsql-2.2.x.jar is just the epsg database packaged as an HSQL
database.

rschulz@Linnaea:~/.maven/repository$ java -classpath
gt2/jars/referencing-2.2.x.jar:gt2/jars/epsg-hsql-2.2.x.jar:geoapi/ja=
rs/geoapi-2.0.jar:hsql/jars/hsqldb-1.8.0.1.jar:vecmath/jars/vecmath-1=
.3.jar:units/jars/units-0.01.jar org.geotools.referencing.factory.eps=
g.DefaultFactory EPSG:31467
EPSG version:    6.7
Database engine: HSQL Database Engine version 1.8.0
Database URL:    jdbc:hsqldb:file:/tmp/Geotools/Databases/HSQL/EPSG

_____________________________________________________________________=
__________

PROJCS["DHDN / Gauss-Kruger zone 3",
  GEOGCS["DHDN",
    DATUM["Deutsches Hauptdreiecksnetz",
      SPHEROID["Bessel 1841", 6377397.155, 299.1528128,
AUTHORITY["EPSG","7004"]],
      TOWGS84[598.1, 73.7, 418.2, 0.202, 0.045, -2.455,
1.3819742018555456],
      AUTHORITY["EPSG","6314"]],
    PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
    UNIT["degree", 0.017453292519943295],
    AXIS["Geodetic latitude", NORTH],
    AXIS["Geodetic longitude", EAST],
    AUTHORITY["EPSG","4314"]],
  PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]],
  PARAMETER["central_meridian", 9.0],
  PARAMETER["latitude_of_origin", 0.0],
  PARAMETER["scale_factor", 1.0],
  PARAMETER["false_easting", 3500000.0],
  PARAMETER["false_northing", 0.0],
  UNIT["m", 1.0],
  AXIS["Northing", NORTH],
  AXIS["Easting", EAST],
  AUTHORITY["EPSG","31467"]]

I really should get some of this stuff into the geotools crs tutorial=
.
Hope it helps.

Rueben

P.S. just for fun, here is the output for EPSG:27700.

EPSG version:    6.7
Database engine: HSQL Database Engine version 1.8.0
Database URL:    jdbc:hsqldb:file:/tmp/Geotools/Databases/HSQL/EPSG

_____________________________________________________________________=
__________

PROJCS["OSGB 1936 / British National Grid",
  GEOGCS["OSGB 1936",
    DATUM["OSGB 1936",
      SPHEROID["Airy 1830", 6377563.396, 299.3249646,
AUTHORITY["EPSG","7001"]],
      TOWGS84[446.448, -125.157, 542.06, 0.15, 0.247, 0.842,
-4.2261596151967575],
      AUTHORITY["EPSG","6277"]],
    PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
    UNIT["degree", 0.017453292519943295],
    AXIS["Geodetic latitude", NORTH],
    AXIS["Geodetic longitude", EAST],
    AUTHORITY["EPSG","4277"]],
  PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]],
  PARAMETER["central_meridian", -2.0],
  PARAMETER["latitude_of_origin", 49.0],
  PARAMETER["scale_factor", 0.9996012717],
  PARAMETER["false_easting", 400000.0],
  PARAMETER["false_northing", -100000.0],
  UNIT["m", 1.0],
  AXIS["Easting", EAST],
  AXIS["Northing", NORTH],
  AUTHORITY["EPSG","27700"]]

On Thu, 2005-15-12 at 12:36 -0800, Justin Deoliveira wrote:
> Hi Lena,
>=20
> Yes I was talking about the bursa wolf parameters issue. And yes yo=
u can=20
> fill in values for TOWGS84. Actually I am not sure how one can get =
an=20
> exact list of what these parameters should be. However I have notic=
ed=20
> that all zero's usually works in most cases so I would try that.
>=20
> I have cc'd Reuben, he is one of the geotools crs experts. He shoul=
d be=20
> able to give us a better answer as to how to obtain the values for =
these=20
> paramters.
>=20
> -Justin
>=20
> Lena Pahl wrote:
> > Hi Justin!
> >=20
> > I'm using GeoServer 1.3RC3. I don't know if I made clear that I n=
ever had that error  from the StreamingRenderer Rob was writing about=
, the WMS is working fine when I put in the WGS84 coords for the BBOX=
. So both of these requests work:
> >=20
> > http://localhost:8080/geoserver1.3RC3/wms?bbox=3D3539209,5803566,=
3541356,5806154&styles=3Dcool&Format=3Dimage/png&request=3DGetMap&lay=
ers=3Dlasy:flurstuecke&width=3D800&height=3D800&srs=3DEPSG:31467
> >=20
> > http://localhost:8080/geoserver1.3RC3/wms?bbox=3D9.57,52.36,9.61,=
52.39&styles=3Dcool&Format=3Dimage/png&request=3DGetMap&layers=3Dlasy=
:flurstuecke&width=3D550&height=3D250&srs=3DEPSG:4326
> >=20
> > I only had the error 'SRS transformation error - missing Bursa pa=
rameters' when creating my feature types. But I would be grateful to =
get this working with the wkt file. Because the way I "generated" the=
 geographic coordinates was a bit dodgy, I put the file in ArcGIS and=
 let the coords display in WGS84 :-)))
> >=20
> > So, stupid question: where do I get those Bursa params from? EPSG=
 database?
> > If I have those parameters all I need to do is to fill the zero v=
alues of this: 'TOWGS84[0,0,0,0,0,0,0]' in the *.prj file, right? And=
 then I go and create a new feature type and the bbox is generated au=
tomatically with the help of the *.prj file?
> >=20
> > Thanks for helping!
> >=20
> > Regards,
> > Lena
> >=20
> >=20
> >=20
> >=20
> > Justin Deoliveira <jde...@un...> schrieb am 15.12.05=
 19:53:20:
> >=20
> >>Hi Lena,
> >>
> >>Did you get a chance to try the wkt parameter I sent out? I have =
looked=20
> >>into making the transforms lenient, but i am now not sure that it=
 will=20
> >>solve the problem. However the fix will be in the next release so=
 you=20
> >>should be able to try it at least.
> >>
> >>-Justin
> >>
> >>Justin Deoliveira wrote:
> >>
> >>>Sorry, forgot to attach the files, rookie mistake :)
> >>>
> >>>Lena Pahl wrote:
> >>>
> >>>
> >>>>Hi Justin!
> >>>>How can I supply the necessary parameters? As for me, I had tha=
t=20
> >>>>'Bursa Wolf params' error message when I tried to create featur=
e types=20
> >>>>for shapefiles and for ArcSDE layers based in Germany (EPSG cod=
e is=20
> >>>>31467, Gauss-Krueger Coordinates)
> >>>>Thanks,
> >>>>Lena
> >>>>
> >>>>
> >>>>Justin Deoliveira <jde...@un...> schrieb am 13.12.=
05=20
> >>>>18:24:05:
> >>>>
> >>>>
> >>>>>Yes, Lena is correct here. Because your data is not in WGS84=
=20
> >>>>>geoserver needs to perform a transform. However in this case t=
he=20
> >>>>>neccessary information (Bursa Wolf Paramters) for your project=
ion is=20
> >>>>>not available which means geoserver cant perform the datum shi=
ft.
> >>>>>
> >>>>>What format is your data in? Depending on the format you may b=
e able=20
> >>>>>to suply the neccesary paramters to have the transformation wo=
rk.
> >>>>>
> >>>>>-Justin
> >>>>>
> >>>>>Lena Pahl wrote:
> >>>>>
> >>>>>
> >>>>>>Hi Rob!
> >>>>>>
> >>>>>>I cannot offer help on your WMS problem. But I had the same=
=20
> >>>>>>experience when I tried to create bounding boxes for my=20
> >>>>>>featuretypes. This is what I found out:
> >>>>>>GeoServer uses geographic coordinates (Lat/Lon) for the bound=
ing box=20
> >>>>>>of a layer. If your data is in WGS84 GeoServer is able to cre=
ate a=20
> >>>>>>bounding box automatically from the data itself. If your data=
 is not=20
> >>>>>>WGS84based you have to put in the min/max values for the boun=
ding=20
> >>>>>>box in geographic coords by hand. As far as I remember you do=
 not=20
> >>>>>>need the exact (transformation) values. Just an approximation=
. I=20
> >>>>>>think this is needed by the WMS but I'm not sure. Give it a t=
ry!=20
> >>>>>>Maybe that solves both of your problems :-)
> >>>>>>
> >>>>>>Cheers,
> >>>>>>Lena
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>><rob...@to...> schrieb am 13.12.05 13:25:11:
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>>Hi
> >>>>>>>
> >>>>>>>I've got some geodata I'd like to be able to view using WMS=
=20
> >>>>>>>provided by Geoserver.
> >>>>>>>
> >>>>>>>I can load the data as a WFS OK, but I cant produce WMS.  Th=
e=20
> >>>>>>>Geoserver console
> >>>>>>>shows the following lines:
> >>>>>>>
> >>>>>>>90330 [FINE] org.geotools.renderer.lite.StreamingRenderer - =
added=20
> >>>>>>>attribute GEOMETRY
> >>>>>>>90330 [SEVERE]
> >>>>>>>org.geotools.renderer.lite.StreamingRenderer$DefaultRenderLi=
stener=20
> >>>>>>>- Error
> >>>>>>>transforming bbox
> >>>>>>>90330 [SEVERE]
> >>>>>>>org.geotools.renderer.lite.StreamingRenderer$DefaultRenderLi=
stener=20
> >>>>>>>- Exception
> >>>>>>>rendering layer org.geotools.map.DefaultMapLayer@1cc1ef4
> >>>>>>>90330 [FINE] org.geotools.renderer.lite.StreamingRenderer - =
Style=20
> >>>>>>>cache hit
> >>>>>>>ratio: NaN , hits 0, requests 0
> >>>>>>>
> >>>>>>>Any clues on how I resolve this?
> >>>>>>>
> >>>>>>>
> >>>>>>>Now this may be unrelated, but I also have a problem creatin=
g=20
> >>>>>>>FeatureTypes for
> >>>>>>>ALL of my data.  I am UK based, and am predominantly using t=
he OS=20
> >>>>>>>BNG (27700). The Geoserver administration page never creates=
=20
> >>>>>>>Bounding Box Lat Longs reporting:
> >>>>>>>
> >>>>>>>Could not construct a SRS transformation - Bursa wolf parame=
ters=20
> >>>>>>>required.
> >>>>>>>
> >>>>>>>340099 [FINE]=20
> >>>>>>>org.geotools.data.jdbc.ConnectionPool$ConnectionListManager =
-
> >>>>>>>Connection closed - adding to available connections.
> >>>>>>>340099 [FINE] DefaultFactory - Unavailable factory: European=
=20
> >>>>>>>Petroleum Survey Group
> >>>>>>>340099 [FINE] DefaultFactory - Unavailable factory: European=
=20
> >>>>>>>Petroleum Survey Group
> >>>>>>>
> >>>>>>>Any clues in resolving this would similarly be gratefully re=
ceived!
> >>>>>>>
> >>>>>>>MTIA
> >>>>>>>
> >>>>>>>Rob
> >>>>>>>
> >>>>>>>---- Message sent via Totalise Webmail - http://www.totalise=
.co.uk/=20
> >>>>>>>MAX Broadband 2Mb now just =C2=A317.99 a month plus 2 months=
 free!
> >>>>>>>
> >>>>>>>
> >>>>>>>-------------------------------------------------------
> >>>>>>>This SF.net email is sponsored by: Splunk Inc. Do you grep t=
hrough=20
> >>>>>>>log files
> >>>>>>>for problems?  Stop!  Download the new AJAX search engine th=
at makes
> >>>>>>>searching your log files as easy as surfing the  web.  DOWNL=
OAD=20
> >>>>>>>SPLUNK!
> >>>>>>>http://ads.osdn.com/?ad_id=3D7637&alloc_id=3D16865&op=3Dclic=
k
> >>>>>>>_______________________________________________
> >>>>>>>Geoserver-users mailing list
> >>>>>>>Geo...@li...
> >>>>>>>https://lists.sourceforge.net/lists/listinfo/geoserver-users
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>_____________________________________________________________=
_________
> >>>>>>XXL-Speicher, PC-Virenschutz, Spartarife & mehr: Nur im WEB.D=
E=20
> >>>>>>Club!      =20
> >>>>>>Jetzt gratis testen! http://freemail.web.de/home/landingpad/?=
mc=3D021130
> >>>>>>
> >>>>>>
> >>>>>>
> >>>>>>-------------------------------------------------------
> >>>>>>This SF.net email is sponsored by: Splunk Inc. Do you grep th=
rough=20
> >>>>>>log files
> >>>>>>for problems?  Stop!  Download the new AJAX search engine tha=
t makes
> >>>>>>searching your log files as easy as surfing the  web.  DOWNLO=
AD SPLUNK!
> >>>>>>http://ads.osdn.com/?ad_idv37&alloc_id=16865&op=3Dclick
> >>>>>>_______________________________________________
> >>>>>>Geoserver-users mailing list
> >>>>>>Geo...@li...
> >>>>>>https://lists.sourceforge.net/lists/listinfo/geoserver-users
> >>>>>>
> >>>>>
> >>>>>
> >>>>>--=20
> >>>>>Justin Deoliveira
> >>>>>The Open Planning Project
> >>>>>http://topp.openplans.org
> >>>>>
> >>>>>
> >>>>>-------------------------------------------------------
> >>>>>This SF.net email is sponsored by: Splunk Inc. Do you grep thr=
ough=20
> >>>>>log files
> >>>>>for problems?  Stop!  Download the new AJAX search engine that=
 makes
> >>>>>searching your log files as easy as surfing the  web.  DOWNLOA=
D SPLUNK!
> >>>>>http://ads.osdn.com/?ad_id=3D7637&alloc_id=3D16865&op=3Dclick
> >>>>>_______________________________________________
> >>>>>Geoserver-users mailing list
> >>>>>Geo...@li...
> >>>>>https://lists.sourceforge.net/lists/listinfo/geoserver-users
> >>>>
> >>>>
> >>>>
> >>>>
> >>>>______________________________________________________________
> >>>>Verschicken Sie romantische, coole und witzige Bilder per SMS!
> >>>>Jetzt bei WEB.DE FreeMail: http://f.web.de/?mc=3D021193
> >>>>
> >>>>
> >>>
> >>>
> >>>
> >>>----------------------------------------------------------------=
--------
> >>>
> >>>PROJCS["DHDN / Gauss-Kruger zone 3",GEOGCS["DHDN",DATUM["Deutsch=
es_Hauptdreiecksnetz",SPHEROID["Bessel 1841",6377397.155,299.1528128,=
AUTHORITY["EPSG","7004"]],AUTHORITY["EPSG","6314"]],PRIMEM["Greenwich=
",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHO=
RITY["EPSG","9122"]],AUTHORITY["EPSG","4314"]],PROJECTION["Transverse=
_Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meri=
dian",9],PARAMETER["scale_factor",1],PARAMETER["false_easting",350000=
0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","900=
1"]],AUTHORITY["EPSG","31467"]]
> >>>
> >>>
> >>>----------------------------------------------------------------=
--------
> >>>
> >>>PROJCS["DHDN / Gauss-Kruger zone 3",GEOGCS["DHDN",DATUM["Deutsch=
es_Hauptdreiecksnetz",SPHEROID["Bessel 1841",6377397.155,299.1528128,=
AUTHORITY["EPSG","7004"]],AUTHORITY["EPSG","6314"]],TOWGS84[0,0,0,0,0=
,0,0],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.=
01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4314"]]=
,PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",0],=
PARAMETER["central_meridian",9],PARAMETER["scale_factor",1],PARAMETER=
["false_easting",3500000],PARAMETER["false_northing",0],UNIT["metre",=
1,AUTHORITY["EPSG","9001"]],AUTHORITY["EPSG","31467"]]
> >>
> >>
> >>--=20
> >>Justin Deoliveira
> >>The Open Planning Project
> >>http://topp.openplans.org
> >>
> >>
> >>-------------------------------------------------------
> >>This SF.net email is sponsored by: Splunk Inc. Do you grep throug=
h log files
> >>for problems?  Stop!  Download the new AJAX search engine that ma=
kes
> >>searching your log files as easy as surfing the  web.  DOWNLOAD S=
PLUNK!
> >>http://ads.osdn.com/?ad_id=3D7637&alloc_id=3D16865&op=3Dclick
> >>_______________________________________________
> >>Geoserver-users mailing list
> >>Geo...@li...
> >>https://lists.sourceforge.net/lists/listinfo/geoserver-users
> >=20
> >=20
> >=20
> > _________________________________________________________________=
_____
> > XXL-Speicher, PC-Virenschutz, Spartarife & mehr: Nur im WEB.DE Cl=
ub!=09=09
> > Jetzt gratis testen! http://freemail.web.de/home/landingpad/?mc=
=3D021130
> >=20
> >=20
> >=20
> > -------------------------------------------------------
> > This SF.net email is sponsored by: Splunk Inc. Do you grep throug=
h log files
> > for problems?  Stop!  Download the new AJAX search engine that ma=
kes
> > searching your log files as easy as surfing the  web.  DOWNLOAD S=
PLUNK!
> > http://ads.osdn.com/?ad_idv37&alloc_id=16865&op=3Dclick
> > _______________________________________________
> > Geoserver-users mailing list
> > Geo...@li...
> > https://lists.sourceforge.net/lists/listinfo/geoserver-users
> >=20
>=20
>=20