thomasaarholt

def P1(P0=(0,0), Pc=(2,2), step=1):
    "Get the point P1, the point to which one moves from P0 with a fixed step in the direction of Pc"
    X0, Y0 = P0
    Xc, Yc = Pc
    
    diffX = Xc - X0
    diffY = Yc - Y0
    
    magnitude = (diffX**2 + diffY**2)**0.5
    

thomasaarholt

from pprint import pprint
import temscript as ts
from math import pi

dwelltime = 1e-5 # seconds

tem = ts.Microscope()
illumination = tem._tem_illumination
instrument = ts.GetInstrument()

thomasaarholt

import hyperspy.api as hs
expr = "where(abs(x - center) < radius, A, 0)"
tophat = hs.model.components1D.Expression(expr, 'Top Hat', position="center", A=1, center=0, radius=1)

thomasaarholt

import hyperspy.api as hs
expr = "where(sqrt((x - xc)**2 + (y - yc)**2) < radius, A, 0)"
tophat = hs.model.components2D.Expression(expr, 'Top Hat', position=("xc", "yc"), A=1, xc=0, yc=0, radius=1)

thomasaarholt

import hyperspy.api as hs
expr = "where(And(left < x, x < right), A, 0)"
tophat = hs.model.components1D.Expression(expr, 'Top Hat', position="left", A=1, left=0, right=1)

thomasaarholt

def loadXRD(filename):
    with open(filename, 'r') as f:
        lines = f.readlines()
    angles = []
    intensity = []
    for line in lines[1:-2]:
        x, y = line.split('\n')[0].split(",")[:-1]
        angles.append(float(x))
        intensity.append(float(y))
    x = np.array(angles)

thomasaarholt

import numpy as np
from scipy.signal import convolve2d, gaussian
from scipy.misc import face
import matplotlib.pyplot as plt

def convolve2d_fft(arr1, arr2):
    s0, s1 = arr1.shape
    
    conv = np.fft.irfft2(
        np.fft.rfft2(arr1) * np.fft.rfft2(arr2), 

thomasaarholt

$ pip install pytetgen
Collecting pytetgen
  Using cached pytetgen-0.1.4.tar.gz (392 kB)
Building wheels for collected packages: pytetgen
  Building wheel for pytetgen (setup.py) ... error
  ERROR: Command errored out with exit status 1:
   command: 'C:\Users\thomasaar\Miniconda3\envs\py38\python.exe' -u -c 'import sys, setuptools, tokenize; sys.argv[0] = '"'"'C:\\Users\\thomasaar\\AppData\\Local\\Temp\\pip-install-xphrb4j7\\pytetgen_dc8134ee876b4f18b1a1373353bd0c09\\setup.py'"'"'; __file__='"'"'C:\\Users\\thomasaar\\AppData\\Local\\Temp\\pip-install-xphrb4j7\\pytetgen_dc8134ee876b4f18b1a1373353bd0c09\\setup.py'"'"';f=getattr(tokenize, '"'"'open'"'"', open)(__file__);code=f.read().replace('"'"'\r\n'"'"', '"'"'\n'"'"');f.close();exec(compile(code, __file__, '"'"'exec'"'"'))' bdist_wheel -d 'C:\Users\thomasaar\AppData\Local\Temp\pip-wheel-igdis3c8'
       cwd: C:\Users\thomasaar\AppData\Local\Temp\pip-install-xphrb4j7\pytetgen_dc8134ee876b4f18b1a1373353bd0c09\
  Complete output (339 lines):
  running bdist_whee

thomasaarholt

import numpy as np
import matplotlib
import matplotlib.pyplot as plt
Arc = matplotlib.patches.Arc

def halfangle(a, b):
    "Gets the middle angle between a and b, when increasing from a to b"
    if b < a:
        b += 360
    return (a + b)/2 % 360

thomasaarholt

def row_roll(arr, shifts, axis=1, fill=np.nan):
    """Apply an independent roll for each dimensions of a single axis.

    Parameters
    ----------
    arr : np.ndarray
        Array of any shape.

    shifts : np.ndarray, dtype int. Shape: `(arr.shape[:axis],)`.
        Amount to roll each row by. Positive shifts row right.