[pymvpa] atlas coordinate system

Lisa S lstorey at gmail.com
Tue Feb 23 20:22:31 UTC 2010


Thank you, that's exactly what I was missing :)

On Tue, Feb 23, 2010 at 2:30 PM, Yaroslav Halchenko
<debian at onerussian.com> wrote:
> Hi Lisa,
>
> wow -- you found a hidden undocumented treasure...
>
> The issue in your case is that atlases are loaded at maximal available
> resolution ... ie in this case, it was 1mm, hence voxel querying provides your
> output for a voxel at 1mm resolution, not 2mm like you've assumed I guess.
>
> In [103]: print atlas._resolution
> 1.0
>
> In [104]: atlas = Atlas(name="HarvardOxford-Cortical", resolution=2.0)
>
> In [105]: atlas.labelVoxel([28,82,55],0)
> Out[105]:
> {'labels': [[{'index': 0, 'label': 'Frontal Pole', 'prob': 50},
>             {'index': 3, 'label': 'Middle Frontal Gyrus', 'prob': 23},
>                 {'index': 2, 'label': 'Superior Frontal Gyrus', 'prob': 2}]],
> 'voxel_queried': [28, 82, 55]}
>
>
> but once again -- all atlas handling is yet to be properly documented and
> exposed.  It is great if it already works for you, but be careful and verify
> the output (like you have already done ;))
>
>
> On Tue, 23 Feb 2010, Lisa S wrote:
>> 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").
>
> --
>                                  .-.
> =------------------------------   /v\  ----------------------------=
> Keep in touch                    // \\     (yoh@|www.)onerussian.com
> Yaroslav Halchenko              /(   )\               ICQ#: 60653192
>                   Linux User    ^^-^^    [175555]
>
>
>
> _______________________________________________
> Pkg-ExpPsy-PyMVPA mailing list
> Pkg-ExpPsy-PyMVPA at lists.alioth.debian.org
> http://lists.alioth.debian.org/mailman/listinfo/pkg-exppsy-pymvpa
>



More information about the Pkg-ExpPsy-PyMVPA mailing list