Waiting
Login processing...

Trial ends in Request Full Access Tell Your Colleague About Jove
Click here for the English version

Behavior

Topographische Schätzung von Visual Bevölkerung rezeptiven Felder von fMRI

Published: February 3, 2015 doi: 10.3791/51811

Abstract

Sehrinde wird retinotop organisiert, so dass benachbarte Populationen von Zellen abzubilden, um benachbarte Teile des Gesichtsfeldes. Funktionelle Magnetresonanz-Bildgebung ermöglicht uns voxelbasierten Population rezeptiven Felder (PRF), dh zu schätzen, der Teil des Sichtfeldes, die die Zellen innerhalb jedes Voxel aktiviert. Vor, direkt, pRF Schätzmethoden 1 leiden unter gewissen Einschränkungen: 1) die pRF Modell ist a priori gewählt und können die tatsächlichen pRF Form nicht vollständig zu erfassen, und 2) pRF Zentren sind anfällig für Fehllokalisation der Nähe der Grenze des Stimulus Raum. Hier eine neue topographische pRF Schätzmethode 2 wird vorgeschlagen, die weitgehend umgeht diese Einschränkungen. Ein lineares Modell verwendet, um den Blutsauerstoffpegel-abhängiger (BOLD) Signal durch Faltung des linearen Antwort des pRF zur visuellen Stimulus mit der kanonischen hämodynamischen Antwortfunktion vorherzusagen. PRF Topographie wird als Gewichtsvektor, dessen Komponenten die str vertretenength der Gesamtreaktion der Voxel-Neuronen zu unterschiedlichen Gesichtsfeld Standorten präsentiert Reize. Die daraus resultierenden linearen Gleichungen können für die pRF Gewichtsvektor mit Firstregressions 3, was das pRF Topographie gelöst werden. Ein pRF Modell, das dem geschätzten Topographie angepassten können dann Post-hoc gewählt werden, so die Schätzungen der pRF Parameter wie pRF Stadtlage, pRF Ausrichtung, Größe, etc. zu verbessern. Mit der pRF Topographie vorhanden ermöglicht auch die visuelle Überprüfung der pRF Parameterschätzungen ermöglicht die Extraktion von verschiedenen pRF Eigenschaften, ohne a priori Annahmen über die pRF Struktur zu machen. Dieser Ansatz verspricht besonders nützlich für die Untersuchung der pRF Organisation von Patienten mit Störungen des visuellen Systems zu sein.

Introduction

Die funktionelle Magnetresonanztomographie (fMRT) misst nicht-invasiv die funktionelle Organisation des visuellen Kortex auf makroskopischer Ebene (in der Regel in der Größenordnung von Millimetern). Frühe fMRI Retinotopie Studien eine Kohärenzmaßes zwischen Reiz Lage und löste BOLD Antworten 4-7. Diese Studien in der Regel nicht zu schätzen Bevölkerung rezeptiven Feld Größe. Später, Dumoulin und Wandell 1 ein Verfahren vorgeschlagen, eine solche Beschränkung durch die explizite Modellierung der pRF Lage und Größe, mit einer linearen Funktion des Modells, das BOLD Antwort vorhersagen zu überwinden. Allerdings ist eine Beschränkung dieser bahnbrechenden Verfahren, dass die parametrische pRF Modell muss a priori gewählt werden und kann zu fehlerhaften pRF schätzt, wenn sich herausstellt, nicht angemessen.

