stggh

[etalon:.+Software/Python]$ git clone https://github.com/PyAbel/PyAbel.git
Cloning into 'PyAbel'...
remote: Enumerating objects: 7858, done.
remote: Counting objects: 100% (675/675), done.
remote: Compressing objects: 100% (372/372), done.
remote: Total 7858 (delta 401), reused 508 (delta 285), pack-reused 7183
Receiving objects: 100% (7858/7858), 22.54 MiB | 9.50 MiB/s, done.
Resolving deltas: 100% (5611/5611), done.
[etalon:.+Software/Python]$ mv PyAbel PyAbel-master
[etalon:.+Software/Python]$ python -V

stggh

import numpy as np
import cse
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit
import scipy.constants as const
from scipy.special import spherical_jn, spherical_yn

au = 0.52917720859
eh = 27.21138386

stggh

Name:       Intel Core i7-6700
uArch:      Skylake
Technology: 14nm
Max Freq:   4.000 GHz
Cores:      4 cores (8 threads)
AVX:        AVX,AVX2
FMA:        FMA3
L1i Size:   32KB (128KB Total)
L1d Size:   32KB (128KB Total)
L2 Size:    256KB (1MB Total)

stggh

import numpy as np
import abel
import matplotlib.pyplot as plt
from scipy.ndimage.interpolation import shift

fig, ax = plt.subplots(2, 2, figsize=(8, 8))

n = 201
n2 = n//2
geo_origin = (n2, n2)

stggh

01-28 19:28:19.603 18040 18040 W tech.ula: JIT profile information will not be recorded: profile file does not exits.
01-28 19:28:19.606 18040 18040 I Perf    : Connecting to perf service.
01-28 19:28:19.609 18040 18040 W ContextImpl: Failed to ensure /data/user/0/tech.ula/cache: mkdir failed: EACCES (Permission denied)
01-28 19:28:19.609 18040 18040 W ContextImpl: Failed to update user.inode_cache: stat failed: EACCES (Permission denied)
01-28 19:28:19.609 18040 18040 W ContextImpl: Failed to ensure /data/user_de/0/tech.ula/code_cache: mkdir failed: ENOENT (No such file or directory)
01-28 19:28:19.609 18040 18040 W ContextImpl: Failed to update user.inode_code_cache: stat failed: ENOENT (No such file or directory)
01-28 19:28:19.619   602   602 I hwservicemanager: getTransport: Cannot find entry vendor.qti.hardware.iop@2.0::IIop/default in either framework or device manifest.
01-28 19:28:19.619   971   971 E ANDR-IOP: IIop:: Iop HAL Service is not available.
01-28 19:28:19.620 18040 18059 E Perf    : Fail t

stggh

import numpy as np
import matplotlib.pyplot as plt
from pbasex import pbasex, loadG
from quadrant import foldQuadrant, resizeFolded
import abel

# Load images
im = np.loadtxt('O2-ANU1024.txt.bz2')
raw = abel.tools.center.center_image(im, square=True)
n = raw.shape[-1]

stggh

import numpy as np
import abel
import bz2
import anisotro

from scipy.signal import find_peaks, peak_widths
import matplotlib.pylab as plt

# O2- PES VM-image
imagefile = bz2.BZ2File('O2-ANU1024.txt.bz2')

stggh

# basex(var)      iabel(ms)           
    5             0.53462
    9             0.53325
   17             0.54541
   33             0.56431
   65             0.67819
  129             0.80625
  257             1.45382
  513             3.34333
 1025            14.87905

stggh

import abel
import matplotlib.pyplot as plt

im = abel.tools.analytical.SampleImage(1001, 'Ominus').image

proj = abel.Transform(im, method='basex', verbose=False,
                      direction='forward').transform

bas = abel.Transform(proj, method='basex', verbose=False,
                     angular_integration=True)

stggh

import numpy as np
import abel
import matplotlib.pyplot as plt

n = 101

fig, axs = plt.subplots(3, 2, sharex=True, sharey=True)

case = '1st-order-diff'
# case = 'gradient'