def testOGCWKT(self): self.map.setWKTProjection('PROJCS["unnamed",PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["standard_parallel_1",65],PARAMETER["standard_parallel_2",55],PARAMETER["latitude_of_center",0],PARAMETER["longitude_of_center",-153],PARAMETER["false_easting",-4943910.68],PARAMETER["false_northing",0]]') proj4 = self.map.getProjection() assert proj4.find( '+proj=aea' ) != -1 assert proj4.find( '+ellps=WGS84' ) != -1 assert (mapscript.projectionObj(proj4)).getUnits() != mapscript.MS_DD
def testESRIWKT(self): self.map.setWKTProjection('ESRI::PROJCS["Pulkovo_1995_GK_Zone_2",GEOGCS["GCS_Pulkovo_1995",DATUM["D_Pulkovo_1995",SPHEROID["Krasovsky_1940",6378245,298.3]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Gauss_Kruger"],PARAMETER["False_Easting",2500000],PARAMETER["False_Northing",0],PARAMETER["Central_Meridian",9],PARAMETER["Scale_Factor",1],PARAMETER["Latitude_Of_Origin",0],UNIT["Meter",1]]') proj4 = self.map.getProjection() assert proj4.find( '+proj=tmerc' ) != -1 assert proj4.find( '+ellps=krass' ) != -1 assert (mapscript.projectionObj(proj4)).getUnits() != mapscript.MS_DD
def testWellKnownGEOGCS(self): self.map.setWKTProjection('WGS84') proj4 = self.map.getProjection() assert proj4.find('+proj=longlat') != -1, proj4 assert proj4.find('+datum=WGS84') != -1, proj4 assert ( mapscript.projectionObj(proj4)).getUnits() != mapscript.MS_METERS
def testWellKnownGEOGCS(self): self.map.setWKTProjection('WGS84') proj4 = self.map.getProjection() assert proj4.find( '+proj=longlat' ) != -1, proj4 assert proj4.find( '+ellps=WGS84' ) != -1, proj4 assert (mapscript.projectionObj(proj4)).getUnits() != mapscript.MS_METERS