ProjectionsΒΆ
code chooser
The proj module provides support for representing spatial reference
systems and transforming data between spatial projections.
The Projection class is used to represent a spatial reference
system. A common way to identify a spatial reference system is with an EPSG code.
GeoScript projections can be created with an EPSG code.
>>> from geoscript.proj import Projection
>>> prj = Projection('epsg:4326')
>>> prj
GEOGCS["WGS 84",
DATUM["World Geodetic System 1984",
SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
UNIT["degree", 0.017453292519943295],
AXIS["Geodetic longitude", EAST],
AXIS["Geodetic latitude", NORTH],
AUTHORITY["EPSG","4326"]]
>> var Projection = require("geoscript/proj").Projection;
>> var prj = new Projection("epsg:4326");
>> prj
<Projection EPSG:4326>
scala> import org.geoscript.projection.Projection
scala> val prj = new Projection("epsg:4326")
proj: org.geoscript.projection.Projection = EPSG:4326
groovy:000> import geoscript.proj.Projection
===> [import geoscript.proj.Projection]
groovy:000> prj = new Projection("epsg:4326")
===> EPSG:4326
The EPSG code can also be accessed for any projection.
>>> from geoscript.proj import Projection
>>> prj = Projection('epsg:4326')
>>> prj.id
'EPSG:4326'
>> var Projection = require("geoscript/proj").Projection;
>> var prj = new Projection("epsg:4326");
>> prj.id
EPSG:4326
scala> import org.geoscript.projection._
scala> val proj = Projection("epsg:4326")
proj: org.geoscript.projection.Projection = EPSG:4326
scala> proj.id
res0: String = EPSG:4326
groovy:000> import geoscript.proj.Projection
===> [import geoscript.proj.Projection]
groovy:000> proj = new Projection("epsg:4326")
===> EPSG:4326
groovy:000> proj.id
===> EPSG:4326
Often a spatial reference system is defined by well-known text
rather than an EPSG code. A Projection object can be created directly
from well known text.
>>> from geoscript.proj import Projection
>>> wkt = 'GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]'
>>> prj = Projection(wkt)
>>> prj
GEOGCS["GCS_WGS_1984",
DATUM["D_WGS_1984",
SPHEROID["WGS_1984", 6378137.0, 298.257223563]],
PRIMEM["Greenwich", 0.0],
UNIT["degree", 0.017453292519943295],
AXIS["Longitude", EAST],
AXIS["Latitude", NORTH]]
>> var Projection = require("geoscript/proj").Projection;
>> var wkt = 'GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]';
>> var prj = new Projection(wkt);
>> prj
<Projection EPSG:4326>
scala> import org.geoscript.projection._
scala> val proj = Projection("""GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]""")
proj: org.geoscript.projection.Projection = GCS_WGS_1984
groovy:000> import geoscript.proj.Projection
===> [import geoscript.proj.Projection]
groovy:000> proj = new Projection("""GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137,298.257223563]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]""")
===> EPSG:4326
The well-known text can also be accessed for any projection, regardless of how it was constructed.
>>> from geoscript.proj import Projection
>>> prj = Projection('epsg:26912')
>>> prj.wkt
'PROJCS["NAD83 / UTM zone 12N", \n GEOGCS["NAD83", \n DATUM["North American Datum 1983", \n SPHEROID["GRS 1980", 6378137.0, 298.257222101, AUTHORITY["EPSG","7019"]], \n TOWGS84[1.0, 1.0, -1.0, 0.0, 0.0, 0.0, 0.0], \n AUTHORITY["EPSG","6269"]], \n PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], \n UNIT["degree", 0.017453292519943295], \n AXIS["Geodetic longitude", EAST], \n AXIS["Geodetic latitude", NORTH], \n AUTHORITY["EPSG","4269"]], \n PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]], \n PARAMETER["central_meridian", -111.0], \n PARAMETER["latitude_of_origin", 0.0], \n PARAMETER["scale_factor", 0.9996], \n PARAMETER["false_easting", 500000.0], \n PARAMETER["false_northing", 0.0], \n UNIT["m", 1.0], \n AXIS["Easting", EAST], \n AXIS["Northing", NORTH], \n AUTHORITY["EPSG","26912"]]'
>> var Projection = require("geoscript/proj").Projection;
>> var prj = new Projection("epsg:26912");
>> prj.wkt
PROJCS["NAD83 / UTM zone 12N",
GEOGCS["NAD83",
DATUM["North American Datum 1983",
SPHEROID["GRS 1980", 6378137.0, 298.257222101, AUTHORITY["EPSG","7019"]],
TOWGS84[1.0, 1.0, -1.0, 0.0, 0.0, 0.0, 0.0],
AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
UNIT["degree", 0.017453292519943295],
AXIS["Geodetic longitude", EAST],
AXIS["Geodetic latitude", NORTH],
AUTHORITY["EPSG","4269"]],
PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]],
PARAMETER["central_meridian", -111.0],
PARAMETER["latitude_of_origin", 0.0],
PARAMETER["scale_factor", 0.9996],
PARAMETER["false_easting", 500000.0],
PARAMETER["false_northing", 0.0],
UNIT["m", 1.0],
AXIS["Easting", EAST],
AXIS["Northing", NORTH],
AUTHORITY["EPSG","26912"]]
scala> import org.geoscript.projection._
scala> val proj = Projection("epsg:26912")
proj: org.geoscript.projection.Projection = EPSG:26912
scala> proj.wkt
res0: String =
PROJCS["NAD83 / UTM zone 12N",
GEOGCS["NAD83",
DATUM["North American Datum 1983",
SPHEROID["GRS 1980", 6378137.0, 298.257222101, AUTHORITY["EPSG","7019"]],
TOWGS84[1.0, 1.0, -1.0, 0.0, 0.0, 0.0, 0.0],
AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
UNIT["degree", 0.017453292519943295],
AXIS["Geo...
groovy:000> import geoscript.proj.Projection
===> [import geoscript.proj.Projection]
groovy:000> proj = new Projection("epsg:26912")
===> EPSG:26912
groovy:000> proj.wkt
===> PROJCS["NAD83 / UTM zone 12N",
GEOGCS["NAD83",
DATUM["North American Datum 1983",
SPHEROID["GRS 1980", 6378137.0, 298.257222101, AUTHORITY["EPSG","7019"]],
TOWGS84[1.0, 1.0, -1.0, 0.0, 0.0, 0.0, 0.0],
AUTHORITY["EPSG","6269"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
UNIT["degree", 0.017453292519943295],
AXIS["Geodetic latitude", NORTH],
AXIS["Geodetic longitude", EAST],
AUTHORITY["EPSG","4269"]],
PROJECTION["Transverse Mercator", AUTHORITY["EPSG","9807"]],
PARAMETER["central_meridian", -111.0],
PARAMETER["latitude_of_origin", 0.0],
PARAMETER["scale_factor", 0.9996],
PARAMETER["false_easting", 500000.0],
PARAMETER["false_northing", 0.0],
UNIT["m", 1.0],
AXIS["Easting", EAST],
AXIS["Northing", NORTH],
AUTHORITY["EPSG","26912"]]
Coordinates in one spatial reference system can be transformed into coordinates in a different spatial reference system.
>>> from geoscript.proj import Projection
>>> src = Projection('epsg:4326')
>>> src.transform((-111, 45.7), 'epsg:26912')
(499999.42501775385, 5060716.092032814)
>>> from geoscript.geom import Point
>>> src.transform(Point(-111, 45.7), 'epsg:26912')
POINT (499999.42501775385 5060716.092032814)
>> var Projection = require("geoscript/proj").Projection;
>> var src = new Projection("epsg:4326");
>> var GEOM = require("geoscript/geom");
>> var point = new GEOM.Point([-111, 45.7]);
>> point.projection = src;
<Projection EPSG:4326>
>> point.transform("epsg:26912")
<Point [500000, 5060716.31816507]>
scala> import org.geoscript.projection._
scala> import org.geoscript.geometry._
scala> val src = Projection("epsg:4326")
src: org.geoscript.projection.Projection = EPSG:4326
scala> val point = Point(-111, 45.7) in src
point: org.geoscript.geometry.Point = POINT (-111 45.7)
scala> point in Projection("epsg:26912")
res0: org.geoscript.geometry.Point = POINT (-412539.9935856778 -7830632.25817681)
groovy:000> import geoscript.proj.Projection
===> [import geoscript.proj.Projection]
groovy:000> import geoscript.geom.*
===> [import geoscript.proj.Projection, import geoscript.geom.*]
groovy:000> src = new Projection("epsg:4326")
===> EPSG:4326
groovy:000> point = new Point(-111, 45.7)
===> POINT (-111 45.7)
groovy:000> src.transform(point,"epsg:26912")
===> POINT (-412539.9935856778 -7830632.25817681)
See also
See the proj module in the Python API reference
for the details on the properties and method available on projection objects.
See also
See the proj module in the JavaScript API reference
for the details on the properties and method available on projection objects.
