def test_reproject_gcps_auto(self): ds = gdal.Open(self.test_file) d = Domain(ds=ds) d.reproject_gcps() gcpproj = NSR(d.vrt.dataset.GetGCPProjection()) self.assertEqual(gcpproj.GetAttrValue('PROJECTION'), 'Stereographic')
def test_reproject_gcps(self): ds = gdal.Open(self.test_file) d = Domain(ds=ds) d.reproject_gcps('+proj=stere +datum=WGS84 +ellps=WGS84 +lat_0=75 +lon_0=10 +no_defs') gcp = d.vrt.dataset.GetGCPs()[0] self.assertTrue(gcp.GCPX > 636161) self.assertTrue(gcp.GCPY < -288344)