Secure your code as it's written. Use Snyk Code to scan source code in minutes - no build needed - and fix issues immediately.
def test_2d_with_time_itransform_original_crs_obs2():
assert_almost_equal(
list(
itransform(4896, 7930, [(3496737.2679, 743254.4507, 2019.0)], time_3rd=True)
),
[(3496737.4105305015, 743254.1014318303, 2019.0)],
)
def test_itransform_no_error():
with pytest.warns(FutureWarning):
pj = Proj(init="epsg:4555")
pjx, pjy = pj(116.366, 39.867)
list(itransform(pj, Proj(4326), [(pjx, pjy)], radians=True, errcheck=True))
def test_itransform_radians():
with pytest.warns(FutureWarning):
WGS84 = pyproj.Proj("+init=EPSG:4326")
ECEF = pyproj.Proj(proj="geocent", ellps="WGS84", datum="WGS84")
assert_almost_equal(
list(
pyproj.itransform(
ECEF, WGS84, [(-2704026.010, -4253051.810, 3895878.820)], radians=True
)
),
[(-2.137113493845668, 0.6613203738996222, -20.531156923621893)],
)
assert_almost_equal(
list(
pyproj.itransform(
WGS84,
ECEF,
[(-2.137113493845668, 0.6613203738996222, -20.531156923621893)],
radians=True,
)
),
[(-2704026.010, -4253051.810, 3895878.820)],
def test_itransform_radians():
with pytest.warns(DeprecationWarning):
WGS84 = pyproj.Proj("+init=EPSG:4326")
ECEF = pyproj.Proj(proj="geocent", ellps="WGS84", datum="WGS84")
assert_almost_equal(
list(
pyproj.itransform(
ECEF, WGS84, [(-2704026.010, -4253051.810, 3895878.820)], radians=True
)
),
[(-2.137113493845668, 0.6613203738996222, -20.531156923621893)],
)
assert_almost_equal(
list(
pyproj.itransform(
WGS84,
ECEF,
[(-2.137113493845668, 0.6613203738996222, -20.531156923621893)],
radians=True,
)
),
[(-2704026.010, -4253051.810, 3895878.820)],
)
def test_always_xy__itransform():
assert_almost_equal(
list(itransform(2193, 4326, [(1625350, 5504853)], always_xy=True)),
[(173.29964730317386, -40.60674802693758)],
)
def test_always_xy__itransform():
assert_almost_equal(
list(itransform(2193, 4326, [(1625350, 5504853)], always_xy=True)),
[(173.29964730317386, -40.60674802693758)],
)
'srs': nativeSrs,
}
# Make sure geographic coordinates do not exceed their limits
if self._proj4Proj(nativeSrs).crs.is_geographic and srs:
try:
self._proj4Proj(srs)(0, 90, errcheck=True)
yBound = 90.0
except RuntimeError:
yBound = 89.999999
for key in ('ll', 'ul', 'lr', 'ur'):
bounds[key]['y'] = max(min(bounds[key]['y'], yBound), -yBound)
if srs and srs != nativeSrs:
inProj = self._proj4Proj(nativeSrs)
outProj = self._proj4Proj(srs)
keys = ('ll', 'ul', 'lr', 'ur')
pts = pyproj.itransform(inProj, outProj, [
(bounds[key]['x'], bounds[key]['y']) for key in keys], always_xy=True)
for idx, pt in enumerate(pts):
key = keys[idx]
bounds[key]['x'] = pt[0]
bounds[key]['y'] = pt[1]
bounds['srs'] = srs.decode('utf8') if isinstance(srs, six.binary_type) else srs
bounds['xmin'] = min(bounds['ll']['x'], bounds['ul']['x'],
bounds['lr']['x'], bounds['ur']['x'])
bounds['xmax'] = max(bounds['ll']['x'], bounds['ul']['x'],
bounds['lr']['x'], bounds['ur']['x'])
bounds['ymin'] = min(bounds['ll']['y'], bounds['ul']['y'],
bounds['lr']['y'], bounds['ur']['y'])
bounds['ymax'] = max(bounds['ll']['y'], bounds['ul']['y'],
bounds['lr']['y'], bounds['ur']['y'])
self._bounds[srs] = bounds
return self._bounds[srs]