Python single_dipole_search 예제들

예제 #1

파일: searchtest.py 프로젝트: CyclotronResearchCentre/forward

import os
import os.path as op
from forward.search import single_dipole_search
import numpy as np
import time

potential_file = op.join(os.environ["FWD_DIR"], "examples", "sphere_cost_datasink", "potential_from_dipole", "potential.npy")
leadfield_file = op.join(os.environ["FWD_DIR"], "examples", "sphere_datasink", "leadfield", "leadfield.hdf5")
mesh_file = op.join(os.environ["FWD_DIR"], "examples", "sphere_datasink", "mesh_file", "4shell_elec.msh")

elements_to_consider = [1001]
electrode_potential = np.load(potential_file)
start_time = time.time()
xopt, qopt, element_data, geo_file = single_dipole_search(electrode_potential, leadfield_file, mesh_file, elements_to_consider)
end_time = time.time()

print("Finished in %f seconds" % (end_time - start_time))
print(xopt)
print(qopt)
print(geo_file)

예제 #2

import os
import os.path as op
from forward.search import single_dipole_search
import numpy as np
import time

potential_file = op.join(os.environ["FWD_DIR"], "examples",
                         "sphere_cost_datasink", "potential_from_dipole",
                         "potential.npy")
leadfield_file = op.join(os.environ["FWD_DIR"], "examples", "sphere_datasink",
                         "leadfield", "leadfield.hdf5")
mesh_file = op.join(os.environ["FWD_DIR"], "examples", "sphere_datasink",
                    "mesh_file", "4shell_elec.msh")

elements_to_consider = [1001]
electrode_potential = np.load(potential_file)
start_time = time.time()
xopt, qopt, element_data, geo_file = single_dipole_search(
    electrode_potential, leadfield_file, mesh_file, elements_to_consider)
end_time = time.time()

print("Finished in %f seconds" % (end_time - start_time))
print(xopt)
print(qopt)
print(geo_file)

예제 #3

파일: brainsearch.py 프로젝트: CyclotronResearchCentre/forward

import os
import os.path as op
from forward.search import single_dipole_search
import numpy as np

from forward.datasets import sample
data_path = sample.data_path(name="leadfield")


conductivity_tensor_included = False
potential_file = op.join(os.environ["FWD_DIR"], "examples", 'forward_datasink','subject','potential_from_dipole','_subject_id_TMS007','potential.npy')

if conductivity_tensor_included:
    leadfield = op.join(data_path, "TMS007_aniso_leadfield.hdf5")
else:
    leadfield = op.join(data_path, "TMS007_iso_leadfield.hdf5")

mesh_file = op.join(os.environ["FWD_DIR"], "examples", "sphere_cost_datasink", "cost_mesh_file", "4shell_elec_dipole_cost.msh")

elements_to_consider = [1002]
electrode_potential = np.load(potential_file)
element_idx, geo_file = single_dipole_search(electrode_potential, leadfield, mesh_file, elements_to_consider)

예제 #4

import os
import os.path as op
from forward.search import single_dipole_search
import numpy as np

from forward.datasets import sample
data_path = sample.data_path(name="leadfield")

conductivity_tensor_included = False
potential_file = op.join(os.environ["FWD_DIR"], "examples", 'forward_datasink',
                         'subject', 'potential_from_dipole',
                         '_subject_id_TMS007', 'potential.npy')

if conductivity_tensor_included:
    leadfield = op.join(data_path, "TMS007_aniso_leadfield.hdf5")
else:
    leadfield = op.join(data_path, "TMS007_iso_leadfield.hdf5")

mesh_file = op.join(os.environ["FWD_DIR"], "examples", "sphere_cost_datasink",
                    "cost_mesh_file", "4shell_elec_dipole_cost.msh")

elements_to_consider = [1002]
electrode_potential = np.load(potential_file)
element_idx, geo_file = single_dipole_search(electrode_potential, leadfield,
                                             mesh_file, elements_to_consider)