Um Einschränkungen des parametrischen pRF-Modell-Methode zu überwinden, wurden neue Methoden vor kurzem entwickelt worden. Diese Methoden direkt sagen die BOLD Antwort auf die stimulus durch die Rekonstruktion der PRF-Topographie. Ein von Greene und Kollegen vorgeschlagene Methode 8 rekonstruiert die pRF Topographie durch Rückprojektion der BOLD-Antworten auf die einzelnen 1D Reiz Flächen und den Bau der pRF Topographie in der 2D-Stimulus Raum wie ein typischer Computer-Tomographie-Technik. Andererseits ist die von uns vorgeschlagene Verfahren 2 direkt schätzt den 2D pRF Topographie unter Verwendung linearer Regression und Anlegen einer Regularisierung Technik. In diesem Verfahren wird die PRF Topographie als ein Satz von Gewichten, die durch den Stimulus, multipliziert wird, um den neuronalen Population Antwort eines gegebenen Voxels Abschätzung dargestellt. Dann wird der endgültige Sauerstoffgehalt des Blutes-Dependent (fett) Reaktion des Stimulus hervorgerufen durch Faltung des neuronalen Population Antwort und die kanonische hämodynamischen Antwortfunktion geschätzt. Um die unter beschränkte Linearsystem zu lösen, zusätzlich Firstregressions Regularisierung wird zur Kargheit durchzusetzen (siehe Abbildung 1unten). Die Regularisierung Technik unterdrückt Rauschen und Artefakte und ermöglicht unsere Methode, um die Topographie pRF robuster zu schätzen.

Die topographischen Methoden nicht pRF Form zu zwingen, eine bestimmte parametrische Form haben, und kann daher die tatsächliche pRF Struktur aufzudecken. Eine geeignete parametrische Modell kann dann auf der Basis der Topographie pRF gewählt werden. Beispielsweise kann die pRF Topographie verwendet, um die pRF Center- und Surround trennen, und dann können die anschließende pRF Zentrum Modellierung genauer zu sein, mit minimalem Einfluss von Surround-Unterdrückung sowie den Einfluss von anderen potenziellen Artefakte, die sich in Bereichen entfernt, um die werden pRF Zentrum. Wir haben vor kurzem durchgeführt, einen quantitativen Vergleich zwischen unseren Verfahren und verschiedene andere Verfahren, die direkt (dh vor der Schätzung der Topographie) fit isotropen Gaußschen 1, anisotrope Gaußsche und Differenz von isotropen Gauß an die pRF 9. Es wurde, dass der topogr gefundenO S-basierte Methode besser als diese Methoden im Hinblick auf pRF Zentrum Modellierung durch das Erreichen höherer erklärte Varianz des BOLD-Signal Zeitreihen.

Genaue Schätzung pRF Eigenschaften in verschiedenen Bereichen zeigt, wie sie das Sichtfeld abdecken, und ist wichtig für die Untersuchung der funktionellen Organisation des visuellen Kortex insbesondere in Bezug auf die visuelle Wahrnehmung. Eigenschaften wie, wie pRF Größenänderungen mit Exzentrizität 1,10 und pRF Surround-Center-Organisation 9 sind in der menschlichen Literatur studierte. Das vorgeschlagene Verfahren zum Schätzen der pRF Topographie führt zu genaueren pRF Parameter Modellierung und eher unbekannten Gesetzmäßigkeiten offenbaren, nicht leicht in den Direkt parametrischer Modelle modelliert a-priori. Dieser Ansatz wird besonders für das Studium pRF Organisation bei Patienten mit Sehbahnläsionen, für die pRF Struktur ist nicht unbedingt vorhersehbar a priori sein. Im Folgenden wird beschrieben, wie man schätzen, the pRF Topographie und wie man die Topographie verwenden, um die pRF Zentrum zu modellieren.

Subscription Required. Please recommend JoVE to your librarian.

Protocol

