Hi Guillaume,

You seem to have discovered a Mapserver bug.  It would
be very nice if you could take the trouble to report it
(so that it can be addressed / fixed).

http://trac.osgeo.org/mapserver/newticket
(everyone can create tickets, you just need to log in -
if you do not have an account, it is easy to create one)

Håvard Tveite

On 11/8/2011 2:58 PM, Guillaume Loiselle wrote:
Hi,

Wow! Adding a projection object on layer level fixed the problem. I was sure 
the map level projection object was used when a layer was lacking it. Now the 
result is the same than that of the WMS GeoTiff output (except for the 
background color but that's not an issue).

Thanks a lot for your help!
Guillaume Loiselle


2011/11/7 Rahkonen Jukka <[email protected] 
<mailto:[email protected]>>

    __
    Hi,
    Mapfile looks good to me. I fear I have only two suggestions left:
    - try adding projection also to the layer even it should use the map level 
projection by default
    - forget HDF5 for a moment and use some geotiff instead to check if there 
are some oddities due to source data format
    -Jukka Rahkonen-

        
-------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------
-----------
        *Lähettäjä:* [email protected] 
<mailto:[email protected]> 
[mailto:[email protected] 
<mailto:[email protected]>] *Puolesta *Guillaume Loiselle
        *Lähetetty:* 7. marraskuuta 2011 23:07
        *Vastaanottaja:* [email protected] 
<mailto:[email protected]>
        *Aihe:* Re: [mapserver-users] Re: WCS Geotiff output issue

        Hi,

        Both mapfile and datasource are available in the first email I sent.

        
http://lists.osgeo.org/pipermail/mapserver-users/2011-November/070670.html

        Filenames appear to be a little bit messed up but the data itself is 
fine. I don't know why it went so wrong.

        Mapfile : 
http://lists.osgeo.org/pipermail/mapserver-users/attachments/20111104/78e5fa24/2011102811-0001.obj

        Data (HDF5 format) : 
http://lists.osgeo.org/pipermail/mapserver-users/attachments/20111104/78e5fa24/data-0001.obj

        I intend to stick with 1.0.0 indeed. At least until it works... I will 
look up 1.1.0's details in due time.




        2011/11/7 Rahkonen Jukka <[email protected] 
<mailto:[email protected]>>

            Pity that is was not just version mismatch.  Could it be possible 
to see your mapfile and info about the datasource?
            Stay with version 1.0.0 while testing. With version 1.1.0 you must 
give EPSG:4326 coordinates in Northing-Easting order which can make things more 
complicated.

            -Jukka Rahkonen-
            ________________________________
            Guillaume Loiselle wrote:

            Hi,

            I'm afraid it isn't. I didn't notice the request version I sent was a 
"mixed" one but I tried Coverage on 1.0.0 and Identifier on 1.1.0 as well and 
the result was basically the same.
            Here is the 1.1.0 version of the error :


            <?xml version="1.0" encoding="ISO-8859-1"?>
            <ows:ExceptionReport xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"; 
xmlns:ows="http://www.opengis.net/ows/1.1"; version="1.1.0" xml:lang="en-US" 
xsi:schemaLocation="http://www.opengis.net/ows/1.1 http://schemas.opengis.net/ows/1.1.0/owsExceptionReport.xsd";>
            <ows:Exception exceptionCode="bbox" locator="NoApplicableCode">
            <ows:ExceptionText>
            msWCSGetCoverage(): WCS server error. Requested BBOX 
(-178.804111261525,2802.84894277903,179.578125552223,2907.18853507175) is outside 
requested coverage BBOX (-70,45,-55,52)</ows:ExceptionText>
            </ows:Exception>
            </ows:ExceptionReport>


            I of course spent a lot of time on the WCS server documentation 
page. I saw no reference to such a problem, though. My request,

            http://localhost/cgi-bin/mapserv.exe?
            map=D:/Guillaume/ProjetCourantsGlaces/WCSproblem/2011102811.map&
            SERVICE=WCS&
            VERSION=1.0.0&
            REQUEST=GetCoverage&
            COVERAGE=thelayer&
            CRS=EPSG:4326&
            BBOX=-67.5,48.922499263758255,-64.6875,50.736455137010665&
            WIDTH=256&
            HEIGHT=256&
            FORMAT=GTiff

            is formatted the same way the example is :

            http://my.host.com/cgi-bin/mapserv?
            map=mywcs.map&
            SERVICE=wcs&
            VERSION=1.0.0&
            REQUEST=GetCoverage&
            coverage=bathymetry&
            CRS=EPSG:42304&
            BBOX=-2200000,-712631,3072800,3840000&
            WIDTH=3199&
            HEIGHT=2833&
            FORMAT=GTiff


            My mapfile contains every required attribute.

            Thank you.
            Guillaume Loiselle


            2011/11/7 Rahkonen Jukka <[email protected] 
<mailto:[email protected]><mailto:[email protected] 
<mailto:[email protected]>>>
            Hi,

            See examples in http://www.mapserver.org/ogc/wcs_server.html and try to avoid mixing 
versions 1.0.0 and 1.1.0. Let's hope using "coverage" instead of "identifier" 
will be enough.

            -Jukka Rahkonen-
            ________________________________
            Lähettäjä: [email protected] 
<mailto:[email protected]><mailto:[email protected] 
<mailto:[email protected]>> [[email protected] 
<mailto:[email protected]><mailto:[email protected] 
<mailto:[email protected]>>] käyttäjän Guillaume Loiselle [[email protected] 
<mailto:[email protected]><mailto:[email protected] <mailto:[email protected]>>] 
puolesta
            Lähetetty: 7. marraskuuta 2011 21:21
            Vastaanottaja: [email protected] 
<mailto:[email protected]><mailto:[email protected] 
<mailto:[email protected]>>
            Aihe: [mapserver-users] Re: WCS Geotiff output issue

            Hi,

            For some reason, the text part of the email didn't make it. It can 
be read below.
            Attachments can still be found in the first message :
            
http://lists.osgeo.org/pipermail/mapserver-users/2011-November/070670.html

            Thanks.
            Guillaume Loiselle


            2011/11/4 Guillaume Loiselle <[email protected] 
<mailto:[email protected]><mailto:[email protected] 
<mailto:[email protected]>><mailto:[email protected] 
<mailto:[email protected]><mailto:[email protected] 
<mailto:[email protected]>>>>
            Hi,

            I installed MapServer 6.0.1 on Windows XP from MS4W. I cannot give you access 
to my MapServer instance but I can show you what "-v" gives me :

            MapServer version 6.0.1 OUTPUT=GIF OUTPUT=PNG OUTPUT=JPEG 
OUTPUT=KML SUPPORTS=PR
            OJ SUPPORTS=AGG SUPPORTS=CAIRO SUPPORTS=FREETYPE SUPPORTS=ICONV 
SUPPORTS=FRIBIDI
              SUPPORTS=WMS_SERVER SUPPORTS=WMS_CLIENT SUPPORTS=WFS_SERVER 
SUPPORTS=WFS_CLIENT
              SUPPORTS=WCS_SERVER SUPPORTS=SOS_SERVER SUPPORTS=FASTCGI 
SUPPORTS=THREADS SUPPO
            RTS=GEOS INPUT=JPEG INPUT=POSTGIS INPUT=OGR INPUT=GDAL 
INPUT=SHAPEFILE

            I am able to generate a GeoTiff file through the use of a WMS 
service (see attachment : wmsoutput.tif). The request is :

            
http://localhost/cgi-bin/mapserv.exe?SERVICE=WMS&map=D:/Guillaume/ProjetCourantsGlaces/WCSproblem/2011102811.map&LAYERS=thelayer&VERSION=1.1.1&REQUEST=GetMap&bbox=-67.5,48.922499263758255,-64.6875,50.736455137010665&SRS=EPSG:4326&WIDTH=256&HEIGHT=256&FORMAT=GTiff
 
<http://localhost/cgi-bin/mapserv.exe?SERVICE=WMS&map=D:/Guillaume/ProjetCourantsGlaces/WCSproblem/2011102811.map&LAYERS=thelayer&VERSION=1.1.1&REQUEST=GetMap&bbox=-67.5,48.922499263758255,-64.6875,50.736455137010665&SRS=EPSG:4326&WIDTH=256&HEIGHT=256&FORMAT=GTiff>

            However, I cannot do so with a WCS service. I first tried with :

            
http://localhost/cgi-bin/mapserv.exe?SERVICE=WCS&map=D:/Guillaume/ProjetCourantsGlaces/WCSproblem/2011102811.map&IDENTIFIER=thelayer&VERSION=1.0.0&REQUEST=GetCoverage&BBOX=-67.5,48.922499263758255,-64.6875,50.736455137010665&CRS=EPSG:4326&WIDTH=256&HEIGHT=256&FORMAT=GTiff
 
<http://localhost/cgi-bin/mapserv.exe?SERVICE=WCS&map=D:/Guillaume/ProjetCourantsGlaces/WCSproblem/2011102811.map&IDENTIFIER=thelayer&VERSION=1.0.0&REQUEST=GetCoverage&BBOX=-67.5,48.922499263758255,-64.6875,50.736455137010665&CRS=EPSG:4326&WIDTH=256&HEIGHT=256&FORMAT=GTiff>

            , to which I can add :

            &RESPONSE_CRS=EPSG:4326

            for the same result (since 4326 is my native projection). Said 
result is not a tiff file, but an error :

            <?xml version='1.0' encoding="ISO-8859-1" ?>
            <ServiceExceptionReport version="1.2.0"
            xmlns="http://www.opengis.net/ogc"; 
xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"; 
xsi:schemaLocation="http://www.opengis.net/ogc 
http://schemas.opengis.net/wcs/1.0.0/OGC-exception.xsd";>
            <ServiceException code="NoApplicableCode" 
locator="bbox">msWCSGetCoverage(): WCS server error. Requested BBOX 
(-178.835708198757,2803.05273104522,179.552848002437,2906.98474680556) is outside requested coverage BBOX 
(-70,45,-55,52)
            </ServiceException>
            </ServiceExceptionReport>

            This error doesn't make sense to me. Which projection is it using 
to get these silly values? Anyway, the BBOX in the request is definitely 
included in the extent... Now what if I change the Response_CRS? I would expect 
it to give me the same error, since it doesn't affect the requested BBOX. Yet, 
this time no error is produced and I obtain a tiff file (see attachment : 
wcsoutput.tif). The data is all wrong (0 everywhere), though. And if I open 
both WMS and WCS obtained tiffs in QuantumGIS, they appear nowhere near each 
other (for the same bbox... yes).

            I attached the results of GetCapabilities and DescribeCoverage 
requests as well as the MapServer logs (modified to map logs with requests). 
The data source is an HDF5 file (see attachment : data.h5). I suspect an issue 
on MapServer's side but I would be happy if you could tell me it's all my fault 
(and how to fix it).

            Oh, a last thing : it may or may not be related. I had issues with 
the KML/KMZ WMS outputformats. Results for epsg:4326 and epsg:900913 were 
completely different (and wrong, of course). I might open a new discussion on 
this subject at some point.

            Thank you all for your help.
            Guillaume Loiselle

            P.S. : If you want to test, know that epsg:900913 is the same as 
epsg:3857 (which is more likely to be in your epsg reference file).


            _______________________________________________
            mapserver-users mailing list
            [email protected] 
<mailto:[email protected]>
            http://lists.osgeo.org/mailman/listinfo/mapserver-users




--
Håvard Tveite
Department of Mathematical Sciences and Technology, UMB
Drøbakveien 31, POBox 5003, N-1432 Ås, NORWAY
Phone: +47 64965483 Fax: +47 64965401 http://www.umb.no/imt/
_______________________________________________
mapserver-users mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/mapserver-users

Reply via email to