def latlon_xyz(latitude, longitude, depth, latitude_0, longitude_0,
    latitude_1, latitude_2):
    """Converts latitude and longitude into x y z

    Uses pyproj abd lambert projection to convert the coordinates using a given origin.

    Args
        latitude: float latitude in degrees
        longitude: float longitude in degrees
        depth: float depth in km
        latitude_0: float latitude of origin in degrees
        longitude_0: float longitude of origin in degrees 
        latitude_1: float latitude of first standard parallel in degrees 
        latitude_2: float latitude of second standard parallel in degrees 

    Returns
        (float,float,float): tuple of x, y, z coordinates in km from latitude_0 and longitude_0 origin
    """
    from pyproj import Proj
    p = Proj('+proj=lcc +lat_0=' + str(latitude_0) + ' +lon_0=' + str(
        longitude_0) + ' +lat_1=' + str(latitude_1) + ' +lat_2=' + str(
        latitude_2), ellps='WGS84')
    [x_0, y_0] = p(float(longitude_0), float(latitude_0))
    [x, y] = p(longitude, latitude)
    return (x - x_0) / 1000.0, (y - y_0) / 1000.0, depth
