problem with wms getfeatureinfo request and time indexed files

[fred]

Hello, I am facing a possible bug whith a WMS GetFeatureInfo? request. The accessed data are time dependant images. As a consequence, I use an index by setting the TILEINDEX attribute in the mapfile. When using such an index whith the above mentionned request, MapServer returns as much results (all the results are the same) as images referenced in this index containing the requested geographic point instead of only one (only one image is asked). I attach to this ticket the mapfile and the index. The images can be downloaded from : -  http://rapidfire.sci.gsfc.nasa.gov/subsets/?subset=AERONET_Skukuza.2010245.terra.2km.tif -  http://rapidfire.sci.gsfc.nasa.gov/subsets/?subset=AERONET_Wits_University.2010244.terra.2km.tif

To reproduce this, please copy these files in the same directory and try the following (what I call the first image is the one of 2010/09/02 ie AERONET_Skukuza.2010245.terra.2km.tif, the second one is the one of 2010/09/01 ie AERONET_Wits_University.2010244.terra.2km.tif) :

GetMap? request with the first image -> OK

http://PathToMapserver/mapserv?map=PathToMapfile/bug.map&TIME=2010-09-02T00:00:00Z&SRS=EPSG:4326&STYLES=&LAYERS=terra&TRANSPARENT=true&FORMAT=image/png&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-35,10,-20,40&WIDTH=512&HEIGHT=256

#======================================

GetMap? request with the second image -> OK

http://PathToMapserver/mapserv?map=PathToMapfile/bug.map&TIME=2010-09-01T00:00:00Z&SRS=EPSG:4326&STYLES=&LAYERS=terra&TRANSPARENT=true&FORMAT=image/png&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&BBOX=-35,10,-20,40&WIDTH=512&HEIGHT=256

#======================================

GetFeatureInfo? whith the first image (this point is in the first image, not in the second) -> 1 result, OK

http://PathToMapserver/mapserv?map=PathToMapfile/bug.map&TIME=2010-09-02T00:00:00Z&SRS=EPSG:4326&STYLES=&LAYERS=terra&TRANSPARENT=true&FORMAT=image/png&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetFeatureInfo&BBOX=-35,10,-20,40&WIDTH=512&HEIGHT=256&X=300&Y=150&INFO_FORMAT=text/plain&QUERY_LAYERS=terra

Output :

GetFeatureInfo results:

Layer 'terra'
  Feature 0: 
    x = '27.599056'
    y = '-28.810878'
    value_0 = '119'
    value_1 = '96'
    value_2 = '55'
    value_list = '119,96,55'
    red = '119'
    green = '96'
    blue = '55'

#======================================

GetFeatureInfo? whith the second image and the same point as above (this point is in the first image, not in the second) -> 1 result, PROBLEM

http://PathToMapserver/mapserv?map=PathToMapfile/bug.map&TIME=2010-09-01T00:00:00Z&SRS=EPSG:4326&STYLES=&LAYERS=terra&TRANSPARENT=true&FORMAT=image/png&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetFeatureInfo&BBOX=-35,10,-20,40&WIDTH=512&HEIGHT=256&X=300&Y=150&INFO_FORMAT=text/plain&QUERY_LAYERS=terra

Output, same as above :

GetFeatureInfo results:

Layer 'terra'
  Feature 0: 
    x = '27.599056'
    y = '-28.810878'
    value_0 = '119'
    value_1 = '96'
    value_2 = '55'
    value_list = '119,96,55'
    red = '119'
    green = '96'
    blue = '55'

#======================================

GetFeatureInfo? whith the first image (this point is in the other image) -> 2 results, PROBLEM

http://PathToMapserver/mapserv?map=PathToMapfile/bug.map&TIME=2010-09-02T00:00:00Z&SRS=EPSG:4326&STYLES=&LAYERS=terra&TRANSPARENT=true&FORMAT=image/png&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetFeatureInfo&BBOX=-35,10,-20,40&WIDTH=512&HEIGHT=256&X=300&Y=130&INFO_FORMAT=text/plain&QUERY_LAYERS=terra

Output :

GetFeatureInfo results:

