code: def lorentzian2d(
      x,
      y,
      amplitude=1.0,
      centerx=0.0,
      centery=0.0,
      sigmax=1.0,
      sigmay=1.0,
      rotation=0,
  ):
      """Return a two dimensional lorentzian.
  
      The maximum of the peak occurs at ``centerx`` and ``centery``
      with widths ``sigmax`` and ``sigmay`` in the x and y directions
      respectively. The peak can be rotated by choosing the value of ``rotation``
      in radians.
      """
      xp = (x - centerx) * np.cos(rotation) - (y - centery) * np.sin(
          rotation
      )
      yp = (x - centerx) * np.sin(rotation) + (y - centery) * np.cos(
          rotation
      )
      R = (xp / sigmax) ** 2 + (yp / sigmay) ** 2
  
      return 2 * amplitude * lorentzian(R) / (np.pi * sigmax * sigmay)
  
imports: import matplotlib.pyplot
  import numpy
  from scipy.interpolate import griddata
  import lmfit
  from lmfit.lineshapes import gaussian2d, lorentzian
label_desc: Use image node_id 1 for calling a global function with example usage: lorentzian2d(x, y, amplitude, centerx, centery, sigmax, sigmay, rotation) and returns: unknown