[pymvpa] atlas coordinate system

[Lisa S]

I was hoping to use the atlas voxel labeler functionality in
mvpa.atlases. I installed lxml, imported mvpa.atlases, and loaded
FSL's Harvard Oxford Cortical Atlas using atlas =
Atlas("HarvardOxford-Cortical.xml").

Then I used did the following:

>>> atlas.labelVoxel([28,82,55],0)
{'labels': [[{'index': 15, 'prob': 58, 'label': 'Inferior Temporal
Gyrus, temporooccipital part'}, {'index': 12, 'prob': 15, 'label':
'Middle Temporal Gyrus, temporooccipital part'}, {'index': 14, 'prob':
4, 'label': 'Inferior Temporal Gyrus, posterior division'}]],
'voxel_queried': [28, 82, 55]}

However, the label for the voxel (28,82,55) does not seem to mach up
with the labels when I look up the same coordinates in the atlas in
FSLView:

Harvard-Oxford Cortical Structural Atlas (28, 82, 55)
50% Frontal Pole, 23% Middle Frontal Gyrus, 2% Superior Frontal Gyrus

In fact, many of the other tuples that I have tried using labelVoxel
have come out with no label. Does anyone have an idea of how the
coordinate system for the labelAtlas call is oriented that might be
the cause of these differences?