def test_vertices_in_crs2(self): point = Point((-123.0, 49.0), crs=LonLatWGS84) BCAlbers = ProjectedCRS("+proj=aea +lat_1=50 +lat_2=58.5 +lat_0=45 " "+lon_0=-126 +x_0=1000000 +y_0=0 +ellps=GRS80 +datum=NAD83 " "+units=m +no_defs", "BC Albers") self.assertTupleAlmostEqual(point.get_vertex(BCAlbers), (1219731.770879303, 447290.49891930853)) return
def test_walk_albers_geodetic(self): AlaskaAlbers = ProjectedCRS("+proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 " "+x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs", "+ellps=GRS80") start = Point((-2658638, 2443580), crs=AlaskaAlbers) dest = start.walk(4500, 195.0, projected=False) self.assertAlmostEqual(dest.x, -2662670.889, places=3) self.assertAlmostEqual(dest.y, 2441551.155, places=3)
def test_walk_albers_projected(self): AlaskaAlbers = ProjectedCRS("+proj=aea +lat_1=55 +lat_2=65 +lat_0=50 +lon_0=-154 " "+x_0=0 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=m +no_defs", "+ellps=GRS80") start = Point((-2658638, 2443580), crs=AlaskaAlbers) dest = start.walk(4500, 195.0) self.assertAlmostEqual(dest.x, -2659802.686, places=3) self.assertAlmostEqual(dest.y, 2439233.334, places=3) return
def test_vertices_in_crs3(self): line = Line([(2.0, 34.0), (2.15, 34.2), (2.7, 34.1)], crs=LonLatWGS84) UTM31N = ProjectedCRS( "+proj=utm +zone=31 +ellps=WGS84 " "+datum=WGS84 +units=m +no_defs", "UTM 31N (WGS 84)") ans = [[407650.39665729, 3762606.65987638], [421687.71905897, 3784658.46708431], [472328.10951276, 3773284.48524179]] for v0, v1 in zip(line.get_vertices(UTM31N), ans): self.assertTupleAlmostEqual(v0, v1, places=6) return