Thanks for your suggestion; I will try this.<div><br></div><div>My end goal, however, is to set up a MapServer WCS that can be consumed by ArcMap. This project is for a client that uses ArcMap, and it's not an option for me to tell them to use something else. So, although writing WCS queries by hand for testing might be instructive, as far as I know there is no way for me to customize the queries that ArcMap generates --- the service itself needs to be set up so that it will just work in ArcMap.</div>
<div><br></div><div>Any ideas on how to do that? Are there any known compatibility issues with ArcMap? Should I force ArcMap to use a specific WCS version? (I've been telling it to use version 1.0.0 so far.)</div><div>
<br></div><div>--Mark<br><br><div class="gmail_quote">On Fri, Jun 3, 2011 at 10:26 AM, Rahkonen Jukka <span dir="ltr"><<a href="mailto:Jukka.Rahkonen@mmmtike.fi">Jukka.Rahkonen@mmmtike.fi</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex;">
Hi,<br>
<br>
Response_crs exists for letting user to select the output CRS. Have you tried it yet?<br>
There are better possibilities for a success and especially for understanding what is really happening if you write WCS requests by hand and semd them with a browser instead of using a GIS client. Request examples at <a href="http://mapserver.org/ogc/wcs_server.html" target="_blank">http://mapserver.org/ogc/wcs_server.html</a> are good to start with.<br>
<br>
-Jukka Rahkonen-<br>
<br>
<br>
-----Alkuperäinen viesti-----<br>
Lähettäjä: <a href="mailto:mapserver-users-bounces@lists.osgeo.org">mapserver-users-bounces@lists.osgeo.org</a> puolesta: Mark Phillips<br>
Lähetetty: to 2.6.2011 22:16<br>
Vastaanottaja: <a href="mailto:mapserver-users@lists.osgeo.org">mapserver-users@lists.osgeo.org</a><br>
Aihe: [mapserver-users] WCS reprojection issue<br>
<div><div></div><div class="h5"><br>
Hi,<br>
<br>
I'm struggling with an issue related to projections in a MapServer WCS<br>
service. I've got a raster data file that is in EPSG:2163, and I'd like to<br>
serve it via WCS in EPSG:4326. The only way I have managed to get this to<br>
work is to reproject the data manually in advance and then reference the<br>
ESPG:4326 copy of the data in my mapfile. When I do that, with the mapfile<br>
below, the service seems to work beautifully. What I really want, though, is<br>
to be able to leave my data in its original projection and have MapServer<br>
reproject on the fly.<br>
<br>
When I change the mapfile to use the original ESPG:2163 copy of the data,<br>
and update the relevant projections in the mapfile from 4326 to 2163,<br>
though, the service breaks, in the sense that the data shows up in (very<br>
much) the wrong place. I first tried just changing the<br>
projections commented in the mapfile below with ###LAYER_PROJECTION and<br>
###LAYER_WCS_SRS from 4326 to 2163 (after changing the DATA line in the<br>
layer to point to the EPSG:2163 copy of the data, of course), but that<br>
didn't work. I then experimented with changing the ###SERVICE_WCS_SRS line<br>
and/or the ###MAP_PROJECTION as well, but no luck.<br>
<br>
What is especially confusing is that I *can* get this to work for a WMS<br>
service. In the case of WMS, MapServer seems quite happy to serve EPSG:2163<br>
data in ESPG:4326, as long as the ###LAYER_PROJECTION and ###LAYER_WCS_SRS<br>
lines say EPSG:2163 (with "wcs" replaced by "wms", of course).<br>
<br>
I'm sure that part of my problem is that I don't completely understand the<br>
purpose served by each of the various projections specified in the mapfile.<br>
Can someone offer some insight into what the correct approach is to getting<br>
MapServer to reproject on the fly in a WCS?<br>
<br>
In case it matters: I'm using MapServer 6.0, and am testing the WCS with<br>
ArcMap 9.3.<br>
<br>
Thanks!<br>
<br>
--Mark<br>
<br>
Here is my working mapfile for ESPG:4326:<br>
<br>
MAP<br>
NAME lucky<br>
STATUS ON<br>
SIZE 600 400<br>
EXTENT -83.938 34.746 -81.340 36.361<br>
IMAGECOLOR 255 255 255<br>
<br>
WEB<br>
METADATA<br>
"ows_enable_request" "*"<br>
"wcs_label" "My WCS Service"<br>
"wcs_srs" "EPSG:4326" ###SERVICE_WCS_SRS<br>
END<br>
END<br>
<br>
PROJECTION<br>
"init=epsg:4326" ###MAP_PROJECTION<br>
END<br>
<br>
LAYER<br>
NAME mylayer<br>
TYPE RASTER<br>
DUMP TRUE<br>
STATUS OFF<br>
DATA sampledata_wgs84/tif_wgs84.tif<br>
INCLUDE "colormap.cmap"<br>
PROJECTION<br>
"init=epsg:4326" ###LAYER_PROJECTION<br>
END<br>
METADATA<br>
"wcs_srs" "EPSG:4326" ###LAYER_WCS_SRS<br>
"wcs_label" "My Layer"<br>
"wcs_rangeset_name" "My Range"<br>
"wcs_rangeset_label" "My Label"<br>
END<br>
END<br>
END<br>
<br>
</div></div></blockquote></div><br></div>