Lawyer-Fujita Algorithm

read 3ds file using MATLAB / Python

Lawler-Fujita Algorithm for STM image process

###################################
# load data from 3ds file
###################################
filename = 'mapping--002.3ds'
Im = 1j
grid = nap.read.Grid(filename)
size_x = grid.header['size_xy'][0]
size_y = grid.header['size_xy'][1]
lines = grid.header['dim_px'][0]
points = grid.header['dim_px'][1]
nm_to_m = 10**(-9)
im = grid.signals['topo']
palette = spiepy.NANOMAP
palette.set_bad('r',1.0)
fig, ax = plt.subplots()
ax.imshow(im, cmap=palette, origin='lower')
plt.show()