Layer 'terra'
  Feature 0: 
    x = '27.599056'
    y = '-27.641763'
    value_0 = '127'
    value_1 = '103'
    value_2 = '59'
    value_list = '127,103,59'
    red = '127'
    green = '103'
    blue = '59'
  Feature 0: 
    x = '27.599056'
    y = '-27.641763'
    value_0 = '127'
    value_1 = '103'
    value_2 = '59'
    value_list = '127,103,59'
    red = '127'
    green = '103'
    blue = '59'

#======================================

GetFeatureInfo? whith the second image (this point is in the other image) -> 2 results, same as above, PROBLEM

http://PathToMapserver/mapserv?map=PathToMapfile/bug.map&TIME=2010-09-01T00:00:00Z&SRS=EPSG:4326&STYLES=&LAYERS=terra&TRANSPARENT=true&FORMAT=image/png&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetFeatureInfo&BBOX=-35,10,-20,40&WIDTH=512&HEIGHT=256&X=300&Y=130&INFO_FORMAT=text/plain&QUERY_LAYERS=terra

Output, same as above :

GetFeatureInfo results:

Layer 'terra'
  Feature 0: 
    x = '27.599056'
    y = '-27.641763'
    value_0 = '127'
    value_1 = '103'
    value_2 = '59'
    value_list = '127,103,59'
    red = '127'
    green = '103'
    blue = '59'
  Feature 0: 
    x = '27.599056'
    y = '-27.641763'
    value_0 = '127'
    value_1 = '103'
    value_2 = '59'
    value_list = '127,103,59'
    red = '127'
    green = '103'
    blue = '59'

[fred]

Replying to sdlime: Hello, The ticket has not evolved since this kind of acknowledgement of receipt. Is it something usual or have I been wrong somewhere ? Thanks for your answer

[dmorissette]

Nothing wrong on your part, just that nobody has had time to get to this ticket yet.

Sounds like the problem is with GetFeatureInfo? on time-based *raster* layers... adding Frank to CC to see if he has an idea.

[aboudreault]

I'm going to take a look at this bug...

[aboudreault]

fred, I think this bug is fixed in svn_trunk. Could you test and confirm? You can see the different results by changing the ΅Y parameter to 130/150 and the TIME parameter. Here's what I obtained with trunk:

{ ...&TIME=2010-09-02T00:00:00Z&CRS=EPSG:4326&STYLES=&LAYERS=terra&TRANSPARENT=true&FORMAT=image/png&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetFeatureInfo&BBOX=-35,10,-20,40&WIDTH=512&HEIGHT=256&X=300&Y=130&INFO_FORMAT=text/plain&QUERY_LAYERS=terra : GetFeatureInfo? results:

Layer 'terra'

Feature 0:

x = '27.61011' y = '-27.646242' value_0 = '126' value_1 = '112' value_2 = '111' value_list = '126,112,111' red = '126' green = '112' blue = '111'

}

{ &TIME=2010-09-01T00:00:00Z&CRS=EPSG:4326&STYLES=&LAYERS=terra&TRANSPARENT=true&FORMAT=image/png&SERVICE=WMS&VERSION=1.3.0&REQUEST=GetFeatureInfo&BBOX=-35,10,-20,40&WIDTH=512&HEIGHT=256&X=300&Y=130&INFO_FORMAT=text/plain&QUERY_LAYERS=terra : GetFeatureInfo? results:

Layer 'terra'

Feature 0:

x = '27.599056' y = '-27.641763' value_0 = '127' value_1 = '103' value_2 = '59' value_list = '127,103,59' red = '127' green = '103' blue = '59'

}

[fred]

Replying to aboudreault:

Hello, I am not sure... For now, I do not have enough time to test. I will try to test before the end of the year, most probably in january. Thanks for your look.

[aboudreault]

This seems to have been fixed in trunk. Please reopen with more details if you can reproduce with a recent release.

[fred]

Indeed, it has been fixed in the trunk, you are right. Though, this fix has not been integrated in the latest release (5.6.6). Can you confirm that it have to wait for a future release ? The 6.0 one ? Have you got an idea of when it will be available ? Thanks a lot.