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: