我正在尝试关注Mouse Connectivity sdk并让他们的两个例子正常运作。
pd
返回无或全部投影信号密度实验。为什么会这样?
from allensdk.api.queries.mouse_connectivity_api import MouseConnectivityApi
mca = MouseConnectivityApi()
# get metadata for all non-Cre experiments
experiments = mca.experiment_source_search(injection_structures='root', transgenic_lines=0)
# download the projection density volume for one of the experiments
#pd = mca.download_projection_density('example.nrrd', experiments[0]['id'], resolution=25)
for exp in range(len(experiments)):
pd = mca.download_projection_density('example.nrrd', experiments[exp]['id'], resolution=25)
print(type(pd))
结果:
C:\Anaconda\python.exe C:/Users/user/PycharmProjects/MyFirstAllenBrain/Test.py
<type 'NoneType'>
<type 'NoneType'>
<type 'NoneType'>
<type 'NoneType'>
<type 'NoneType'>
<type 'NoneType'>
... etc
事情是,experiments
确实收到了一个值,因此似乎是MouseConnectivityApi和nrrd(我根据以下post安装)正常工作的情况。< / p>
是的,现在是第二个例子
from allensdk.core.mouse_connectivity_cache import MouseConnectivityCache
# tell the cache class what resolution (in microns) of data you want to download
mcc = MouseConnectivityCache(resolution=25)
# use the ontology class to get the id of the isocortex structure
ontology = mcc.get_ontology()
isocortex = ontology['Isocortex']
# a list of dictionaries containing metadata for non-Cre experiments
experiments = mcc.get_experiments(injection_structure_ids=isocortex['id'])
# download the projection density volume for one of the experiments
pd = mcc.get_projection_density(experiments[0]['id'])
这是从艾伦一字不差地复制的,但这是我得到的错误信息:
C:\Anaconda\python.exe C:/Users/user/PycharmProjects/MyFirstAllenBrain/Test.py
Traceback (most recent call last):
File "C:/Users/user/PycharmProjects/MyFirstAllenBrain/Test.py", line 14, in <module>
pd = mcc.get_projection_density(experiments[0]['id'])
File "C:\Users\user\AppData\Roaming\Python\Python27\site-packages\allensdk\core\mouse_connectivity_cache.py", line 170, in get_projection_density
return nrrd.read(file_name)
File "C:\Anaconda\lib\site-packages\nrrd.py", line 384, in read
data = read_data(header, filehandle, filename)
File "C:\Anaconda\lib\site-packages\nrrd.py", line 235, in read_data
mmap.MAP_PRIVATE, mmap.PROT_READ)
AttributeError: 'module' object has no attribute 'MAP_PRIVATE'
Process finished with exit code 1
为什么会发生这种情况?
再次像以前一样(不需要我假设的另一张图片),experiments
变量确实收到了每个实验的值。
答案 0 :(得分:0)
不幸的是,这是pynrrd / master github存储库的Windows特定问题。我知道一个特定的修订工作:
https://github.com/mhe/pynrrd/commit/3c0f3d577b0b435fb4825c14820322a574311af0
要从Windows命令提示符安装此修订版,您可以:
> git clone https://github.com/mhe/pynrrd.git
> cd pynrrd
> git checkout 3c0f3d5
> cd ..
> pip install --upgrade pynrrd\
最后的反斜杠非常重要。它告诉pip从本地路径安装而不是检查PyPI。