def osr_compd_6(): if not osr_proj4.have_proj480(): return 'skip' srs = osr.SpatialReference() srs.SetFromUserInput( '+proj=utm +zone=11 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +geoidgrids=g2003conus.gtx,g2003alaska.gtx,g2003h01.gtx,g2003p01.gtx +vunits=us-ft +no_defs ' ) if srs.Validate() != 0: gdaltest.post_reason('Does not validate') return 'fail' exp_wkt = """COMPD_CS["UTM Zone 11, Northern Hemisphere + Unnamed Vertical Datum", PROJCS["UTM Zone 11, Northern Hemisphere", GEOGCS["GRS 1980(IUGG, 1980)", DATUM["unknown", SPHEROID["GRS80",6378137,298.257222101], TOWGS84[0,0,0,0,0,0,0]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",-117], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",0], UNIT["Meter",1]], VERT_CS["Unnamed", VERT_DATUM["Unnamed",2005, EXTENSION["PROJ4_GRIDS","g2003conus.gtx,g2003alaska.gtx,g2003h01.gtx,g2003p01.gtx"]], UNIT["Foot_US",0.3048006096012192], AXIS["Up",UP]]]""" wkt = srs.ExportToPrettyWkt() if gdaltest.equal_srs_from_wkt(exp_wkt, wkt) == 0: gdaltest.post_reason('fail') return 'fail' elif exp_wkt != wkt: print('warning they are equivalent, but not completely the same') print(wkt) if wkt.find('g2003conus.gtx') == -1: gdaltest.post_reason('Did not get PROJ4_GRIDS EXTENSION node') return 'fail' exp_proj4 = '+proj=utm +zone=11 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +geoidgrids=g2003conus.gtx,g2003alaska.gtx,g2003h01.gtx,g2003p01.gtx +vunits=us-ft +no_defs ' proj4 = srs.ExportToProj4() if proj4 != exp_proj4: gdaltest.post_reason('Did not get expected proj.4 string, got:' + proj4) return 'fail' return 'success'
def osr_compd_6(): if not osr_proj4.have_proj480(): return 'skip' srs = osr.SpatialReference() srs.SetFromUserInput( '+proj=utm +zone=11 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +geoidgrids=g2003conus.gtx,g2003alaska.gtx,g2003h01.gtx,g2003p01.gtx +vunits=us-ft +no_defs ' ) if srs.Validate() != 0: gdaltest.post_reason( 'Does not validate' ) return 'fail' exp_wkt = """COMPD_CS["UTM Zone 11, Northern Hemisphere + Unnamed Vertical Datum", PROJCS["UTM Zone 11, Northern Hemisphere", GEOGCS["GRS 1980(IUGG, 1980)", DATUM["unknown", SPHEROID["GRS80",6378137,298.257222101], TOWGS84[0,0,0,0,0,0,0]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",-117], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",500000], PARAMETER["false_northing",0], UNIT["Meter",1]], VERT_CS["Unnamed", VERT_DATUM["Unnamed",2005, EXTENSION["PROJ4_GRIDS","g2003conus.gtx,g2003alaska.gtx,g2003h01.gtx,g2003p01.gtx"]], UNIT["Foot_US",0.3048006096012192], AXIS["Up",UP]]]""" wkt = srs.ExportToPrettyWkt() if gdaltest.equal_srs_from_wkt( exp_wkt, wkt ) == 0: gdaltest.post_reason('fail') return 'fail' elif exp_wkt != wkt: print('warning they are equivalent, but not completely the same') print(wkt) if wkt.find('g2003conus.gtx') == -1: gdaltest.post_reason( 'Did not get PROJ4_GRIDS EXTENSION node' ) return 'fail' exp_proj4 = '+proj=utm +zone=11 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +geoidgrids=g2003conus.gtx,g2003alaska.gtx,g2003h01.gtx,g2003p01.gtx +vunits=us-ft +no_defs ' proj4 = srs.ExportToProj4() if proj4 != exp_proj4: gdaltest.post_reason( 'Did not get expected proj.4 string, got:' + proj4 ) return 'fail' return 'success'