From pe_lord at protonmail.com Thu Oct 2 06:20:24 2025 From: pe_lord at protonmail.com (PE Lord) Date: Thu, 02 Oct 2025 13:20:24 +0000 Subject: [MapServer-users] Verbosity of PROJ_DEBUG ON Message-ID: <2_hvmcRXxkquQQ6NfLa2B9pmXoHfKsJE7me4-KjAcAFX0YJ85BYt6IbnU72l3WqRrp1rtsLdsIjOEQJuoRLMheT1pJZAlwKcuPWz43-yGRg=@protonmail.com> I have Mapserver 7.6.2 compiled with proj 7.2.1 and gdal version (?). Sorry. I have solved an PROJ issue on my current server. I have modified the PROJ_DATA (PROJ_LIB) sqlite database to handle properly my local projection, including a Canadian NTv2 grid transformation. My ...proj/bin/projinfo, on my server seem to use the right Canadian nadgrid file (ca_nrc_NA83SCRS.tif). Here my question, how to debug PROJ, through Mapserver will a high level of verbose? Here my mapfile CONFIG "CPL_DEBUG" "ON" CONFIG "PROJ_DEBUG" "ON" CONFIG "CPL_CURL_VERBOSE" "ON" CONFIG "CPL_TIMESTAMP" "ON" CONFIG "MS_ERRORFILE" "debug_proj.log" ########################### LAYER DEBUG 5 NAME "MX9-01_20240730_000784-001_0030" DATA "[...path...]/MX9-01_20240730_000784-001_0030.tif" TYPE RASTER PROJECTION "init=epsg:3799" END END ########################### Here the log output: [Wed Aug 27 10:13:25 2025].225353 msDrawRasterLayerLow(MX9-01_20240730_000784-001_0030): entering. [Wed Aug 27 10:13:25 2025].225880 msDrawRasterLayerGDAL(): Entering transform. [Wed Aug 27 10:13:25 2025].225894 msDrawRasterLayerGDAL(): src=1534,1715,419,296, dst=0,0,1300,919 [Wed Aug 27 10:13:25 2025].225899 msDrawRasterLayerGDAL(): source raster PL (1533.977,1715.259) for dst PL (0,0). [Wed Aug 27 10:13:25 2025].225910 msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,0,0,0 [Wed Aug 27 10:13:25 2025].251631 msDrawMap(): Layer 0 (MX9-01_20240730_000784-001_0030), 0.026s [Wed Aug 27 10:13:25 2025].327641 freeLayer(): freeing layer at 0x2727a80. [Wed Aug 27 10:13:25 2025].327746 freeLayer(): freeing layer at 0x27a4850. [Wed Aug 27 10:13:27 2025].569803 freeLayer(): freeing layer at 0x3e92d90.[Wed Aug 27 10:13:27 2025].569862 freeLayer(): freeing layer at 0x2660f60. -------------- next part -------------- An HTML attachment was scrubbed... URL: From alistaire at met.no Fri Oct 3 02:59:45 2025 From: alistaire at met.no (Alistair Everett) Date: Fri, 3 Oct 2025 11:59:45 +0200 Subject: [MapServer-users] UVRASTER reprojection Message-ID: Hi all, I'm trying to use the UVRASTER connection to plot wind barbs. I noticed if I request the layer in a different projection the location of the arrows are reprojected but the angle, and presumably vector components, are not - rather they remain in the same direction relative to my screen. It's a little unclear from the documentation if this is the expected behaviour for UVRASTER or not? I included a pasteboard link to a screenshot to show what I mean, compare for example vectors in the corners of the data and you will see they point in the same direction relative to my screen, rather than being rotated to match the projection. If this is expected, does anyone have any suggested ways of working around this? If this is not expected, does anyone have any ideas what I've done wrong? I can provide some sample data if someone wants to try to reproduce. I'm using: MapServer version 8.4.0 PROJ version 9.4 GDAL version 3.11 OUTPUT=PNG OUTPUT=JPEG SUPPORTS=PROJ SUPPORTS=AGG SUPPORTS=FREETYPE SUPPORTS=CAIRO SUPPORTS=ICONV SUPPORTS=FRIBIDI SUPPORTS=WMS_SERVER SUPPORTS=WMS_CLIENT SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT SUPPORTS=WCS_SERVER SUPPORTS=SOS_SERVER SUPPORTS=OGCAPI_SERVER SUPPORTS=FASTCGI SUPPORTS=THREADS SUPPORTS=GEOS SUPPORTS=PBF INPUT=JPEG INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE INPUT=FLATGEOBUF Screenshot: https://pasteboard.co/WTgbexDmQaAe.png Example mapfile: https://pastebin.com/ibPc6gVK Thanks, Alistair ----------------------------- Dr Alistair Everett Senioringeni?r / Senior Software Engineer Avdeling for Geoutvikling Phone: +47 939 68 985 Email: alistair.everett at met.no Meteorologisk Institutt, Henrik Mohns plass 1, 0313 Oslo ----------------------------- *Dr Alistair Everett* Senioringeni?r / Senior Software Engineer Avdeling for Geoutvikling *Phone:* +47 939 68 985 *Email:* alistair.everett at met.no *Meteorologisk Institutt*, Henrik Mohns plass 1, 0313 Oslo -------------- next part -------------- An HTML attachment was scrubbed... URL: From sethg at geographika.co.uk Sun Oct 5 01:00:35 2025 From: sethg at geographika.co.uk (Seth G) Date: Sun, 05 Oct 2025 10:00:35 +0200 Subject: [MapServer-users] UVRASTER reprojection In-Reply-To: References: Message-ID: Hi Alistair, Could you give an example request (I presume it is WMS?)? Are you rotating your map at all, or is it just the projection? There are some tests which reproject a UVRASTER LAYER at https://github.com/MapServer/MapServer/blob/29e4ae81620a2c50201101c763df0b9d8d540b67/msautotest/wxs/wms_uvraster_map_reprojection.map It sounds like an oversight, so if you are able to create a standalone test case, it would be worth opening an issue. Seth -- web:https://geographika.net & https://mapserverstudio.net mastodon: @geographika at mastodon.social On Fri, Oct 3, 2025, at 11:59 AM, Alistair Everett via MapServer-users wrote: > Hi all, > > I'm trying to use the UVRASTER connection to plot wind barbs. I noticed if I request the layer in a different projection the location of the arrows are reprojected but the angle, and presumably vector components, are not - rather they remain in the same direction relative to my screen. It's a little unclear from the documentation if this is the expected behaviour for UVRASTER or not? I included a pasteboard link to a screenshot to show what I mean, compare for example vectors in the corners of the data and you will see they point in the same direction relative to my screen, rather than being rotated to match the projection. > > If this is expected, does anyone have any suggested ways of working around this? > > If this is not expected, does anyone have any ideas what I've done wrong? I can provide some sample data if someone wants to try to reproduce. > > I'm using: MapServer version 8.4.0 PROJ version 9.4 GDAL version 3.11 OUTPUT=PNG OUTPUT=JPEG SUPPORTS=PROJ SUPPORTS=AGG SUPPORTS=FREETYPE SUPPORTS=CAIRO SUPPORTS=ICONV SUPPORTS=FRIBIDI SUPPORTS=WMS_SERVER SUPPORTS=WMS_CLIENT SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT SUPPORTS=WCS_SERVER SUPPORTS=SOS_SERVER SUPPORTS=OGCAPI_SERVER SUPPORTS=FASTCGI SUPPORTS=THREADS SUPPORTS=GEOS SUPPORTS=PBF INPUT=JPEG INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE INPUT=FLATGEOBUF > > Screenshot: https://pasteboard.co/WTgbexDmQaAe.png > > Example mapfile: https://pastebin.com/ibPc6gVK > > Thanks, > Alistair > > ----------------------------- > Dr Alistair Everett > Senioringeni?r / Senior Software Engineer > Avdeling for Geoutvikling > > Phone: +47 939 68 985 > Email: alistair.everett at met.no > > Meteorologisk Institutt, > Henrik Mohns plass 1, > 0313 Oslo > ----------------------------- > *Dr Alistair Everett* > Senioringeni?r / Senior Software Engineer > Avdeling for Geoutvikling > > *Phone:* +47 939 68 985 > *Email:* alistair.everett at met.no > > *Meteorologisk Institutt*, > Henrik Mohns plass 1, > 0313 Oslo > > _______________________________________________ > 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: From sethg at geographika.co.uk Sun Oct 5 00:53:06 2025 From: sethg at geographika.co.uk (Seth G) Date: Sun, 05 Oct 2025 09:53:06 +0200 Subject: [MapServer-users] Verbosity of PROJ_DEBUG ON In-Reply-To: <2_hvmcRXxkquQQ6NfLa2B9pmXoHfKsJE7me4-KjAcAFX0YJ85BYt6IbnU72l3WqRrp1rtsLdsIjOEQJuoRLMheT1pJZAlwKcuPWz43-yGRg=@protonmail.com> References: <2_hvmcRXxkquQQ6NfLa2B9pmXoHfKsJE7me4-KjAcAFX0YJ85BYt6IbnU72l3WqRrp1rtsLdsIjOEQJuoRLMheT1pJZAlwKcuPWz43-yGRg=@protonmail.com> Message-ID: Hi, You have all the debug options on for MapServer that I am aware of. Maybe you can try CONFIG "PROJ_DEBUG" "5" to see if that is more verbose? https://proj.org/en/stable/usage/environmentvars.html#envvar-PROJ_DEBUG MapServer passes this value on as an environment variable. Maybe something has changed as the MapServer docs always use "ON" and "OFF" whereas PROJ states: > Starting with PROJ 9.3, ON can be used as an alias for 2, and OFF as an alias for 1. Seth -- web:https://geographika.net & https://mapserverstudio.net mastodon: @geographika at mastodon.social On Thu, Oct 2, 2025, at 3:20 PM, PE Lord via MapServer-users wrote: > I have Mapserver 7.6.2 compiled with proj 7.2.1 and gdal version (?). Sorry. > > I have solved an PROJ issue on my current server. I have modified the PROJ_DATA (PROJ_LIB) sqlite database to handle properly my local projection, including a Canadian NTv2 grid transformation. > > My ...proj/bin/projinfo, on my server seem to use the right Canadian nadgrid file (ca_nrc_NA83SCRS.tif). > > Here my question, how to debug PROJ, through Mapserver will a high level of verbose? > > Here my mapfile > > CONFIG "CPL_DEBUG" "ON" > CONFIG "PROJ_DEBUG" "ON" > CONFIG "CPL_CURL_VERBOSE" "ON" > CONFIG "CPL_TIMESTAMP" "ON" > CONFIG "MS_ERRORFILE" "debug_proj.log" > > > > ########################### > > LAYER > DEBUG 5 > NAME "MX9-01_20240730_000784-001_0030" > DATA "[...path...]/MX9-01_20240730_000784-001_0030.tif" > TYPE RASTER > PROJECTION > "init=epsg:3799" > END > END > > ########################### > Here the log output: > > [Wed Aug 27 10:13:25 2025].225353 msDrawRasterLayerLow(MX9-01_20240730_000784-001_0030): entering. > [Wed Aug 27 10:13:25 2025].225880 msDrawRasterLayerGDAL(): Entering transform. > [Wed Aug 27 10:13:25 2025].225894 msDrawRasterLayerGDAL(): src=1534,1715,419,296, dst=0,0,1300,919 > [Wed Aug 27 10:13:25 2025].225899 msDrawRasterLayerGDAL(): source raster PL (1533.977,1715.259) for dst PL (0,0). > [Wed Aug 27 10:13:25 2025].225910 msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,0,0,0 > [Wed Aug 27 10:13:25 2025].251631 msDrawMap(): Layer 0 (MX9-01_20240730_000784-001_0030), 0.026s > [Wed Aug 27 10:13:25 2025].327641 freeLayer(): freeing layer at 0x2727a80. > [Wed Aug 27 10:13:25 2025].327746 freeLayer(): freeing layer at 0x27a4850. > [Wed Aug 27 10:13:27 2025].569803 freeLayer(): freeing layer at 0x3e92d90. > [Wed Aug 27 10:13:27 2025].569862 freeLayer(): freeing layer at 0x2660f60. > > _______________________________________________ > 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: From sethg at geographika.co.uk Sun Oct 5 01:06:34 2025 From: sethg at geographika.co.uk (Seth G) Date: Sun, 05 Oct 2025 10:06:34 +0200 Subject: [MapServer-users] WMS requests over the antimeridian : inconsistent responses, depending on east or west side In-Reply-To: References: Message-ID: <3c453f4c-615d-4eb5-923a-02558cc2339f@app.fastmail.com> Hi Philippe, I'd agree it seems strange this isn't consistent. There is a lot of code in the MapServer codebase to handle data crossing the antemeridian, so I guess there is an oversight somewhere (and fixing it is likely to be tricky). I think this should be added as an issue, with the examples you provided. Seth -- web:https://geographika.net & https://mapserverstudio.net mastodon: @geographika at mastodon.social On Fri, Sep 19, 2025, at 4:13 PM, Philippe Ghesquiere via MapServer-users wrote: > Dear all, > > I noticed a weird behaviour with a WMS request based on EPSG:4326, when the BBOX exeeds the [-180 ; +180] longitude intervalle. > We bump into this problem when we try to display images over the antimeridian line. > > I could observe this problem with the "MapServer Demonstration Server". > > *_Observed _*Mapserver response is not consistent whether the BBOX overlaps eastward or westward : > 1) Westward > BBOX=-90,-200,90,100 > The image is clipped on both east and west sides, that is in the [-180? ; +100?] longitude intervalle. > On the west side, pixels in the [-200? ; -180?] longitude intervalle are transparent. > > The complete request URL is : > https://demo.mapserver.org/cgi-bin/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-90,-200,90,100&CRS=EPSG:4326&WIDTH=953&HEIGHT=480&LAYERS=bluemarble&STYLES=&FORMAT=image/png&TRANSPARENT=true&STYLES= > > 2) Eastward > BBOX=-90,-100,90,200 > The image is clipped on the west side, which is OK. > On the east side, pixels in the [+180?; +200?] longitude intervalle are not correct. > > The complete request URL is : > https://demo.mapserver.org/cgi-bin/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-90,-100,90,200&CRS=EPSG:4326&WIDTH=953&HEIGHT=480&LAYERS=bluemarble&STYLES=&FORMAT=image/png&TRANSPARENT=true&STYLES= > > > *_Expected _*At least, I would expect to have the same behaviour on both sides of the antimeridian. > > Furthermore, couldn't we get the "repeated world" on the [-200? ; -180?] (respectively [+180? ; +200?]) longitude intervalle ? > > Thanks for your help > Philippe > The information in this e-mail is confidential. The contents may not be disclosed or used by anyone other than the addressee. Access to this e-mail by anyone else is unauthorised. > If you are not the intended recipient, please notify Airbus immediately and delete this e-mail. > Airbus cannot accept any responsibility for the accuracy or completeness of this e-mail as it has been sent over public networks. If you have any concerns over the content of this message or its Accuracy or Integrity, please contact Airbus immediately. > All outgoing e-mails from Airbus are checked using regularly updated virus scanning software but you should take whatever measures you deem to be appropriate to ensure that this message and any attachments are virus free. > _______________________________________________ > 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: From philippe.ghesquiere at airbus.com Mon Oct 6 00:26:35 2025 From: philippe.ghesquiere at airbus.com (Philippe Ghesquiere) Date: Mon, 6 Oct 2025 09:26:35 +0200 Subject: [MapServer-users] WMS requests over the antimeridian : inconsistent responses, depending on east or west side In-Reply-To: <3c453f4c-615d-4eb5-923a-02558cc2339f@app.fastmail.com> References: <3c453f4c-615d-4eb5-923a-02558cc2339f@app.fastmail.com> Message-ID: Hi Seth, Thanks for your answer. I just created a new issue: https://github.com/MapServer/MapServer/issues/7357 I think this new issue is related to #7212 and #5162 Philippe On Sun, Oct 5, 2025 at 10:07?AM Seth G wrote: > Hi Philippe, > > I'd agree it seems strange this isn't consistent. There is a lot of code > in the MapServer codebase to handle data crossing the antemeridian, so I > guess there is an oversight somewhere (and fixing it is likely to be > tricky). I think this should be added as an issue, with the examples you > provided. > > Seth > > -- > web:https://geographika.net & https://mapserverstudio.net > mastodon: @geographika at mastodon.social > > On Fri, Sep 19, 2025, at 4:13 PM, Philippe Ghesquiere via MapServer-users > wrote: > > Dear all, > > I noticed a weird behaviour with a WMS request based on EPSG:4326, when > the BBOX exeeds the [-180 ; +180] longitude intervalle. > We bump into this problem when we try to display images over the > antimeridian line. > > I could observe this problem with the "MapServer Demonstration Server". > > > *Observed*Mapserver response is not consistent whether the BBOX overlaps > eastward or westward : > 1) Westward > BBOX=-90,-200,90,100 > The image is clipped on both east and west sides, that is in the [-180? ; > +100?] longitude intervalle. > On the west side, pixels in the [-200? ; -180?] longitude intervalle are > transparent. > > The complete request URL is : > > https://demo.mapserver.org/cgi-bin/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-90,-200,90,100&CRS=EPSG:4326&WIDTH=953&HEIGHT=480&LAYERS=bluemarble&STYLES=&FORMAT=image/png&TRANSPARENT=true&STYLES= > > 2) Eastward > BBOX=-90,-100,90,200 > The image is clipped on the west side, which is OK. > On the east side, pixels in the [+180?; +200?] longitude intervalle are > not correct. > > The complete request URL is : > > https://demo.mapserver.org/cgi-bin/wms?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-90,-100,90,200&CRS=EPSG:4326&WIDTH=953&HEIGHT=480&LAYERS=bluemarble&STYLES=&FORMAT=image/png&TRANSPARENT=true&STYLES= > > > > *Expected*At least, I would expect to have the same behaviour on both > sides of the antimeridian. > > Furthermore, couldn't we get the "repeated world" on the [-200? ; -180?] > (respectively [+180? ; +200?]) longitude intervalle ? > > Thanks for your help > Philippe > The information in this e-mail is confidential. The contents may not be > disclosed or used by anyone other than the addressee. Access to this e-mail > by anyone else is unauthorised. > If you are not the intended recipient, please notify Airbus immediately > and delete this e-mail. > Airbus cannot accept any responsibility for the accuracy or completeness > of this e-mail as it has been sent over public networks. If you have any > concerns over the content of this message or its Accuracy or Integrity, > please contact Airbus immediately. > All outgoing e-mails from Airbus are checked using regularly updated virus > scanning software but you should take whatever measures you deem to be > appropriate to ensure that this message and any attachments are virus free. > _______________________________________________ > MapServer-users mailing list > MapServer-users at lists.osgeo.org > https://lists.osgeo.org/mailman/listinfo/mapserver-users > > > The information in this e-mail is confidential. The contents may not be disclosed or used by anyone other than the addressee. Access to this e-mail by anyone else is unauthorised. If you are not the intended recipient, please notify Airbus immediately and delete this e-mail. Airbus cannot accept any responsibility for the accuracy or completeness of this e-mail as it has been sent over public networks. If you have any concerns over the content of this message or its Accuracy or Integrity, please contact Airbus immediately. All outgoing e-mails from Airbus are checked using regularly updated virus scanning software but you should take whatever measures you deem to be appropriate to ensure that this message and any attachments are virus free. -------------- next part -------------- An HTML attachment was scrubbed... URL: From alistaire at met.no Mon Oct 6 03:16:42 2025 From: alistaire at met.no (Alistair Everett) Date: Mon, 6 Oct 2025 12:16:42 +0200 Subject: [MapServer-users] UVRASTER reprojection In-Reply-To: References: Message-ID: Hi Seth, Thanks for your reply, we're just reprojecting the data. I'll take a closer look at the test case too, since it seems like that should be reprojecting from stereographic to geographic. I should also be able to pull together some test data in the next few days too, I can open an issue with it when I've put it together. Here are a couple of example requests: One in EPSG:4326 localhost:8000/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.3.0&FORMAT=image%2Fpng&STYLES=wind_barb&TRANSPARENT=TRUE&TIME=2025-10-06T07%3A00%3A00Z&DIM_PRESSURE=50.0&DIM_REFERENCE_TIME=2025-10-06T03%3A00%3A00Z&LAYERS=wind_pl_vector_meps_det_vdiv_2_5km_calculations&WIDTH=1304&HEIGHT=1504&CRS=EPSG%3A3857&BBOX=-2460529.2703191666%2C4596419.518351211%2C6809846.057963322%2C15288631.55342549 https://pasteboard.co/jI1eMKaj3PT1.png This is one in EPSG:3411: localhost:8000/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.3.0&FORMAT=image%2Fpng&STYLES=wind_barb&TRANSPARENT=TRUE&TIME=2025-10-06T16%3A00%3A00Z&DIM_PRESSURE=50.0&DIM_REFERENCE_TIME=2025-10-06T03%3A00%3A00Z&LAYERS=wind_pl_vector_meps_det_vdiv_2_5km_calculations&WIDTH=1438&HEIGHT=1504&CRS=EPSG%3A3411&BBOX=404935.5094684181%2C-3769428.4764078585%2C4977748.7994076125%2C1013263.5876175575 https://pasteboard.co/IHZjgXibtLyJ.png Thanks, Alistair On Sun, 5 Oct 2025 at 10:00, Seth G wrote: > Hi Alistair, > > Could you give an example request (I presume it is WMS?)? Are you rotating > your map at all, or is it just the projection? > > There are some tests which reproject a UVRASTER LAYER at > https://github.com/MapServer/MapServer/blob/29e4ae81620a2c50201101c763df0b9d8d540b67/msautotest/wxs/wms_uvraster_map_reprojection.map > It sounds like an oversight, so if you are able to create a standalone > test case, it would be worth opening an issue. > > Seth > > -- > web:https://geographika.net & https://mapserverstudio.net > mastodon: @geographika at mastodon.social > > On Fri, Oct 3, 2025, at 11:59 AM, Alistair Everett via MapServer-users > wrote: > > Hi all, > > I'm trying to use the UVRASTER connection to plot wind barbs. I noticed if > I request the layer in a different projection the location of the arrows > are reprojected but the angle, and presumably vector components, are not - > rather they remain in the same direction relative to my screen. It's a > little unclear from the documentation if this is the expected behaviour for > UVRASTER or not? I included a pasteboard link to a screenshot to show what > I mean, compare for example vectors in the corners of the data and you will > see they point in the same direction relative to my screen, rather than > being rotated to match the projection. > > If this is expected, does anyone have any suggested ways of working around > this? > > If this is not expected, does anyone have any ideas what I've done wrong? > I can provide some sample data if someone wants to try to reproduce. > > I'm using: MapServer version 8.4.0 PROJ version 9.4 GDAL version 3.11 > OUTPUT=PNG OUTPUT=JPEG SUPPORTS=PROJ SUPPORTS=AGG SUPPORTS=FREETYPE > SUPPORTS=CAIRO SUPPORTS=ICONV SUPPORTS=FRIBIDI SUPPORTS=WMS_SERVER > SUPPORTS=WMS_CLIENT SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT > SUPPORTS=WCS_SERVER SUPPORTS=SOS_SERVER SUPPORTS=OGCAPI_SERVER > SUPPORTS=FASTCGI SUPPORTS=THREADS SUPPORTS=GEOS SUPPORTS=PBF INPUT=JPEG > INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE INPUT=FLATGEOBUF > > Screenshot: https://pasteboard.co/WTgbexDmQaAe.png > > Example mapfile: https://pastebin.com/ibPc6gVK > > Thanks, > Alistair > > ----------------------------- > Dr Alistair Everett > Senioringeni?r / Senior Software Engineer > Avdeling for Geoutvikling > > Phone: +47 939 68 985 > Email: alistair.everett at met.no > > Meteorologisk Institutt, > Henrik Mohns plass 1, > 0313 Oslo > ----------------------------- > *Dr Alistair Everett* > Senioringeni?r / Senior Software Engineer > Avdeling for Geoutvikling > > *Phone:* +47 939 68 985 > *Email:* alistair.everett at met.no > > *Meteorologisk Institutt*, > Henrik Mohns plass 1, > 0313 Oslo > > _______________________________________________ > 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: From Louis-Philippe.RousseauLambert at ec.gc.ca Mon Oct 6 04:04:45 2025 From: Louis-Philippe.RousseauLambert at ec.gc.ca (Rousseau Lambert,Louis-Philippe (ECCC)) Date: Mon, 6 Oct 2025 11:04:45 +0000 Subject: [MapServer-users] UVRASTER reprojection In-Reply-To: References: Message-ID: Hi, I am tagging along this thread as I?m also very interested to see if there is a simple answer to fix this behaviour ? This is also an issue we had with using UVRaster connection in MapServer with non-global weather models on our side. Luckily for us, we also have access to the wind direction and wind speed variables for our non-global models and ending up using those variables instead, because the wind speed and direction are corrected properly to be relative to the North instead of using the wind vector components (U and V) which are always relative to the model?s grid axis. For this we needed to update the styles, we were using though as it is not the regular use for the UVRaster connection type. Because we are passing the wind direction and speed directly and variables, we are using the U and V values accordingly for the styles. In our case U becomes the speed and V the direction instead of using uv_length and uv_angle. MapServer class example: CLASS NAME "0 - 5 kts" EXPRESSION ([v] > 0.00 AND [v] <= 5) STYLE COLOR 0 0 127 SYMBOL "arrow_wind" ANGLE [u] SIZE 12 END END For your curiosity, here is a layer example using this method in our viewer MSC AniMet Maybe this could also be useful on your side ? LP -- Louis-Philippe Rousseau-Lambert, B.Sc. Sp?cialiste en g?omatique / Geomatics Specialist Syst?mes de donn?es g?ospatiales et ouvertes Geospatial and Open Data Systems Service M?t?orologique du Canada Meteorological Service of Canada Environnement et Changement Climatique Canada Environment and Climate Change Canada louis-philippe.rousseaulambert at ec.gc.ca (438) 801-0599 From: MapServer-users On Behalf Of Alistair Everett via MapServer-users Sent: Monday, October 6, 2025 6:17 AM To: Seth G Cc: MapServer Users Subject: Re: [MapServer-users] UVRASTER reprojection Hi Seth, Thanks for your reply, we're just reprojecting the data. I'll take a closer look at the test case too, since it seems like that should be reprojecting from stereographic to geographic. I should also be able to pull together some test data in the next few days too, I can open an issue with it when I've put it together. Here are a couple of example requests: One in EPSG:4326 localhost:8000/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.3.0&FORMAT=image%2Fpng&STYLES=wind_barb&TRANSPARENT=TRUE&TIME=2025-10-06T07%3A00%3A00Z&DIM_PRESSURE=50.0&DIM_REFERENCE_TIME=2025-10-06T03%3A00%3A00Z&LAYERS=wind_pl_vector_meps_det_vdiv_2_5km_calculations&WIDTH=1304&HEIGHT=1504&CRS=EPSG%3A3857&BBOX=-2460529.2703191666%2C4596419.518351211%2C6809846.057963322%2C15288631.55342549 https://pasteboard.co/jI1eMKaj3PT1.png This is one in EPSG:3411: localhost:8000/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.3.0&FORMAT=image%2Fpng&STYLES=wind_barb&TRANSPARENT=TRUE&TIME=2025-10-06T16%3A00%3A00Z&DIM_PRESSURE=50.0&DIM_REFERENCE_TIME=2025-10-06T03%3A00%3A00Z&LAYERS=wind_pl_vector_meps_det_vdiv_2_5km_calculations&WIDTH=1438&HEIGHT=1504&CRS=EPSG%3A3411&BBOX=404935.5094684181%2C-3769428.4764078585%2C4977748.7994076125%2C1013263.5876175575 https://pasteboard.co/IHZjgXibtLyJ.png Thanks, Alistair On Sun, 5 Oct 2025 at 10:00, Seth G > wrote: Hi Alistair, Could you give an example request (I presume it is WMS?)? Are you rotating your map at all, or is it just the projection? There are some tests which reproject a UVRASTER LAYER at https://github.com/MapServer/MapServer/blob/29e4ae81620a2c50201101c763df0b9d8d540b67/msautotest/wxs/wms_uvraster_map_reprojection.map It sounds like an oversight, so if you are able to create a standalone test case, it would be worth opening an issue. Seth -- web:https://geographika.net & https://mapserverstudio.net mastodon: @geographika at mastodon.social On Fri, Oct 3, 2025, at 11:59 AM, Alistair Everett via MapServer-users wrote: Hi all, I'm trying to use the UVRASTER connection to plot wind barbs. I noticed if I request the layer in a different projection the location of the arrows are reprojected but the angle, and presumably vector components, are not - rather they remain in the same direction relative to my screen. It's a little unclear from the documentation if this is the expected behaviour for UVRASTER or not? I included a pasteboard link to a screenshot to show what I mean, compare for example vectors in the corners of the data and you will see they point in the same direction relative to my screen, rather than being rotated to match the projection. If this is expected, does anyone have any suggested ways of working around this? If this is not expected, does anyone have any ideas what I've done wrong? I can provide some sample data if someone wants to try to reproduce. I'm using: MapServer version 8.4.0 PROJ version 9.4 GDAL version 3.11 OUTPUT=PNG OUTPUT=JPEG SUPPORTS=PROJ SUPPORTS=AGG SUPPORTS=FREETYPE SUPPORTS=CAIRO SUPPORTS=ICONV SUPPORTS=FRIBIDI SUPPORTS=WMS_SERVER SUPPORTS=WMS_CLIENT SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT SUPPORTS=WCS_SERVER SUPPORTS=SOS_SERVER SUPPORTS=OGCAPI_SERVER SUPPORTS=FASTCGI SUPPORTS=THREADS SUPPORTS=GEOS SUPPORTS=PBF INPUT=JPEG INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE INPUT=FLATGEOBUF Screenshot: https://pasteboard.co/WTgbexDmQaAe.png Example mapfile: https://pastebin.com/ibPc6gVK Thanks, Alistair ----------------------------- Dr Alistair Everett Senioringeni?r / Senior Software Engineer Avdeling for Geoutvikling Phone: +47 939 68 985 Email: alistair.everett at met.no Meteorologisk Institutt, Henrik Mohns plass 1, 0313 Oslo ----------------------------- Dr Alistair Everett Senioringeni?r / Senior Software Engineer Avdeling for Geoutvikling Phone: +47 939 68 985 Email: alistair.everett at met.no Meteorologisk Institutt, Henrik Mohns plass 1, 0313 Oslo _______________________________________________ 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: From Robert.Sanson at asurequality.com Wed Oct 8 17:47:15 2025 From: Robert.Sanson at asurequality.com (Robert Sanson) Date: Thu, 9 Oct 2025 00:47:15 +0000 Subject: [MapServer-users] How to reference an ArcGIS FeatureService REST endpoint in the mapfile? In-Reply-To: References: Message-ID: Hi All Does anyone have an example of setting up a layer in a Mapfile to cascade an ArcGIS FeatureService REST endpoint? Many thanks Robert Sanson New Zealand This email, including any attachments, may contain information which is confidential or subject to copyright. If you are not the intended recipient, please notify us immediately and delete this email from your system. Email communications are not secure and are not guaranteed by AsureQuality Limited to be free of unauthorised interference, error or virus. Anyone who communicates with us by email is taken to accept this risk. Anything in this email which does not relate to the official business of AsureQuality Limited is neither given nor endorsed by AsureQuality Limited. Please refer to www.asurequality.com or contact AsureQuality Limited for more information, including about our terms of business. -------------- next part -------------- An HTML attachment was scrubbed... URL: From sethg at geographika.co.uk Sat Oct 11 02:57:02 2025 From: sethg at geographika.co.uk (Seth G) Date: Sat, 11 Oct 2025 11:57:02 +0200 Subject: [MapServer-users] How to reference an ArcGIS FeatureService REST endpoint in the mapfile? In-Reply-To: References: Message-ID: <0b1e88be-0082-4cd9-bc58-3c5876a66310@app.fastmail.com> Hi Robert, This works, but renders very slowly. It uses https://gdal.org/en/stable/drivers/vector/esrijson.html underneath. map2img -m arcgis.map -o map.png MAP NAME "PoolPermits" EXTENT -13075816 4014771 -13073005 4016869 SIZE 400 400 PROJECTION "init=epsg:3857" END LAYER NAME "buildings" TYPE POLYGON STATUS DEFAULT CONNECTIONTYPE OGR CONNECTION "https://sampleserver6.arcgisonline.com/arcgis/rest/services/PoolPermits/FeatureServer/0/query?resultRecordCount=10&f=pjson" CLASS STYLE COLOR 246 241 223 OUTLINECOLOR 0 0 0 END END END END Seth -- web:https://geographika.net & https://mapserverstudio.net mastodon: @geographika at mastodon.social On Thu, Oct 9, 2025, at 2:47 AM, Robert Sanson via MapServer-users wrote: > Hi All > > Does anyone have an example of setting up a layer in a Mapfile to cascade an ArcGIS FeatureService REST endpoint? > > Many thanks > > Robert Sanson > New Zealand > > This email, including any attachments, may contain information which is confidential or subject to copyright. If you are not the intended recipient, please notify us immediately and delete this email from your system. Email communications are not secure and are not guaranteed by AsureQuality Limited to be free of unauthorised interference, error or virus. Anyone who communicates with us by email is taken to accept this risk. Anything in this email which does not relate to the official business of AsureQuality Limited is neither given nor endorsed by AsureQuality Limited. Please refer to www.asurequality.com or contact AsureQuality Limited for more information, including about our terms of business. > _______________________________________________ > 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: From sethg at geographika.co.uk Sat Oct 11 03:00:38 2025 From: sethg at geographika.co.uk (Seth G) Date: Sat, 11 Oct 2025 12:00:38 +0200 Subject: [MapServer-users] UVRASTER reprojection In-Reply-To: References: Message-ID: <74b6adef-6a57-47c3-9b62-e6e11f17aa41@app.fastmail.com> Hi Alistair, The Mapfile and requests look fine to me, so this is likely an oversight in the implementation. It would be good to get this logged as an issue, ideally with a recreatable test case at https://github.com/MapServer/MapServer/issues. Seth -- web:https://geographika.net & https://mapserverstudio.net mastodon: @geographika at mastodon.social On Mon, Oct 6, 2025, at 12:16 PM, Alistair Everett wrote: > Hi Seth, > > Thanks for your reply, we're just reprojecting the data. I'll take a closer look at the test case too, since it seems like that should be reprojecting from stereographic to geographic. I should also be able to pull together some test data in the next few days too, I can open an issue with it when I've put it together. > > Here are a couple of example requests: > > One in EPSG:4326 > localhost:8000/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.3.0&FORMAT=image%2Fpng&STYLES=wind_barb&TRANSPARENT=TRUE&TIME=2025-10-06T07%3A00%3A00Z&DIM_PRESSURE=50.0&DIM_REFERENCE_TIME=2025-10-06T03%3A00%3A00Z&LAYERS=wind_pl_vector_meps_det_vdiv_2_5km_calculations&WIDTH=1304&HEIGHT=1504&CRS=EPSG%3A3857&BBOX=-2460529.2703191666%2C4596419.518351211%2C6809846.057963322%2C15288631.55342549 > > https://pasteboard.co/jI1eMKaj3PT1.png > > This is one in EPSG:3411: > localhost:8000/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.3.0&FORMAT=image%2Fpng&STYLES=wind_barb&TRANSPARENT=TRUE&TIME=2025-10-06T16%3A00%3A00Z&DIM_PRESSURE=50.0&DIM_REFERENCE_TIME=2025-10-06T03%3A00%3A00Z&LAYERS=wind_pl_vector_meps_det_vdiv_2_5km_calculations&WIDTH=1438&HEIGHT=1504&CRS=EPSG%3A3411&BBOX=404935.5094684181%2C-3769428.4764078585%2C4977748.7994076125%2C1013263.5876175575 > > https://pasteboard.co/IHZjgXibtLyJ.png > > Thanks, > Alistair > > > On Sun, 5 Oct 2025 at 10:00, Seth G wrote: >> __ >> Hi Alistair, >> >> Could you give an example request (I presume it is WMS?)? Are you rotating your map at all, or is it just the projection? >> >> There are some tests which reproject a UVRASTER LAYER at https://github.com/MapServer/MapServer/blob/29e4ae81620a2c50201101c763df0b9d8d540b67/msautotest/wxs/wms_uvraster_map_reprojection.map >> It sounds like an oversight, so if you are able to create a standalone test case, it would be worth opening an issue. >> >> Seth >> >> -- >> web:https://geographika.net & https://mapserverstudio.net >> mastodon: @geographika at mastodon.social >> >> On Fri, Oct 3, 2025, at 11:59 AM, Alistair Everett via MapServer-users wrote: >>> Hi all, >>> >>> I'm trying to use the UVRASTER connection to plot wind barbs. I noticed if I request the layer in a different projection the location of the arrows are reprojected but the angle, and presumably vector components, are not - rather they remain in the same direction relative to my screen. It's a little unclear from the documentation if this is the expected behaviour for UVRASTER or not? I included a pasteboard link to a screenshot to show what I mean, compare for example vectors in the corners of the data and you will see they point in the same direction relative to my screen, rather than being rotated to match the projection. >>> >>> If this is expected, does anyone have any suggested ways of working around this? >>> >>> If this is not expected, does anyone have any ideas what I've done wrong? I can provide some sample data if someone wants to try to reproduce. >>> >>> I'm using: MapServer version 8.4.0 PROJ version 9.4 GDAL version 3.11 OUTPUT=PNG OUTPUT=JPEG SUPPORTS=PROJ SUPPORTS=AGG SUPPORTS=FREETYPE SUPPORTS=CAIRO SUPPORTS=ICONV SUPPORTS=FRIBIDI SUPPORTS=WMS_SERVER SUPPORTS=WMS_CLIENT SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT SUPPORTS=WCS_SERVER SUPPORTS=SOS_SERVER SUPPORTS=OGCAPI_SERVER SUPPORTS=FASTCGI SUPPORTS=THREADS SUPPORTS=GEOS SUPPORTS=PBF INPUT=JPEG INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE INPUT=FLATGEOBUF >>> >>> Screenshot: https://pasteboard.co/WTgbexDmQaAe.png >>> >>> Example mapfile: https://pastebin.com/ibPc6gVK >>> >>> Thanks, >>> Alistair >>> >>> ----------------------------- >>> Dr Alistair Everett >>> Senioringeni?r / Senior Software Engineer >>> Avdeling for Geoutvikling >>> >>> Phone: +47 939 68 985 >>> Email: alistair.everett at met.no >>> >>> Meteorologisk Institutt, >>> Henrik Mohns plass 1, >>> 0313 Oslo >>> ----------------------------- >>> *Dr Alistair Everett* >>> Senioringeni?r / Senior Software Engineer >>> Avdeling for Geoutvikling >>> >>> *Phone:* +47 939 68 985 >>> *Email:* alistair.everett at met.no >>> >>> *Meteorologisk Institutt*, >>> Henrik Mohns plass 1, >>> 0313 Oslo >>> >>> _______________________________________________ >>> 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: From Robert.Sanson at asurequality.com Sun Oct 12 17:43:46 2025 From: Robert.Sanson at asurequality.com (Robert Sanson) Date: Mon, 13 Oct 2025 00:43:46 +0000 Subject: [MapServer-users] [EXT] Re: How to reference an ArcGIS FeatureService REST endpoint in the mapfile? In-Reply-To: <0b1e88be-0082-4cd9-bc58-3c5876a66310@app.fastmail.com> References: <0b1e88be-0082-4cd9-bc58-3c5876a66310@app.fastmail.com> Message-ID: Thanks Seth! Good to know it can work. Still trying to get it to show anything. Here is my layer: LAYER NAME "NZ_Property" TYPE POLYGON STATUS DEFAULT CONNECTIONTYPE OGR CONNECTION https://services.arcgis.com/xdsHIIxuCWByZiCB/arcgis/rest/services/LINZ_NZ_Property_Boundaries/FeatureServer/0/query?resultRecordCount=100&f=geoJSON PROJECTION "init=epsg:2193" END CLASS STYLE WIDTH 2 COLOR 255 255 0 OUTLINECOLOR 0 0 0 END END END Mapserver v. 8.2.2 is not showing any error messages even though I have DEBUG 5 on my map. The layer projection matches the MAP projection. According to the documentation for the site, Supported query Formats: JSON, geoJSON, PBF Thanks Robert From: Seth G Sent: Saturday, 11 October 2025 10:57 p.m. To: Robert Sanson ; MapServer Users Subject: [EXT] Re: [MapServer-users] How to reference an ArcGIS FeatureService REST endpoint in the mapfile? CAUTION: This email originated from outside of the organisation. Do not click links or open attachments unless you recognise the sender and know the content is safe. Hi Robert, This works, but renders very slowly. It uses https://gdal.org/en/stable/drivers/vector/esrijson.html underneath. map2img -m arcgis.map -o map.png MAP NAME "PoolPermits" EXTENT -13075816 4014771 -13073005 4016869 SIZE 400 400 PROJECTION "init=epsg:3857" END LAYER NAME "buildings" TYPE POLYGON STATUS DEFAULT CONNECTIONTYPE OGR CONNECTION "https://sampleserver6.arcgisonline.com/arcgis/rest/services/PoolPermits/FeatureServer/0/query?resultRecordCount=10&f=pjson" CLASS STYLE COLOR 246 241 223 OUTLINECOLOR 0 0 0 END END END END Seth -- web:https://geographika.net & https://mapserverstudio.net mastodon: @geographika at mastodon.social On Thu, Oct 9, 2025, at 2:47 AM, Robert Sanson via MapServer-users wrote: Hi All Does anyone have an example of setting up a layer in a Mapfile to cascade an ArcGIS FeatureService REST endpoint? Many thanks Robert Sanson New Zealand This email, including any attachments, may contain information which is confidential or subject to copyright. If you are not the intended recipient, please notify us immediately and delete this email from your system. Email communications are not secure and are not guaranteed by AsureQuality Limited to be free of unauthorised interference, error or virus. Anyone who communicates with us by email is taken to accept this risk. Anything in this email which does not relate to the official business of AsureQuality Limited is neither given nor endorsed by AsureQuality Limited. Please refer to www.asurequality.com or contact AsureQuality Limited for more information, including about our terms of business. _______________________________________________ MapServer-users mailing list MapServer-users at lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/mapserver-users This email, including any attachments, may contain information which is confidential or subject to copyright. If you are not the intended recipient, please notify us immediately and delete this email from your system. Email communications are not secure and are not guaranteed by AsureQuality Limited to be free of unauthorised interference, error or virus. Anyone who communicates with us by email is taken to accept this risk. Anything in this email which does not relate to the official business of AsureQuality Limited is neither given nor endorsed by AsureQuality Limited. Please refer to www.asurequality.com or contact AsureQuality Limited for more information, including about our terms of business. -------------- next part -------------- An HTML attachment was scrubbed... URL: From michael.smith.erdc at gmail.com Mon Oct 13 02:45:29 2025 From: michael.smith.erdc at gmail.com (Michael Smith) Date: Mon, 13 Oct 2025 05:45:29 -0400 Subject: [MapServer-users] [EXT] Re: How to reference an ArcGIS FeatureService REST endpoint in the mapfile? Message-ID: <2E398F76-3CB7-44EA-8B57-F3BC60D47E70@gmail.com> Robert, your request is not returning data if you try it independent of mapserver. You need to add a where option it seems. Try gdal vector info --if esrijson https://services.arcgis.com/xdsHIIxuCWByZiCB/ArcGIS/rest/services/LINZ_NZ_Property_Boundaries/FeatureServer/0/query\?where\=objectID%3E%3D0\&f\=pgeojson which returns: "featureCount":2732041, Or keep the resultRecordCount to limit Mike -- Michael Smith RSGIS Center ? ERDC CRREL NH US Army Corps From: MapServer-users on behalf of Robert Sanson via MapServer-users Reply-To: Robert Sanson Date: Sunday, October 12, 2025 at 8:43?PM To: Seth G Cc: MapServer Users Subject: Re: [MapServer-users] [EXT] Re: How to reference an ArcGIS FeatureService REST endpoint in the mapfile? Thanks Seth! Good to know it can work. Still trying to get it to show anything. Here is my layer: LAYER NAME "NZ_Property" TYPE POLYGON STATUS DEFAULT CONNECTIONTYPE OGR CONNECTION https://services.arcgis.com/xdsHIIxuCWByZiCB/arcgis/rest/services/LINZ_NZ_Property_Boundaries/FeatureServer/0/query?resultRecordCount=100&f=geoJSON PROJECTION "init=epsg:2193" END CLASS STYLE WIDTH 2 COLOR 255 255 0 OUTLINECOLOR 0 0 0 END END END Mapserver v. 8.2.2 is not showing any error messages even though I have DEBUG 5 on my map. The layer projection matches the MAP projection. According to the documentation for the site, Supported query Formats: JSON, geoJSON, PBF Thanks Robert From: Seth G Sent: Saturday, 11 October 2025 10:57 p.m. To: Robert Sanson ; MapServer Users Subject: [EXT] Re: [MapServer-users] How to reference an ArcGIS FeatureService REST endpoint in the mapfile? CAUTION: This email originated from outside of the organisation. Do not click links or open attachments unless you recognise the sender and know the content is safe. Hi Robert, This works, but renders very slowly. It uses https://gdal.org/en/stable/drivers/vector/esrijson.html underneath. map2img -m arcgis.map -o map.png MAP NAME "PoolPermits" EXTENT -13075816 4014771 -13073005 4016869 SIZE 400 400 PROJECTION "init=epsg:3857" END LAYER NAME "buildings" TYPE POLYGON STATUS DEFAULT CONNECTIONTYPE OGR CONNECTION "https://sampleserver6.arcgisonline.com/arcgis/rest/services/PoolPermits/FeatureServer/0/query?resultRecordCount=10&f=pjson" CLASS STYLE COLOR 246 241 223 OUTLINECOLOR 0 0 0 END END END END Seth -- web:https://geographika.net & https://mapserverstudio.net mastodon: @geographika at mastodon.social On Thu, Oct 9, 2025, at 2:47 AM, Robert Sanson via MapServer-users wrote: Hi All Does anyone have an example of setting up a layer in a Mapfile to cascade an ArcGIS FeatureService REST endpoint? Many thanks Robert Sanson New Zealand This email, including any attachments, may contain information which is confidential or subject to copyright. If you are not the intended recipient, please notify us immediately and delete this email from your system. Email communications are not secure and are not guaranteed by AsureQuality Limited to be free of unauthorised interference, error or virus. Anyone who communicates with us by email is taken to accept this risk. Anything in this email which does not relate to the official business of AsureQuality Limited is neither given nor endorsed by AsureQuality Limited. Please refer to www.asurequality.com or contact AsureQuality Limited for more information, including about our terms of business. _______________________________________________ MapServer-users mailing list MapServer-users at lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/mapserver-users This email, including any attachments, may contain information which is confidential or subject to copyright. If you are not the intended recipient, please notify us immediately and delete this email from your system. Email communications are not secure and are not guaranteed by AsureQuality Limited to be free of unauthorised interference, error or virus. Anyone who communicates with us by email is taken to accept this risk. Anything in this email which does not relate to the official business of AsureQuality Limited is neither given nor endorsed by AsureQuality Limited. Please refer to www.asurequality.com or contact AsureQuality Limited for more information, including about our terms of business. _______________________________________________ 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: From alistaire at met.no Tue Oct 14 05:28:49 2025 From: alistaire at met.no (Alistair Everett) Date: Tue, 14 Oct 2025 14:28:49 +0200 Subject: [MapServer-users] UVRASTER reprojection In-Reply-To: <74b6adef-6a57-47c3-9b62-e6e11f17aa41@app.fastmail.com> References: <74b6adef-6a57-47c3-9b62-e6e11f17aa41@app.fastmail.com> Message-ID: Hi Seth, Thanks for taking a look. I've had a look at some of the existing test cases and managed to make some tests which show the issue more clearly. I can open an issue and link to a branch with these updated tests, I'll try to do this later today or tomorrow. I've written an algorithm in python which we will use in the mean time to pre-calculate the corrected angles. I'm more than happy to share that code if someone would like to try to port it to mapserver. I can also take a look if I have time, but have to admit my C++ isn't very up to scratch. To avoid taking a lot of time calculating on the fly it would also be smart to cache the point-wise rotation angle used to convert between different projections, and I'm unsure how this would be handled in mapserver. Thanks, Alistair On Sat, 11 Oct 2025 at 12:01, Seth G wrote: > Hi Alistair, > > The Mapfile and requests look fine to me, so this is likely an oversight > in the implementation. > It would be good to get this logged as an issue, ideally with a > recreatable test case at https://github.com/MapServer/MapServer/issues. > > Seth > > -- > web:https://geographika.net & https://mapserverstudio.net > mastodon: @geographika at mastodon.social > > On Mon, Oct 6, 2025, at 12:16 PM, Alistair Everett wrote: > > Hi Seth, > > Thanks for your reply, we're just reprojecting the data. I'll take a > closer look at the test case too, since it seems like that should be > reprojecting from stereographic to geographic. I should also be able to > pull together some test data in the next few days too, I can open an issue > with it when I've put it together. > > Here are a couple of example requests: > > One in EPSG:4326 > > localhost:8000/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.3.0&FORMAT=image%2Fpng&STYLES=wind_barb&TRANSPARENT=TRUE&TIME=2025-10-06T07%3A00%3A00Z&DIM_PRESSURE=50.0&DIM_REFERENCE_TIME=2025-10-06T03%3A00%3A00Z&LAYERS=wind_pl_vector_meps_det_vdiv_2_5km_calculations&WIDTH=1304&HEIGHT=1504&CRS=EPSG%3A3857&BBOX=-2460529.2703191666%2C4596419.518351211%2C6809846.057963322%2C15288631.55342549 > > https://pasteboard.co/jI1eMKaj3PT1.png > > This is one in EPSG:3411: > > localhost:8000/wms?REQUEST=GetMap&SERVICE=WMS&VERSION=1.3.0&FORMAT=image%2Fpng&STYLES=wind_barb&TRANSPARENT=TRUE&TIME=2025-10-06T16%3A00%3A00Z&DIM_PRESSURE=50.0&DIM_REFERENCE_TIME=2025-10-06T03%3A00%3A00Z&LAYERS=wind_pl_vector_meps_det_vdiv_2_5km_calculations&WIDTH=1438&HEIGHT=1504&CRS=EPSG%3A3411&BBOX=404935.5094684181%2C-3769428.4764078585%2C4977748.7994076125%2C1013263.5876175575 > > https://pasteboard.co/IHZjgXibtLyJ.png > > Thanks, > Alistair > > > On Sun, 5 Oct 2025 at 10:00, Seth G wrote: > > > Hi Alistair, > > Could you give an example request (I presume it is WMS?)? Are you rotating > your map at all, or is it just the projection? > > There are some tests which reproject a UVRASTER LAYER at > https://github.com/MapServer/MapServer/blob/29e4ae81620a2c50201101c763df0b9d8d540b67/msautotest/wxs/wms_uvraster_map_reprojection.map > It sounds like an oversight, so if you are able to create a standalone > test case, it would be worth opening an issue. > > Seth > > -- > web:https://geographika.net & https://mapserverstudio.net > mastodon: @geographika at mastodon.social > > On Fri, Oct 3, 2025, at 11:59 AM, Alistair Everett via MapServer-users > wrote: > > Hi all, > > I'm trying to use the UVRASTER connection to plot wind barbs. I noticed if > I request the layer in a different projection the location of the arrows > are reprojected but the angle, and presumably vector components, are not - > rather they remain in the same direction relative to my screen. It's a > little unclear from the documentation if this is the expected behaviour for > UVRASTER or not? I included a pasteboard link to a screenshot to show what > I mean, compare for example vectors in the corners of the data and you will > see they point in the same direction relative to my screen, rather than > being rotated to match the projection. > > If this is expected, does anyone have any suggested ways of working around > this? > > If this is not expected, does anyone have any ideas what I've done wrong? > I can provide some sample data if someone wants to try to reproduce. > > I'm using: MapServer version 8.4.0 PROJ version 9.4 GDAL version 3.11 > OUTPUT=PNG OUTPUT=JPEG SUPPORTS=PROJ SUPPORTS=AGG SUPPORTS=FREETYPE > SUPPORTS=CAIRO SUPPORTS=ICONV SUPPORTS=FRIBIDI SUPPORTS=WMS_SERVER > SUPPORTS=WMS_CLIENT SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT > SUPPORTS=WCS_SERVER SUPPORTS=SOS_SERVER SUPPORTS=OGCAPI_SERVER > SUPPORTS=FASTCGI SUPPORTS=THREADS SUPPORTS=GEOS SUPPORTS=PBF INPUT=JPEG > INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE INPUT=FLATGEOBUF > > Screenshot: https://pasteboard.co/WTgbexDmQaAe.png > > Example mapfile: https://pastebin.com/ibPc6gVK > > Thanks, > Alistair > > ----------------------------- > Dr Alistair Everett > Senioringeni?r / Senior Software Engineer > Avdeling for Geoutvikling > > Phone: +47 939 68 985 > Email: alistair.everett at met.no > > Meteorologisk Institutt, > Henrik Mohns plass 1, > 0313 Oslo > ----------------------------- > *Dr Alistair Everett* > Senioringeni?r / Senior Software Engineer > Avdeling for Geoutvikling > > *Phone:* +47 939 68 985 > *Email:* alistair.everett at met.no > > *Meteorologisk Institutt*, > Henrik Mohns plass 1, > 0313 Oslo > > _______________________________________________ > 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: From pe_lord at protonmail.com Wed Oct 15 13:28:50 2025 From: pe_lord at protonmail.com (PE Lord) Date: Wed, 15 Oct 2025 20:28:50 +0000 Subject: [MapServer-users] Verbosity of PROJ_DEBUG ON In-Reply-To: References: <2_hvmcRXxkquQQ6NfLa2B9pmXoHfKsJE7me4-KjAcAFX0YJ85BYt6IbnU72l3WqRrp1rtsLdsIjOEQJuoRLMheT1pJZAlwKcuPWz43-yGRg=@protonmail.com> Message-ID: I'm migrating to Mapserver 8.4.1 + proj 9.6 I've tried your solution without any success. Here my debug config: CONFIG "CPL_DEBUG" "ON" CONFIG "PROJ_DEBUG" "5" CONFIG "CPL_CURL_VERBOSE" "ON" CONFIG "CPL_TIMESTAMP" "ON" CONFIG "MS_ERRORFILE" "stderr" Here the output mod_fcgid: stderr: msDrawRasterLayerLow(MX9-01_20240730_000784-001_0030): entering. mod_fcgid: stderr: msDrawRasterLayerLow(MX9-01_20240730_000784-001_0030): Filename is: .../MX9-01_20240730_000784-001_0030.tif mod_fcgid: stderr: msDrawRasterLayerLow(MX9-01_20240730_000784-001_0030): Path is: .../MX9-01_20240730_000784-001_0030.tif mod_fcgid: stderr: msResampleGDALToMap in effect: cellsize = 0.223749 mod_fcgid: stderr: msDrawGDAL(MX9-01_20240730_000784-001_0030): using RAW_WINDOW=0 0 2101 2851, dst=0,0,939,1274 mod_fcgid: stderr: msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,0,0,0 mod_fcgid: stderr: msDrawMap(): Layer 0 (MX9-01_20240730_000784-001_0030), 0.027s mod_fcgid: stderr: freeLayer(): freeing layer at 0x13d42d0. mod_fcgid: stderr: freeLayer(): freeing layer at 0x1abd450.mod_fcgid: stderr: freeLayer(): freeing layer at 0x1b12e40. Le dimanche 5 octobre 2025 ? 03:53, Seth G a ?crit : > Hi, > > You have all the debug options on for MapServer that I am aware of. Maybe you can try CONFIG "PROJ_DEBUG" "5" to see if that is more verbose? > > https://proj.org/en/stable/usage/environmentvars.html#envvar-PROJ_DEBUG > > MapServer passes this value on as an environment variable. Maybe something has changed as the MapServer docs always use "ON" and "OFF" whereas PROJ states: > >> Starting with PROJ 9.3, ON can be used as an alias for 2, and OFF as an alias for 1. > > Seth > > -- > web:https://geographika.net & https://mapserverstudio.net > mastodon: @geographika at mastodon.social > > On Thu, Oct 2, 2025, at 3:20 PM, PE Lord via MapServer-users wrote: > >> I have Mapserver 7.6.2 compiled with proj 7.2.1 and gdal version (?). Sorry. >> >> I have solved an PROJ issue on my current server. I have modified the PROJ_DATA (PROJ_LIB) sqlite database to handle properly my local projection, including a Canadian NTv2 grid transformation. >> >> My ...proj/bin/projinfo, on my server seem to use the right Canadian nadgrid file (ca_nrc_NA83SCRS.tif). >> >> Here my question, how to debug PROJ, through Mapserver will a high level of verbose? >> >> Here my mapfile >> >> CONFIG "CPL_DEBUG" "ON" >> CONFIG "PROJ_DEBUG" "ON" >> CONFIG "CPL_CURL_VERBOSE" "ON" >> CONFIG "CPL_TIMESTAMP" "ON" >> CONFIG "MS_ERRORFILE" "debug_proj.log" >> >> ########################### >> >> LAYER >> DEBUG 5 >> NAME "MX9-01_20240730_000784-001_0030" >> DATA "[...path...]/MX9-01_20240730_000784-001_0030.tif" >> TYPE RASTER >> PROJECTION >> "init=epsg:3799" >> END >> END >> >> ########################### >> Here the log output: >> >> [Wed Aug 27 10:13:25 2025].225353 msDrawRasterLayerLow(MX9-01_20240730_000784-001_0030): entering. >> [Wed Aug 27 10:13:25 2025].225880 msDrawRasterLayerGDAL(): Entering transform. >> [Wed Aug 27 10:13:25 2025].225894 msDrawRasterLayerGDAL(): src=1534,1715,419,296, dst=0,0,1300,919 >> [Wed Aug 27 10:13:25 2025].225899 msDrawRasterLayerGDAL(): source raster PL (1533.977,1715.259) for dst PL (0,0). >> [Wed Aug 27 10:13:25 2025].225910 msDrawRasterLayerGDAL(): red,green,blue,alpha bands = 1,0,0,0 >> [Wed Aug 27 10:13:25 2025].251631 msDrawMap(): Layer 0 (MX9-01_20240730_000784-001_0030), 0.026s >> [Wed Aug 27 10:13:25 2025].327641 freeLayer(): freeing layer at 0x2727a80. >> [Wed Aug 27 10:13:25 2025].327746 freeLayer(): freeing layer at 0x27a4850. >> [Wed Aug 27 10:13:27 2025].569803 freeLayer(): freeing layer at 0x3e92d90. >> [Wed Aug 27 10:13:27 2025].569862 freeLayer(): freeing layer at 0x2660f60. >> >> _______________________________________________ >> 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: From karsten at terragis.net Thu Oct 23 12:11:30 2025 From: karsten at terragis.net (Karsten Vennemann) Date: Thu, 23 Oct 2025 12:11:30 -0700 Subject: [MapServer-users] CUGOS Fall Fling - registration is now open! In-Reply-To: References: Message-ID: <30fda3d8-9ef2-484d-83d3-ad5e2759e124@terragis.net> Hi FOSS geospatial aficionados, The wait is finally over ? registration is now open for the *CUGOS Fall Fling 2025!* ? Come join us for a full day of geospatial fun, talks, and community connection in Seattle: ? *Date:* Friday, November 14, 2025 ? 9:00 AM ? 5:00 PM ? *Location:* Maple Hall Great Room, University of Washington, Seattle, WA 98105 ? *Details:* cugos.org/2025-fall-fling ? *Registration:* Eventbrite ? CUGOS Fall Fling 2025 We?re looking forward to seeing you there for a great day of open geospatial exchange and community spirit! Cheers, *Karsten Vennemann* CUGOS Board Member ___________________________________________ Subject: CUGOS Fall Fling - reminder about call for presentations ! Date: Mon, 22 Sep 2025 11:59:05 +0200 From: Karsten Vennemann To: osgeo-bc Hi GIS and FOSS enthusiasts, This is a friendly reminder about the *CUGOS Fall Fling 2025!* ?, details see below (more to come). The deadline for submitting a talk is approaching and is only a little over a week away, so don't delay and submit your favorite talk in regards to open source geospatial! *Call for Proposals* Submit your talk ideas by *Tuesday, September 30th, 5:00 pm*. We look forward to your contributions and to another day filled with geospatial energy, learning, and collaboration! ? *Date:* Friday, November 14, 2025 ?*Location:* University of Washington, Seattle, WA 98105 ? *Details:* https://cugos.org/2025-fall-fling Registration will open once the agenda is finalized. Cheers, Karsten Vennemann CUGOS Board Member ------------------- /(For those new to CUGOS or curious about what the ?Fall Fling? is, read on below.)/ The *Cascadia Users of Geospatial Open Source (CUGOS)* was founded in 2007 and serves as the Cascadia Chapter of the OSGeo Foundation .?We are a non-profit community of open source GIS enthusiasts ? from students to professionals, from academia to industry, and including active OSGeo contributors (board members, charter members, and project participants). CUGOS provides a space to share work, exchange ideas, and celebrate the power of open source geospatial tools. Learn more: https://cugos.org/about/ This year?s *Fall Fling* will be a one-day, low-cost event bringing together academics, professionals, students, and enthusiasts for a full day of presentations on open source geospatial tools and open data in and beyond the Puget Sound region. Talks will cover projects like QGIS, PostGIS, MapServer, GDAL/OGR, OpenLayers, and more (see also: OSGeo projects ?) While we?re still building the program, you can get a taste of what to expect by checking out schedules from our previous Fling events: * 2015 https://cugos.org/2015-spring-fling/ * 2016 https://cugos.org/2016-spring-fling/ * 2017 Spring Fling * 2019 Fall Fling * 2023 Spring Fling -- Karsten Vennemann *Terra GIS LTD * 2119 Boyer Ave E Seattle, WA 98112 Tel ++1 206 905 1711 _www.terragis.net _ -- Karsten Vennemann *Terra GIS LTD * 2119 Boyer Ave E Seattle, WA 98112 Tel ++1 206 905 1711 _www.terragis.net _ -------------- next part -------------- An HTML attachment was scrubbed... URL: From trygve at aspenes.priv.no Wed Oct 29 06:04:55 2025 From: trygve at aspenes.priv.no (Trygve Aspenes) Date: Wed, 29 Oct 2025 14:04:55 +0100 Subject: [MapServer-users] uvraster using vrt files from tileitem fails: GDALRasterIO() failed In-Reply-To: References: Message-ID: Sorry Somehow I managed to send the email before I had completed it. Here is all again > Hello > > I have a layer using uvraster as connectiontype, with the data comming > from a helper layer via tileitem. > > So the helper layer is like this: > > LAYER > PROJECTION > "init=epsg:3857" > END > STATUS ON > NAME "time_idx_vrt" > TYPE POLYGON > PROCESSING "CLOSE_CONNECTION=DEFER" > DATA "geom from (select * from ) as foo using unique id using > srid=3857" > METADATA > "wms_title" "x_timeindex" > "wms_timeextent" "
) as foo using unique id using srid=3857" METADATA "wms_title" "x_timeindex" "wms_timeextent" "