[mapserver-users] TinyOWS: Axis order fix wanted

Even Rouault even.rouault at mines-paris.org
Sat Jun 28 14:38:47 PDT 2014


Hi,

The following pull request should fix Jukka's issue :

https://github.com/mapserver/tinyows/pull/65 

This was indeed quite severe and affected very popular projected SRS like UTM 
WGS84 SRS (EPSG:326XX), that are similar to EPSG:3067 as far as their WKT AXIS 
definition is concerned, i.e. ending with :

AXIS["Easting",EAST],AXIS["Northing",NORTH]]

I guess this went unnoticed since the test	/wfs_110 tests only tests EPSG:4326, 
which was unaffected.
And test/wfs_100 tests EPSG:32615 which was affected, but in WFS 1.0, there's 
no axis swapping...

Even

Le jeudi 19 juin 2014 07:05:19, Rahkonen Jukka (Tike) a écrit :
> Hi,
> 
> 
> 
> EPSG:2393 is Northing-Easting
> EPSG:3047 is Northing-Easting
> EPSG:3067 is Easting-Northing
> 
> How TinyOWS is planned to work it that axis should be flipped, if needed,
> when srsName is given as urn. Short EPSG:xxxx should return always E-N
> just like WFS 1.0.0 always.
> 
> The old tinyows is still available at hip.latuviitta.org/cgi-bin/tinyows
> and you can compare it with this release version installed from Ubuntugis
> 
> 2393 is OK
> http://XXXXXXX/cgi-bin/tinyows?service=wfs&version=1.1.0&request=getf
> eature&typename=lv:municipalities&maxfeatures=1&srsname=EPSG:2393 Comes as
> E-N
> http://XXXXXXX/cgi-bin/tinyows?service=wfs&version=1.1.0&request=getf
> eature&typename=lv:municipalities&maxfeatures=1&srsname=urn:ogc:def:crs:EPS
> G::2393 Comes as N-E
> 
> 3047 is wrong
> http://XXXXXXX/cgi-bin/tinyows?service=wfs&version=1.1.0&request=getf
> eature&typename=lv:municipalities&maxfeatures=1&srsname=EPSG:3047 Comes as
> N-E, should be E-N
> http://XXXXXXX/cgi-bin/tinyows?service=wfs&version=1.1.0&request=getf
> eature&typename=lv:municipalities&maxfeatures=1&srsname=urn:ogc:def:crs:EPS
> G::3047 Comes as  E-N, should be N-E
> 
> 3067 is wrong
> http://XXXXXXX/cgi-bin/tinyows?service=wfs&version=1.1.0&request=getf
> eature&typename=lv:municipalities&maxfeatures=1&srsname=EPSG:3067 Comes as
> E-N, which is correct
> http://XXXXXXX/cgi-bin/tinyows?service=wfs&version=1.1.0&request=getf
> eature&typename=lv:municipalities&maxfeatures=1&srsname=urn:ogc:def:crs:EPS
> G::3067 Comes as N-E but should come as E-N
> This is my default projection so this should also come as E-N but it does
> not
> http://XXXXXXX/cgi-bin/tinyows?service=wfs&version=1.1.0&request=get
> feature&typename=lv:municipalities&maxfeatures=1
> 
> Here comes srtext from the database for those three and everything is OK.
> 
> gis93=> select srid,srtext from spatial_ref_sys where srid in
> (2393,3047,3067);
> UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Transverse_Mercator"]
> ,PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",27],PARAMET
> ER["scale_factor",1],PARAMETER["false_easting",3500000],PARAMETER["false_no
> rthing",0],AUTHORITY["EPSG","2393"],AXIS["X",NORTH],AXIS["Y",EAST]] 3047 |
> PROJCS["ETRS89 /
> ETRS-TM35",GEOGCS["ETRS89",DATUM["European_Terrestrial_Reference_System_19
> 89",SPHEROID["GRS
> 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0,0,0,0,0,0,
> 0],AUTHORITY["EPSG","6258"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]]
> ,UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG
> ","4258"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Transverse_
> Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",2
> 7],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAM
> ETER["false_northing",0],AUTHORITY["EPSG","3047"],AXIS["Northing",NORTH],AX
> IS["Easting",EAST]] 3067 | PROJCS["ETRS89 /
> ETRS-TM35FIN",GEOGCS["ETRS89",DATUM["European_Terrestrial_Reference_System
> _1989",SPHEROID["GRS
> 1980",6378137,298.257222101,AUTHORITY["EPSG","7019"]],TOWGS84[0,0,0,0,0,0,
> 0],AUTHORITY["EPSG","6258"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]]
> ,UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG
> ","4258"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Transverse_
> Mercator"],PARAMETER["latitude_of_origin",0],PARAMETER["central_meridian",2
> 7],PARAMETER["scale_factor",0.9996],PARAMETER["false_easting",500000],PARAM
> ETER["false_northing",0],AUTHORITY["EPSG","3067"],AXIS["Easting",EAST],AXIS
> ["Northing",NORTH]]
> 

-- 
Geospatial professional services
http://even.rouault.free.fr/services.html



More information about the MapServer-users mailing list