Link
WLCD
Search

Perform Omnipole Voltage Mapping

Table of contents

  1. Background
  2. Dependencies
  3. Quickstart
  4. Process
    1. Create omnipole dataset
    2. Visualise the data
    3. Convert the data to VTK Format
  5. Publications using this code

Background

Omnipole mapping refers to … This page describes how to create an omnipole voltage map

Dependencies

This process depends on unipole data from Carto. Specifically, the following fields must be present in the userdata structure: tbc In addition the following functions are required:

RBFConductionVelocity.m % for interpolation
connectaxes % for visualisation
drawMap % for visualisation
carto2vtk % for saving data

Quickstart

[interpOtmax, interpBip, interpUni] = run_maxOmnipoleVoltage(userdata);

Process

Create omnipole dataset

datapath = '/Users/Steven/Documents/Dropbox/Matlab_Code/workingdir/OmnipoleMapping/data/485LAiain.mat';
load(datapath)
[interpOtmax, interpBip, interpUni] = run_maxOmnipoleVoltage(userdata);

Each of the above output variables has numel() equal to size(userdata.surface.triRep.X,2). In other words, this represents point data, which can be used directly with carto2vtk.

Visualise the data

ax(1) = subplot(1,3,1);
drawMap(userdata, 'data', interpBip', 'coloraxis', [0.05 0.5]); title('Bipolar');
ax(2)= subplot(1,3,2);
drawMap(userdata, 'data', interpUni', 'coloraxis', [0.5 5]); title('Unipolar');
ax(3) = subplot(1,3,3);
drawMap(userdata, 'data', interpOtmax', 'coloraxis', [0.5 5]); title('Omnipolar');
connectaxes(ax)

Convert the data to VTK Format

See carto2vtk

datapath = '/Users/Steven/Documents/Dropbox/Matlab_Code/workingdir/OmnipoleMapping/data/485LAiain.mat';
[d,f] = fileparts(datapath);
carto2vtk(userdata, [d filesep() f 'omni.vtk'], 'manual', 'pointdata', interpOtmax);

Publications using this code

tbc