(************** Content-type: application/mathematica ************** CreatedBy='Mathematica 5.2' Mathematica-Compatible Notebook This notebook can be used with any Mathematica-compatible application, such as Mathematica, MathReader or Publicon. The data for the notebook starts with the line containing stars above. To get the notebook into a Mathematica-compatible application, do one of the following: * Save the data starting with the line of stars above into a file with a name ending in .nb, then open the file inside the application; * Copy the data starting with the line of stars above to the clipboard, then use the Paste menu command inside the application. Data for notebooks contains only printable 7-bit ASCII and can be sent directly in email or through ftp in text mode. Newlines can be CR, LF or CRLF (Unix, Macintosh or MS-DOS style). NOTE: If you modify the data for this notebook not in a Mathematica- compatible application, you must delete the line below containing the word CacheID, otherwise Mathematica-compatible applications may try to use invalid cache data. For more information on notebooks and Mathematica-compatible applications, contact Wolfram Research: web: http://www.wolfram.com email: info@wolfram.com phone: +1-217-398-0700 (U.S.) Notebook reader applications are available free of charge from Wolfram Research. *******************************************************************) (*CacheID: 232*) (*NotebookFileLineBreakTest NotebookFileLineBreakTest*) (*NotebookOptionsPosition[ 10905128, 274291]*) (*NotebookOutlinePosition[ 10905771, 274313]*) (* CellTagsIndexPosition[ 10905727, 274309]*) (*WindowFrame->Normal*) Notebook[{ Cell[CellGroupData[{ Cell["Spatial localization of a 2D Gaussian signal", "Title"], Cell[CellGroupData[{ Cell["PACKAGES", "Section"], Cell["\<\ First, we call two packages that we will use for our analysis, namely \ Graphic3D and NonlinearFit for making movies and curve fitting, respectively.\ \ \>", "Text"], Cell[BoxData[{ \(<< Graphics`Graphics3D`\), "\[IndentingNewLine]", \(<< Statistics`NonlinearFit`\)}], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["INPUTS", "Section"], Cell["\<\ The routine requires several manual inputs from the user. First, the \ approximate center of a centroid needs to be declared. The easiest way to \ obtain this is using WASABI program from HAMAMATSU.\ \>", "Text"], Cell[BoxData[{ \(\(x\_o = 64;\)\), "\[IndentingNewLine]", \(\(y\_o = 67;\)\)}], "Input"], Cell["\<\ Next, the user specify the size of the size of a square region that will be \ used in the movies and curve fitting. \ \>", "Text"], Cell[BoxData[ \(\(\[CapitalDelta] = 6;\)\)], "Input"], Cell["and the first and last files to use", "Text"], Cell[BoxData[{ \(\(start = 10000;\)\), "\[IndentingNewLine]", \(\(finish = 10060;\)\)}], "Input"] }, Open ]], Cell[CellGroupData[{ Cell["MOVIE", "Section"], Cell[CellGroupData[{ Cell[BoxData[ \(For[\[IndentingNewLine]i = start, i < finish + 1, \[IndentingNewLine]\[IndentingNewLine]data = \ \[IndentingNewLine]BinaryReadList["\" <> ToString[ i] <> "\<.RBF\>", "\"]; \[IndentingNewLine]\ \[IndentingNewLine]data = Partition[Drop[data, 32], 144]; \[IndentingNewLine]\[IndentingNewLine]datafit = \ \[IndentingNewLine]Table[{i, j, Extract[data, {i, j}]}, {i, 1, 144}, {j, 1, 144}]; \[IndentingNewLine]\[IndentingNewLine]crop = Table[datafit[\([i, j]\)], {i, x\_o - \[CapitalDelta], x\_o + \[CapitalDelta]}, {j, y\_o - \[CapitalDelta], y\_o + \[CapitalDelta]}]; \ \[IndentingNewLine]\[IndentingNewLine]cropfit = Flatten[crop, 1]; \[IndentingNewLine]\[IndentingNewLine]p2 = ListContourPlot[\[IndentingNewLine]Partition[ Table[Extract[cropfit[\([i]\)], 3], {i, 1, \((2\ \[CapitalDelta] + 1)\)\^2}], \((2\ \[CapitalDelta] + 1)\)], \[IndentingNewLine]ImageSize \[Rule] 300, \[IndentingNewLine]ContourLines \[Rule] False, \[IndentingNewLine]Contours \[Rule] 30, \[IndentingNewLine]FrameLabel \[Rule] TraditionalForm /@ {"\", "\", "\