def import_geotiff(file_path):
    """
    Import a GeoTIFF to a numpy array and prints
    information of the Coordinate Reference System (CRS).

    **Requires `osgeo`.**

    Args:
        file_path : str
            path to the GeoTIFF

    Returns:
        data : 2D array
        extent : tuple
            bounding box in the projection of the GeoTIFF
            e.g. [xmin, xmax, ymin, ymax]
    """
    from osgeo import gdal, osr
    gtiff = gdal.Open(file_path)
    data = gtiff.ReadAsArray()
    gt = gtiff.GetGeoTransform()
    gtproj = gtiff.GetProjection()
    inproj = osr.SpatialReference()
    inproj.ImportFromWkt(gtproj)
    gtextent = gt[0], gt[0] + gtiff.RasterXSize * gt[1], gt[3], gt[3
        ] + gtiff.RasterYSize * gt[5]
    print(inproj)
    gtiff = None
    return data, gtextent