1. Datenerfassung

  1. Bereiten Sie einen Stimulus-Protokoll, das wirksam bei der Auslösung einer zuverlässigen retinotope visuelle Antwort wie zuvor in Dumoulin und Wandell 1 und Lee et al ist. 2. Jedoch sind auch andere etablierte Paradigmen auch anwendbar in Abhängigkeit von der spezifischen experimentellen Fragen angesprochen werden.
  2. Präsentieren bar Reize driften über den Bildschirm nacheinander entlang 8 Richtungen des Raumes, in Schritten von 45 Grad. Stellen Sie sicher, dass die Bewegung synchron mit Scanner Rahmen Erwerb (TR ~ 2 sec), so dass der Balken bewegt sich einen Schritt einmal eine fMRI Rahmen beginnt und bleibt an der neuen Position, bis der Rahmen endet.
  3. Um eine korrekte Ausgangssignal zu messen, fügen Epochen ohne Stimulation 1 bar.
    1. Definiert ein Sichtfeld (10 bis 15 ° Radius) in Sichtwinkel, über den der dargeboten wird. Präsentieren bewegenden oder flackerndes Schachbrettmuster (checker size = 0,94 x 0,94 Grad 2, pattern Aktualisierungsrate = 250 ms / Muster) in der Bar zu starken visuellen Reaktionen hervorrufen.
    2. Geben Sie die folgenden spezifischen Parameter: 8 gleichmäßig verteilte Bewegungsrichtungen, Breite bar gleich 1,875 ° und Bars zu bewegen um die halbe Strichbreite pro Bild (2 sec). Weitere Details finden sich in Lee et al. 2.
    3. Generieren Sie einen Ort (~ 0,25 °) in die Mitte des Bildschirms, auf dem die Augen der Person während des Experiments zu fixieren. Farbe des Spots ändern zufällig in der Zeit.
  4. Scannen Sie das Gehirn eines Subjekts in einem MRI-Scanner mit einem typischen Echo-Planar-Bildgebung (EPI) Scan, 192 Frames Dauer (24 Bilder in jede Bewegungsrichtung) hat. Wiederholen Sie die Scans 4-8 mal, um das Signal-Rausch-Verhältnis zu erhöhen.
  5. Stellen Sie die Parameter für den EPI-Sequenz wie folgt: TR = 2 sec, TE = 40 ms, Matrixgröße = 64 x 64, 28 Scheiben, Voxelgröße = 3 x 3 x 3 mm 3, Flip-Winkel = 90 °, alternativ gelten Sequenzen mit einer feineren Auflösung (z.g., 2 x 2 x 2 mm 3) oder ein Kurz TR (beispielsweise 1-1,5 sec), die ausschließlich die Sehrinde 2.
  6. Spur Augenbewegungen mit einem Augenverfolgungssystem während funktionellen Scans Fixierung sicherzustellen innerhalb 1-1,5 ° vom Fixationspunkt gehalten.
    Anmerkung: Hier wird ein Kopf-Koordinate basierend Augenverfolger in einer Schutzbrille System verwendet wird, aber andere geeignete Augenverfolgungssysteme können stattdessen verwendet werden.
  7. Weisen Sie die Themen, um die Stelle auf die Mitte des Bildschirms in Schritt 1.3.2 erzeugt fixieren. Um sicherzustellen, dass die Probanden fixiert, sie anweisen, die Farbveränderungen der Befestigungspunkt zu melden.
  8. Erhalten anatomischen Scans auf 1 x 1 x 1 mm 3 Auflösung (zB T1-MPRAGE; TR = 1900 ms, TE = 2,26 ms, TI = 900 ms, Flipwinkel = 9 °, 176 Partitionen).
    ANMERKUNG: Diese anatomischen Scans zur Segmentierung sowie zum Ausrichten der funktionellen Bilder der Anatomie sowohl innerhalb der Abtastungen verwendet werden. Für eine bessere Ausrichtung zwischen Funktional (EPI) von Bildern und die Anatomie, erhalten auch eine inplane Anatomie Scans mit einer Auflösung identisch mit der EPI, mit T1-gewichteten schnell verdorben Gradientenecho (SPGR) Folge 1.

2. Datenvorverarbeitung

