sij/mapnik
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

Merge branch 'master' of github.com:mapnik/mapnik into lazy-proj4

Browse source
This commit is contained in:
Dane Springmeyer 2013-01-28 15:27:18 -05:00
commit 6c911dc43c
1 changed files with 26 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

35
tests/python_tests/projection_test.py
View file

@ -5,6 +5,7 @@ from nose.tools import *
import mapnik import mapnik
import random import random
import math import math
from utilities import run_all
# Tests that exercise map projections. # Tests that exercise map projections.
@ -77,25 +78,41 @@ def test_proj_transform_between_init_and_literal():
one = mapnik.Projection('+init=epsg:4326') one = mapnik.Projection('+init=epsg:4326')
two = mapnik.Projection('+init=epsg:3857') two = mapnik.Projection('+init=epsg:3857')
tr1 = mapnik.ProjTransform(one,two) tr1 = mapnik.ProjTransform(one,two)
tr1b = mapnik.ProjTransform(two,one)
wgs84 = '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs' wgs84 = '+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs'
merc = '+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs +over' merc = '+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext +no_defs +over'
src = mapnik.Projection(wgs84) src = mapnik.Projection(wgs84)
dest = mapnik.Projection(merc) dest = mapnik.Projection(merc)
tr2 = mapnik.ProjTransform(src,dest) tr2 = mapnik.ProjTransform(src,dest)
tr2b = mapnik.ProjTransform(dest,src)
for x in xrange(-180,180,10): for x in xrange(-180,180,10):
for y in xrange(-60,60,10): for y in xrange(-60,60,10):
coord = mapnik.Coord(x,y) coord = mapnik.Coord(x,y)
merc_coord1 = tr1.forward(coord) merc_coord1 = tr1.forward(coord)
merc_coord2 = tr2.forward(coord) merc_coord2 = tr1b.backward(coord)
assert_almost_equal(merc_coord1.x,merc_coord2.x) merc_coord3 = tr2.forward(coord)
assert_almost_equal(merc_coord1.y,merc_coord2.y) merc_coord4 = tr2b.backward(coord)
eq_(math.fabs(merc_coord1.x - merc_coord1.x) < 1,True)
eq_(math.fabs(merc_coord1.x - merc_coord2.x) < 1,True)
eq_(math.fabs(merc_coord1.x - merc_coord3.x) < 1,True)
eq_(math.fabs(merc_coord1.x - merc_coord4.x) < 1,True)
eq_(math.fabs(merc_coord1.y - merc_coord1.y) < 1,True)
eq_(math.fabs(merc_coord1.y - merc_coord2.y) < 1,True)
eq_(math.fabs(merc_coord1.y - merc_coord3.y) < 1,True)
eq_(math.fabs(merc_coord1.y - merc_coord4.y) < 1,True)
lon_lat_coord1 = tr1.backward(merc_coord1) lon_lat_coord1 = tr1.backward(merc_coord1)
lon_lat_coord2 = tr2.backward(merc_coord2) lon_lat_coord2 = tr1b.forward(merc_coord2)
assert_almost_equal(coord.x,lon_lat_coord1.x) lon_lat_coord3 = tr2.backward(merc_coord3)
assert_almost_equal(coord.y,lon_lat_coord1.y) lon_lat_coord4 = tr2b.forward(merc_coord4)
assert_almost_equal(coord.x,lon_lat_coord2.x) eq_(math.fabs(coord.x - lon_lat_coord1.x) < 1,True)
assert_almost_equal(coord.y,lon_lat_coord2.y) eq_(math.fabs(coord.x - lon_lat_coord2.x) < 1,True)
eq_(math.fabs(coord.x - lon_lat_coord3.x) < 1,True)
eq_(math.fabs(coord.x - lon_lat_coord4.x) < 1,True)
eq_(math.fabs(coord.y - lon_lat_coord1.y) < 1,True)
eq_(math.fabs(coord.y - lon_lat_coord2.y) < 1,True)
eq_(math.fabs(coord.y - lon_lat_coord3.y) < 1,True)
eq_(math.fabs(coord.y - lon_lat_coord4.y) < 1,True)
if __name__ == "__main__": if __name__ == "__main__":
[eval(run)() for run in dir() if 'test_' in run] run_all(eval(x) for x in dir() if x.startswith("test_"))