Hi Felix,
Thanks for including all of the files.
I ran this locally, with MS4W 5.0 beta2 (MapServer 8.0-dev, GDAL 3.4.0,
PROJ 8.2.0) using the map2img commandline utility (or "shp2img" for
earlier MapServer versions) and received what looks like a valid map
image each time:
- map image: https://pasteboard.co/2mTyHAzprtQk.png
- command:
map2img -m uvraster.map -e -154 26 -42 70 -l "HRDPS.CONTINENTAL_UU"
-o ttt.png -map_debug 3
- command's debug output: https://pastebin.com/qPsHr4wJ
My output appears valid, I think. I wonder if you are having a PROJ
issue (try upgrading to MapServer-main and PROJ 8, and re-test)
Hope this somewhat helps,
Have a nice weekend, and hello from the east coast.
-jeff
--
Jeff McKenna
GatewayGeo: Developers of MS4W, MapServer Consulting and Training
co-founder of FOSS4G
http://gatewaygeo.com/
On 2022-01-21 3:08 p.m., Laframboise,Felix (ECCC) wrote:
Hi,
I am currently using GDAL 3.1.3 and Mapserver 7.6.2 on a Ubuntu machine.
I created a UVRaster layer with a style using arrow symbols. The source
data is GRIB2 files using a Rotated Latitude-Longitude grid data (RLatLon).
When I tested the layer, I noticed that it does not output correctly.
Sometimes, no image outputs, and sometimes all the arrows/data seem to
be clustered somewhere on the extent. The expected extent should be the
one of North America (approximatively [-154, 26, -42, 70] in the
/[min_x, min_y, max_x, max_y]/ format), where all the arrows are spread
out on this extent.
Following this output, I did investigate the issue, and the mapfile is
working fine as I tested it with other data that is not using a RLatLon
grid.
I created a test case to show you the issue here:
·Sample mapfile: https://pastebin.com/h06AsJzz
<https://pastebin.com/h06AsJzz>
·Sample VRT: https://pastebin.com/ERViMRpR
·2 UV GRIB2 files (available here:
https://dd.weather.gc.ca/model_hrdps/continental/2.5km/grib2/00/000/
<https://dd.weather.gc.ca/model_hrdps/continental/2.5km/grib2/00/000/> )
for the source data (using following patterns) :
o/*_MSC_HRDPS_UGRD_AGL-40m_RLatLon0.0225_*.grib2/
o/*_MSC_HRDPS_VGRD_AGL-40m_RLatLon0.0225_*.grib2/
·Output image when arrows are clustered:
https://pasteboard.co/WMoun7bhLZSl.png
<https://pasteboard.co/WMoun7bhLZSl.png>
We serve the mapfiles using Python mapscript but I tested the mapfiles
with a mapserv command and got the same result.
Do you have an idea what could be the issue and how to solve it? RLatLon
data is working fine with Mapserver, but it seems to only affect
UVraster layers.
Thanks a lot in advance for your answer(s)!
Félix Laframboise
_______________________________________________
MapServer-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/mapserver-dev
_______________________________________________
MapServer-dev mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/mapserver-dev