HINWEIS: Vor der Schätzung pRF Eigenschaften, einige typische fMRT-Daten vor, eine Bearbeitungsschritte benötigt werden, wie zum Beispiel Kopfbewegungskorrektur und die Ausrichtung der Funktionsvolumen der anatomischen Scan. In diesem Artikel werden alle Vorverarbeitung, Schätzung, Analyse und Präsentation der Ergebnisse erhalten werden unter Verwendung des Open-Source-MATLAB-basierte Software-Toolbox VISTA LAB auf der VISTA-Software zur Verfügung. http://white.stanford.edu/newlm/index.php/Main_Page.

  1. Laden Sie die anatomischen Scan in MATLAB und bereiten ein Volumen Anatomie mit einer Funktion namens createVolAnat.
  2. Segment graue Substanz, weiße Substanz und CSF mit der Funktion "ItkGray".
  3. Bereiten Funktionsdaten durch Umwandlung DICOM (dh (dh Standard-funktionellen MRT-Format) -Dateien, und Laden von Daten in VISTA mit einer Funktion namens mrInit.
  4. Correct Kopf-Bewegung, und richten Funktionsaufnahmen der Anatomie in Schritt 2.1 mit rxAlign basierend auf einer affinen Matrix-Transformation geladen.
  5. Durchschnittliche funktionelle Bewegungskorrigierten Scans für die Verbesserung der Signal-zu-Rausch-Verhältnis, indem Sie auf mrVISTA Analyse Timeseries Durchschnittliche tseries. Ausschließen von durchschnittlich Scans, während der Augenbewegungen weicht von Fixierung mehr als 1 bis 1,5 °. Wenn Signale von verschiedenen Läufen unterschiedliche DC-Drifts, durchschnittliche Funktions Scans nach dem Entfernen der DC-Drifts.
  6. Berechnen Sie die Mapping-Koordinaten zwischen funktionalen Scans und graue Substanz und Bestimmung der entsprechenden Grau Sache Voxel in den funktionalen Scans, indem Sie die folgenden Menüs: mrVISTA Fenster öffnen Grau Fenster 3-View. Vergeben BOLD-Signale in der grauen Substanz Voxel durch Interpolation, die Wahl einer der Optionen einvailable in mrVISTA.

3. Die Einschätzung der pRF Topographie und Parametric Modeling

  1. Laden Sie sich die Code-Dateien über den folgenden Link: https://sites.google.com/site/leesangkyun/prf/codes.zip, extrahieren Sie die komprimierte Datei und legen Sie sie in einem bevorzugten Ort für den lokalen Computer. Fügen Sie den Pfad des Ordners, in MATLAB.
  2. Stellen Sie die Stimulationsparameter in das Experiment, indem Sie die folgenden Menüs: mrVISTA Analyse retinotope Modell Parameter einstellen. Geben Sie die folgenden Parameter wie Stimulus Bilder, Stimulus Größe, die kanonische hämodynamischen Funktion, die Bildrate der fMRI-Scanner.
  3. Vor der pRF Schätzung, bereiten die Anfangsparametersätze (Abbildung 1B).
    1. Stellen Sie die Kreuzvalidierung Sets in "tprf_set_params.m" aus den Codedateien. Teilen Sie Zeitreihen in mindestens zwei Teilmengen (ein Satz für die Prüfung und die verbleibenden Sätze für die Ausbildung), die lang genug für die Bar zu Swee sindp die ganze Reiz Raum. Alternativ ohne Mittelung Scans in Schritt 2.4, Validierung Scans durch das Weglassen von einem Scan zum Testen und mit den restlichen Scans für die Ausbildung.
    2. Stellen Sie eine grobe Parametersatz (λ in 1; λ = [10 -2 10 -1 1 10 1 10 2]) in "tprf_set_params.m". Stellen Sie dann eine feine Skalenbereich ([0,1 0,3 0,5 0,7 0,9 1 3 5 7 9]) in "tprf_set_params.m".
      HINWEIS: Das Programm verwendet die Grob gesetzt, um den λ was zu den höchsten erklärte Varianz zu wählen. Dann sucht das Programm den Raum um den ausgewählten λ mit dem Feinmessbereich, weitere Verfeinerung der Auswahl von λ, die den höchsten erklärte Varianz ergibt.
    3. Einen Schwellenwert (0,2) der erklärten Varianz für visuell ansprechende Voxel in "tprf_set_params.m".
      HINWEIS: Diese Schwelle wird als Referenz für die Auswahl von visuell ansprechenden Voxel verwendet. Alternativ können Sie einROI für eine nicht-optisch ansprechenden Bereich (beispielsweise durch Ziehen einer Kugel mit einem Radius von 1 cm in einer nicht optisch ansprechenden Bereich des Gehirns), wobei der Schwellenwert automatisch berechnet werden.
    4. Festlegen eines Satzes von Schwellenwerten ([0,3, 0,5, 0,7]) zum Definieren der pRF Mittelbereich in dem normalisierten Topographie "tprf_set_params.m" (dh, [0-1] oder [-1 bis 1] mit Epochen ohne bar Stimulations in Schritt 1.3.1).
      HINWEIS: Aus dem Satz von Schwellenwerten das Programm vorgesehen wählt den "besten" Schwelle, dh die Schwelle, die eine pRF zentralen Region, für die das Zentrum pRF Modell erklärt die größte Signalvarianz definiert. Alternativ wählen Sie eine andere Menge von Schwellenwerten in Abhängigkeit von den Eigenschaften der Topographie.
  4. Ausführen "tprf_runpRFest.m" berechnen die pRF Topographie (1) entfernen und 2D anisotrope Gaußsche. Nach Angabe aller in diesem Protokoll beschriebenen Parameter, undAusführen des Codes, erhalten die endgültigen Schätzergebnisse.

Figur 1
Abbildung 1: PRF Schätzverfahren. (A) Schematische Darstellung des Prozesses folgte für pRF Topographie Schätzung h (t):. Hämodynamischen Antwortfunktion A (t): Reiz, m: pRF, Reg. L2-Norm Regularisierung (B) Spezifische Schritte für pRF Topographie Schätzung und pRF Zentrum Modellierung. Die Menge der für die Bestimmung erforderlichen Parameter in jedem Schritt aufgelistet. Eine eindimensionale Abschnitt Topographie und sein Modell veranschaulicht. Unter "Modell Fitting", schwarze und rote Kurven stellen die Topographie und seiner pRF Zentrum-Modell mit einem Zentrum Schwelle von 0,5. Die blau gestrichelte Linie gibt einen Schwellenwert für die pRFZentralregion.

Subscription Required. Please recommend JoVE to your librarian.

Representative Results

Genaue pRF Modellierung erfordert korrekt erfassen pRF Formen. Ohne Kenntnis der pRF Topographie, die Auswahl der kreissymmetrische Modelle in früheren Studien verwendet 1,9-11 ist eine vernünftige Wahl. Dies ist, weil, wenn die lokale retinotopic Organisation ist in allen Richtungen des Gesichtsfeldes homogen, könnte eine lokale Bevölkerung Antwort als eine kreissymmetrische kumulativen Summe von neuronalen Antworten dargestellt werden. Jedoch haben unsere Beobachtungen zeigen, dass dies nicht unbedingt der Fall ist (Abbildung 2). Daher kann die Beobachtung des pRF Topographie entscheidend für die Auswahl eines geeigneten parametrische Funktion eine pRF Modell sein. Dies ist ein Vorteil der pRF Topographie, so dass die Topographie basierte Modelle treffen die Direkt fit isotropen Gaußschen Modellen in pRF Zentrum Modellierung, was typischerweise in höheren erklärte Varianz (Fig. 2; siehe Lee et al 2 für weitere Vergleiche mit anderen Modelle). Diese Beispiele zeigen den Vorteil der Schätzung der pRF Topographie vor der Modellanpassung.

Abbildung 2
Abbildung 2: Beispiele für pRF Topographie Schätzung und Passung pRF Zentrum Modelle (A) Ein typisches pRF Topographie.. In der Topographie, rote Farbe zeigt die meisten reagieren Gebiet, das zeigt die pRF Zentrum liegt auf der Mitte der rechten horizontalen Meridian. Im pRF Topographie, sind Balkenmuster auf der pRF Center-Struktur mit geringen Gewichten auch manchmal beobachtet. Dies bezieht sich auf die Tatsache, dass der Bereich entlang der Balkenöffnung durch den pRF Zentrum wird auch gleichzeitig mit dem pRF Zentrum stimuliert. . Sie sind leicht in der Schwellenwertschritt eliminiert (B) Vergleich zwischen einem früheren Verfahren (DIG, Direkt fit isotropen Gaußschen) 1 und Topographie-basierte pRF Zentrum model (T-Modell). Die entsprechende Prozent erklärten Varianz ist über jedes Modell gezeigt. T-Modelle zeigen höhere Varianzaufklärung in allen Beispielen mit genaueren pRF Formerfassung. Siehe Lee et al. 2 für weitere Informationen und weitere Beispiele.

Eine wichtige Anforderung ist, dass das verwendete fMRI Paradigma bietet gute Retinotopie Daten. Dann wird die pRF Topographie Methode kann retinotope Exzentrizität und Azimut Karten (Abbildung 3) zu schätzen. Diese Karten zeigen ähnliche Grund retinotope Architektur als bisherige Methoden 1,4-7, aber sie sind genauer, da die Beobachtung der pRF Topographie ermöglicht es uns, besser zu trennen die pRF Zentrum von den hinteren Surround und von Potentialrauschen oder Artefakte entfernt, um die PRF-Center. Diese, unter anderem zu einer besseren Abschätzung der retinotope Karten bei hohen Exzentrizitäten (eine detaillierte Darstellung der beobachteten Unterschiede in Lee et al. 2).

Figur 3
Abbildung 3: retinotope Karten und pRF Größe (A) Exzentrizität und Polarwinkel Karten in der linken Hemisphäre eines Subjekts.. CS gibt die Sulkus. Im rechten Teil der Abbildung A, der schwarze Kreis zeigt eine Region-of-Interest (ROI), von dem die Voxel, deren pRF ist in Abbildung 4 dargestellt wird genommen. (B) Zusammenhang zwischen pRF Größe und Exzentrizität. Die pRF Größe mit Exzentrizität erhöht in visuellen Arealen V1-3. Dieses Grundstück befindet sich von (A) gezogen.

Die Topographie-basiertes Modell (T-Modell) Methode kann verwendet werden, um verschiedene pRF Eigenschaften wie pRF Größe, Dehnung, Ausrichtung und Surround-Unterdrückung effizient zu schätzen, ohne dass viele verschiedene parametrische Modelle zu testen. Zur Visualisierung von solchen Eigenschaften zu helfen, ein MATLAB Funktion (tprf_plotpRF.m) vorgesehen, die den pRF Topographie, die entsprechende pRF Zentrum Modell und ihre Anpassung an die rohen BOLD-Signal (Abbildung 4) aufgetragen ist. Man beachte, dass in einigen Fällen pRF Eigenschaften können auch direkt von der Topographie geschätzt werden, wodurch die Notwendigkeit für pRF Modellierung.

4
Abbildung 4: Demonstration der MATLAB-Toolbox von den Autoren entwickelt. Dieser Plot zeigt die pRF Topographie und entsprechende pRF Modellanpassung eines Voxels von einem Benutzer ausgewählt. Mit dem pRF Zentrum parametrische Modellvorhersage: Die dargestellte Voxel aus dem in 3A gezeigten ROI ausgewählten Rohstoffen. Tatsächlichen BOLD Antwort pred t: Vorhersage mit dem pRF Topographie, pred m.1fig4highres.jpg "target =" _ blank "> Bitte klicken Sie hier, um eine größere Version dieses Bild anzuzeigen.

Subscription Required. Please recommend JoVE to your librarian.

Discussion

Dieser Artikel beschreibt, wie die Topographie der visuellen Bevölkerung rezeptiven Felder im menschlichen visuellen Kortex und wie man es benutzt, um eine geeignete parametrisches Modell für die rezeptiven Feld auswählen zu schätzen. Für eine erfolgreiche Retinotopie sollte eine geeignete Stimulationsprotokoll und eine effiziente Analyseverfahren ausgewählt werden, und Versuchsparameter des Subjekts (Bewegung und Fixierung) optimiert werden. Bar Reize Ablauf eines ganzen Gesichtsfeld sind eine effiziente Anregung Paradigma für pRF Schätzung, da es deutliche BOLD Antworten von verschiedenen Stimuluspositionen erzeugt. Die vorgesehenen Verfahren erstellt die pRF Topographie. Da das Problem der pRF Schätzung ist in der Regel unterbestimmt wird ein mathematisches Werkzeug namens Firstregressions 3 verwendet, um die angemessene Einschränkung der Kargheit auf dem pRF igen Lösung durchzusetzen. Diese Bereinigung Technik ist sehr effektiv bei der Schätzung der pRF Modell, wenn die Anzahl der Beobachtungen (Zeitpunkts der BOLD-Signal) ist wesentlich kleiner als die Anzahl der Pixel für den räumlichen Dimension des Stimulus.

Diese Methode bietet robustere Schätzung des pRF Zentrum als bisherige Verfahren. Es gibt mehrere Gründe dafür: 1) es ersten Segmente der pRF zentralen Region von der pRF Topographie und dann ein geeignetes Modell passt, Vermeidung potenzieller Verzerrungen, pRF Modell beeinflussen können, passt in direktem Modelle (zB Surround-Unterdrückung oder Rauschartefakte vom pRF Mitte). 2) Mit der Fähigkeit, zu inspizieren die Topographie bietet eine visuell die Möglichkeit, die Leistung des fertigen Modells fit Aufdeckung systematischer Fehler sowie 3) die Möglichkeit, Merkmale des pRF Struktur, die ansonsten unentdeckt bleiben, kann erkennen, zu bestätigen. 4) Durch Einschränken der Passbereich, dieses Modell ist weniger wahrscheinlich, dass der PRF-Karte innerhalb der Grenze von Stimuluspräsentation falsch gegenüber der direkten für Modelle (siehe Abbildung 2B). Keinerdings ein Benutzer braucht sich bewusst sein, dass das vorgeschlagene Verfahren auch Einschränkungen für die genaue Erfassung von pRF Form in der Nähe des Stimulus Grenze. Dies ist aufgrund der Tatsache, dass in der Nähe der Grenze der Bar Reize aktivieren Teil rezeptiven Felder, um Voxel, deren pRF Zentrum würde normalerweise außerhalb des Stimuluspräsentation Region gehören. Jede rezeptiven Feld Mapping-Verfahren unterliegen würden dieses Problem und zeigt eine relative Peak an der Grenze, wenn es perfekt aus dem Teil des rezeptiven Feldes Zentrum, das für die gesamte zugeordnet ist, zu extrapolieren. Abgesehen davon, unsere Methode ist genauer als direkte sitz Methoden 1,9, die deutlich überschätzen die Entfernung zum Zentrum von PRFs, die in der Nähe der Reizdarbietung Grenze liegen neigen (siehe Abbildungen 5 und 6 von Lee et al. 2 für weitere Details).

