def setUp(self): v = peaks(500)[:100, :] utm7 = karta.crs.ProjectedCRS("+proj=utm +zone=7 +north +datum=WGS84", "UTM 7N (WGS 84)") g = karta.RegularGrid([15.0, 15.0, 30.0, 30.0, 0.0, 0.0], v, crs=utm7) fpath = os.path.join(TMPDATA, "test.tif") g.to_geotiff(fpath, compress=None) self.grid = karta.read_geotiff(fpath, in_memory=False)
def setUp(self): v = peaks(500)[:100,:] utm7 = karta.crs.ProjectedCRS("+proj=utm +zone=7 +north +datum=WGS84", "UTM 7N (WGS 84)") g = karta.RegularGrid([15.0, 15.0, 30.0, 30.0, 0.0, 0.0], v, crs=utm7) fpath = os.path.join(TMPDATA, "test.tif") g.to_geotiff(fpath, compress=None) self.grid = karta.read_geotiff(fpath, in_memory=False)
def test_write_read(self): # try writing a file, then read it back in and verify that it matches v = peaks(500)[:100, :] utm7 = karta.crs.ProjectedCRS("+proj=utm +zone=7 +north +datum=WGS84", "UTM 7N (WGS 84)") g = karta.RegularGrid([15.0, 15.0, 30.0, 30.0, 0.0, 0.0], v, crs=utm7) fpath = os.path.join(TMPDATA, "test.tif") g.to_geotiff(fpath, compress=None) gnew = karta.read_geotiff(fpath) self.assertTrue("+proj=utm" in gnew.crs.get_proj4()) self.assertTrue("+zone=7" in gnew.crs.get_proj4()) self.assertEqual(g.transform, gnew.transform) self.assertEqual(g.values.dtype, gnew.values.dtype) self.assertTrue(np.all(g[:, :] == gnew[:, :])) return
def test_write_read_disk(self): # try writing a file, then open it without loading into memory and verify v = peaks(500)[:100,:] utm7 = karta.crs.ProjectedCRS("+proj=utm +zone=7 +north +datum=WGS84", "UTM 7N (WGS 84)") g = karta.RegularGrid([15.0, 15.0, 30.0, 30.0, 0.0, 0.0], v, crs=utm7) fpath = os.path.join(TMPDATA, "test.tif") g.to_geotiff(fpath, compress=None) gnew = karta.read_geotiff(fpath, in_memory=False) self.assertEqual(g.transform, gnew.transform) self.assertEqual(g.values.dtype, gnew.values.dtype) self.assertEqual(g.size, gnew.size) self.assertTrue(np.all(g[10:50:3, 15:45:2] == gnew[10:50:3, 15:45:2])) self.assertTrue(np.all(g[10:50:3, 45:15:-2] == gnew[10:50:3, 45:15:-2])) self.assertTrue(np.all(g[:,:] == gnew[:,:])) return
def test_write_read(self): # try writing a file, then read it back in and verify that it matches v = peaks(500)[:100,:] utm7 = karta.crs.ProjectedCRS("+proj=utm +zone=7 +north +datum=WGS84", "UTM 7N (WGS 84)") g = karta.RegularGrid([15.0, 15.0, 30.0, 30.0, 0.0, 0.0], v, crs=utm7) fpath = os.path.join(TMPDATA, "test.tif") ret = g.to_geotiff(fpath, compress=None) self.assertEqual(ret, g) gnew = karta.read_geotiff(fpath) self.assertTrue("+proj=utm" in gnew.crs.get_proj4()) self.assertTrue("+zone=7" in gnew.crs.get_proj4()) self.assertEqual(g.transform, gnew.transform) self.assertEqual(g.values.dtype, gnew.values.dtype) self.assertTrue(np.all(g[:,:] == gnew[:,:])) return
def test_write_read_disk(self): # try writing a file, then open it without loading into memory and verify v = peaks(500)[:100, :] utm7 = karta.crs.ProjectedCRS("+proj=utm +zone=7 +north +datum=WGS84", "UTM 7N (WGS 84)") g = karta.RegularGrid([15.0, 15.0, 30.0, 30.0, 0.0, 0.0], v, crs=utm7) fpath = os.path.join(TMPDATA, "test.tif") g.to_geotiff(fpath, compress=None) gnew = karta.read_geotiff(fpath, in_memory=False) self.assertEqual(g.transform, gnew.transform) self.assertEqual(g.values.dtype, gnew.values.dtype) self.assertEqual(g.size, gnew.size) self.assertTrue(np.all(g[10:50:3, 15:45:2] == gnew[10:50:3, 15:45:2])) self.assertTrue(np.all(g[10:50:3, 45:15:-2] == gnew[10:50:3, 45:15:-2])) self.assertTrue(np.all(g[:, :] == gnew[:, :])) return