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