Wie bereits erwähnt, zu konstruieren eine robuste pRF Topographie hängt von der freien Regularisierungsparameter, λ (Abbildung 1), die Be optimiert separat einzelner Voxel oder als gemeinsame Parameter für alle Voxel. Der Parameter Regularisierung beeinflusst pRF Topographie durch Einstellen des Ausmaßes der Armatur (Überanpassung oder Unter Fitting) auf die Daten. Während ein kleiner λ führt zu laut pRF Topographien (dh Überanpassung) im Vergleich zu dem eigentlichen pRF, ein großes λ drückt visuellen Reaktionen und damit in mehr Verbreitung Topographien führen als durch die tatsächliche Größe pRF gerechtfertigt (dh unter sitz). Auswahl des optimalen Lambda ist entscheidend für eine erfolgreiche pRF Schätzung. Wir schätzten, λ die in verschiedenen Teilmengen von Daten und bewertet diese Schätzungen mit Hilfe eines Kreuzvalidierungsstrategie. Dies minimiert Verzerrungen in pRF Topographie Schätzung. Mögliche Rest Verzerrungen werden in der Mitte pRF Modellierungsschritt, in dem verschiedene Topographie Schwellen werden untersucht, um eine, die in der höchsten erklärte Varianz ergibt wählen reduziert (siehe Lee et al. 2).

