WCS and Reprojection with TileItem

[Norman Barker]

Hi Jerome,

 

I know Jeremy Morley doesn't watch these lists too much, but his ICEDS service

 

http://iceds.ge.ucl.ac.uk/cgi-bin/icedswcs?SERVICE=WCS&REQUEST=DescribeCoverage&Coverage=srtm

 

offers srtm data in different projections.

 

His website also gives full config details for MapServer configuration, a subset of which is on the MapServer website under WCS use cases.

 

Thanks,

 

Norman

 

________________________________

From: UMN MapServer Users List [mailto:MAPSERVER-USERS at LISTS.UMN.EDU] On Behalf Of Jérôme Martin
Sent: 01 February 2007 01:06
To: MAPSERVER-USERS at LISTS.UMN.EDU
Subject: [UMN_MAPSERVER-USERS] WCS and Reprojection with TileItem

 

Hi,

I'm using mapserv like WMS and now, I would like use the WCS service to get DEM, from DTED files. 

To do this, i used gdaltindex software to create a shapefile to have only one entry on my mapfile with the "TILEITEM" keyword. 

I precise my dted file are LatLon, in EPSG:4326. 

So when i did a request Getmap from the server, i had the result i wanted for, in lat-lon. 
But my problem was when i tried to send the request, but with CRS=epsg:32631 (UTM Zone N31) by example (BBOX in 32631 too), nothing work and the server send me a file with the good boundingbox envelloppe, but with no data, only 0. 

When I tried to do a resquest with CRS=epsg:4326 (BBOX in 4326)  and RESPONSE_CRS=epsg:32631, i have the same empty 0 file, and the server doesn't warn me any error. Just a file with 0m data elevation. 

I'm really surprised, because the OGC documents about the WCS Service tell the response_crs can be used to have a reprojection, and the how-to user case from the mapserver website tell this option can be used. 

So when i try to put only one file with the DATA keyword, i was really surprised to see a request in another crs gave me the good reprojection and the right results, so the reprojection works on only one file, but not in shapefile tileitem. So Did I do something wrong, or this is it a bug from mapserv ?

Moreover, i am already using mapserver with the WMS service, and i have only shapefile tileitem entries on my wms mapfile : the reprojection works perfectly ! (the shapefile are created from the gdaltindex tool). 

I can't put each file on my mapfile, because i thousands files and this will be affect the performances of mapserv (moreover, i think it exist a limitation of the maximum number of layers in a mapfile). 

So do you would have any idea about this problem, or do you know an easy exemple of a wcs service server who delivers DEM data in differents reprojections ?

thanks you for the help,

sincerly, 

Jérôme Martin




J'ai donc configuré un server wcs, et j'ai également crée un fichier shp via gdaltindex afin de n'avoir qu'une seule entrée dans le mapfile regroupant plusieurs fichiers DTED. 

Je précise que les fichiers sont en latlon EPSG:4326

En faisant une requête avec CRS=epsg:4326, j'ai donc le résultat attendu, soit un fichier tiff (ce que je voulais) en latlon. 
Par contre, lorsque je fais une requête avec CRS=epsg:32631 (UTM zone N31) par exemple, là, tout se complique puisque la réponse que me renvoye le server semble bonne au niveau des dimensions, mais je n'ai uniquement que des 0, les données d'élévation disparaissent. 

Je précise que j'ai déjà un server WMS utilisant exclusivement des tileindex (via gdaltindex) et que la reprojection se passe parfaitement. Mais en WCS, rien ne marche plus pour la reprojection. 

Par contre, si je prend un fichier seul (fichier.dt1) que je rajoute dans mon mapfile wcs, alors dans ce cas là, la reprojection marche bien et ne pose pas de problème. Je ne peux cependant pas résoudre mon problème comme cela car je ne peux inclure des milliers de fichiers dans mon mapfile, et je ne peux pas non plus modifier les fichiers dted en entrée (d'autant plus que je ne veux pas distinguer séparement chaque fichier dted dnas mes requêtes). 

Comment puis-je résoudre mon problème ? quelqu'un a-t-il déjà été confronté à ce type de problème sous WCS ? 

cordialement,

Jérôme   <http://www.forumsig.org/editpost.php?do=editpost&p=93083> 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20070201/c40e49a4/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image001.gif
Type: image/gif
Size: 1834 bytes
Desc: image001.gif
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20070201/c40e49a4/attachment.gif>