docs: def get_profile(self, img_type, coordinate, num_points):
          ''' Extract a profile from (lat1,lon1) to (lat2,lon2)
  
          Args:
              img_type (str): Either lola or wac.
              coordinate (float,float,float,flaot): A tupple
                  ``(lon0,lon1,lat0,lat1)`` with:
  
                  - lon0: First point longitude
                  - lat0: First point latitude
                  - lon1: Second point longitude
                  - lat1: Second point latitude
  
              num_points (int): Number of points to use in the
                  interpolation process.
  
          Note:
              Be carefull, longitude has to be in between 0-360 !
          '''
  
          lon0, lon1, lat0, lat1 = coordinate
          X, Y, Z = self.get_arrays(img_type)
          y0, x0 = np.argmin(np.abs(X[0, :] - lon0)
                             ), np.argmin(np.abs(Y[:, 0] - lat0))
          y1, x1 = np.argmin(np.abs(X[0, :] - lon1)
                             ), np.argmin(np.abs(Y[:, 0] - lat1))
          x, y = np.linspace(x0, x1, num_points), np.linspace(y0, y1, num_points)
          zi = scipy.ndimage.map_coordinates(Z, np.vstack((x, y)))
  
          return zi