Fischließlich die vorgeschlagene Topographie Ansatz ist rechnerisch effizient. Die Schätzung der pRF Topographien über alle Voxel, einschließlich der Suche nach der optimalen Regelungsparameter λ, dauert nur wenige Minuten in einer PC-Umgebung. Identifizieren visuell reagierenden Voxel in diesem Schritt schließt sie von der mehr rechenintensive Schritt der pRF Zentrum Modellierung, die Effizienz weiter zu verbessern. Vielleicht noch wichtiger ist, die Ermittler müssen nicht mehr mehrere verschiedene pRF Modelle zu testen, um eine, die gut passt, da sie bei der Auswahl des geeigneten Modells durch den pRF Topographie geführt zu finden.

Die in diesem Protokoll demonstriert Verfahren misst Bevölkerung rezeptiven Feld Topographie und verwendet sie, um Bevölkerung empfänglich Feldmodellierung zu führen. Dieser Ansatz reduziert die in direktem Bevölkerung rezeptiven Feld Zuordnungsmethoden vorhanden Bias, was zu mehr robust und genau pRF Schätzungen. Es minimiert auch systematische Fehler und ermöglicht es uns, die funktionelle Organ studierensierung des visuellen Kortex mit höherer Empfindlichkeit. Es ist im Falle von Patienten mit Läsionen der Sehbahn besonders anwendbar, in denen pRF Struktur kann nicht einfach, ein priori antizipieren.

