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