From sethg at geographika.co.uk  Fri Oct  1 10:12:00 2021
From: sethg at geographika.co.uk (Seth G)
Date: Fri, 01 Oct 2021 19:12:00 +0200
Subject: [mapserver-users] State of MapServer 2021 Slides
Message-ID: <72339445-14b7-4553-b7de-218cebf65a54@www.fastmail.com>

Hi all,

For those who weren't able to attend FOSS4G this year the slides of the State of MapServer 2021 presentation can be found at https://geographika.github.io/mapserver-state-2021/

Also of interest to MapServer users are:

Mapserver ? Labeling and Circles ? my best tips and tricks [1]
Geostyler Mapfile Parser [2]

I believe the talks will all be published online in around a month.

Seth

[1] https://2021-foss4g.venueless.events/schedule/talks/M8JCVH
[2] https://2021-foss4g.venueless.events/schedule/talks/FG7PBN
--
web:http://geographika.co.uk
twitter: @geographika


From atlantageek at gmail.com  Tue Oct  5 11:51:40 2021
From: atlantageek at gmail.com (Atlanta Geek)
Date: Tue, 5 Oct 2021 14:51:40 -0400
Subject: [mapserver-users] (no subject)
Message-ID: <CACnRHCRWGLj1Lo2iAS40+mdcB3gas7aK7WZ4WwcVY9S_zu_FiQ@mail.gmail.com>

Im attempting to configure mapcache.
I assume that this is covered in this mailing list as well.  Mapcache is
returning images of the map for me but the cache files are not being
written.  Its strange because I can use the seeding command with the same
xml file and its populating the expected directory but when I go through
the mapcache script its not writing to that cache directory.

Any guidance would be appreciated.

The URL Im testing with is this.
http://192.168.112.225/mapcache?service=WMS&request=GetMap&layers=default&styles=&format=image%2Fpng&transparent=true&version=1.1.1&width=2048&height=2048&srs=EPSG%3A3857&bbox=-10018754.171394622,0,0,10018754.171394628



mapcache.xml is as follows:
<mapcache mode="combined_mirror">
   <cache name="disk" type="disk">
           <base>/tmp/mapcache</base>
      <symlink_blank/>
   </cache>
<format name="png" type="PNG">
   <compression>fast</compression>
   <colors>256</colors>
</format>
   <source name="vmap0" type="wms">
      <getmap>
         <params>
            <FORMAT>image/jpeg</FORMAT>
            <LAYERS>default</LAYERS>
            <map>/usr/lib/cgi-bin/osm-bw.map</map>
            <transparent>true</transparent>
         </params>
      </getmap>
      <http>
         <url>http://192.168.112.225/cgi-bin/mapserv?</url>
      </http>
   </source>
   <tileset name="mirror">
      <source>vmap0</source>
      <cache>disk</cache>
      <grid>WGS84</grid>
      <grid>g</grid>
      <format>png</format>
              <expires>3600</expires>
   </tileset>
   <default_format>png</default_format>

   <service type="wms" enabled="true">
      <format allow_client_override="true/false">png</format>
      <full_wms>assemble</full_wms>
      <resample_mode>bilinear</resample_mode>
      <maxsize>4096</maxsize>
      <forwarding_rule name="catch all">
            <http>
              <url>
http://192.168.112.225/cgi-bin/mapserv?transparent=true&map=/usr/lib/cgi-bin/osm-bw.map
</url>
            </http>
      </forwarding_rule>
   </service>
   <service type="wmts" enabled="true"/>
   <service type="tms" enabled="true"/>
   <service type="kml" enabled="true"/>
   <service type="gmaps" enabled="true"/>
   <service type="ve" enabled="true"/>
   <service type="demo" enabled="true"/>
   <errors>report</errors>
   <lock_dir>/tmp</lock_dir>
   <log_level>debug</log_level>
</mapcache>
-- 
http://www.atlantageek.com
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211005/34bf0939/attachment.htm>

From sdlime at gmail.com  Thu Oct  7 05:57:13 2021
From: sdlime at gmail.com (Steve Lime)
Date: Thu, 7 Oct 2021 07:57:13 -0500
Subject: [mapserver-users] (no subject)
In-Reply-To: <CACnRHCRWGLj1Lo2iAS40+mdcB3gas7aK7WZ4WwcVY9S_zu_FiQ@mail.gmail.com>
References: <CACnRHCRWGLj1Lo2iAS40+mdcB3gas7aK7WZ4WwcVY9S_zu_FiQ@mail.gmail.com>
Message-ID: <CAMrKZ98-U0+Lo4+xZbA83BhN9utzRj3TUzzZzn6c01nZi01+bA@mail.gmail.com>

I've not used a similar setup but it seems like there's no need to create
tiles if you're passing WMS calls through to another WMS server. Do you get
tiles if you make tiled service requests?

On Tue, Oct 5, 2021 at 1:52 PM Atlanta Geek <atlantageek at gmail.com> wrote:

> Im attempting to configure mapcache.
> I assume that this is covered in this mailing list as well.  Mapcache is
> returning images of the map for me but the cache files are not being
> written.  Its strange because I can use the seeding command with the same
> xml file and its populating the expected directory but when I go through
> the mapcache script its not writing to that cache directory.
>
> Any guidance would be appreciated.
>
> The URL Im testing with is this.
>
> http://192.168.112.225/mapcache?service=WMS&request=GetMap&layers=default&styles=&format=image%2Fpng&transparent=true&version=1.1.1&width=2048&height=2048&srs=EPSG%3A3857&bbox=-10018754.171394622,0,0,10018754.171394628
>
>
>
> mapcache.xml is as follows:
> <mapcache mode="combined_mirror">
>    <cache name="disk" type="disk">
>            <base>/tmp/mapcache</base>
>       <symlink_blank/>
>    </cache>
> <format name="png" type="PNG">
>    <compression>fast</compression>
>    <colors>256</colors>
> </format>
>    <source name="vmap0" type="wms">
>       <getmap>
>          <params>
>             <FORMAT>image/jpeg</FORMAT>
>             <LAYERS>default</LAYERS>
>             <map>/usr/lib/cgi-bin/osm-bw.map</map>
>             <transparent>true</transparent>
>          </params>
>       </getmap>
>       <http>
>          <url>http://192.168.112.225/cgi-bin/mapserv?</url>
>       </http>
>    </source>
>    <tileset name="mirror">
>       <source>vmap0</source>
>       <cache>disk</cache>
>       <grid>WGS84</grid>
>       <grid>g</grid>
>       <format>png</format>
>               <expires>3600</expires>
>    </tileset>
>    <default_format>png</default_format>
>
>    <service type="wms" enabled="true">
>       <format allow_client_override="true/false">png</format>
>       <full_wms>assemble</full_wms>
>       <resample_mode>bilinear</resample_mode>
>       <maxsize>4096</maxsize>
>       <forwarding_rule name="catch all">
>             <http>
>               <url>
> http://192.168.112.225/cgi-bin/mapserv?transparent=true&map=/usr/lib/cgi-bin/osm-bw.map
> </url>
>             </http>
>       </forwarding_rule>
>    </service>
>    <service type="wmts" enabled="true"/>
>    <service type="tms" enabled="true"/>
>    <service type="kml" enabled="true"/>
>    <service type="gmaps" enabled="true"/>
>    <service type="ve" enabled="true"/>
>    <service type="demo" enabled="true"/>
>    <errors>report</errors>
>    <lock_dir>/tmp</lock_dir>
>    <log_level>debug</log_level>
> </mapcache>
> --
> http://www.atlantageek.com
> _______________________________________________
> MapServer-users mailing list
> MapServer-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/mapserver-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211007/c814ca3b/attachment.htm>

From bob.basques at ci.stpaul.mn.us  Thu Oct  7 15:43:01 2021
From: bob.basques at ci.stpaul.mn.us (Basques, Bob (CI-StPaul))
Date: Thu, 7 Oct 2021 22:43:01 +0000
Subject: [mapserver-users] OSGeo Twin Cities (aka TCMUG) Local Chapter
 Meeting - Oct 13th
Message-ID: <D6083A53-F6AB-4D04-ACD6-E57AE20EE193@ci.stpaul.mn.us>

All,

Last month?s meeting, while it technically happened, didn?t amount to too much cross pollination, since there were only two of us that showed up.  :c).  We did get rained out at the end with a short lived sun shower.

No matter.  This month I have a virtual Presentation for the group related to the EPC/ UUMPT<https://dp-pro.com/minnesota-underground-utilities-mapping-project-team/> group I?m a part of.  Some more info available here<https://www.youtube.com/channel/UC3hwp5_9t3BkiTt-hyALArg> in Video meeting recordings.  We?re just starting back up from the Summer  hiatus and the work going on in this group I think will be interesting to a number of, if not most of you.

I?m going to try and get Damon Nelton, another member of the EPC/UUMPT group, to tag team on this presentation with me as a rehearsal for the GIS/LIS meeting where we?ll be presenting something similar.   Not sure yet if he can be ready before next Wednesday though, in which case I?ll tackle it on my own.



Online Presentation:  https://meet.jit.si/osgeo_tcmug

When Oct, 13th, 4:30 PM



bobb

OSGeo, Twin Cities (aka TCMUG), MN, USA Local Chapter Page.<https://www.osgeo.org/local-chapters/twin-cities-mn-usa-chapter/>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211007/13188826/attachment.htm>

From cgodkin at gmail.com  Thu Oct  7 16:12:31 2021
From: cgodkin at gmail.com (Carl Godkin)
Date: Thu, 7 Oct 2021 16:12:31 -0700
Subject: [mapserver-users] Mapserver + WMS
Message-ID: <CAE=8kmmYT8d-BUXeStKAzAJzA25apvABGwbAVq4PiNYOhkhbJA@mail.gmail.com>

Hi,

I'm new to mapserver but fairly experienced with GDAL.  I hope that someone
can explain where I've gone wrong with setting up WMS.

