WCS and Reprojection with TileItem

Jérôme Martin gooogle at FREE.FR
Wed Jan 31 20:05:32 EST 2007


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 Modifier/Supprimer le message 
<http://www.forumsig.org/editpost.php?do=editpost&p=93083>
-------------- next part --------------
Skipped content of type multipart/related


More information about the mapserver-users mailing list