Subscription Required. Please recommend JoVE to your librarian.

Materials

Name Company Catalog Number Comments
MRI scanner Siemens/Philips/GE
MATLAB The Mathworks, Inc. http://www.mathworks.com 
VISTA software VISTA software group http://white.stanford.edu/newlm/index.php/Software
PsychoToolbox PsychoToolbox  http://psychtoolbox.org 
Eye Tracker (VisuaStimDigital) Resonance Technology Inc http://mrivideo.com/

DOWNLOAD MATERIALS LIST

References

  1. Dumoulin, S. O., Wandell, B. A. Population receptive field estimates in human visual cortex. Neuroimage. 39, 647-660 (2008).
  2. Lee, S., Papanikolaou, A., Logothetis, N. K., Smirnakis, S. M., Keliris, G. A. A new method for estimating population receptive field topography in visual cortex. Neuroimage. 81, 144-157 (2013).
  3. Hastie, T., Tibshirani, R., Friedman, J. H. The elements of statistical learning : data mining, inference, and prediction. , 2nd edn, Springer. (2009).
  4. Sereno, M. I., et al. Borders of multiple visual areas in humans revealed by functional magnetic resonance imaging. Science. 268, 889-893 (1995).
  5. Engel, S. A., Glover, G. H., Wandell, B. A. Retinotopic organization in human visual cortex and the spatial precision of functional MRI. Cereb Cortex. 7, 181-192 (1997).
  6. Engel, S. A., et al. fMRI of human visual cortex. Nature. 369, 525 (1994).
  7. DeYoe, E. A., et al. Mapping striate and extrastriate visual areas in human cerebral cortex. Proc Natl Acad Sci U S A. 93, 2382-2386 (1996).
  8. Greene, C. A., Dumoulin, S. O., Harvey, B. M., Ress, D. Measurement of population receptive fields in human early visual cortex using back-projection tomography. J Vis. , (2014).
  9. Zuiderbaan, W., Harvey, B. M., Dumoulin, S. O. Modeling center-surround configurations in population receptive fields using fMRI. J Vis. , (2012).
  10. Harvey, B. M., Dumoulin, S. O. The relationship between cortical magnification factor and population receptive field size in human visual cortex: constancies in cortical architecture. J Neurosci. 31, 13604-13612 (2011).
  11. Haak, K. V., Cornelissen, F. W., Morland, A. B. Population receptive field dynamics in human visual cortex. PLoS One. 7, e37686 (2012).

Tags

Verhalten Bevölkerungs rezeptiven Feld Vision die funktionelle Magnetresonanztomographie Retinotopie
Topographische Schätzung von Visual Bevölkerung rezeptiven Felder von fMRI
Play Video
PDF DOI DOWNLOAD MATERIALS LIST

Cite this Article

Lee, S., Papanikolaou, A., Keliris,More

Lee, S., Papanikolaou, A., Keliris, G. A., Smirnakis, S. M. Topographical Estimation of Visual Population Receptive Fields by fMRI. J. Vis. Exp. (96), e51811, doi:10.3791/51811 (2015).

Less
Copy Citation Download Citation Reprints and Permissions
View Video

Get cutting-edge science videos from JoVE sent straight to your inbox every month.

Waiting X
Simple Hit Counter