def test_middle_point_with_GRS80(self):
# Same as test_middle_point(), but with a different reference system, GRS80; the
# result should be the same.
point = hspatial.coordinates2point(325077, 4205177, srid=2100)
self.assertAlmostEqual(hspatial.extract_point_from_raster(
point, self.fp),
2.2,
places=2)
def test_with_prefix(self):
# Same as test_with_list_of_files(), but with prefix.
point = hspatial.coordinates2point(22.01001, 37.98001)
prefix = os.path.join(self.tempdir, "test")
ts = hspatial.PointTimeseries(point,
prefix=prefix,
default_time=dt.time(23, 58)).get()
self._check_against_expected(ts)
def test_bottom_left_point_with_GRS80(self):
# Same as test_bottom_left_point(), but with a different reference system,
# GRS80; the result should be the same.
point = hspatial.coordinates2point(324076, 4205176, srid=2100)
self.assertAlmostEqual(hspatial.extract_point_from_raster(
point, self.fp),
3.1,
places=2)
def test_bottom_left_point(self):
# Use almost exactly same point as test_middle_point(), only slightly altered
# so that we get bottom left point instead.
point = hspatial.coordinates2point(22.00999, 37.97999)
self.assertAlmostEqual(hspatial.extract_point_from_raster(
point, self.fp),
3.1,
places=2)
def test_middle_point(self):
# We use co-ordinates almost to the common corner of the four lower left points,
# only a little bit towards the center.
point = hspatial.coordinates2point(22.01001, 37.98001)
self.assertAlmostEqual(hspatial.extract_point_from_raster(
point, self.fp),
2.2,
places=2)
def test_with_list_of_files(self):
# Use co-ordinates almost to the common corner of the four lower left points,
# and only a little bit towards the center.
point = hspatial.coordinates2point(22.01001, 37.98001)
filenames = [
os.path.join(self.tempdir, "test-2014-11-22-16-1.tif"),
os.path.join(self.tempdir, "test-2014-11-21-16-1.tif"),
os.path.join(self.tempdir, "test-2014-11-23-16-1.tif"),
]
ts = hspatial.PointTimeseries(point, filenames=filenames).get()
self._check_against_expected(ts)
def test_transform_does_not_modify_srid_of_gdal_point(self):
pppoint = hspatial.PassepartoutPoint(
hspatial.coordinates2point(324651, 4205742, srid=2100))
original_spatial_reference = pppoint.point.GetSpatialReference(
).ExportToWkt()
sr = osr.SpatialReference()
sr.ImportFromEPSG(4326)
pppoint.transform_to(sr.ExportToWkt())
self.assertEqual(
pppoint.point.GetSpatialReference().ExportToWkt(),
original_spatial_reference,
)
def setUp(self):
super().setUp()
self.point = hspatial.coordinates2point(22.01001, 37.98001)
self.prefix = os.path.join(self.tempdir, "test")
self.dest = os.path.join(self.tempdir, "dest.hts")
def test_unit_of_measurement(self):
point = hspatial.coordinates2point(22.01001, 37.98001)
prefix = os.path.join(self.tempdir, "test")
ts = hspatial.PointTimeseries(point, prefix=prefix).get()
self.assertEqual(ts.unit, "microkernels")
def test_with_prefix_and_geodjango(self):
point = hspatial.coordinates2point(22.01001, 37.98001)
prefix = os.path.join(self.tempdir, "test")
ts = hspatial.PointTimeseries(point, prefix=prefix).get()
self._check_against_expected(ts)
def test_fails_gracefully_when_osr_point_is_really_outside_crs_limits(
self):
point = hspatial.coordinates2point(125.0, 85.0)
with self.assertRaises(RuntimeError):
hspatial.extract_point_from_raster(point, self.fp)
def test_point_outside_raster(self):
point = hspatial.coordinates2point(21.0, 38.0)
with self.assertRaises(RuntimeError):
hspatial.extract_point_from_raster(point, self.fp)
def test_does_not_modify_srid_of_point(self):
point = hspatial.coordinates2point(325077, 4205177, srid=2100)
original_spatial_reference = point.GetSpatialReference().ExportToWkt()
hspatial.extract_point_from_raster(point, self.fp)
self.assertEqual(point.GetSpatialReference().ExportToWkt(),
original_spatial_reference)
def test_top_middle_point(self):
point = hspatial.coordinates2point(22.015, 37.995)
self.assertTrue(
math.isnan(hspatial.extract_point_from_raster(point, self.fp)))
def test_top_left_point(self):
point = hspatial.coordinates2point(22.005, 37.995)
self.assertAlmostEqual(hspatial.extract_point_from_raster(
point, self.fp),
1.1,
places=2)