Skip to content

Commit daf3608

Browse files
committed
Use DiracBasis3D for jsb examples
1 parent 61a25a8 commit daf3608

2 files changed

Lines changed: 8 additions & 2 deletions

File tree

gallery/experiments/experimental_abinitio_pipeline_10028_jsb.py

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -29,6 +29,7 @@
2929
import logging
3030
from pathlib import Path
3131

32+
from aspire.basis import DiracBasis3D
3233
from aspire.denoising import LegacyClassAvgSource
3334
from aspire.reconstruction import MeanEstimator
3435
from aspire.source import OrientedSource, RelionSource
@@ -156,7 +157,9 @@
156157
logger.info("Begin Volume reconstruction")
157158

158159
# Set up an estimator to perform the backprojection.
159-
estimator = MeanEstimator(oriented_src)
160+
# Legacy MATLAB FIRM used Dirac basis.
161+
basis3d = DiracBasis3D(oriented.L)
162+
estimator = MeanEstimator(oriented_src, basis=basis3d)
160163

161164
# Perform the estimation and save the volume.
162165
estimated_volume = estimator.estimate()

gallery/experiments/experimental_abinitio_pipeline_10073_jsb.py

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -32,6 +32,7 @@
3232
import numpy as np
3333

3434
from aspire.abinitio import CLSync3N
35+
from aspire.basis import DiracBasis3D
3536
from aspire.denoising import LegacyClassAvgSource
3637
from aspire.reconstruction import MeanEstimator
3738
from aspire.source import ArrayImageSource, OrientedSource, RelionSource
@@ -160,7 +161,9 @@
160161
logger.info("Begin Volume reconstruction")
161162

162163
# Set up an estimator to perform the backprojection.
163-
estimator = MeanEstimator(oriented_src)
164+
# Legacy MATLAB FIRM used Dirac basis.
165+
basis3d = DiracBasis3D(oriented.L)
166+
estimator = MeanEstimator(oriented_src, basis=basis3d)
164167

165168
# Perform the estimation and save the volume.
166169
estimated_volume = estimator.estimate()

0 commit comments

Comments
 (0)