I followed the instructions here [
https://worldwind.earth/elevation-server.html] almost exactly to set up a
WMS server and it sort of works.  (I am using MS4W actually, so I didn't
install Apache & Mapserver as that page described and of course I made a
few other adaptations for Windows.)

I downloaded the GEBCO tiles from here [
https://www.gebco.net/data_and_products/gridded_bathymetry_data/] without
any problems.

What I don't understand is that while the GEBCO GeoTIFF files that I am
using have one band that gdalinfo describes as

 Band 1 Block=512x512 Type=Int16, ColorInterp=Gray

when I look at the WMS bands using gdalinfo, I see three byte-sized bands:

Band 1 Block=1024x1024 Type=Byte, ColorInterp=Red
Band 2 Block=1024x1024 Type=Byte, ColorInterp=Green
Band 3 Block=1024x1024 Type=Byte, ColorInterp=Blue

I'd really like to get 16-bit integer elevations from my WMS server but I
don't know why this didn't just happen or how I mis-configured things to
make it not happen.

I know these are not complete details, but rather than write a novel for my
first post here, I'm hoping someone can suggest where the problem might be
or suggest some further documentation to read.

Thanks very much,
carl
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211007/589cdd0f/attachment.htm>

From jukka.rahkonen at maanmittauslaitos.fi  Fri Oct  8 09:50:05 2021
From: jukka.rahkonen at maanmittauslaitos.fi (Rahkonen Jukka (MML))
Date: Fri, 8 Oct 2021 16:50:05 +0000
Subject: [mapserver-users] WCS: Write error at scanline xxxx
Message-ID: <78a286a963bb429eb86e7d924d258643@maanmittauslaitos.fi>

Hi,

When I try to get any bigger output from WCS (in this case 12000x12000 pixels RGB) with a request like

http://localhost:8060/cgi-bin/mapserv.exe?map=c:/ms4w_data/wcs_vv_jr.map&service=WCS&version=2.0.0&request=GetCoverage&CoverageID=orto_vv&format=image/tiff&geotiff:tiling=true&geotiff:tileheight=256&geotiff:tilewidth=256<http://localhost:8060/cgi-bin/mapserv.exe?map=c:/ms4w_data/wcs_vv_jr.map&service=WCS&version=2.0.0&request=GetCoverage&CoverageID=orto_vv&format=image/tiff&geotiff:tiling=true&geotiff:tileheight=256&geotiff:tilewidth=256&geotiff:compression=jpeg>
I get errors like this

<ows:ExceptionReport version="2.0.0" xml:lang="en-US" xsi:schemaLocation=http://www.opengis.net/ows/2.0 http://schemas.opengis.net/ows/2.0/owsExceptionReport.xsd<http://www.opengis.net/ows/2.0%20http:/schemas.opengis.net/ows/2.0/owsExceptionReport.xsd>>
<ows:Exception exceptionCode="NoApplicableCode" locator="mapserv">
<ows:ExceptionText>
msWCSWriteFile20(): General error message. msSaveImage() failed msSaveImageGDAL(): General error message. Failed to create output GTiff file. TIFFAppendToStrip:Write error at scanline 1024
</ows:ExceptionText>
</ows:Exception>
</ows:ExceptionReport>

I managed to get bigger images = more pixels by asking for compressed output with &geotiff:compression=jpeg
I wonder what is the limiting resource and if there is something I could improve with some configuration. Because the output is a tiled GeoTIFF but the error is about writing a strip I guess that Mapserver is writing first an interim TIFF.

-Jukka Rahkonen-
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211008/ed4b8347/attachment.htm>

From ycharieau at sma-netagis.fr  Mon Oct 11 00:56:49 2021
From: ycharieau at sma-netagis.fr (Yoann CHARIEAU)
Date: Mon, 11 Oct 2021 09:56:49 +0200 (CEST)
Subject: [mapserver-users] (no subject)
In-Reply-To: <CAMrKZ98-U0+Lo4+xZbA83BhN9utzRj3TUzzZzn6c01nZi01+bA@mail.gmail.com>
References: <CACnRHCRWGLj1Lo2iAS40+mdcB3gas7aK7WZ4WwcVY9S_zu_FiQ@mail.gmail.com>
 <CAMrKZ98-U0+Lo4+xZbA83BhN9utzRj3TUzzZzn6c01nZi01+bA@mail.gmail.com>
Message-ID: <20817046.118470.1633939009493.JavaMail.zimbra@sma-netagis.fr>

Did you check permissions on your /tmp/mapcache directory ?
It should be writeable by your http server user (usually www-data).

Yoann

----- Mail original -----
De: "Steve Lime" <sdlime at gmail.com>
?: "Atlanta Geek" <atlantageek at gmail.com>
Cc: "mapserver-users" <mapserver-users at lists.osgeo.org>
Envoy?: Jeudi 7 Octobre 2021 14:57:13
Objet: Re: [mapserver-users] (no subject)

I've not used a similar setup but it seems like there's no need to create tiles if you're passing WMS calls through to another WMS server. Do you get tiles if you make tiled service requests? 

On Tue, Oct 5, 2021 at 1:52 PM Atlanta Geek < [ mailto:atlantageek at gmail.com | atlantageek at gmail.com ] > wrote: 



Im attempting to configure mapcache. 
I assume that this is covered in this mailing list as well. Mapcache is returning images of the map for me but the cache files are not being written. Its strange because I can use the seeding command with the same xml file and its populating the expected directory but when I go through the mapcache script its not writing to that cache directory. 

Any guidance would be appreciated. 

The URL Im testing with is this. 
[ http://192.168.112.225/mapcache?service=WMS&request=GetMap&layers=default&styles=&format=image%2Fpng&transparent=true&version=1.1.1&width=2048&height=2048&srs=EPSG%3A3857&bbox=-10018754.171394622,0,0,10018754.171394628 | http://192.168.112.225/mapcache?service=WMS&request=GetMap&layers=default&styles=&format=image%2Fpng&transparent=true&version=1.1.1&width=2048&height=2048&srs=EPSG%3A3857&bbox=-10018754.171394622,0,0,10018754.171394628 ] 



mapcache.xml is as follows: 
<mapcache mode="combined_mirror"> 
<cache name="disk" type="disk"> 
<base>/tmp/mapcache</base> 
<symlink_blank/> 
</cache> 
<format name="png" type="PNG"> 
<compression>fast</compression> 
<colors>256</colors> 
</format> 
<source name="vmap0" type="wms"> 
<getmap> 
<params> 
<FORMAT>image/jpeg</FORMAT> 
<LAYERS>default</LAYERS> 
<map>/usr/lib/cgi-bin/osm-bw.map</map> 
<transparent>true</transparent> 
</params> 
</getmap> 
<http> 
<url> [ http://192.168.112.225/cgi-bin/mapserv | http://192.168.112.225/cgi-bin/mapserv ] ?</url> 
</http> 
</source> 
<tileset name="mirror"> 
<source>vmap0</source> 
<cache>disk</cache> 
<grid>WGS84</grid> 
<grid>g</grid> 
<format>png</format> 
<expires>3600</expires> 
</tileset> 
<default_format>png</default_format> 

<service type="wms" enabled="true"> 
<format allow_client_override="true/false">png</format> 
<full_wms>assemble</full_wms> 
<resample_mode>bilinear</resample_mode> 
<maxsize>4096</maxsize> 
<forwarding_rule name="catch all"> 
<http> 
<url> [ http://192.168.112.225/cgi-bin/mapserv?transparent=true&map=/usr/lib/cgi-bin/osm-bw.map | http://192.168.112.225/cgi-bin/mapserv?transparent=true&map=/usr/lib/cgi-bin/osm-bw.map ] </url> 
</http> 
</forwarding_rule> 
</service> 
<service type="wmts" enabled="true"/> 
<service type="tms" enabled="true"/> 
<service type="kml" enabled="true"/> 
<service type="gmaps" enabled="true"/> 
<service type="ve" enabled="true"/> 
<service type="demo" enabled="true"/> 
<errors>report</errors> 
<lock_dir>/tmp</lock_dir> 
<log_level>debug</log_level> 
</mapcache> 
-- 
[ http://www.atlantageek.com/ | http://www.atlantageek.com ] 
_______________________________________________ 
MapServer-users mailing list 
[ mailto:MapServer-users at lists.osgeo.org | MapServer-users at lists.osgeo.org ] 
[ https://lists.osgeo.org/mailman/listinfo/mapserver-users | https://lists.osgeo.org/mailman/listinfo/mapserver-users ] 




_______________________________________________ 
MapServer-users mailing list 
MapServer-users at lists.osgeo.org 
https://lists.osgeo.org/mailman/listinfo/mapserver-users 



From aimdev at mail.com  Mon Oct 11 11:04:31 2021
From: aimdev at mail.com (aimdev at mail.com)
Date: Mon, 11 Oct 2021 19:04:31 +0100
Subject: [mapserver-users] Xml support in Mapserver v8
Message-ID: <66EF6C84-562A-4862-B882-FBE05439D3C5@mail.com>

Hi

I maintain my maps in xml, converting to map files utilising xsltproc.

1. will there be an xsltproc in v8?
2. will there be an updated xsd file for v8?
3  if 1 & 2 are yes, is it possible to get an early release of the xsd file?

Thanks

Aimee


From jmckenna at gatewaygeomatics.com  Mon Oct 11 11:35:41 2021
From: jmckenna at gatewaygeomatics.com (Jeff McKenna)
Date: Mon, 11 Oct 2021 15:35:41 -0300
Subject: [mapserver-users] Mapserver + WMS
In-Reply-To: <CAE=8kmmYT8d-BUXeStKAzAJzA25apvABGwbAVq4PiNYOhkhbJA@mail.gmail.com>
References: <CAE=8kmmYT8d-BUXeStKAzAJzA25apvABGwbAVq4PiNYOhkhbJA@mail.gmail.com>
Message-ID: <449649f1-4a96-597c-f2e7-155967281c1d@gatewaygeomatics.com>

Hi Carl,

First, welcome to the MapServer community!

Thanks for the detailed question, and pointing to the source data, this 
always helps to give a proper answer.

You've done well and are very close.  Here are some points to consider:

- the OGC specification for serving raw raster data is actually WCS 
(which MapServer does well, see https://mapserver.org/ogc/wcs_server.html )

- but you're right that a cool trick is to use the (more) popular WMS 
specification and set an OUTPUTFORMAT in your server's mapfile, to 
generate a GeoTIFF through a WMS GetMap request.

- I see that you've already setup an OUTPUTFORMAT such as:

OUTPUTFORMAT
   NAME "GEOTIFF_16"
   DRIVER "GDAL/GTiff"
   MIMETYPE "image/tiff"
   IMAGEMODE INT16
   EXTENSION "tif"
END

- now the trick is to call that specific outputformat.  I recommend 
always testing this at the commandline first, with shp2img (which will 
be renamed to 'map2img' shortly, for the MapServer 8.0 release).  It 
takes a mapfile and generates an output (or gives exact error).

- since you are using MS4W, just open a CMD window and cd into /ms4w. 
Then execute 'setenv.bat' to set necessary paths.  Then, in that same 
windows, cd to where your mapfile is, and execute the following (notice 
the "-i" switch, which requests the named outputformat set in the mapfile) :

   shp2img -m local.map -o ttt16.tif -map_debug 3 -i GEOTIFF_16

- that command should return a draw speed for each layer, and generate 
your 16-bit GeoTIFF.  And look closely at that response, the first words 
will be something like:

    rendering using outputformat named GEOTIFF_16 (GDAL/GTiff)

Great! That's what we want.

- once you're satisfied, now we can turn to WMS requests.  But always 
first start with removing any "WARNING" messages in your GetCapabilities 
request, and then also verifying that the <Format> section of that 
response includes a listed "image/tiff" format option.  (example 
GetCapabilities request) :

 
http://127.0.0.1/cgi-bin/mapserv.exe?map=/ms4w/apps/ticket-wms-gebco-outputformat/local.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetCapabilities

- once you're happy with the GetCapabilities response, then you can try 
a GetMap request, that includes "&FORMAT=image/tiff"


 
http://127.0.0.1/cgi-bin/mapserv.exe?map=/ms4w/apps/ticket-wms-gebco-outputformat/local.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-90,-180,90,180&CRS=EPSG:4326&WIDTH=346&HEIGHT=173&LAYERS=gebco&STYLES=&FORMAT=image/tiff&DPI=96&MAP_RESOLUTION=96&FORMAT_OPTIONS=dpi:96&TRANSPARENT=TRUE

- or you can stay right there at your commandline, and make that GetMap 
request there with the "mapserv -nh QUERY_STRING="" trick, and pipe that 
to a new image file, such as:

   mapserv -nh 
QUERY_STRING="map=/ms4w/apps/ticket-wms-gebco-outputformat/local.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-90,-180,90,180&CRS=EPSG:4326&WIDTH=346&HEIGHT=173&LAYERS=gebco&STYLES=&FORMAT=image/tiff&DPI=96&MAP_RESOLUTION=96&FORMAT_OPTIONS=dpi:96&TRANSPARENT=TRUE" 
 > getmap.tif

Hope that helps.

Also thanks for the positive feedback on MS4W, its goal is to serve big 
data fast for organizations, on enterprise Windows servers, and it 
sounds like it has met your needs well.  Look for a big MS4W release soon ;)

And a happy (Canadian) Thanksgiving holiday to you.

-jeff



-- 
Jeff McKenna
GatewayGeo: Developers of MS4W, MapServer Consulting and Training
co-founder of FOSS4G
http://gatewaygeo.com/



On 2021-10-07 8:12 p.m., Carl Godkin wrote:
> Hi,
> 
> I'm new to mapserver but fairly experienced with GDAL.? I hope that 
> someone can explain where I've gone wrong with setting up WMS.
> 
> I followed the instructions here 
> [https://worldwind.earth/elevation-server.html 
> <https://worldwind.earth/elevation-server.html>] almost exactly to set 
> up a WMS server and it sort of works.? (I am using MS4W actually, so I 
> didn't install Apache & Mapserver as that page described and of course I 
> made a few other adaptations for Windows.)
> 
> I downloaded the GEBCO tiles from here 
> [https://www.gebco.net/data_and_products/gridded_bathymetry_data/ 
> <https://www.gebco.net/data_and_products/gridded_bathymetry_data/>] 
> without any problems.
> 
> What I don't understand is that while the GEBCO GeoTIFF files that I am 
> using have one band that gdalinfo describes as
> 
>  ?Band 1 Block=512x512 Type=Int16, ColorInterp=Gray
> 
> when I look at the WMS bands using gdalinfo, I see three byte-sized bands:
> 
> Band 1 Block=1024x1024 Type=Byte, ColorInterp=Red
> Band 2 Block=1024x1024 Type=Byte, ColorInterp=Green
> Band 3 Block=1024x1024 Type=Byte, ColorInterp=Blue
> 
> I'd really like to get 16-bit integer elevations from my WMS server but 
> I don't know why this didn't just happen or how I mis-configured things 
> to make it not happen.
> 
> I know these are not complete details, but rather than write a novel for 
> my first post here, I'm hoping someone can suggest where the problem 
> might be or suggest some further documentation to read.
> 
> Thanks very much,
> carl
> 
> 
> 


From jmckenna at gatewaygeomatics.com  Mon Oct 11 11:47:12 2021
From: jmckenna at gatewaygeomatics.com (Jeff McKenna)
Date: Mon, 11 Oct 2021 15:47:12 -0300
Subject: [mapserver-users] Mapserver + WMS
In-Reply-To: <449649f1-4a96-597c-f2e7-155967281c1d@gatewaygeomatics.com>
References: <CAE=8kmmYT8d-BUXeStKAzAJzA25apvABGwbAVq4PiNYOhkhbJA@mail.gmail.com>
 <449649f1-4a96-597c-f2e7-155967281c1d@gatewaygeomatics.com>
Message-ID: <f3d77dee-1977-3370-b128-5532df156716@gatewaygeomatics.com>

Hi Carl,

Also, as you travel down this path, please aim to adding your upcoming 
elevation WMS service into the MapServer Site Gallery : 
https://github.com/MapServer/MapServer/wiki/MapServer-Site-Gallery

This really helps new users see the powers of MapServer.

Thanks!

-jeff



-- 
Jeff McKenna
GatewayGeo: Developers of MS4W, MapServer Consulting and Training
co-founder of FOSS4G
http://gatewaygeo.com/


On 2021-10-11 3:35 p.m., Jeff McKenna wrote:
> Hi Carl,
> 
> First, welcome to the MapServer community!
> 
> Thanks for the detailed question, and pointing to the source data, this 
> always helps to give a proper answer.
> 
> You've done well and are very close.? Here are some points to consider:
> 
> - the OGC specification for serving raw raster data is actually WCS 
> (which MapServer does well, see https://mapserver.org/ogc/wcs_server.html )
> 
> - but you're right that a cool trick is to use the (more) popular WMS 
> specification and set an OUTPUTFORMAT in your server's mapfile, to 
> generate a GeoTIFF through a WMS GetMap request.
> 
> - I see that you've already setup an OUTPUTFORMAT such as:
> 
> OUTPUTFORMAT
>  ? NAME "GEOTIFF_16"
>  ? DRIVER "GDAL/GTiff"
>  ? MIMETYPE "image/tiff"
>  ? IMAGEMODE INT16
>  ? EXTENSION "tif"
> END
> 
> - now the trick is to call that specific outputformat.? I recommend 
> always testing this at the commandline first, with shp2img (which will 
> be renamed to 'map2img' shortly, for the MapServer 8.0 release).? It 
> takes a mapfile and generates an output (or gives exact error).
> 
> - since you are using MS4W, just open a CMD window and cd into /ms4w. 
> Then execute 'setenv.bat' to set necessary paths.? Then, in that same 
> windows, cd to where your mapfile is, and execute the following (notice 
> the "-i" switch, which requests the named outputformat set in the 
> mapfile) :
> 
>  ? shp2img -m local.map -o ttt16.tif -map_debug 3 -i GEOTIFF_16
> 
> - that command should return a draw speed for each layer, and generate 
> your 16-bit GeoTIFF.? And look closely at that response, the first words 
> will be something like:
> 
>  ?? rendering using outputformat named GEOTIFF_16 (GDAL/GTiff)
> 
> Great! That's what we want.
> 
> - once you're satisfied, now we can turn to WMS requests.? But always 
> first start with removing any "WARNING" messages in your GetCapabilities 
> request, and then also verifying that the <Format> section of that 
> response includes a listed "image/tiff" format option.? (example 
> GetCapabilities request) :
> 
> 
> http://127.0.0.1/cgi-bin/mapserv.exe?map=/ms4w/apps/ticket-wms-gebco-outputformat/local.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetCapabilities 
> 
> 
> - once you're happy with the GetCapabilities response, then you can try 
> a GetMap request, that includes "&FORMAT=image/tiff"
> 
> 
> 
> http://127.0.0.1/cgi-bin/mapserv.exe?map=/ms4w/apps/ticket-wms-gebco-outputformat/local.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-90,-180,90,180&CRS=EPSG:4326&WIDTH=346&HEIGHT=173&LAYERS=gebco&STYLES=&FORMAT=image/tiff&DPI=96&MAP_RESOLUTION=96&FORMAT_OPTIONS=dpi:96&TRANSPARENT=TRUE 
> 
> 
> - or you can stay right there at your commandline, and make that GetMap 
> request there with the "mapserv -nh QUERY_STRING="" trick, and pipe that 
> to a new image file, such as:
> 
>  ? mapserv -nh 
> QUERY_STRING="map=/ms4w/apps/ticket-wms-gebco-outputformat/local.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-90,-180,90,180&CRS=EPSG:4326&WIDTH=346&HEIGHT=173&LAYERS=gebco&STYLES=&FORMAT=image/tiff&DPI=96&MAP_RESOLUTION=96&FORMAT_OPTIONS=dpi:96&TRANSPARENT=TRUE" 
>  > getmap.tif
> 
> Hope that helps.
> 
> Also thanks for the positive feedback on MS4W, its goal is to serve big 
> data fast for organizations, on enterprise Windows servers, and it 
> sounds like it has met your needs well.? Look for a big MS4W release 
> soon ;)
> 
> And a happy (Canadian) Thanksgiving holiday to you.
> 
> -jeff
> 
> 
> 



From cgodkin at gmail.com  Mon Oct 11 15:35:20 2021
From: cgodkin at gmail.com (Carl Godkin)
Date: Mon, 11 Oct 2021 15:35:20 -0700
Subject: [mapserver-users] Mapserver + WMS
In-Reply-To: <449649f1-4a96-597c-f2e7-155967281c1d@gatewaygeomatics.com>
References: <CAE=8kmmYT8d-BUXeStKAzAJzA25apvABGwbAVq4PiNYOhkhbJA@mail.gmail.com>
 <449649f1-4a96-597c-f2e7-155967281c1d@gatewaygeomatics.com>
Message-ID: <CAE=8kmmS-y9jx-iw13tj8n4axgsWXJDZvsO2YQYiP-5PDxQu4w@mail.gmail.com>

Hi Jeff,

Thank you very much for the detailed response.  When I run the commands you
provided, everything matches and the output GeoTIFF file definitely has the
type "Type=Int16."

I was originally asking why I didn't see that type when querying the server
with the sub-datasets returned by this command

> gdalinfo "WMS:http://localhost:8080//elev"

but I'm probably just coloring outside the lines of what WMS is meant to
support.

I think that the next thing I will do is use WCS for this instead.  I
confess that I went with WMS since I'm more familiar with it.  I expected a
single-band 16-bit "image" to act like a 2D grid which is what happens when
I import "BIL" files or other single-band images.

Thanks for stepping me through this and pointing me at WCS.

carl


On Mon, Oct 11, 2021 at 11:35 AM Jeff McKenna <jmckenna at gatewaygeomatics.com>
wrote:

> Hi Carl,
>
> First, welcome to the MapServer community!
>
> Thanks for the detailed question, and pointing to the source data, this
> always helps to give a proper answer.
>
> You've done well and are very close.  Here are some points to consider:
>
> - the OGC specification for serving raw raster data is actually WCS
> (which MapServer does well, see https://mapserver.org/ogc/wcs_server.html
> )
>
> - but you're right that a cool trick is to use the (more) popular WMS
> specification and set an OUTPUTFORMAT in your server's mapfile, to
> generate a GeoTIFF through a WMS GetMap request.
>
> - I see that you've already setup an OUTPUTFORMAT such as:
>
> OUTPUTFORMAT
>    NAME "GEOTIFF_16"
>    DRIVER "GDAL/GTiff"
>    MIMETYPE "image/tiff"
>    IMAGEMODE INT16
>    EXTENSION "tif"
> END
>
> - now the trick is to call that specific outputformat.  I recommend
> always testing this at the commandline first, with shp2img (which will
> be renamed to 'map2img' shortly, for the MapServer 8.0 release).  It
> takes a mapfile and generates an output (or gives exact error).
>
> - since you are using MS4W, just open a CMD window and cd into /ms4w.
> Then execute 'setenv.bat' to set necessary paths.  Then, in that same
> windows, cd to where your mapfile is, and execute the following (notice
> the "-i" switch, which requests the named outputformat set in the mapfile)
> :
>
>    shp2img -m local.map -o ttt16.tif -map_debug 3 -i GEOTIFF_16
>
> - that command should return a draw speed for each layer, and generate
> your 16-bit GeoTIFF.  And look closely at that response, the first words
> will be something like:
>
>     rendering using outputformat named GEOTIFF_16 (GDAL/GTiff)
>
> Great! That's what we want.
>
> - once you're satisfied, now we can turn to WMS requests.  But always
> first start with removing any "WARNING" messages in your GetCapabilities
> request, and then also verifying that the <Format> section of that
> response includes a listed "image/tiff" format option.  (example
> GetCapabilities request) :
>
>
>
> http://127.0.0.1/cgi-bin/mapserv.exe?map=/ms4w/apps/ticket-wms-gebco-outputformat/local.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetCapabilities
>
> - once you're happy with the GetCapabilities response, then you can try
> a GetMap request, that includes "&FORMAT=image/tiff"
>
>
>
>
> http://127.0.0.1/cgi-bin/mapserv.exe?map=/ms4w/apps/ticket-wms-gebco-outputformat/local.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-90,-180,90,180&CRS=EPSG:4326&WIDTH=346&HEIGHT=173&LAYERS=gebco&STYLES=&FORMAT=image/tiff&DPI=96&MAP_RESOLUTION=96&FORMAT_OPTIONS=dpi:96&TRANSPARENT=TRUE
>
> - or you can stay right there at your commandline, and make that GetMap
> request there with the "mapserv -nh QUERY_STRING="" trick, and pipe that
> to a new image file, such as:
>
>    mapserv -nh
> QUERY_STRING="map=/ms4w/apps/ticket-wms-gebco-outputformat/local.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-90,-180,90,180&CRS=EPSG:4326&WIDTH=346&HEIGHT=173&LAYERS=gebco&STYLES=&FORMAT=image/tiff&DPI=96&MAP_RESOLUTION=96&FORMAT_OPTIONS=dpi:96&TRANSPARENT=TRUE"
>
>  > getmap.tif
>
> Hope that helps.
>
> Also thanks for the positive feedback on MS4W, its goal is to serve big
> data fast for organizations, on enterprise Windows servers, and it
> sounds like it has met your needs well.  Look for a big MS4W release soon
> ;)
>
> And a happy (Canadian) Thanksgiving holiday to you.
>
> -jeff
>
>
>
> --
> Jeff McKenna
> GatewayGeo: Developers of MS4W, MapServer Consulting and Training
> co-founder of FOSS4G
> http://gatewaygeo.com/
>
>
>
> On 2021-10-07 8:12 p.m., Carl Godkin wrote:
> > Hi,
> >
> > I'm new to mapserver but fairly experienced with GDAL.  I hope that
> > someone can explain where I've gone wrong with setting up WMS.
> >
> > I followed the instructions here
> > [https://worldwind.earth/elevation-server.html
> > <https://worldwind.earth/elevation-server.html>] almost exactly to set
> > up a WMS server and it sort of works.  (I am using MS4W actually, so I
> > didn't install Apache & Mapserver as that page described and of course I
> > made a few other adaptations for Windows.)
> >
> > I downloaded the GEBCO tiles from here
> > [https://www.gebco.net/data_and_products/gridded_bathymetry_data/
> > <https://www.gebco.net/data_and_products/gridded_bathymetry_data/>]
> > without any problems.
> >
> > What I don't understand is that while the GEBCO GeoTIFF files that I am
> > using have one band that gdalinfo describes as
> >
> >   Band 1 Block=512x512 Type=Int16, ColorInterp=Gray
> >
> > when I look at the WMS bands using gdalinfo, I see three byte-sized
> bands:
> >
> > Band 1 Block=1024x1024 Type=Byte, ColorInterp=Red
> > Band 2 Block=1024x1024 Type=Byte, ColorInterp=Green
> > Band 3 Block=1024x1024 Type=Byte, ColorInterp=Blue
> >
> > I'd really like to get 16-bit integer elevations from my WMS server but
> > I don't know why this didn't just happen or how I mis-configured things
> > to make it not happen.
> >
> > I know these are not complete details, but rather than write a novel for
> > my first post here, I'm hoping someone can suggest where the problem
> > might be or suggest some further documentation to read.
> >
> > Thanks very much,
> > carl
> >
> >
> >
> _______________________________________________
> MapServer-users mailing list
> MapServer-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/mapserver-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211011/404d44b7/attachment.htm>

From bob.basques at ci.stpaul.mn.us  Tue Oct 12 10:10:58 2021
From: bob.basques at ci.stpaul.mn.us (Basques, Bob (CI-StPaul))
Date: Tue, 12 Oct 2021 17:10:58 +0000
Subject: [mapserver-users] [REMINDER] OSGeo Twin Cities (aka TCMUG) Local
 Chapter Meeting - Oct 13th
Message-ID: <E3FE2619-F46D-4820-99EE-1FD906CFB6D5@ci.stpaul.mn.us>

UPDATE: Damon will be sharing the presentation duties.



All,

Last month?s meeting, while it technically happened, didn?t amount to too much cross pollination, since there were only two of us that showed up.  :c).  We did get rained out at the end with a short lived sun shower.

No matter.  This month I have a virtual Presentation for the group related to the EPC/ UUMPT<https://dp-pro.com/minnesota-underground-utilities-mapping-project-team/> group I?m a part of.  Some more info available here<https://www.youtube.com/channel/UC3hwp5_9t3BkiTt-hyALArg> in Video meeting recordings.  We?re just starting back up from the Summer  hiatus and the work going on in this group I think will be interesting to a number of, if not most of you.

I?m going to try and get Damon Nelton, another member of the EPC/UUMPT group, to tag team on this presentation with me as a rehearsal for the GIS/LIS meeting where we?ll be presenting something similar.   Not sure yet if he can be ready before next Wednesday though, in which case I?ll tackle it on my own.



Online Presentation:  https://meet.jit.si/osgeo_tcmug

When Oct, 13th, 4:30 PM



bobb

OSGeo, Twin Cities (aka TCMUG), MN, USA Local Chapter Page.<https://www.osgeo.org/local-chapters/twin-cities-mn-usa-chapter/>

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211012/583fdd67/attachment.htm>

From n.schipkov at gmail.com  Thu Oct 14 03:43:22 2021
From: n.schipkov at gmail.com (=?UTF-8?B?0J3QuNC60LjRgtCwINCp0LjQv9C60L7Qsg==?=)
Date: Thu, 14 Oct 2021 13:43:22 +0300
Subject: [mapserver-users] localhost refused to connect
Message-ID: <CAB6P0bM=fH1JEZqgF8UoXgyTu-xpSyfWAifehUOAuMD65rAsGA@mail.gmail.com>

Hello everyone!

Yesterday I opened demo apps from MapCache. Today in the morning the
localhost is not connecting to the server. How to diagnose the
problem?

Thank you in advance!


From bahbouba.909 at gmail.com  Sat Oct 16 04:12:55 2021
From: bahbouba.909 at gmail.com (Boubacar Bah)
Date: Sat, 16 Oct 2021 12:12:55 +0100
Subject: [mapserver-users] asking help
Message-ID: <CALOg8kGKN=hQ1P+tSNnZ47L2iFBj1yJOb0Hf4iR9iQM5YQn5Yw@mail.gmail.com>

hello everyone ; recently l tried to use openlayers to create my own map
with mapserver. my program work well if try just to display a open layer
map in a browser *but if l try to add the interaction tools like modify,
snap , draw* the program don't work. this is my source code bellow if
someone can hem me with the interaction library of openlayer.
<!DOCTYPE html>
<html lang="en">
  <head>
    <meta charset="utf-8">
    <link rel="stylesheet" href="
https://cdn.jsdelivr.net/gh/openlayers/openlayers.github.io at master/en/v6.8.1/css/ol.css
" type="text/css">
    <style>
      .map {
        height: 400px;
        width: 100%;
      }
    </style>
    <script src="
https://cdn.jsdelivr.net/gh/openlayers/openlayers.github.io at master/en/v6.8.1/build/ol.js
"></script>
    <script src="https://unpkg.com/elm-pep"></script>
    <title>OpenLayers example</title>
  </head>
  <body>
    <h2>My bah</h2>
    <div id="map" class="map"></div>
    <div id="mouse-position"></div>
    <div id="raster-button">
      <button onclick= "raster_switch()">Raster</button>
    </div>
    <div id="dessin-button">
      <button onclick="add_creer()">Dessiner</button>
    </div>
    <div id="stop-button">
      <button onclick="stop()">Stop</button>
    </div>
    <div id="modify-button">
      <button onclick="add_modify()">Modifier</button>
    </div>
    <script type="text/javascript">
      var raster = new ol.layer.Tile({
           source: new ol.source.OSM()
      });
      var vector = new ol.layer.Vector({
        source: new ol.source.Vector({
        url: '
https://raw.githubusercontent.com/openlayers/openlayers/main/examples/data/geojson/countries.geojson
',
        format: new ol.format.GeoJSON(),
         })
      });
      var map = new ol.Map({
        controls: ol.control.defaults().extend([
          new ol.control.MousePosition({
            coordinateFormat: ol.coordinate.createStringXY(4),
            projection: 'EPSG:4326',
            className: 'custom-mouse-position',
            target: document.getElementById('mouse-position'),
            undefinedHTML: '&nbsp;'
          }),
          new ol.control.ScaleLine()
        ]),
        target: 'map',
        layers: [raster, vector],
        view: new ol.View({
          center: [0, 0],
          zoom: 2
        })

      });
      function raster_switch(){
        if (raster.getVisible()){   //si le raster est visible
        raster.setVisible(false)    //Nous l'?teignons
        }
        else{
          raster.setVisible(true)   //Sinon nous l'affichons
        }
      }
this is my code witch work but if l try to add this code below
let draw; // global so we can remove it later
      function add_creer(){
          draw= new.ol.interaction.Draw({
            source: vector.getSource(),
            type: "Polygon",
          });
          map.addInteraction(draw);
          //appeler snap apres draw
          add_snap()
      }
the browser send a blank page if some one can help me to correct this code
thanks
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211016/fcbb4393/attachment.htm>

From cgodkin at gmail.com  Mon Oct 18 11:23:39 2021
From: cgodkin at gmail.com (Carl Godkin)
Date: Mon, 18 Oct 2021 11:23:39 -0700
Subject: [mapserver-users] Getting started with WCS on mapserver...
Message-ID: <CAE=8kmkz-=m02eQ8tKVAd65VuB1whLpGuSCw7+f2oTFRS2YRjw@mail.gmail.com>

Hi,

A week and a half ago I asked about using Mapserver +WMS to serve elevation
data since I was following a tutorial that discussed WMS.

In one of the replies I received, Jeff suggested that I actually wanted WCS
for my use case which was a forehead slap + d'oh! moment for me.  Of course
I do!

Anyway, I have followed the Mapserver WCS docs and everything seems to work
well.

However, I have one question that may be kind of niche or will perhaps
betray my overall ignorance.

A while ago I wrote a WCS client for another project and it expects WCS
version 2.0 requests.

The map file that I made up (based closely on the one in the docs:
https://mapserver.org/ogc/wcs_server.html#example-wcs-server-mapfile )
seems to work fine for WCS 1.0, 1.1, and 2.0.  However, I don't seem to get
a bounding box from GetCapabilities in the CoverageSummary from my 2.0
requests whereas I certainly do from the 1.1 request:

I get this from VERSION=1.1.0:

<Contents>
<CoverageSummary>
<ows:Title>bathymetry</ows:Title>
<ows:Abstract/>
<ows:WGS84BoundingBox dimensions="2">
<ows:LowerCorner>-180 -90</ows:LowerCorner>
<ows:UpperCorner>180 90</ows:UpperCorner>
</ows:WGS84BoundingBox>
<SupportedCRS>urn:ogc:def:crs:EPSG::4326</SupportedCRS>
<SupportedFormat>image/tiff</SupportedFormat>
<Identifier>bathymetry</Identifier>
</CoverageSummary>
</Contents>

whereas I get just this from VERSION=2.0.1:

<wcs:Contents>
<wcs:CoverageSummary>
<wcs:CoverageId>bathymetry</wcs:CoverageId>
<wcs:CoverageSubtype>RectifiedGridCoverage</wcs:CoverageSubtype>
</wcs:CoverageSummary>
</wcs:Contents>

Is there something I can add to my WCS map file to get the
"ows:WGS84BoundingBox"
element returned in my GetVersion?  I believe I was working with Geoserver
sources when I wrote that client but I would really like to get the
geographic extents.

Thank you very much for the great documentation!

carl
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211018/123495da/attachment.htm>

From jukka.rahkonen at maanmittauslaitos.fi  Mon Oct 18 11:49:54 2021
From: jukka.rahkonen at maanmittauslaitos.fi (Rahkonen Jukka (MML))
Date: Mon, 18 Oct 2021 18:49:54 +0000
Subject: [mapserver-users] Getting started with WCS on mapserver...
Message-ID: <d5f98bb2a5344e3c908f407a434cb759@maanmittauslaitos.fi>

Hi,

The wgs84-BoundingBox is an optional element in the CoverageSummary by the WCS 2.0.1 standard. There may be no way to configure Mapserver to add that information into GetCapabilities. You must read extents coverage by coverage with DescribeCoverage requests. But the boundedBy in that place is using the native CRS of the coverage like here:

<gml:boundedBy>
<gml:Envelope srsName="http://www.opengis.net/def/crs/EPSG/0/3067" axisLabels="x y" uomLabels="m m" srsDimension="2">
<gml:lowerCorner>212000 6762000</gml:lowerCorner>
<gml:upperCorner>218000 6768000</gml:upperCorner>
</gml:Envelope>
</gml:boundedBy>

The extent that appears in DescribeCoverage can be configured with the ?wcs_extent? metadata if it does not get generated automatically right. Having wgs84-BoundingBox in CoverageSummary could be useful for WCS clients but even the GetCapabilities example in the WCS Core 2.0.1 standard does not contain it.

-Jukka Rahkonen-

L?hett?j?: MapServer-users <mapserver-users-bounces at lists.osgeo.org> Puolesta Carl Godkin
L?hetetty: maanantai 18. lokakuuta 2021 21.24
Vastaanottaja: mapserver-users at lists.osgeo.org
Aihe: [mapserver-users] Getting started with WCS on mapserver...

Hi,

A week and a half ago I asked about using Mapserver +WMS to serve elevation data since I was following a tutorial that discussed WMS.

In one of the replies I received, Jeff suggested that I actually wanted WCS for my use case which was a forehead slap + d'oh! moment for me.  Of course I do!

Anyway, I have followed the Mapserver WCS docs and everything seems to work well.

However, I have one question that may be kind of niche or will perhaps betray my overall ignorance.

A while ago I wrote a WCS client for another project and it expects WCS version 2.0 requests.

The map file that I made up (based closely on the one in the docs: https://mapserver.org/ogc/wcs_server.html#example-wcs-server-mapfile ) seems to work fine for WCS 1.0, 1.1, and 2.0.  However, I don't seem to get a bounding box from GetCapabilities in the CoverageSummary from my 2.0 requests whereas I certainly do from the 1.1 request:

I get this from VERSION=1.1.0:

<Contents>
<CoverageSummary>
<ows:Title>bathymetry</ows:Title>
<ows:Abstract/>
<ows:WGS84BoundingBox dimensions="2">
<ows:LowerCorner>-180 -90</ows:LowerCorner>
<ows:UpperCorner>180 90</ows:UpperCorner>
</ows:WGS84BoundingBox>
<SupportedCRS>urn:ogc:def:crs:EPSG::4326</SupportedCRS>
<SupportedFormat>image/tiff</SupportedFormat>
<Identifier>bathymetry</Identifier>
</CoverageSummary>
</Contents>

whereas I get just this from VERSION=2.0.1:

<wcs:Contents>
<wcs:CoverageSummary>
<wcs:CoverageId>bathymetry</wcs:CoverageId>
<wcs:CoverageSubtype>RectifiedGridCoverage</wcs:CoverageSubtype>
</wcs:CoverageSummary>
</wcs:Contents>

Is there something I can add to my WCS map file to get the "ows:WGS84BoundingBox" element returned in my GetVersion?  I believe I was working with Geoserver sources when I wrote that client but I would really like to get the geographic extents.

Thank you very much for the great documentation!

carl

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211018/2901d45d/attachment.htm>

From cgodkin at gmail.com  Mon Oct 18 12:39:57 2021
From: cgodkin at gmail.com (Carl Godkin)
Date: Mon, 18 Oct 2021 12:39:57 -0700
Subject: [mapserver-users] Getting started with WCS on mapserver...
In-Reply-To: <d5f98bb2a5344e3c908f407a434cb759@maanmittauslaitos.fi>
References: <d5f98bb2a5344e3c908f407a434cb759@maanmittauslaitos.fi>
Message-ID: <CAE=8kmkKr0q-8RRg+DVZUKARkwxM2Jg5fw0soybJquX4+hQScA@mail.gmail.com>

Hi Jukka,

Thank you for the explanation.  I will see about using DescribeCoverage
instead.  I just tried it with my server and the gml:Envelope is indeed
filled in properly and I should be able to use it.

carl

On Mon, Oct 18, 2021 at 11:49 AM Rahkonen Jukka (MML) <
jukka.rahkonen at maanmittauslaitos.fi> wrote:

> Hi,
>
>
>
> The wgs84-BoundingBox is an optional element in the CoverageSummary by the
> WCS 2.0.1 standard. There may be no way to configure Mapserver to add that
> information into GetCapabilities. You must read extents coverage by
> coverage with DescribeCoverage requests. But the boundedBy in that place is
> using the native CRS of the coverage like here:
>
>
>
> <gml:boundedBy>
>
> <gml:Envelope srsName="http://www.opengis.net/def/crs/EPSG/0/3067"
> axisLabels="x y" uomLabels="m m" srsDimension="2">
>
> <gml:lowerCorner>212000 6762000</gml:lowerCorner>
>
> <gml:upperCorner>218000 6768000</gml:upperCorner>
>
> </gml:Envelope>
>
> </gml:boundedBy>
>
>
>
> The extent that appears in DescribeCoverage can be configured with the
> ?wcs_extent? metadata if it does not get generated automatically right.
> Having wgs84-BoundingBox in CoverageSummary could be useful for WCS clients
> but even the GetCapabilities example in the WCS Core 2.0.1 standard does
> not contain it.
>
>
>
> -Jukka Rahkonen-
>
>
>
> *L?hett?j?:* MapServer-users <mapserver-users-bounces at lists.osgeo.org> *Puolesta
> *Carl Godkin
> *L?hetetty:* maanantai 18. lokakuuta 2021 21.24
> *Vastaanottaja:* mapserver-users at lists.osgeo.org
> *Aihe:* [mapserver-users] Getting started with WCS on mapserver...
>
>
>
> Hi,
>
>
>
> A week and a half ago I asked about using Mapserver +WMS to serve
> elevation data since I was following a tutorial that discussed WMS.
>
>
>
> In one of the replies I received, Jeff suggested that I actually wanted
> WCS for my use case which was a forehead slap + d'oh! moment for me.  Of
> course I do!
>
>
>
> Anyway, I have followed the Mapserver WCS docs and everything seems to
> work well.
>
>
>
> However, I have one question that may be kind of niche or will perhaps
> betray my overall ignorance.
>
>
>
> A while ago I wrote a WCS client for another project and it expects WCS
> version 2.0 requests.
>
>
>
> The map file that I made up (based closely on the one in the docs:
> https://mapserver.org/ogc/wcs_server.html#example-wcs-server-mapfile )
> seems to work fine for WCS 1.0, 1.1, and 2.0.  However, I don't seem to get
> a bounding box from GetCapabilities in the CoverageSummary from my 2.0
> requests whereas I certainly do from the 1.1 request:
>
>
>
> I get this from VERSION=1.1.0:
>
>
>
> <Contents>
>
> <CoverageSummary>
>
> <ows:Title>bathymetry</ows:Title>
>
> <ows:Abstract/>
>
> <ows:WGS84BoundingBox dimensions="2">
>
> <ows:LowerCorner>-180 -90</ows:LowerCorner>
>
> <ows:UpperCorner>180 90</ows:UpperCorner>
>
> </ows:WGS84BoundingBox>
>
> <SupportedCRS>urn:ogc:def:crs:EPSG::4326</SupportedCRS>
>
> <SupportedFormat>image/tiff</SupportedFormat>
>
> <Identifier>bathymetry</Identifier>
>
> </CoverageSummary>
>
> </Contents>
>
>
>
> whereas I get just this from VERSION=2.0.1:
>
>
>
> <wcs:Contents>
>
> <wcs:CoverageSummary>
>
> <wcs:CoverageId>bathymetry</wcs:CoverageId>
>
> <wcs:CoverageSubtype>RectifiedGridCoverage</wcs:CoverageSubtype>
>
> </wcs:CoverageSummary>
>
> </wcs:Contents>
>
>
>
> Is there something I can add to my WCS map file to get the "ows:WGS84BoundingBox"
> element returned in my GetVersion?  I believe I was working with Geoserver
> sources when I wrote that client but I would really like to get the
> geographic extents.
>
>
>
> Thank you very much for the great documentation!
>
>
>
> carl
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211018/7cae2fb8/attachment.htm>

From cgodkin at gmail.com  Tue Oct 19 15:16:40 2021
From: cgodkin at gmail.com (Carl Godkin)
Date: Tue, 19 Oct 2021 15:16:40 -0700
Subject: [mapserver-users] WCS Tile index question
Message-ID: <CAE=8kmnnOgbzj4F2G4EfKpPLZ8takDH7RHOBwbSzKhA8tcyFkw@mail.gmail.com>

Hi again,

I am working with Mapserver trying to set up a WCS layer based on a large
collection of SRTM .hgt files.

Since I have many input files instead of a single file, I want to use a
tile index.  I know how to make and use a tile index for WMS, but there's a
twist for WCS that I don't understand.

According to this section of the docs:
 https://mapserver.org/optimization/tileindex.html?highlight=wcs_resolution#note-for-wcs-layers
<https://mapserver.org/optimization/tileindex.html?highlight=wcs_resolution#note-for-wcs-layers>

I need to set ?wcs_extent?, ?wcs_size?, and ?wcs_resolution? but I don't
see where they are defined or documented.

Can you point me to documentation?  The WCS Use Cases page gives some hints
so I think I understand wcs_extent (though it's "ows_extent" in some
examples there so I guess these are interchangeable?)

How are the others defined?

Thank you for your help,
carl
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211019/461c232a/attachment.htm>

From jukka.rahkonen at maanmittauslaitos.fi  Tue Oct 19 17:37:07 2021
From: jukka.rahkonen at maanmittauslaitos.fi (Rahkonen Jukka (MML))
Date: Wed, 20 Oct 2021 00:37:07 +0000
Subject: [mapserver-users] WCS Tile index question
Message-ID: <ca6bc3836c7741caa11076402322bab9@maanmittauslaitos.fi>

Hi,

Here is more accurate documentation but only in the text https://www.mapserver.org/ogc/wcs_server.html

?The convention is that once (wcs|ows)_extent and one of (wcs|ows)_size and (wcs|ows)_resolution is set in the layer metadata, all the coverage specific metadata will be retrieved from there. Otherwise the source image is queried via GDAL, if possible.?

It seems that only wcs_extent is documented in the list of metadata elements but this is what they do:

  *   wcs_extent defines the bounding box or your coverage -> that will go into DescribeCoverage
  *   wcs_size is the size of the coverage as pixels -> pixel size can be computed by extent and size
  *   wcs_resolution tells the pixel size explicitly -> size of the coverage in pixels can be computed

We have in the LAYER-METADATA:
"wcs_resolution" "0.5 0.5"
"wcs_extent"          "19970 6597278 749676 7785300"

That should be enough and Describe Coverage will have correct limits, origin, and offset vectors in RectifiedGrid.

-Jukka Rahkonen-


L?hett?j?: MapServer-users <mapserver-users-bounces at lists.osgeo.org> Puolesta Carl Godkin
L?hetetty: keskiviikko 20. lokakuuta 2021 1.17
Vastaanottaja: mapserver-users at lists.osgeo.org
Aihe: [mapserver-users] WCS Tile index question

Hi again,

I am working with Mapserver trying to set up a WCS layer based on a large collection of SRTM .hgt files.

Since I have many input files instead of a single file, I want to use a tile index.  I know how to make and use a tile index for WMS, but there's a twist for WCS that I don't understand.

According to this section of the docs:
 https://mapserver.org/optimization/tileindex.html?highlight=wcs_resolution#note-for-wcs-layers
<https://mapserver.org/optimization/tileindex.html?highlight=wcs_resolution#note-for-wcs-layers>

I need to set ?wcs_extent?, ?wcs_size?, and ?wcs_resolution? but I don't see where they are defined or documented.

Can you point me to documentation?  The WCS Use Cases page gives some hints so I think I understand wcs_extent (though it's "ows_extent" in some examples there so I guess these are interchangeable?)

How are the others defined?

Thank you for your help,
carl


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211020/130e2633/attachment.htm>

From cgodkin at gmail.com  Wed Oct 20 13:19:12 2021
From: cgodkin at gmail.com (Carl Godkin)
Date: Wed, 20 Oct 2021 13:19:12 -0700
Subject: [mapserver-users] WCS Tile index question
In-Reply-To: <ca6bc3836c7741caa11076402322bab9@maanmittauslaitos.fi>
References: <ca6bc3836c7741caa11076402322bab9@maanmittauslaitos.fi>
Message-ID: <CAE=8kmnVm8B3_2n=xL-JrHYnp7t0YiJTCfO7NYb4JuPFSVGzxg@mail.gmail.com>

On Tue, Oct 19, 2021 at 5:37 PM Rahkonen Jukka (MML) <
jukka.rahkonen at maanmittauslaitos.fi> wrote:

> Hi,
>
>
>
> Here is more accurate documentation but only in the text
> https://www.mapserver.org/ogc/wcs_server.html
>
>
>
> ?The convention is that once (wcs|ows)_extent and one of (wcs|ows)_size
> and (wcs|ows)_resolution is set in the layer metadata, all the coverage
> specific metadata will be retrieved from there. Otherwise the source image
> is queried via GDAL, if possible.?
>
>
>
> It seems that only wcs_extent is documented in the list of metadata
> elements but this is what they do:
>
>    - wcs_extent defines the bounding box or your coverage -> that will go
>    into DescribeCoverage
>    - wcs_size is the size of the coverage as pixels -> pixel size can be
>    computed by extent and size
>    - wcs_resolution tells the pixel size explicitly -> size of the
>    coverage in pixels can be computed
>
> Thanks a lot, Jukka.  That all makes some sense as far as it goes, but I'm
still not clear on a few details.  I have this range of data: 24N to 50N
and 125W to 66W which is 26 degrees by 59 degrees.

Since the SRTM data has 3 arc-second spacing, that's 1200 pixels per degree
plus one for the extra edge so I have tried a number of things without
success.

1. Is wcs_size the size FULL range of the tiles that I have?  And what is
the order?
2. Again, what is the order for wcs_resolution?  Are the units degrees
since the coordinate system of the original data is EPSG:4326?

I've tried all of these, one at a time, and get various errors.  (1 / 1200
= 0.00083333...)

        #"wcs_resolution" "0.00083333 0.00083333"
        #"wcs_size" "70801 31201"
        "wcs_size" "31201 70801"

In case it's helpful, here's my LAYER definition (with probably extra
things I've tried based on my searching and trial & error):

    LAYER
      NAME SRTM_3_arc-second_grid
      METADATA
        "wcs_label"           "SRTM - 3 arc-second grid"  ### required
        "wcs_rangeset_name"   "Range 1"  ### required to support
DescribeCoverage request
        "wcs_rangeset_label"  "Lower 48" ### required to support
DescribeCoverage request
        "ows_extent" "-125 24 -66 50"
        #"wcs_resolution" "0.00083333 0.00083333" # 1/1200 for 3" spacing?
        #"wcs_size" "70801 31201" # Not sure of order
        "wcs_size" "31201 70801"
        "wcs_bandcount" "1"
        "wcs_native_format" "SRTMHGT" # What should this be? I can't find
examples...
        "wcs_srs" "EPSG:4326"
      END

      TYPE RASTER ### required
      STATUS ON
      TILEINDEX "../srtm_3_hgt/srtm_3_hgt-index.shp" # Path is relative to
SHAPEPATH
      TILEITEM "Location"
      PROJECTION
        "init=epsg:4326"
      END
    END
I feel like I must be close. Thanks a lot for any further pointers,

carl
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211020/3fc1d351/attachment.htm>

From jukka.rahkonen at maanmittauslaitos.fi  Wed Oct 20 15:23:09 2021
From: jukka.rahkonen at maanmittauslaitos.fi (Rahkonen Jukka (MML))
Date: Wed, 20 Oct 2021 22:23:09 +0000
Subject: [mapserver-users] WCS Tile index question
Message-ID: <4bd9fd33781142dfbd37e14a2bcc179b@maanmittauslaitos.fi>



L?hett?j?: MapServer-users <mapserver-users-bounces at lists.osgeo.org> Puolesta Carl Godkin



>> On Tue, Oct 19, 2021 at 5:37 PM Rahkonen Jukka (MML) <mailto:jukka.rahkonen at maanmittauslaitos.fi> wrote:
>> Hi,
?
>> Here is more accurate documentation but only in the text 
>> https://www.mapserver.org/ogc/wcs_server.html
?
>> ?The convention is that once (wcs|ows)_extent and one of 
>> (wcs|ows)_size and (wcs|ows)_resolution is set in the layer metadata, all the coverage specific metadata will be retrieved from there. Otherwise the source image is queried via GDAL, if possible.?
?
>> It seems that only wcs_extent is documented in the list of metadata elements but this is what they do:
>> ? wcs_extent defines the bounding box or your coverage -> that will 
>> go into DescribeCoverage ? wcs_size is the size of the coverage as 
>> pixels -> pixel size can be computed by extent and size ? 
>> wcs_resolution tells the pixel size explicitly -> size of the 
>> coverage in pixels can be computed

> Thanks a lot, Jukka.? That all makes some sense as far as it goes, but 
> I'm still not clear on a few details. I have this range of data: 24N to 50N and 125W to 66W which is 26 degrees by 59 degrees.

>  Since the SRTM data has 3 arc-second spacing, that's 1200 pixels per 
> degree plus one for the extra edge so I have tried a number of things without success.

What fails? Don't you get anything into DescribeCoverage or is all that you get wrong?

> 1. Is wcs_size the size FULL range of the tiles that I have?? And what is the order?

You do not need to bother about the size if you set the resolution, but yes, it is the full range of the coverage. See this https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?service=WCS&REQUEST=DescribeCoverage&VERSION=2.0.1&COVERAGEID=ortokuva__ortokuva

The size is converted into the pixel space of the coverage gml:GridEnvelope>
<gml:low>0 0</gml:low>
<gml:high>1348356 2316499</gml:high>
</gml:GridEnvelope

I would say that the order is left-to-right - down-to-up but I fear that by the GML standard it can be something else. You can read about rectified grid from the OpenGIS Geography Markup Language (GML) Encoding Standard 3.2.1. I have never understood it but there are some images which may clarify this: "rectified grid: grid for which there is an affine transformation between the grid coordinates and the coordinates of an external coordinate reference system"

> 2. Again, what is the order for wcs_resolution?? Are the units degrees since the coordinate system of the original data is EPSG:4326?
With Mapserver the order is easting/longitude - northing/latitude

> I've tried all of these, one at a time, and get various errors.? (1 /
> 1200 = 0.00083333...)
What errors?

>? ? ? ? #"wcs_resolution" "0.00083333 0.00083333"
Looks right
>? ? ? ? #"wcs_size" "70801 31201"
Looks right if you prefer to give size instead of resolution
>? ? ? ? "wcs_size" "31201 70801"
This is wrong.

> In case it's helpful, here's my LAYER definition (with probably extra things I've tried based on my searching and trial & error):

? ? LAYER
? ? ? NAME SRTM_3_arc-second_grid
? ? ? METADATA
? ? ? ? "wcs_label" ? ? ? ? ? "SRTM - 3 arc-second grid" ?### required
? ? ? ? "wcs_rangeset_name" ? "Range 1" ?### required to support DescribeCoverage request
? ? ? ? "wcs_rangeset_label" ?"Lower 48" ### required to support DescribeCoverage request
? ? ? ? "ows_extent" "-125 24 -66 50"
? ? ? ? #"wcs_resolution" "0.00083333 0.00083333"    # 1/1200 for 3" spacing?
? ? ? ? #"wcs_size" "70801 31201"                    # Not sure of order
? ? ? ? "wcs_size" "31201 70801"
? ? ? ? "wcs_bandcount" "1"
? ? ? ? "wcs_native_format" "SRTMHGT"    # What should this be?  I can't find examples...

It is supposed to be the native format of the layer data, format that theoretically requires no processing, just selecting the right pixels. In practice the format that you users get by default with GetCoverage. It goes into DescribeCoverage about this way:
<wcs:nativeFormat>image/tiff</wcs:nativeFormat> 
and it is documented in the WCS 2.0.1 Core standard:
"4.7 Native Format
encoding format where, in a GetCoverage request, the range set values can be obtained unaltered"

? ? ? ? "wcs_srs" "EPSG:4326"
? ? ? END

? ? ? TYPE RASTER ### required
? ? ? STATUS ON
? ? ? TILEINDEX "../srtm_3_hgt/srtm_3_hgt-index.shp" # Path is relative to SHAPEPATH
? ? ? TILEITEM "Location"
? ? ? PROJECTION
? ? ? ? "init=epsg:4326"
? ? ? END
? ? END

> I feel like I must be close.  Thanks a lot for any further pointers,

-Jukka-

> carl


From cgodkin at gmail.com  Wed Oct 20 15:46:05 2021
From: cgodkin at gmail.com (Carl Godkin)
Date: Wed, 20 Oct 2021 15:46:05 -0700
Subject: [mapserver-users] WCS Tile index question
In-Reply-To: <4bd9fd33781142dfbd37e14a2bcc179b@maanmittauslaitos.fi>
References: <4bd9fd33781142dfbd37e14a2bcc179b@maanmittauslaitos.fi>
Message-ID: <CAE=8kmnbTfEoMDtcx+8uAubkOmDiWniERR07DOkzRtMzXDnHsA@mail.gmail.com>

Hi again Jukka,

Thanks for the speedy response.  I now think that I have the resolution and
extent correct with these.  (I was getting various errors trying other
things but this
set seems okay for GetCoverage, GetCapabilities and DescribeCoverage):

        "ows_extent" "-125 24 -66 50"
        "wcs_resolution" "0.0008333333 0.0008333333"

I actually retrieve a correctly located GeoTIFF (yay!!) out with this query
for instance:

http://localhost:8080/wcs?service=WCS&version=2.0.1& request=GetCoverage&
coverageid=SRTM_3_arc-second_grid& format=geotiff_16&
subset=Lat(43.256274,43.721318)& subset=Long(-112.357549,-111.719135)

but it's all zeros.

My data files are all native SRTM ".hgt" files. (And yes, I did "gdalinfo
-stats" on one tile for the above query and there are good values.)

I could gdal_translate them all to GeoTIFF, but since GDAL supports HGT
files I hoped I could use them as is.

I have tried everything I can think of for the "wcs_native_format" metadata:

      "wcs_native_format" "application/srtmhgt"

and application/hgt, image/hgt, srtmhgt, etc. I don't see errors in
mapserver log though no matter what I try. Is this my problem?

What do you suggest?

Thanks again,
carl









On Wed, Oct 20, 2021 at 3:23 PM Rahkonen Jukka (MML) <
jukka.rahkonen at maanmittauslaitos.fi> wrote:

>
>
> L?hett?j?: MapServer-users <mapserver-users-bounces at lists.osgeo.org>
> Puolesta Carl Godkin
>
>
>
> >> On Tue, Oct 19, 2021 at 5:37 PM Rahkonen Jukka (MML) <mailto:
> jukka.rahkonen at maanmittauslaitos.fi> wrote:
> >> Hi,
>
> >> Here is more accurate documentation but only in the text
> >> https://www.mapserver.org/ogc/wcs_server.html
>
> >> ?The convention is that once (wcs|ows)_extent and one of
> >> (wcs|ows)_size and (wcs|ows)_resolution is set in the layer metadata,
> all the coverage specific metadata will be retrieved from there. Otherwise
> the source image is queried via GDAL, if possible.?
>
> >> It seems that only wcs_extent is documented in the list of metadata
> elements but this is what they do:
> >> ? wcs_extent defines the bounding box or your coverage -> that will
> >> go into DescribeCoverage ? wcs_size is the size of the coverage as
> >> pixels -> pixel size can be computed by extent and size ?
> >> wcs_resolution tells the pixel size explicitly -> size of the
> >> coverage in pixels can be computed
>
> > Thanks a lot, Jukka.  That all makes some sense as far as it goes, but
> > I'm still not clear on a few details. I have this range of data: 24N to
> 50N and 125W to 66W which is 26 degrees by 59 degrees.
>
> >  Since the SRTM data has 3 arc-second spacing, that's 1200 pixels per
> > degree plus one for the extra edge so I have tried a number of things
> without success.
>
> What fails? Don't you get anything into DescribeCoverage or is all that
> you get wrong?
>
> > 1. Is wcs_size the size FULL range of the tiles that I have?  And what
> is the order?
>
> You do not need to bother about the size if you set the resolution, but
> yes, it is the full range of the coverage. See this
> https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?service=WCS&REQUEST=DescribeCoverage&VERSION=2.0.1&COVERAGEID=ortokuva__ortokuva
>
> The size is converted into the pixel space of the coverage
> gml:GridEnvelope>
> <gml:low>0 0</gml:low>
> <gml:high>1348356 2316499</gml:high>
> </gml:GridEnvelope
>
> I would say that the order is left-to-right - down-to-up but I fear that
> by the GML standard it can be something else. You can read about rectified
> grid from the OpenGIS Geography Markup Language (GML) Encoding Standard
> 3.2.1. I have never understood it but there are some images which may
> clarify this: "rectified grid: grid for which there is an affine
> transformation between the grid coordinates and the coordinates of an
> external coordinate reference system"
>
> > 2. Again, what is the order for wcs_resolution?  Are the units degrees
> since the coordinate system of the original data is EPSG:4326?
> With Mapserver the order is easting/longitude - northing/latitude
>
> > I've tried all of these, one at a time, and get various errors.  (1 /
> > 1200 = 0.00083333...)
> What errors?
>
> >        #"wcs_resolution" "0.00083333 0.00083333"
> Looks right
> >        #"wcs_size" "70801 31201"
> Looks right if you prefer to give size instead of resolution
> >        "wcs_size" "31201 70801"
> This is wrong.
>
> > In case it's helpful, here's my LAYER definition (with probably extra
> things I've tried based on my searching and trial & error):
>
>     LAYER
>       NAME SRTM_3_arc-second_grid
>       METADATA
>         "wcs_label"           "SRTM - 3 arc-second grid"  ### required
>         "wcs_rangeset_name"   "Range 1"  ### required to support
> DescribeCoverage request
>         "wcs_rangeset_label"  "Lower 48" ### required to support
> DescribeCoverage request
>         "ows_extent" "-125 24 -66 50"
>         #"wcs_resolution" "0.00083333 0.00083333"    # 1/1200 for 3"
> spacing?
>         #"wcs_size" "70801 31201"                    # Not sure of order
>         "wcs_size" "31201 70801"
>         "wcs_bandcount" "1"
>         "wcs_native_format" "SRTMHGT"    # What should this be?  I can't
> find examples...
>
> It is supposed to be the native format of the layer data, format that
> theoretically requires no processing, just selecting the right pixels. In
> practice the format that you users get by default with GetCoverage. It goes
> into DescribeCoverage about this way:
> <wcs:nativeFormat>image/tiff</wcs:nativeFormat>
> and it is documented in the WCS 2.0.1 Core standard:
> "4.7 Native Format
> encoding format where, in a GetCoverage request, the range set values can
> be obtained unaltered"
>
>         "wcs_srs" "EPSG:4326"
>       END
>
>       TYPE RASTER ### required
>       STATUS ON
>       TILEINDEX "../srtm_3_hgt/srtm_3_hgt-index.shp" # Path is relative to
> SHAPEPATH
>       TILEITEM "Location"
>       PROJECTION
>         "init=epsg:4326"
>       END
>     END
>
> > I feel like I must be close.  Thanks a lot for any further pointers,
>
> -Jukka-
>
> > carl
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211020/e175c59a/attachment.htm>

From jukka.rahkonen at maanmittauslaitos.fi  Wed Oct 20 15:55:06 2021
From: jukka.rahkonen at maanmittauslaitos.fi (Rahkonen Jukka (MML))
Date: Wed, 20 Oct 2021 22:55:06 +0000
Subject: [mapserver-users] WCS Tile index question
Message-ID: <f64373ad8ad643308c93bb72ee31a134@maanmittauslaitos.fi>

Hi,

I was just about to suggest to use DEBUG. GDAL does not use ?wcs_native_format?, it is used for setting what will be written into DescribeCoverage.

-Jukka Rahkonen-

L?hett?j?: Carl Godkin <cgodkin at gmail.com>
L?hetetty: torstai 21. lokakuuta 2021 1.50
Vastaanottaja: Rahkonen Jukka (MML) <jukka.rahkonen at maanmittauslaitos.fi>
Aihe: Re: [mapserver-users] WCS Tile index question

Sorry, right after I sent that last reply to the list, I thought of turning up the DEBUG value in the mapfile to 5 and now I see it's not finding my .hgt files for some reason.

I will dig into this; maybe mapserver + GDAL don't actually needs wcs_native_format.

More later to the list if I sort this out.  I really appreciate the help you and Jeff and others on the list have given me!

carl

On Tue, Oct 19, 2021 at 5:37 PM Rahkonen Jukka (MML) <jukka.rahkonen at maanmittauslaitos.fi<mailto:jukka.rahkonen at maanmittauslaitos.fi>> wrote:
Hi,

Here is more accurate documentation but only in the text https://www.mapserver.org/ogc/wcs_server.html

?The convention is that once (wcs|ows)_extent and one of (wcs|ows)_size and (wcs|ows)_resolution is set in the layer metadata, all the coverage specific metadata will be retrieved from there. Otherwise the source image is queried via GDAL, if possible.?

It seems that only wcs_extent is documented in the list of metadata elements but this is what they do:

  *   wcs_extent defines the bounding box or your coverage -> that will go into DescribeCoverage
  *   wcs_size is the size of the coverage as pixels -> pixel size can be computed by extent and size
  *   wcs_resolution tells the pixel size explicitly -> size of the coverage in pixels can be computed

We have in the LAYER-METADATA:
"wcs_resolution" "0.5 0.5"
"wcs_extent"          "19970 6597278 749676 7785300"

That should be enough and Describe Coverage will have correct limits, origin, and offset vectors in RectifiedGrid.

-Jukka Rahkonen-


L?hett?j?: MapServer-users <mapserver-users-bounces at lists.osgeo.org<mailto:mapserver-users-bounces at lists.osgeo.org>> Puolesta Carl Godkin
L?hetetty: keskiviikko 20. lokakuuta 2021 1.17
Vastaanottaja: mapserver-users at lists.osgeo.org<mailto:mapserver-users at lists.osgeo.org>
Aihe: [mapserver-users] WCS Tile index question

Hi again,

I am working with Mapserver trying to set up a WCS layer based on a large collection of SRTM .hgt files.

Since I have many input files instead of a single file, I want to use a tile index.  I know how to make and use a tile index for WMS, but there's a twist for WCS that I don't understand.

According to this section of the docs:
 https://mapserver.org/optimization/tileindex.html?highlight=wcs_resolution#note-for-wcs-layers
<https://mapserver.org/optimization/tileindex.html?highlight=wcs_resolution#note-for-wcs-layers>

I need to set ?wcs_extent?, ?wcs_size?, and ?wcs_resolution? but I don't see where they are defined or documented.

Can you point me to documentation?  The WCS Use Cases page gives some hints so I think I understand wcs_extent (though it's "ows_extent" in some examples there so I guess these are interchangeable?)

How are the others defined?

Thank you for your help,
carl


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211020/a13bd1e4/attachment.htm>

From cgodkin at gmail.com  Wed Oct 20 16:04:14 2021
From: cgodkin at gmail.com (Carl Godkin)
Date: Wed, 20 Oct 2021 16:04:14 -0700
Subject: [mapserver-users] WCS Tile index question
In-Reply-To: <f64373ad8ad643308c93bb72ee31a134@maanmittauslaitos.fi>
References: <f64373ad8ad643308c93bb72ee31a134@maanmittauslaitos.fi>
Message-ID: <CAE=8kmm0t-8BhjFFtmNd6sDqFvy8zF8=+HT168mV44SrSKsg9Q@mail.gmail.com>

Thanks again; I had messed up the gdaltindex paths so (as the debug log
told me), mapserver couldn't find the .hgt files.

I just fixed that and I'm thrilled to report I've got a WCS serving terrain
now.

Thanks so much!!
carl

On Wed, Oct 20, 2021 at 3:55 PM Rahkonen Jukka (MML) <
jukka.rahkonen at maanmittauslaitos.fi> wrote:

> Hi,
>
>
>
> I was just about to suggest to use DEBUG. GDAL does not use ?wcs_native_format?,
> it is used for setting what will be written into DescribeCoverage.
>
>
>
> -Jukka Rahkonen-
>
>
>
> *L?hett?j?:* Carl Godkin <cgodkin at gmail.com>
> *L?hetetty:* torstai 21. lokakuuta 2021 1.50
> *Vastaanottaja:* Rahkonen Jukka (MML) <jukka.rahkonen at maanmittauslaitos.fi
> >
> *Aihe:* Re: [mapserver-users] WCS Tile index question
>
>
>
> Sorry, right after I sent that last reply to the list, I thought of
> turning up the DEBUG value in the mapfile to 5 and now I see it's not
> finding my .hgt files for some reason.
>
>
>
> I will dig into this; maybe mapserver + GDAL don't actually needs
> wcs_native_format.
>
>
>
> More later to the list if I sort this out.  I really appreciate the help
> you and Jeff and others on the list have given me!
>
>
>
> carl
>
>
>
> On Tue, Oct 19, 2021 at 5:37 PM Rahkonen Jukka (MML) <
> jukka.rahkonen at maanmittauslaitos.fi> wrote:
>
> Hi,
>
>
>
> Here is more accurate documentation but only in the text
> https://www.mapserver.org/ogc/wcs_server.html
>
>
>
> ?The convention is that once (wcs|ows)_extent and one of (wcs|ows)_size
> and (wcs|ows)_resolution is set in the layer metadata, all the coverage
> specific metadata will be retrieved from there. Otherwise the source image
> is queried via GDAL, if possible.?
>
>
>
> It seems that only wcs_extent is documented in the list of metadata
> elements but this is what they do:
>
>    - wcs_extent defines the bounding box or your coverage -> that will go
>    into DescribeCoverage
>    - wcs_size is the size of the coverage as pixels -> pixel size can be
>    computed by extent and size
>    - wcs_resolution tells the pixel size explicitly -> size of the
>    coverage in pixels can be computed
>
>
>
> We have in the LAYER-METADATA:
>
> "wcs_resolution" "0.5 0.5"
>
> "wcs_extent"          "19970 6597278 749676 7785300"
>
>
>
> That should be enough and Describe Coverage will have correct limits,
> origin, and offset vectors in RectifiedGrid.
>
>
>
> -Jukka Rahkonen-
>
>
>
>
>
> *L?hett?j?:* MapServer-users <mapserver-users-bounces at lists.osgeo.org> *Puolesta
> *Carl Godkin
> *L?hetetty:* keskiviikko 20. lokakuuta 2021 1.17
> *Vastaanottaja:* mapserver-users at lists.osgeo.org
> *Aihe:* [mapserver-users] WCS Tile index question
>
>
>
> Hi again,
>
>
>
> I am working with Mapserver trying to set up a WCS layer based on a large
> collection of SRTM .hgt files.
>
>
>
> Since I have many input files instead of a single file, I want to use a
> tile index.  I know how to make and use a tile index for WMS, but there's a
> twist for WCS that I don't understand.
>
>
>
> According to this section of the docs:
>
>
>  https://mapserver.org/optimization/tileindex.html?highlight=wcs_resolution#note-for-wcs-layers
>
> <https://mapserver.org/optimization/tileindex.html?highlight=wcs_resolution#note-for-wcs-layers>
>
>
>
> I need to set ?wcs_extent?, ?wcs_size?, and ?wcs_resolution? but I don't
> see where they are defined or documented.
>
>
>
> Can you point me to documentation?  The WCS Use Cases page gives some
> hints so I think I understand wcs_extent (though it's "ows_extent" in some
> examples there so I guess these are interchangeable?)
>
>
>
> How are the others defined?
>
>
>
> Thank you for your help,
>
> carl
>
>
>
>
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211020/5ca5819c/attachment.htm>

From MarkVolz at co.lyon.mn.us  Fri Oct 22 09:14:32 2021
From: MarkVolz at co.lyon.mn.us (Mark Volz)
Date: Fri, 22 Oct 2021 16:14:32 +0000
Subject: [mapserver-users] creating a Mapcache that has the same tile schema
 as ArcGIS vector tiles
Message-ID: <be25ab7efe6c4425906d44444f762bfc@co.lyon.mn.us>

Hello,

I would to create an ArcGIS Server cache that uses an equivalent schema as MapServer Mapcache.  I would like to know if anyone has a tips to make sure that the schema will work between the products.  I am specifically concerned about the tile origin.  So here is what I did:


*         I am using the same projection on both caches.

*         I set the pixel size in ArcGIS to use the same resolutions as MapCache.

*         I am using the same tile width and size.

*         I am using the same units.

*         Of particular concern ArcGIS wants to put the tile origin at -119278400, 141351400 by default while Mapcache uses 436450, 80000.

Anyways, I just wanted to see if I could make the tile cache more compatible between the two systems.  Let me know if you have any tips.

Thanks!

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211022/85deb865/attachment.htm>

From cgodkin at gmail.com  Mon Oct 25 14:54:49 2021
From: cgodkin at gmail.com (Carl Godkin)
Date: Mon, 25 Oct 2021 14:54:49 -0700
Subject: [mapserver-users] Getting started with mapcache for WMTS...
Message-ID: <CAE=8kmnspDr+VOGQaHYHn5H3-9Opd5DB0mZmbBLOatjO5PBviw@mail.gmail.com>

Hi again,

Is there a way to search through the mapserver-users mailing list archive
except month-by-month?  If so, I can't find it.  (I usually use my Gmail
archives for GDAL and PROJ since I've subscribed forever to those but I'm a
mapserver newbie.)  Anyway...

I am trying to configure mapcache to serve from my (local and AFAICT fully
working) WMS server on my MS4W 4.0.5 installation.  I am following the
directions at https://mapserver.org/mapcache/index.html as well as the
steps listed from the local MS4W page (#mapcache-apache-module).

I think I understand the basic principles and have added a new source that
looks like this:

   <source name="faa" type="wms">
      <getmap>
         <params>
            <FORMAT>image/geotiff</FORMAT>
            <LAYERS>FAA_Sectional_Charts</LAYERS>
            <TRANSPARENT>TRUE</TRANSPARENT>
            <SRS>EPSG:4326</SRS>
            <BBOX>-125,24,-61,49</BBOX>
         </params>
      </getmap>
      <http>
         <url>http://localhost:8080/wms</url>
      </http>
   </source>

and a new tileset that looks like this:

   <tileset name="faa">
      <source>faa</source>
      <cache>disk</cache>
      <grid>WGS84</grid>
      <format>PNG</format>
      <metatile>5 5</metatile>
      <metabuffer>10</metabuffer>
      <expires>3600</expires>
   </tileset>
(This is the same as the "test" tileset with a different source.)

But how do I test this? The query

http://localhost:8080/mapcache/wmts?service=wmts&request=getcapabilities&version=1.0.0

lists my new source as a layer.  The docs say

   1.

   Test your MapCache configuration by following steps in the
associated MapCache
   testing docs
   <https://mapserver.org/mapcache/install.html#test-your-mapcache-module>

but that seems to be about testing the mapcache/demo.  I don't see how to
add something to the mapcache/demo if that's possible.  What am I
overlooking?

Thanks a lot,
carl
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211025/e300f448/attachment.htm>

From public at postholer.com  Mon Oct 25 15:24:50 2021
From: public at postholer.com (Scott)
Date: Mon, 25 Oct 2021 15:24:50 -0700
Subject: [mapserver-users] Getting started with mapcache for WMTS...
In-Reply-To: <CAE=8kmnspDr+VOGQaHYHn5H3-9Opd5DB0mZmbBLOatjO5PBviw@mail.gmail.com>
References: <CAE=8kmnspDr+VOGQaHYHn5H3-9Opd5DB0mZmbBLOatjO5PBviw@mail.gmail.com>
Message-ID: <ca447442-bb09-4ff6-0bbd-130bd89fc697@postholer.com>

Try using mapcache seeder. That will create some cached images to work 
with. It will also test your mapcache config:

https://mapserver.org/mapcache/seed.html

How you use those tiles is a whole different subject.



On 10/25/21 2:54 PM, Carl Godkin wrote:
> Hi again,
> 
> Is there a way to search through the mapserver-users mailing list 
> archive except month-by-month?? If so, I can't find it.? (I usually use 
> my Gmail archives for GDAL and PROJ since I've subscribed forever to 
> those but I'm a mapserver newbie.)? Anyway...
> 
> I am trying to configure mapcache to serve from my (local and AFAICT 
> fully working) WMS server on?my MS4W 4.0.5 installation.? I am following 
> the? directions at https://mapserver.org/mapcache/index.html 
> <https://mapserver.org/mapcache/index.html> as well as the steps listed 
> from the local MS4W page?(#mapcache-apache-module).
> 
> I think I understand the basic principles and have added a new source 
> that looks like this:
> 
> <source name="faa" type="wms">
> <getmap>
> <params>
> <FORMAT>image/geotiff</FORMAT>
> <LAYERS>FAA_Sectional_Charts</LAYERS>
> <TRANSPARENT>TRUE</TRANSPARENT>
> <SRS>EPSG:4326</SRS>
> <BBOX>-125,24,-61,49</BBOX>
> </params>
> </getmap>
> <http>
> <url>http://localhost:8080/wms <http://localhost:8080/wms></url>
> </http>
> </source>
> 
> and a new tileset that looks like this:
> 
> <tileset name="faa">
> <source>faa</source>
> <cache>disk</cache>
> <grid>WGS84</grid>
> <format>PNG</format>
> <metatile>5 5</metatile>
> <metabuffer>10</metabuffer>
> <expires>3600</expires>
> </tileset>
> (This is the same as the "test" tileset with a different source.)
> 
> But how do I test this? The query
> 
> http://localhost:8080/mapcache/wmts?service=wmts&request=getcapabilities&version=1.0.0 
> <http://localhost:8080/mapcache/wmts?service=wmts&request=getcapabilities&version=1.0.0>
> 
> lists my new source as a layer. The docs say
> 
>  5.
> 
>     Test your MapCache configuration by following steps in the
>     associated MapCache testing docs
>     <https://mapserver.org/mapcache/install.html#test-your-mapcache-module>
> 
> but that seems to be about testing the mapcache/demo.? I don't see how 
> to add something to the mapcache/demo if that's possible.? What am I 
> overlooking?
> 
> Thanks a lot,
> carl
> 
> 
> 
> 
> 
> 
> _______________________________________________
> MapServer-users mailing list
> MapServer-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/mapserver-users
> 


From cgodkin at gmail.com  Mon Oct 25 16:46:49 2021
From: cgodkin at gmail.com (Carl Godkin)
Date: Mon, 25 Oct 2021 16:46:49 -0700
Subject: [mapserver-users] Getting started with mapcache for WMTS...
In-Reply-To: <ca447442-bb09-4ff6-0bbd-130bd89fc697@postholer.com>
References: <CAE=8kmnspDr+VOGQaHYHn5H3-9Opd5DB0mZmbBLOatjO5PBviw@mail.gmail.com>
 <ca447442-bb09-4ff6-0bbd-130bd89fc697@postholer.com>
Message-ID: <CAE=8kmmxyjf61rzQSrEWErPxnuyu_j2_jPQhY6dThtQ3c9_9DQ@mail.gmail.com>

Thanks for the suggestion, Scott.  I think the seeder is optional, right?

Anyway, I ended up figuring out that my <source> was wrong in a couple of
ways in mapcache.xml.  I misunderstood the
<getmap>  parameters.  Now I have this:

   <source name="faa" type="wms">
         <getmap>
         <params>
            <FORMAT>image/png</FORMAT>
            <LAYERS>FAA_Sectional_Charts</LAYERS>
            <SRS>EPSG:3857</SRS>
         </params>
      </getmap>
      <http>
         <url>http://localhost:8080/wms</url>
      </http>
   </source>

I wrote a tiny little leaflet test and, after making every possible
mistake, finally
got it all to work so I'm all set. Thanks,

carl


On Mon, Oct 25, 2021 at 3:24 PM Scott <public at postholer.com> wrote:

> Try using mapcache seeder. That will create some cached images to work
> with. It will also test your mapcache config:
>
> https://mapserver.org/mapcache/seed.html
>
> How you use those tiles is a whole different subject.
>
>
>
> On 10/25/21 2:54 PM, Carl Godkin wrote:
> > Hi again,
> >
> > Is there a way to search through the mapserver-users mailing list
> > archive except month-by-month?  If so, I can't find it.  (I usually use
> > my Gmail archives for GDAL and PROJ since I've subscribed forever to
> > those but I'm a mapserver newbie.)  Anyway...
> >
> > I am trying to configure mapcache to serve from my (local and AFAICT
> > fully working) WMS server on my MS4W 4.0.5 installation.  I am following
> > the  directions at https://mapserver.org/mapcache/index.html
> > <https://mapserver.org/mapcache/index.html> as well as the steps listed
> > from the local MS4W page (#mapcache-apache-module).
> >
> > I think I understand the basic principles and have added a new source
> > that looks like this:
> >
> > <source name="faa" type="wms">
> > <getmap>
> > <params>
> > <FORMAT>image/geotiff</FORMAT>
> > <LAYERS>FAA_Sectional_Charts</LAYERS>
> > <TRANSPARENT>TRUE</TRANSPARENT>
> > <SRS>EPSG:4326</SRS>
> > <BBOX>-125,24,-61,49</BBOX>
> > </params>
> > </getmap>
> > <http>
> > <url>http://localhost:8080/wms <http://localhost:8080/wms></url>
> > </http>
> > </source>
> >
> > and a new tileset that looks like this:
> >
> > <tileset name="faa">
> > <source>faa</source>
> > <cache>disk</cache>
> > <grid>WGS84</grid>
> > <format>PNG</format>
> > <metatile>5 5</metatile>
> > <metabuffer>10</metabuffer>
> > <expires>3600</expires>
> > </tileset>
> > (This is the same as the "test" tileset with a different source.)
> >
> > But how do I test this? The query
> >
> >
> http://localhost:8080/mapcache/wmts?service=wmts&request=getcapabilities&version=1.0.0
> > <
> http://localhost:8080/mapcache/wmts?service=wmts&request=getcapabilities&version=1.0.0
> >
> >
> > lists my new source as a layer. The docs say
> >
> >  5.
> >
> >     Test your MapCache configuration by following steps in the
> >     associated MapCache testing docs
> >     <
> https://mapserver.org/mapcache/install.html#test-your-mapcache-module>
> >
> > but that seems to be about testing the mapcache/demo.  I don't see how
> > to add something to the mapcache/demo if that's possible.  What am I
> > overlooking?
> >
> > Thanks a lot,
> > carl
> >
> >
> >
> >
> >
> >
> > _______________________________________________
> > MapServer-users mailing list
> > MapServer-users at lists.osgeo.org
> > https://lists.osgeo.org/mailman/listinfo/mapserver-users
> >
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211025/8704b1ac/attachment.htm>

From public at postholer.com  Mon Oct 25 18:05:11 2021
From: public at postholer.com (Scott)
Date: Mon, 25 Oct 2021 18:05:11 -0700
Subject: [mapserver-users] Getting started with mapcache for WMTS...
In-Reply-To: <CAE=8kmmxyjf61rzQSrEWErPxnuyu_j2_jPQhY6dThtQ3c9_9DQ@mail.gmail.com>
References: <CAE=8kmnspDr+VOGQaHYHn5H3-9Opd5DB0mZmbBLOatjO5PBviw@mail.gmail.com>
 <ca447442-bb09-4ff6-0bbd-130bd89fc697@postholer.com>
 <CAE=8kmmxyjf61rzQSrEWErPxnuyu_j2_jPQhY6dThtQ3c9_9DQ@mail.gmail.com>
Message-ID: <2ffc19b8-215f-b66a-74a4-638bfed39815@postholer.com>

Yeah, the seeder is optional. If you have tiles that won't update often 
it might be handy to seed your entire tileset. Your server will love you.

Here's another take on mapcache config. This is how I set it up.
Note the absence of the SRS tag and the use of grid, format, cache tags. 
All 3 are defined elsewhere, not in the tileset itself, which means you 
can reuse them in any new tileset. See docs for more info:

    <source name="snowlevel2day-source" type="wms">
       <http>
 
<url>https://somedomain.com/cgi-bin/mapserverBin?map=/path/to/mapfile/mainmap.map</url>
          <connection_timeout>10</connection_timeout>
          <timeout>12</timeout>
       </http>
       <getmap>
          <params>
             <TRANSPARENT>true</TRANSPARENT>
             <FORMAT>image/png</FORMAT>
 
<LAYERS>ndfdsnowlevel,countylines,statelines,ndfdpttemp</LAYERS>
             <ndfdsnowlevelband>33</ndfdsnowlevelband>
             <ndfdpttempband>2</ndfdpttempband>
          </params>
       </getmap>
    </source>

    <tileset name="snowlevel2day">
       <metadata>
          <title>Day 2 Snow Levels</title>
       </metadata>
       <source>snowlevel2day-source</source>
       <cache>disk</cache>
       <format>mixed</format>
       <grid>google-900913</grid>
       <expires>14400</expires>
       <auto_expire>14400</auto_expire>
    </tileset>


On 10/25/21 4:46 PM, Carl Godkin wrote:
> Thanks for the suggestion, Scott.? I think the seeder is optional, right?
> 
> Anyway, I ended up figuring?out that my <source> was wrong in a couple 
> of ways in mapcache.xml.? I misunderstood the
> <getmap>? parameters.? Now I have this:
> 
> <source name="faa" type="wms">
> <getmap>
> <params>
> <FORMAT>image/png</FORMAT>
> <LAYERS>FAA_Sectional_Charts</LAYERS>
> <SRS>EPSG:3857</SRS>
> </params>
> </getmap>
> <http>
> <url>http://localhost:8080/wms <http://localhost:8080/wms></url>
> </http>
> </source>
> 
> I wrote a tiny little leaflet test and, after making every possible 
> mistake, finally
> got it all to work so I'm all set. Thanks,
> 
> carl
> 
> 
> On Mon, Oct 25, 2021 at 3:24 PM Scott <public at postholer.com 
> <mailto:public at postholer.com>> wrote:
> 
>     Try using mapcache seeder. That will create some cached images to work
>     with. It will also test your mapcache config:
> 
>     https://mapserver.org/mapcache/seed.html
>     <https://mapserver.org/mapcache/seed.html>
> 
>     How you use those tiles is a whole different subject.
> 
> 
> 
>     On 10/25/21 2:54 PM, Carl Godkin wrote:
>      > Hi again,
>      >
>      > Is there a way to search through the mapserver-users mailing list
>      > archive except month-by-month?? If so, I can't find it.? (I
>     usually use
>      > my Gmail archives for GDAL and PROJ since I've subscribed forever to
>      > those but I'm a mapserver newbie.)? Anyway...
>      >
>      > I am trying to configure mapcache to serve from my (local and AFAICT
>      > fully working) WMS server on?my MS4W 4.0.5 installation.? I am
>     following
>      > the? directions at https://mapserver.org/mapcache/index.html
>     <https://mapserver.org/mapcache/index.html>
>      > <https://mapserver.org/mapcache/index.html
>     <https://mapserver.org/mapcache/index.html>> as well as the steps
>     listed
>      > from the local MS4W page?(#mapcache-apache-module).
>      >
>      > I think I understand the basic principles and have added a new
>     source
>      > that looks like this:
>      >
>      > <source name="faa" type="wms">
>      > <getmap>
>      > <params>
>      > <FORMAT>image/geotiff</FORMAT>
>      > <LAYERS>FAA_Sectional_Charts</LAYERS>
>      > <TRANSPARENT>TRUE</TRANSPARENT>
>      > <SRS>EPSG:4326</SRS>
>      > <BBOX>-125,24,-61,49</BBOX>
>      > </params>
>      > </getmap>
>      > <http>
>      > <url>http://localhost:8080/wms <http://localhost:8080/wms>
>     <http://localhost:8080/wms <http://localhost:8080/wms>></url>
>      > </http>
>      > </source>
>      >
>      > and a new tileset that looks like this:
>      >
>      > <tileset name="faa">
>      > <source>faa</source>
>      > <cache>disk</cache>
>      > <grid>WGS84</grid>
>      > <format>PNG</format>
>      > <metatile>5 5</metatile>
>      > <metabuffer>10</metabuffer>
>      > <expires>3600</expires>
>      > </tileset>
>      > (This is the same as the "test" tileset with a different source.)
>      >
>      > But how do I test this? The query
>      >
>      >
>     http://localhost:8080/mapcache/wmts?service=wmts&request=getcapabilities&version=1.0.0
>     <http://localhost:8080/mapcache/wmts?service=wmts&request=getcapabilities&version=1.0.0>
> 
>      >
>     <http://localhost:8080/mapcache/wmts?service=wmts&request=getcapabilities&version=1.0.0
>     <http://localhost:8080/mapcache/wmts?service=wmts&request=getcapabilities&version=1.0.0>>
>      >
>      > lists my new source as a layer. The docs say
>      >
>      >? 5.
>      >
>      >? ? ?Test your MapCache configuration by following steps in the
>      >? ? ?associated MapCache testing docs
>      >   
>      ?<https://mapserver.org/mapcache/install.html#test-your-mapcache-module <https://mapserver.org/mapcache/install.html#test-your-mapcache-module>>
>      >
>      > but that seems to be about testing the mapcache/demo.? I don't
>     see how
>      > to add something to the mapcache/demo if that's possible.? What am I
>      > overlooking?
>      >
>      > Thanks a lot,
>      > carl
>      >
>      >
>      >
>      >
>      >
>      >
>      > _______________________________________________
>      > MapServer-users mailing list
>      > MapServer-users at lists.osgeo.org
>     <mailto:MapServer-users at lists.osgeo.org>
>      > https://lists.osgeo.org/mailman/listinfo/mapserver-users
>     <https://lists.osgeo.org/mailman/listinfo/mapserver-users>
>      >
> 
> 
> _______________________________________________
> MapServer-users mailing list
> MapServer-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/mapserver-users
> 


From jmckenna at gatewaygeomatics.com  Mon Oct 25 18:59:20 2021
From: jmckenna at gatewaygeomatics.com (Jeff McKenna)
Date: Mon, 25 Oct 2021 22:59:20 -0300
Subject: [mapserver-users] Getting started with mapcache for WMTS...
In-Reply-To: <CAE=8kmnspDr+VOGQaHYHn5H3-9Opd5DB0mZmbBLOatjO5PBviw@mail.gmail.com>
References: <CAE=8kmnspDr+VOGQaHYHn5H3-9Opd5DB0mZmbBLOatjO5PBviw@mail.gmail.com>
Message-ID: <efc37a9c-7d57-f6e7-dc06-8df19011917f@gatewaygeomatics.com>

Hi Carl,

You can also find those same local MapCache configuration steps for MS4W 
online at https://ms4w.com/README_INSTALL.html#mapcache-apache-module 
(but I'm very glad to see you following the steps in your local README, 
great!)  Thanks again for the positive feedback for MS4W.

PS.  someone asked recently the same question about the mailing list 
search, here was my response:


I just use Google search to query the archives:
   - goto https://google.ca
   - in the search box enter: 
site:lists.osgeo.org/pipermail/mapserver-users/ yoursearchterm

This may not always contain the latest archives (depending on the speed 
of Google-mothership bots), but I find it very helpful.




--

Thank-you for using MS4W.
"MS4W: open doors as well as windows"


-jeff




-- 
Jeff McKenna
GatewayGeo: Developers of MS4W, MapServer Consulting and Training
co-founder of FOSS4G
http://gatewaygeo.com/

On 2021-10-25 6:54 p.m., Carl Godkin wrote:
> Hi again,
> 
> Is there a way to search through the mapserver-users mailing list 
> archive except month-by-month?? If so, I can't find it.? (I usually use 
> my Gmail archives for GDAL and PROJ since I've subscribed forever to 
> those but I'm a mapserver newbie.)? Anyway...
> 
> I am trying to configure mapcache to serve from my (local and AFAICT 
> fully working) WMS server on?my MS4W 4.0.5 installation.? I am following 
> the? directions at https://mapserver.org/mapcache/index.html 
> <https://mapserver.org/mapcache/index.html> as well as the steps listed 
> from the local MS4W page?(#mapcache-apache-module).
> 
> I think I understand the basic principles and have added a new source 
> that looks like this:
> 
> <source name="faa" type="wms">
> <getmap>
> <params>
> <FORMAT>image/geotiff</FORMAT>
> <LAYERS>FAA_Sectional_Charts</LAYERS>
> <TRANSPARENT>TRUE</TRANSPARENT>
> <SRS>EPSG:4326</SRS>
> <BBOX>-125,24,-61,49</BBOX>
> </params>
> </getmap>
> <http>
> <url>http://localhost:8080/wms <http://localhost:8080/wms></url>
> </http>
> </source>
> 
> and a new tileset that looks like this:
> 
> <tileset name="faa">
> <source>faa</source>
> <cache>disk</cache>
> <grid>WGS84</grid>
> <format>PNG</format>
> <metatile>5 5</metatile>
> <metabuffer>10</metabuffer>
> <expires>3600</expires>
> </tileset>
> (This is the same as the "test" tileset with a different source.)
> 
> But how do I test this? The query
> 
> http://localhost:8080/mapcache/wmts?service=wmts&request=getcapabilities&version=1.0.0 
> <http://localhost:8080/mapcache/wmts?service=wmts&request=getcapabilities&version=1.0.0>
> 
> lists my new source as a layer. The docs say
> 
>  5.
> 
>     Test your MapCache configuration by following steps in the
>     associated MapCache testing docs
>     <https://mapserver.org/mapcache/install.html#test-your-mapcache-module>
> 
> but that seems to be about testing the mapcache/demo.? I don't see how 
> to add something to the mapcache/demo if that's possible.? What am I 
> overlooking?
> 
> Thanks a lot,
> carl
> 
> 
> 
> 
> 
> 
> _______________________________________________
> MapServer-users mailing list
> MapServer-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/mapserver-users
> 


From jukka.rahkonen at maanmittauslaitos.fi  Tue Oct 26 01:32:27 2021
From: jukka.rahkonen at maanmittauslaitos.fi (Rahkonen Jukka (MML))
Date: Tue, 26 Oct 2021 08:32:27 +0000
Subject: [mapserver-users] Times and axis labels in WCS 2.0.1
Message-ID: <cdc207f4574b401da9b9b1681a618487@maanmittauslaitos.fi>

Hi,

Is there any way to get the time positions of a layer advertised in WCS 2.0.1 DescribeCoverage? They are listed with WCS 1.0.0 but as can be checked with the demo service:
WCS 1.0.0 with time positions https://demo.mapserver.org/cgi-bin/wcs?SERVICE=wcs&VERSION=1.0.0&REQUEST=DescribeCoverage&COVERAGE=modis
WCS 2.0.1 that is unaware of times
https://demo.mapserver.org/cgi-bin/wcs?SERVICE=wcs&VERSION=2.0.1&REQUEST=DescribeCoverage&COVERAGEID=modis

Another service https://beta-karttakuva.maanmittauslaitos.fi/wcs/service/ows?service=WCS&version=2.0.1&request=DescribeCoverage&coverageID=ortokuva__ortokuva
puts timePositions into WCS 2.0.1 like this
<gmlcov:metadata>
<gmlcov:Extension>
<wcsgs:TimeDomain default="2021-12-31T00:00:00.000Z">
<gml:TimeInstant gml:id="ortokuva__ortokuva_td_0">
<gml:timePosition>1985-12-31T00:00:00.000Z</gml:timePosition>
</gml:TimeInstant>

I have already learned that Mapserver cannot handle time as a dimension like Geoserver does https://lists.osgeo.org/pipermail/mapserver-users/2019-January/080917.html but it would be nice at least to have some way to show the users which time periods are available.

Another little issue with Mapserver WCS 2.0.1 is that it shows the coordinate axis labels of EPSG:3067 this way
<gml:Envelope srsName=http://www.opengis.net/def/crs/EPSG/0/3067 axisLabels="x y" uomLabels="m m" srsDimension="2">

Obviously there is no way to change axisLabels into "E N" but I hope I am wrong.

-Jukka Rahkonen-
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211026/9b08f3e8/attachment.htm>

From jmckenna at gatewaygeomatics.com  Tue Oct 26 08:03:44 2021
From: jmckenna at gatewaygeomatics.com (Jeff McKenna)
Date: Tue, 26 Oct 2021 12:03:44 -0300
Subject: [mapserver-users] Getting started with mapcache for WMTS...
In-Reply-To: <efc37a9c-7d57-f6e7-dc06-8df19011917f@gatewaygeomatics.com>
References: <CAE=8kmnspDr+VOGQaHYHn5H3-9Opd5DB0mZmbBLOatjO5PBviw@mail.gmail.com>
 <efc37a9c-7d57-f6e7-dc06-8df19011917f@gatewaygeomatics.com>
Message-ID: <b5396afa-e20b-5c7a-4b67-5701d4dae5a7@gatewaygeomatics.com>

Hi Carl,

I forgot to mention that to keep an eye on the MS4W wiki also, that 
users sometimes post tips etc.  In this case I remember a user asking 
how to configure MS4W + MapCache + WMTS + OpenLayers, so I created the 
wiki page at https://ms4w.com/trac/wiki/MapCacheOpenLayers 
(edits/updates welcomed, as it is indeed a 'wiki' and depends on your 
updates as new software versions etc often change those steps)

As you travel down this page you could also create a wiki page and share 
your tips (or pitfalls) for other users: 
https://ms4w.com/trac#AdditionalNotesTips

Have a nice week,

-jeff




-- 
Jeff McKenna
GatewayGeo: Developers of MS4W, MapServer Consulting and Training
co-founder of FOSS4G
http://gatewaygeo.com/



On 2021-10-25 10:59 p.m., Jeff McKenna wrote:
> Hi Carl,
> 
> You can also find those same local MapCache configuration steps for MS4W 
> online at https://ms4w.com/README_INSTALL.html#mapcache-apache-module 
> (but I'm very glad to see you following the steps in your local README, 
> great!)? Thanks again for the positive feedback for MS4W.
> 
> PS.? someone asked recently the same question about the mailing list 
> search, here was my response:
> 
> 
> I just use Google search to query the archives:
>  ? - goto https://google.ca
>  ? - in the search box enter: 
> site:lists.osgeo.org/pipermail/mapserver-users/ yoursearchterm
> 
> This may not always contain the latest archives (depending on the speed 
> of Google-mothership bots), but I find it very helpful.
> 
> 
> 
> 
> -- 
> 
> Thank-you for using MS4W.
> "MS4W: open doors as well as windows"
> 
> 
> -jeff
> 
> 
> 
> 


From public at postholer.com  Tue Oct 26 16:45:15 2021
From: public at postholer.com (Scott)
Date: Tue, 26 Oct 2021 16:45:15 -0700
Subject: [mapserver-users] Raster pixel alpha
Message-ID: <d89dc572-6870-1bff-3180-e519e703ead6@postholer.com>

I'm trying this but it doesn't seem to work, note the alpha byte on the 
end of colorrange values, IMAGEMODE for png is RGBA:


CLASS
    EXPRESSION ([pixel] >= 0 && [pixel] <= 100)
    STYLE
       datarange 0 100
       colorrange '#ffffffff' '#ffffff00'
    END
END


-- 
www.postholer.com


From cgodkin at gmail.com  Fri Oct 29 16:09:57 2021
From: cgodkin at gmail.com (Carl Godkin)
Date: Fri, 29 Oct 2021 16:09:57 -0700
Subject: [mapserver-users] Point query from Mapserver raster
Message-ID: <CAE=8kmkK08btu-k1tRx1Q0gW9DopxHxt76eyyU7_sUY7rsQz5A@mail.gmail.com>

Hi,
Now that I have my WMS and WCS servers working well (thanks, list!!), I
have one more thing I would like to figure out how to do.

I want to do point queries on [Lat,Lon] points for their elevation from the
terrain layer I have.

I have done some searching and come across various ideas such as Raster
Query (https://mapserver.org/input/raster.html#raster-query) and a query
template  (
https://gis.stackexchange.com/questions/167239/mapserver-raster-picking-elevation-data)
but haven't found enough information for either to implement it.

I came up with a solution that works fine, but might be considered as going
"All Around the Barn"   What works is to add an ASCII output format to my
map file that serves terrain like this:

    OUTPUTFORMAT
      NAME XYZ
      DRIVER "GDAL/XYZ"
      MIMETYPE "text/plain"
      IMAGEMODE FLOAT32
      EXTENSION "txt"
    END

and then use GetCoverage on a small rectangle whose lower-left corner
is the point whose elevation I want:

http://localhost:8080/wcs?SERVICE=WCS&VERSION=2.0.1&REQUEST=GetCoverage&coverageid=SRTM_3_arc-second_grid&format=XYZ&subset=Lat(42.948,42.958)&subset=Lon(-122.086,-122.076)

and then take the first line of the returned output file.

What would Mapserver best practice be for a problem like this?

Thanks very much,

carl
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211029/bee09db7/attachment.htm>

From public at postholer.com  Fri Oct 29 16:24:06 2021
From: public at postholer.com (Scott)
Date: Fri, 29 Oct 2021 16:24:06 -0700
Subject: [mapserver-users] Point query from Mapserver raster
In-Reply-To: <CAE=8kmkK08btu-k1tRx1Q0gW9DopxHxt76eyyU7_sUY7rsQz5A@mail.gmail.com>
References: <CAE=8kmkK08btu-k1tRx1Q0gW9DopxHxt76eyyU7_sUY7rsQz5A@mail.gmail.com>
Message-ID: <6db76431-1878-9e60-9e6f-7e2cfa0340ae@postholer.com>

If you're using postgres/postgis it's easy using the raster version of 
st_intersection. It returns lat/lon/value or more precisely 
geometry/value. See more here:

https://postgis.net/docs/RT_ST_Intersection.html

Then it's just WFS.


On 10/29/21 4:09 PM, Carl Godkin wrote:
> Hi,
> Now that I have my WMS and WCS servers working well (thanks, list!!), I 
> have one more thing I would like to figure out how to do.
> 
> I want to do point queries on [Lat,Lon] points for their elevation from 
> the terrain layer I have.
> 
> I have done some searching and come across various ideas such as Raster 
> Query (https://mapserver.org/input/raster.html#raster-query 
> <https://mapserver.org/input/raster.html#raster-query>) and a query 
> template  
> (https://gis.stackexchange.com/questions/167239/mapserver-raster-picking-elevation-data 
> <https://gis.stackexchange.com/questions/167239/mapserver-raster-picking-elevation-data>) 
> but haven't found enough information for either to implement it.
> 
> I came up with a solution that works fine, but?might be considered as 
> going "All Around the Barn"? ?What works is to add an ASCII output 
> format to my map file that serves terrain like this:
> 
> OUTPUTFORMAT
> NAME XYZ
> DRIVER "GDAL/XYZ"
> MIMETYPE "text/plain"
> IMAGEMODE FLOAT32
> EXTENSION "txt"
> END
> 
> and then use GetCoverage on a small rectangle whose lower-left corner
> is the point whose elevation I want:
> 
> http://localhost:8080/wcs?SERVICE=WCS&VERSION=2.0.1&REQUEST=GetCoverage&coverageid=SRTM_3_arc-second_grid&format=XYZ&subset=Lat(42.948,42.958)&subset=Lon(-122.086,-122.076) 
> <http://localhost:8080/wcs?SERVICE=WCS&VERSION=2.0.1&REQUEST=GetCoverage&coverageid=SRTM_3_arc-second_grid&format=XYZ&subset=Lat(42.948,42.958)&subset=Lon(-122.086,-122.076)>
> 
> and then take the first line of the returned output file.
> 
> What would Mapserver best practice be for a problem like this?
> 
> Thanks very much,
> 
> carl
> 
> 
> 
> 
> 
> 
> 
> 
> 
> _______________________________________________
> MapServer-users mailing list
> MapServer-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/mapserver-users
> 


From cneves at lrec.pt  Sat Oct 30 13:08:01 2021
From: cneves at lrec.pt (Carlos Neves (LREC))
Date: Sat, 30 Oct 2021 21:08:01 +0100
Subject: [mapserver-users] msCGILoadMap(): Web application error. CGI
 variable "map" fails to validate
Message-ID: <75ab939e-0081-a2b4-3912-48a4ffd38fae@lrec.pt>

Hi List,

I have upgrade from 7.6.1 to 7.6.3 and now I have an unexpected error 
message.

"msCGILoadMap(): Web application error. CGI variable "map" fails to 
validate"

I'm running OpenSuSE 15.3. and I'm still using "../" path? in main map 
file path

I have insert in Apache2 and export in OpenSuSE the environment 
MS_MAP_BAD_PATTERN to "false" with no results.

What I have missed?

Any help would be appreciated

Carlos


-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211030/2cdb5872/attachment.htm>

From jukka.rahkonen at maanmittauslaitos.fi  Sun Oct 31 11:48:03 2021
From: jukka.rahkonen at maanmittauslaitos.fi (Rahkonen Jukka (MML))
Date: Sun, 31 Oct 2021 18:48:03 +0000
Subject: [mapserver-users] msCGILoadMap(): Web application error. CGI
 variable "map" fails to validate
Message-ID: <2df3c3ab1b444a859c5fc051c39b6a4d@maanmittauslaitos.fi>

Hi,

Read https://mapserver.org/optimization/limit_mapfile_access.html and https://mapserver.org/development/rfc/ms-rfc-56.html and set MS_MAP_PATTERN.

-Jukka Rahkonen-

L?hett?j?: MapServer-users <mapserver-users-bounces at lists.osgeo.org> Puolesta Carlos Neves (LREC)
L?hetetty: lauantai 30. lokakuuta 2021 23.0
Vastaanottaja: mapserver-users at lists.osgeo.org
Aihe: [mapserver-users] msCGILoadMap(): Web application error. CGI variable "map" fails to validate


Hi List,

I have upgrade from 7.6.1 to 7.6.3 and now I have an unexpected error message.

"msCGILoadMap(): Web application error. CGI variable "map" fails to validate"

I'm running OpenSuSE 15.3. and I'm still using "../" path  in main map file path


I have insert in Apache2 and export in OpenSuSE the environment  MS_MAP_BAD_PATTERN to "false" with no results.

What I have missed?

Any help would be appreciated

Carlos




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/mapserver-users/attachments/20211031/99efd315/attachment.htm>