[mapserver-users] GetFetureInfo with rotated grid

[Jeff McKenna]

Hi LP,

Can you also confirm through a shp2img command, that you generate a 
valid map image, in your output EPSG:4326 projection, such as:

   shp2img -m test.map -o ttt.png -map_debug 3

(also, please add "STATUS ON" to your layer)

I would always start with a shp2img command, make sure that MapServer 
(and GDAL+PROJ in this case) can access your data, reproject it, and 
return a nice valid image.  Then, after that works, move onto the fun of 
WMS and OGC service testing.


-jeff



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



On 2021-07-09 12:13 p.m., Rousseau Lambert,Louis-Philippe (ECCC) wrote:
> Hi all,
> 
> I have an issue with WMS GetFeatureInfo request on a rotated GRIB2 file. 
> When I try to do a GetFeatureInfo request on a valid location on my map, 
> I get: “Search returned no results”. I tried the GetFeatureInfo request, 
> using mapscript and the mapserv command, both give me the “Search 
> returned no results”.
> 
> The file projection is “+proj=ob_tran +o_proj=longlat +o_lon_p=0 
> +o_lat_p=31.758312 +lon_0=-92.402969 +R=6371229 +no_defs”
> 
> I’m using MapServer version 7.6.2, GDAL 3.1.3 (released 2020/09/01) and 
> proj 7.2.0.
> 
> But If I try to get the value with gdallocationinfo using the –wsg84 
> flag, I get a valid result, so I guess this issue is not with gdal ?:
> 
> ```
> 
> gdallocationinfo -wgs84 rotated.grib2 -75 45
> 
> Report:
> 
>    Location: (578P,373L)
> 
>    Band 1:
> 
>      Value: 3.1022446478346e-07
> 
> ```
> 
> I can get a GetFeatureInfo request to work, but I have to define a 
> custom EPSG code (for example: EPSG:102987), add this EPSG code (102987) 
> to my epsg file and reference it in my mapfile. This way, the 
> GetFeatureInfo only works if requested in this custom EPSG using the 
> rotated grid BBOX, if I do the request in another CRS, it still return 
> no results:
> 
> ```
> 
> mapserv -nh "QUERY_STRING=map= 
> test.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetFeatureInfo&BBOX=0,0,1,1&CRS=EPSG:102987&WIDTH=1000&HEIGHT=1000&LAYERS=TEST&INFO_FORMAT=text/plain&QUERY_LAYERS=TEST&I=5&J=5&FORMAT=image/png"
> 
> ```
> 
> Here is a GetFeatureInfo request that returns no results using EPSG:4326:
> 
> ```
> 
> mapserv -nh 
> "QUERY_STRING=map=test.map&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetFeatureInfo&BBOX=45,-74,46,-73&CRS=EPSG:4326&WIDTH=1000&HEIGHT=1000&LAYERS=TEST&INFO_FORMAT=text/plain&QUERY_LAYERS=TEST&I=5&J=5"
> 
> ```
> 
> I have included a small mapfile that shows the previous behavior, you 
> can use this curl command to download a sample GRIB2 data:
> 
> ```
> 
> curl -o rotated.grib2 
> 'https://dd.weather.gc.ca/model_raqdps/10km/grib2/00/013/20210709T00Z_MSC_RAQDPS_PM2.5_EATM_RLatLon0.09x0.09_P013.grib2'
> 
> ```
> 
> Do you have any ideas if I’m doing something wrong or if the issue is in 
> MapServer ?
> 
> Thanks
> 
> LP
> 
> 
> _______________________________________________
> mapserver-users mailing list
> mapserver-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/mapserver-users
>