Interferometric ultra-high resolution 3D imaging through brain sections.
The 8 matches · 2 of them tie a paragraph to a whole file, not to given lines: weak matches, whose lines are not tinted
- [1] § Methods › Systematic configuration and 4Pi-BRAINSPOT algorithm › 4Pi in-situ PSF retrieval ↔ Support/PSF Toolbox_4pi/@PSF_4pi/PSF_4pi.m, lines 2–91 · score 0.85 · coherent pupil functions, Pi PSF model, objective misalignment, phase retrieval, cavity phase, depth
- [2] § Results › Design of 4Pi-BRAINSPOT ↔ Support/PSF Toolbox_4pi/@PSF_4pi/PSF_4pi.m, lines 2–91 · score 0.81 · coherent pupil functions, Pi PSF model, PSFs generated, objective misalignments, phase retrieval, cavity phase
- [3] § Methods › Systematic configuration and 4Pi-BRAINSPOT algorithm › 4Pi in-situ PSF retrieval ↔ 4Pi-BRAINSPOT toolbox/4Pi_insitu_PSF_retrieval/INSPR4Pi_model_generation.m, lines 12–104 · score 0.76 · interferometric patterns, Pi PSF, PSF retrieval, pupil functions, XYZ, refined
- [4] § Results › Design of 4Pi-BRAINSPOT ↔ 4Pi-BRAINSPOT toolbox/4Pi_insitu_PSF_retrieval/INSPR4Pi_model_generation.m, lines 12–104 · score 0.69 · Pi interference, Interferometric pattern, Pi PSF, pupil functions, refined, single molecule
- [5] § Results › Characterization of 4Pi-BRAINSPOT performance ↔ 4Pi-BRAINSPOT toolbox/brainspot_4pi_GUI.m, lines 140–224 · score 0.63 · vertical astigmatism, Zernike modes, axial directions, PSF model, axial positions, lateral
- [6] § Methods › Quantification and statistical analysis › Analysis pipeline of 3D structural morphology and circumference ↔ 4Pi-BRAINSPOT Supplementary analysis/Function/Analysis_spine_area_3D.m, lines 126–173 · score 0.58 · alphaShape algorithm, boundaries, sliced, fitting, positions
- [7] § Results › Nanoscale visualization and quantitative analysis of subcellular organisms in mammalian cells ↔ 4Pi-BRAINSPOT toolbox/4Pi_insitu_PSF_retrieval/gen_init_4PiPupil_v2.m, the whole file · a weak match · score 0.58 · immersion oil, cavity phase, refractive, depths, emission, objective
- [8] § Results › Characterization of 4Pi-BRAINSPOT performance ↔ 4Pi-BRAINSPOT toolbox/Dynamic_model_update/caliphi0_and_misObj_NCC_v3.m, the whole file · a weak match · score 0.50 · objective misalignment, model update, pupil functions, NCC, guess, median
Paper
Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
MATLAB · 359 lines · 16 KB · other · 2 matches
- classdef PSF_4pi < handle
- % (C) Copyright Huang Lab, Weldon School of Biomedical Engineering,
- % All rights reserved Purdue University, West Lafayette, IN, USA
- %
- %
- %
- % PSF_4pi class for generating 4PI PSF models
- % create object: obj = PSF_4pi(PRstruct)
- %
- % PSF_4pi Methods:
- % precomputeParam - generate images for k space operation
- % gen2Pupil - generate pupil functions for top and bottom emission paths from user defined or phase retrieved Zernike coefficients
- % genPupil - generate pupil functions from user defined or phase retrieved Zernike coefficients
- % genPupil_4pi - generate coherent pupil functions for the four detection channels
- % genPSF - generate PSF from a given pupil function
- % genPSF_4pi - generate interferometric PSFs for the four detection channels assuming complete interference where the modulation depth is at the maximum
- % genPSF_4pi_md - generate interferometric PSFs for the four detection channels assuming partial interference
- % scalePSF - generate OTF rescaled PSFs
- % addMisalign - add additional objective misalignment into pupil function
- properties
- % PRstruct - define necessary parameters for a PSF model
- % NA
- % Lambda
- % RefractiveIndex
- % Pupil: phase retrieved pupil function
- % phase: phase image
- % mag: magnitude image
- % Zernike_phase: coefficient of zernike polynomials representing the pupil phase
- % Zernike_mag: coefficient of zernike polynomials representing the pupil phase
- % SigmaX: sigmax of Gaussian filter for OTF rescale, unit is 1/micron in k space, the conversion to real space is 1/(2*pi*SigmaX), unit is micron
- % SigmaY: sigmay of Gaussian filter for OTF rescale, unit is 1/micron in k space, the conversion to real space is 1/(2*pi*SigmaY), unit is micron
- PRstruct;
- PRstruct1;
- PRstruct2;
- Xpos;% x positions of simulated emitters, a vector of N elements, unit is pixel
- Ypos;% y positions of simulated emitters, a vector of N elements, unit is pixel
- Zpos;% z positions of simulated emitters, a vector of N elements, unit is micron
- nMed;% refractive index of the sample medium
- PSFsize; % image size of the pupil function
- Boxsize; % image size of out put PSF
- Pixelsize;% pixel size at sample plane, unit is micron
- Z;% object from Zernike_Polynomials class
- StagePosUp;% used for PSF model with index mismatch aberration, the distance the top objective moved to focus on the bead away from the coverslip, unit: micron
- StagePosDown;% used for PSF model with index mismatch aberration, the distance the sample stage (relative to bottom objective) moved to focus on the bead away from the coverslip, unit: micron
- Phasediff;% phase difference between s- and p-polarizations
- Iratio = 1;% transmission ratio between top and bottom emission path, from 0 to 1
- Phi0;% cavity phase
- Zoffset = 0; % z offset caused by index mismatch aberration, set to zero when using 'mvbead'
- ChamberH; % the height of the sample chamber
- PlaneDis = 0; % parameter used for dual focal plane method in 4pi
- ModulationDepth = 1;% modulation strength of interferometric PSFs,from 0 to 1
- % PSFs - out put PSFs from Fourier transform of the pupil function,
- % it's a 3D matrix of Boxsize x Boxsize x N, N is the number of
- % elements in Xpos.
- PSFs;
- end
- properties (SetAccess = private, GetAccess = private)
- % precompute images for k space operation
- Zo;% r coordinates of out put PSF, it's a image of PSFsize x PSFsize
- k_r;% k_r coordinates of out put OTF, it's a image of PSFsize x PSFsize
- k_z;% k_z coordinates of out put OTF, it's a image of PSFsize x PSFsize
- Phi;% phi coordinates out put PSF, it's a image of PSFsize x PSFsize
- NA_constrain;% a circle function defining the limit of k_r, it's a image of PSFsize x PSFsize
- Cos1;% cos(theta1), theta1 is the angle of between the k vector and the optical axis in the sample medium
- Cos3;% cos(theta3), theta3 is the angle of between the k vector and the optical axis in the immersion medium
- end
- properties (SetAccess = private, GetAccess = public)
- % IMMPSFs - out put PSFs from Fourier transform of the pupil
- % function modified with the index mismatch aberration, it's a 3D
- % matrix of Boxsize x Boxsize x N, N is the number of elements in
- % Xpos.
- IMMPSFs;
- % ScaledPSFs - out put PSFs after OTF rescale, it's a 3D matrix of
- % Boxsize x Boxsize x N, N is the number of elements in Xpos.
- ScaledPSFs;
- % Pupil - pupil function generated from a set of zernike polynomials
- % phase: phase image of PSFsize x PSFsize
- % mag: magnitude image of PSFsize x PSFsize
- Pupila; % Pupil from top objective lens
- Pupilb; % Pupil from bottom objective lens
- Pupil;
- Pupil4pi;% coherent pupil function for each quadrant
- PSF4pi;% 4PiPSF model, containing PSFs from the four quadrants
- end
- methods
- function obj=PSF_4pi(PRstruct)
- obj.PRstruct=PRstruct;
- end
- function precomputeParam(obj)
- % precomputeParam - generate images for k space operation, and saved in
- % precomputed parameters.
- [X,Y]=meshgrid(-obj.PSFsize/2:obj.PSFsize/2-1,-obj.PSFsize/2:obj.PSFsize/2-1);
- obj.Zo=sqrt(X.^2+Y.^2);
- scale=obj.PSFsize*obj.Pixelsize;
- obj.k_r=obj.Zo./scale;
- obj.Phi=atan2(Y,X);
- n=obj.PRstruct.RefractiveIndex;
- Freq_max=obj.PRstruct.NA/obj.PRstruct.Lambda;
- obj.NA_constrain=obj.k_r<Freq_max;
- obj.k_z=sqrt((n/obj.PRstruct.Lambda)^2-obj.k_r.^2).*obj.NA_constrain;
- sin_theta3=obj.k_r.*obj.PRstruct.Lambda./n;
- sin_theta1=n./obj.nMed.*sin_theta3;
- obj.Cos1=sqrt(1-sin_theta1.^2);
- obj.Cos3=sqrt(1-sin_theta3.^2);
- end
- function genZernike(obj)
- % create Zernike_Polynomials object
- zk = Zernike_Polynomials();
- zk.Ordering = 'Wyant';
- %zk.Ordering = 'Noll';
- ZN=sqrt(numel(obj.PRstruct.Zernike_phase))-1;
- zk.setN(ZN);
- zk.initialize();
- [Zrho, Ztheta, Zinit] = ...
- zk.params3_Zernike(obj.Phi, obj.k_r, obj.PRstruct.NA, obj.PRstruct.Lambda);
- zk.matrix_Z(Zrho, Ztheta, Zinit);
- obj.Z = zk;
- end
- function gen2Pupil(obj,PRstruct1,PRstruct2)
- obj.PRstruct = PRstruct1;
- obj.precomputeParam();
- obj.genZernike();
- obj.genPupil();
- obj.Pupila.phase = obj.Pupil.phase;
- obj.Pupila.mag = obj.Pupil.mag;
- obj.PRstruct1.Pupil.phase = exp(1i.*obj.Pupila.phase);
- obj.PRstruct1.Pupil.mag = obj.Pupila.mag;
- obj.PRstruct = PRstruct2;
- obj.precomputeParam();
- obj.genZernike();
- obj.genPupil();
- obj.Pupilb.phase = obj.Pupil.phase;
- obj.Pupilb.mag = obj.Pupil.mag;
- obj.PRstruct2.Pupil.phase = exp(1i.*obj.Pupilb.phase);
- obj.PRstruct2.Pupil.mag = obj.Pupilb.mag;
- end
- function genPupil(obj)
- % genPupil - generate pupil function from Zernike polynomials
- % Zernike polynomials are a set of images generated by using
- % Zernike_Polynomials class. The coefficients of the Zernike
- % polynomials are given from the 'PRstruct'. The resulting
- % pupil function includes a phase image and a magnitude image
- %
- % see also Zernike_Polynomials
- R=obj.PSFsize;
- ceffp=obj.PRstruct.Zernike_phase;
- ceffm=obj.PRstruct.Zernike_mag;
- pupil_phase=zeros(R,R);
- pupil_mag=zeros(R,R);
- N=numel(ceffp);
- for k = 1 : N
- pupil_phase = pupil_phase + obj.Z.ZM(:, :, k) .* ceffp(k);
- end
- for k = 1 : N
- pupil_mag = pupil_mag + obj.Z.ZM(:, :, k) .* ceffm(k);
- end
- % Normalize Zernike coefficients.
- % tmp = pupil_mag .* (1/R); changed by FX
- tmp = pupil_mag;
- normF = sqrt(sum(sum(tmp .* conj(tmp))));
- Ceffnorm = ceffm ./ normF;
- pupil_magnorm = zeros(R, R);
- for k = 1 : N
- pupil_magnorm = pupil_magnorm + obj.Z.ZM(:, :, k) .* Ceffnorm(k);
- end
- obj.Pupil.phase=pupil_phase;
- obj.Pupil.mag=pupil_magnorm;
- end
- function set2Pupil(obj,PRstruct1,PRstruct2) % set pupil function from reterived model
- obj.PRstruct1.Pupil.phase = PRstruct1.Pupil.phase;
- obj.PRstruct1.Pupil.mag = PRstruct1.Pupil.mag;
- obj.PRstruct2.Pupil.phase = PRstruct2.Pupil.phase;
- obj.PRstruct2.Pupil.mag = PRstruct2.Pupil.mag;
- obj.Pupila.phase = PRstruct1.Pupil.phase;
- obj.Pupila.mag = PRstruct1.Pupil.mag;
- obj.Pupilb.phase = PRstruct2.Pupil.phase;
- obj.Pupilb.mag = PRstruct2.Pupil.mag;
- end
- function pad2Pupil(obj)
- pad_size = (obj.PSFsize - size(obj.Pupila.phase,1)) / 2;
- obj.Pupila.phase = padarray(obj.Pupila.phase, [pad_size pad_size]);
- obj.Pupila.mag = padarray(obj.Pupila.mag, [pad_size pad_size]);
- obj.Pupilb.phase = padarray(obj.Pupilb.phase, [pad_size pad_size]);
- obj.Pupilb.mag = padarray(obj.Pupilb.mag, [pad_size pad_size]);
- end
- function addMisalign(obj,aber_misalign)
- R=obj.PSFsize;
- pupila_phase=zeros(R,R);
- pupilb_phase=zeros(R,R);
- N=numel(aber_misalign);
- for k = 1 : N
- if k == 3
- sign = 1;
- else
- sign = -1;
- end
- pupila_phase = pupila_phase + sign*obj.Z.ZM(:, :, k+1) .* aber_misalign(k);
- pupilb_phase = pupilb_phase + obj.Z.ZM(:, :, k+1) .* aber_misalign(k);
- end
- obj.Pupila.phase = obj.PRstruct1.Pupil.phase .* exp(1i*pupila_phase);
- obj.Pupilb.phase = obj.PRstruct2.Pupil.phase .* exp(1i*pupilb_phase);
- end
- function setUnifMag(obj) % set uniform magnitude for pupil function
- mag=obj.Pupila.mag;
- mag(mag>0)=1;
- mag = mag.^2;
- mag = mag./sum(sum(mag));
- mag = sqrt(mag);
- obj.Pupila.mag = mag;
- mag=obj.Pupilb.mag;
- mag(mag>0)=1;
- mag = mag.^2;
- mag = mag./sum(sum(mag));
- mag = sqrt(mag);
- obj.Pupilb.mag = mag;
- end
- function genPupil_4pi_2(obj)
- pupilA = obj.Pupila.mag.*obj.Pupila.phase;%top
- pupilB = obj.Pupilb.mag.*obj.Pupilb.phase.*exp(1i.*obj.Phi0).*obj.Iratio;%bottom
- kz = obj.k_z;
- zpos = obj.Zpos;
- phia = 0;% no effect on 4pi PSF
- phib = obj.Phasediff; % changed by FX
- N = numel(obj.Xpos);
- R = obj.PSFsize;
- obj.Pupil4pi.s1 = zeros(R,R,N);
- obj.Pupil4pi.s2 = zeros(R,R,N);
- obj.Pupil4pi.p1 = zeros(R,R,N);
- obj.Pupil4pi.p2 = zeros(R,R,N);
- obj.Pupil.top = zeros(R,R,N);
- obj.Pupil.bot = zeros(R,R,N);
- for ii=1:N
- defocusphaseA = exp(-2.*pi.*1i.*(zpos(ii)+obj.Zoffset).*kz);%top
- defocusphaseB = exp(2.*pi.*1i.*zpos(ii).*kz); %bottom
- obj.Pupil4pi.s1(:,:,ii) = pupilA.*exp(1i*pi).*exp(1i*phia).*defocusphaseA+pupilB.*exp(1i*phib).*defocusphaseB;
- obj.Pupil4pi.s2(:,:,ii) = pupilA.*exp(1i*phia).*defocusphaseA+pupilB.*exp(1i*phib).*defocusphaseB;
- obj.Pupil4pi.p1(:,:,ii) = pupilA.*exp(1i*pi).*defocusphaseA+pupilB.*defocusphaseB; %pi delay in ha
- obj.Pupil4pi.p2(:,:,ii) = pupilA.*defocusphaseA+pupilB.*defocusphaseB;
- obj.Pupil.top(:,:,ii) = pupilA.*defocusphaseA;
- obj.Pupil.bot(:,:,ii) = pupilB.*defocusphaseB;
- end
- end
- function psfs=genPSF(obj,pupil)
- % genPSF - generate PSFs from the given pupil function.
- % The PSFs are directly calculated from the Fourier transform
- % of pupil functions modified by shift phase in x, y and
- % defocus phase in z. The out put is 'PSFs'
- N=numel(obj.Xpos);
- R=obj.PSFsize;
- Ri=obj.Boxsize;
- psfs=zeros(Ri,Ri,N);
- for ii=1:N
- shiftphase=-obj.k_r.*cos(obj.Phi).*obj.Xpos(ii).*obj.Pixelsize-obj.k_r.*sin(obj.Phi).*obj.Ypos(ii).*obj.Pixelsize;
- shiftphaseE=exp(-1i.*2.*pi.*shiftphase);
- defocusphaseDual = exp(2.*pi.*1i.*obj.PlaneDis.*obj.k_z);
- if nargin>1
- pupil_complex=pupil(:,:,ii).*shiftphaseE.*defocusphaseDual;
- else
- defocusphaseE=exp(2.*pi.*1i.*obj.Zpos(ii).*obj.k_z);
- pupil_complex=obj.Pupil.mag.*exp(obj.Pupil.phase.*1i).*shiftphaseE.*defocusphaseE;
- end
- psfA=abs(fftshift(fft2(pupil_complex)));
- Fig2=psfA.^2;
- realsize0=floor(Ri/2);
- realsize1=ceil(Ri/2);
- startx=-realsize0+R/2+1;endx=realsize1+R/2;
- starty=-realsize0+R/2+1;endy=realsize1+R/2;
- psfs(:,:,ii)=Fig2(startx:endx,starty:endy)./R^2;
- end
- obj.PSFs=psfs;
- end
- function genPSF_4pi(obj)
- obj.PSF4pi.s1=obj.genPSF(obj.Pupil4pi.s1)./4;
- obj.PSF4pi.s2=obj.genPSF(obj.Pupil4pi.s2)./4;
- obj.PSF4pi.p1=obj.genPSF(obj.Pupil4pi.p1)./4;
- obj.PSF4pi.p2=obj.genPSF(obj.Pupil4pi.p2)./4;
- end
- function genPSF_4pi_md(obj)
- obj.PSF4pi.s1 = obj.ModulationDepth.*obj.genPSF(obj.Pupil4pi.s1)./4 + (1-obj.ModulationDepth).*(obj.genPSF(obj.Pupil.top)+obj.genPSF(obj.Pupil.bot))./4;
- obj.PSF4pi.s2 = obj.ModulationDepth.*obj.genPSF(obj.Pupil4pi.s2)./4 + (1-obj.ModulationDepth).*(obj.genPSF(obj.Pupil.top)+obj.genPSF(obj.Pupil.bot))./4;
- obj.PSF4pi.p1 = obj.ModulationDepth.*obj.genPSF(obj.Pupil4pi.p1)./4 + (1-obj.ModulationDepth).*(obj.genPSF(obj.Pupil.top)+obj.genPSF(obj.Pupil.bot))./4;
- obj.PSF4pi.p2 = obj.ModulationDepth.*obj.genPSF(obj.Pupil4pi.p2)./4 + (1-obj.ModulationDepth).*(obj.genPSF(obj.Pupil.top)+obj.genPSF(obj.Pupil.bot))./4;
- % FX
- obj.PSF4pi.top = obj.genPSF(obj.Pupil.top);
- obj.PSF4pi.bot = obj.genPSF(obj.Pupil.bot);
- end
- function scalePSF(obj)
- % scalePSF - generate OTF rescaled PSFs
- % It operates 'PSFs' using the OTFrescale class. The OTF
- % rescale acts as a 2D Gaussian filter, the resulting PSFs
- % are smoother than the orignal PSFs.
- %
- % see also OTFrescale
- otfobj=OTFrescale;
- otfobj.SigmaX=obj.PRstruct.SigmaX;
- otfobj.SigmaY=obj.PRstruct.SigmaY;
- otfobj.Pixelsize=obj.Pixelsize;
- otfobj.PSFs=obj.PSFs;
- otfobj.scaleRspace();
- obj.ScaledPSFs=otfobj.Modpsfs;
- end
- end
- end
PSF_4pi.m at commit 0d9e7e7, under other · at the source
Overview
- Weldon School of Biomedical Engineering, Purdue University, West Lafayette, IN USA
- School of Optics and Photonics, Beijing Institute of Technology, Beijing, China
- School of Medical Engineering, Beijing Institute of Technology, Zhuhai, China
- Department of Biological Science, Purdue University, West Lafayette, IN USA
- Purdue Institute for Integrative Neuroscience, Purdue University, West Lafayette, IN USA
- Institute for Cancer Research, Purdue University, West Lafayette, IN USA
Abstract
The abstract is not reproduced here: the paper's license (CC BY-NC-ND) does not allow it. Read it in the paper, at the publisher or on Europe PMC.
Repository
Its files are read in the Code ↔ Paper reader above, with 8 matches between paragraphs and lines of code.
HuanglabPurdue/4Pi-BRAINSPOT
0d9e7e78021da4b143614701ec8a9430554c8dbf, 18 December 2025Availability: 1 check, the latest on 29 September 2026: the link answers
- 29 September 2026: the link answers
79 files
- 4Pi-BRAINSPOT Supplementary analysis/
Analysis3D_main.m , MATLAB, 44 lines - 4Pi-BRAINSPOT Supplementary analysis/
Function/ , MATLAB, 322 lines, 1 matchAnalysis_spine_area_3D.m - 4Pi-BRAINSPOT Supplementary analysis/
Function/ , MATLAB, 319 linesAnalysis_spine_section_3 D.m - 4Pi-BRAINSPOT Supplementary analysis/
Function/ , MATLAB, 254 linesfit_ellipse.m - 4Pi-BRAINSPOT Supplementary analysis/
Function/ , MATLAB, 532 linesinterparc.m - 4Pi-BRAINSPOT Supplementary analysis/
Function/ , MATLAB, 293 linesmyginput.m - 4Pi-BRAINSPOT Supplementary analysis/
Function/ , MATLAB, 22 linespickspot.m - 4Pi-BRAINSPOT Supplementary analysis/
Function/ , MATLAB, 42 linessrhist_color_ca.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 44 lines/ FourierShift2D.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 109 lines, 2 matches/ INSPR4Pi_model_generatio n.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 1,190 lines/ PR4Pi_4channels/ Zernike_Polynomials.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 22 lines/ PR4Pi_4channels/ datapreprocess.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 132 lines/ PR4Pi_4channels/ findOTFparam_fixedsigma. m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 15 lines/ PR4Pi_4channels/ fitdefocus.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 66 lines/ PR4Pi_4channels/ genMpsf.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 143 lines/ PR4Pi_4channels/ genPRfigs_4ch.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 59 lines/ PR4Pi_4channels/ genZKresult.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 215 lines/ PR4Pi_4channels/ phaseretrieve4pi_4channe ls_extraDefocus_v3.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 30 lines/ PR4Pi_4channels/ precomputeParam.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 122 lines/ PR4Pi_fromAveZ_4channels _Zshift.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 9 lines/ cc2.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 190 lines/ classify_4PiPSF_4channel s_seperately.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 226 lines/ dftregistration.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 52 lines/ genIniguess.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 94 lines/ gen_aber4PiPSF_directlyf romPR.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 54 lines, 1 match/ gen_init_4PiPupil_v2.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 21 lines/ registration_in_each_cha nnel.m - 4Pi-BRAINSPOT toolbox/
4Pi_insitu_PSF_retrieval , MATLAB, 25 lines/ subregion_normalization_ 4pi.m - 4Pi-BRAINSPOT toolbox/
4Pi_localization/ , MATLAB, 207 linesanalysis3D_from4pi_pupil _v2.m - 4Pi-BRAINSPOT toolbox/
4Pi_localization/ , MATLAB, 31 linescal_model_affine_4pi.m - 4Pi-BRAINSPOT toolbox/
4Pi_localization/ , MATLAB, 8 linescatstruct.m - 4Pi-BRAINSPOT toolbox/
4Pi_localization/ , MATLAB, 217 linescrop_4pi_subregion_trans Model_sCMOS.m - 4Pi-BRAINSPOT toolbox/
4Pi_localization/ , MATLAB, 52 linesgenIniguess.m - 4Pi-BRAINSPOT toolbox/
4Pi_localization/ , MATLAB, 98 linesgenini4pi_z_parfor_affin e.m - 4Pi-BRAINSPOT toolbox/
4Pi_localization/ , MATLAB, 38 linesgenpsf_4pi_real.m - 4Pi-BRAINSPOT toolbox/
4Pi_localization/ , MATLAB, 75 linesgenpsfstruct.m - 4Pi-BRAINSPOT toolbox/
4Pi_localization/ , MATLAB, 52 linesgensamplepsf_4pi_affine. m - 4Pi-BRAINSPOT toolbox/
4Pi_localization/ , MATLAB, 165 linesloc_channel_specific_4Pi model.m - 4Pi-BRAINSPOT toolbox/
Channel_alignment/ , MATLAB, 43 linesbiplane_registration.m - 4Pi-BRAINSPOT toolbox/
Channel_alignment/ , MATLAB, 35 linesregistration_4pi.m - 4Pi-BRAINSPOT toolbox/
Dynamic_model_update/ , MATLAB, 82 linesanalysis_cavityPhase.m - 4Pi-BRAINSPOT toolbox/
Dynamic_model_update/ , MATLAB, 88 linesanalysis_cavityPhase_and _misObj_v2.m - 4Pi-BRAINSPOT toolbox/
Dynamic_model_update/ , MATLAB, 121 linescaliphi0_NCC_v3.m - 4Pi-BRAINSPOT toolbox/
Dynamic_model_update/ , MATLAB, 126 lines, 1 matchcaliphi0_and_misObj_NCC_ v3.m - 4Pi-BRAINSPOT toolbox/
Dynamic_model_update/ , MATLAB, 50 linesest_misobj_from_hist_v2. m - 4Pi-BRAINSPOT toolbox/
Dynamic_model_update/ , MATLAB, 47 linesest_phi0_from_hist_v2.m - 4Pi-BRAINSPOT toolbox/
Dynamic_model_update/ , MATLAB, 143 linesfind_4pi_subregion.m - 4Pi-BRAINSPOT toolbox/
Dynamic_model_update/ , MATLAB, 52 linesgenIniguess.m - 4Pi-BRAINSPOT toolbox/
Dynamic_model_update/ , MATLAB, 130 linesgenini4pi_misalign_parfo r_v4.m - 4Pi-BRAINSPOT toolbox/
Dynamic_model_update/ , MATLAB, 112 linesgenini4pi_phi0_parfor_v3 .m - 4Pi-BRAINSPOT toolbox/
Segmentation/ , MATLAB, 119 linescrop_subregion_4pi.m - 4Pi-BRAINSPOT toolbox/
brainspot_4pi_GUI.m , MATLAB, 3,576 lines, 1 match - 4Pi-BRAINSPOT toolbox/
export2csv.m , MATLAB, 59 lines - 4Pi-BRAINSPOT toolbox/
genPupilfigs_4ch.m , MATLAB, 160 lines - 4Pi-BRAINSPOT toolbox/
main.m , MATLAB, 25 lines - 4Pi-BRAINSPOT toolbox/
srhist_color.m , MATLAB, 46 lines - Support/
Helpers/ , MATLAB, 13 linesfindcoord_seg.m - Support/
Helpers/ , MATLAB, 4 linesfindmax.m - Support/
Helpers/ , MATLAB, 6 linesfindmax3d.m - Support/
Helpers/ , MATLAB, 5 linesimstretch_linear.m - Support/
Helpers/ , MATLAB, 138 linesnormxcorr3_sparse.m - Support/
Helpers/ , MATLAB, 6 linesunif_img.m - Support/
PSF Toolbox_4pi/ , MATLAB, 105 lines@OTFrescale/ OTFrescale.m - Support/
PSF Toolbox_4pi/ , MATLAB, 359 lines, 2 matches@PSF_4pi/ PSF_4pi.m - Support/
PSF Toolbox_4pi/ , MATLAB, 1,192 lines@Zernike_Polynomials/ Zernike_Polynomials.m - Support/
SRsCMOS/ , MATLAB, 28 linesSRreconstructhist.m - Support/
SRsCMOS/ , MATLAB, 166 linesSRscmos.m - Support/
SRsCMOS/ , MATLAB, 32 linesdriftcorrection_Redun3D. m - Support/
SRsCMOS/ , MATLAB, 124 linesdriftcorrection_Redun_co re2D_parfor.m - Support/
SRsCMOS/ , MATLAB, 162 linesdriftcorrection_Redun_co re3D_parfor.m - Support/
SRsCMOS/ , MATLAB, 43 linesdriftcorrection_core2D.m - Support/
SRsCMOS/ , MATLAB, 65 linesdriftcorrection_core3D.m - Support/
SRsCMOS/ , MATLAB, 16 linesgencormask.m - Support/
SRsCMOS/ , MATLAB, 29 linesshiftcoords.m - Support/
SRsCMOS/ , MATLAB, 8 linesshiftcoords_stack.m - Support/
SRsCMOS/ , MATLAB, 149 linesstackalignment_Redun3D.m - Support/
SRsCMOS/ , MATLAB, 62 linesstackalignment_core3D.m - LICENSE.md, License, 1 line
- README.md, Text, 96 lines
Code availability statement
The paper has a code availability statement. Its license (CC BY-NC-ND) does not allow reproducing it here; in short, from what the harvester recognized in it:
- it points to the authors' code: HuanglabPurdue/
4Pi-BRAINSPOT - it says that the code is available on request
Read it in the paper: doi.org/10.1038/s41467-026-71614-6.
Tracing map
Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.
What the map holds:
- 1 repository of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 77 scripts, each with its path and the digest of its content;
- 8 matches between paragraphs of the paper and lines of the code (method lexical-v1);
- neither the text of the paper nor the code itself.
Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.
Data
Datasets cited
- doi:10.6084/
m9.figshare.c.7721774 , at figshare; found in “Data availability” - figshare:28603676, at figshare; found in DataCite
Data availability statement
The paper has a data availability statement. Its license (CC BY-NC-ND) does not allow reproducing it here; in short, from what the harvester recognized in it:
- it points to a dataset: figshare 10.6084/
m9.figshare.c.7721774 - it says that the data are available on request
Read it in the paper: doi.org/10.1038/s41467-026-71614-6.
Versions
The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.
Version 1, 29 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 10 authors, 3 keywords, 7 MeSH terms, 4 funders, 67 references.
Cite
This paper
Gao, H.-C., Xu, F., Cheng, X., Chen, T., Bi, C., Zheng, Y., Li, Y., Li, Y., Chubykin, A. A., & Huang, F. (2026). Interferometric ultra-high resolution 3D imaging through brain sections. Nature communications, 17(1), 5550. https://
BibTeX
@article{gao2026interfer
author = {Gao, Hao-Cheng and Xu, Fan and Cheng, Xi and Chen, Tailong and Bi, Cheng and Zheng, Yue and Li, Yilun and Li, Yumian and Chubykin, Alexander A and Huang, Fang},
title = {{Interferometric ultra-high resolution 3D imaging through brain sections}},
journal = {Nature communications},
year = {2026},
month = apr,
volume = {17},
number = {1},
pages = {5550},
publisher = {Nature Publishing Group},
issn = {2041-1723},
doi = {10.1038/
url = {https://
pmid = {42020380},
pmcid = {PMC13287808}
}
RIS
TY - JOUR
AU - Gao, Hao-Cheng
AU - Xu, Fan
AU - Cheng, Xi
AU - Chen, Tailong
AU - Bi, Cheng
AU - Zheng, Yue
AU - Li, Yilun
AU - Li, Yumian
AU - Chubykin, Alexander A
AU - Huang, Fang
TI - Interferometric ultra-high resolution 3D imaging through brain sections
T2 - Nature communications
J2 - Nat Commun
PY - 2026
DA - 2026/
VL - 17
IS - 1
SP - 5550
SN - 2041-1723
PB - Nature Publishing Group
DO - 10.1038/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1038/
"type": "article-journal",
"title": "Interferometric ultra-high resolution 3D imaging through brain sections",
"container-title": "Nature communications",
"author": [
{
"family": "Gao",
"given": "Hao-Cheng"
},
{
"family": "Xu",
"given": "Fan"
},
{
"family": "Cheng",
"given": "Xi"
},
{
"family": "Chen",
"given": "Tailong"
},
{
"family": "Bi",
"given": "Cheng"
},
{
"family": "Zheng",
"given": "Yue"
},
{
"family": "Li",
"given": "Yilun"
},
{
"family": "Li",
"given": "Yumian"
},
{
"family": "Chubykin",
"given": "Alexander A"
},
{
"family": "Huang",
"given": "Fang"
}
],
"container-title-short":
"volume": "17",
"issue": "1",
"page": "5550",
"DOI": "10.1038/
"PMID": "42020380",
"PMCID": "PMC13287808",
"ISSN": "2041-1723",
"publisher": "Nature Publishing Group",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
4,
22
]
]
}
}
The tracing map gets a citation of its own once an author has validated it and it has a DOI.
Similar papers
The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.
- [1] doi:10.1038/s41467-026-73045-9 [code]
- Aberration-aware 3D localization microscopy via self-supervised neural-physics learning.Journal: Nature communicationsIn common: histology / microscopy, cellular / molecular, 7 references
- [2] doi:10.1038/s41467-026-69798-y [code]
- Dual deconvolution in multiphoton structured illumination microscopy for deep-tissue super-resolution imaging.Journal: Nature communicationsIn common: Image Processing Toolbox, histology / microscopy, mouse, 5 references
- [3] doi:10.1016/j.celrep.2026.117590 [code]
- Impaired behavioral inhibition in Fmr1 KO mice is linked to disrupted visual cortex theta oscillations.Journal: Cell reportsIn common: Image Processing Toolbox, Statistics and Machine Learning Toolbox, mouse, 1 reference, author Alexander A Chubykin
- [4] doi:10.7554/elife.109083
- Spine nanostructure profiling of cultured neurons from mouse models reveals a schizophrenia-linked role for Ecrg4.Journal: eLifeIn common: histology / microscopy, mouse, cellular / molecular, 4 references
- [5] doi:10.1016/j.crmeth.2026.101472 [code]
- Rapid neuronal labeling and functional imaging in the developing mouse brain with AAV-PHP.eB.Journal: Cell reports methodsIn common: Curve Fitting Toolbox, Image Processing Toolbox, Statistics and Machine Learning Toolbox, histology / microscopy, mouse, 1 reference
- [6] doi:10.1016/j.isci.2026.116769 [code]
- Graph-based modeling of optical system enables adaptive optics on dynamic samples with self-calibration.Journal: iScienceIn common: Curve Fitting Toolbox, Image Processing Toolbox, 2 references
- [7] doi:10.1038/s41467-026-76146-7 [code]
- Resolving endogenous protein organization in cells with nanometer resolution.Journal: Nature communicationsIn common: histology / microscopy, cellular / molecular, 3 references
- [8] doi:10.1016/j.isci.2026.117187 [code]
- Functional and structural characterization of dendritic spine pathology in a mouse model of tauopathy.Journal: iScienceIn common: Curve Fitting Toolbox, Image Processing Toolbox, Statistics and Machine Learning Toolbox, histology / microscopy, mouse
- [9] doi:10.1002/glia.70141 [code]
- Conservation of Neuron-Astrocyte Correlated Activity in Developing Sensory Pathways.Journal: GliaIn common: Curve Fitting Toolbox, Image Processing Toolbox, Statistics and Machine Learning Toolbox, mouse, cellular / molecular
- [10] doi:10.1113/jp289290 [code]
- Robust activity-dependent mitochondrial calcium dynamics at the AIS is dispensable for action potential generation.Journal: The Journal of physiologyIn common: Curve Fitting Toolbox, Image Processing Toolbox, Statistics and Machine Learning Toolbox, mouse, cellular / molecular
Contribute
The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.
Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 1 repository of the authors' code, each at its verified commit and with its license, 77 scripts, and 8 matches between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:4240acc64e0779a4…
Add the badge to its README
The badge links the code to this page. Copy one of these into the README of the paper's code: only you decide where it goes, and nothing is changed for you.
Markdown
[, paste the snippet at the top, then “Commit changes…” and, to review it first, “Create a new branch and start a pull request”. You open the pull request; OSCR asks for no permission.
Request its removal
To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).
Discussion, reproductions, activity
Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.
Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
