def METHOD_NAME(self):
    cube = self._cube(GeogCS(6377000))
    with self.temp_filename(".tif") as temp_filename:
        export_geotiff(cube, temp_filename)
        dataset = gdal.Open(temp_filename, gdal.GA_ReadOnly)
        projection_string = dataset.GetProjection()
        # String has embedded floating point values,
        # Test with values to N decimal places, using a regular expression.
        re_pattern = (
            r'GEOGCS\["unknown",DATUM\["unknown",'
            r'SPHEROID\["unknown",637....,0\]\],PRIMEM\["Greenwich",0\],'
            r'UNIT\["degree",0.01745[0-9]*,AUTHORITY\["EPSG","9122"\]\],'
            r'AXIS\["Latitude",NORTH\],AXIS\["Longitude",EAST\]\]'
        )
        re_exp = re.compile(re_pattern)
        self.assertIsNotNone(
            re_exp.match(projection_string),
            "projection string {!r} does not match {!r}".format(
                projection_string, re_pattern
            ),
        )