PROJ

High-level Rust bindings for the latest stable version of PROJ (7.0.x), compatible with the Georust ecosystem.

Requirements

By default, this crate depends on a pre-built library, so PROJ v7.0.x must be present on your system. While this crate may be backwards-compatible with older PROJ 6 versions, this is neither tested nor supported.

Two features are available:

proj = { version = "0.16.1", features = ["pkg_config"] }
proj = = { version = "0.16.1", features = ["bundled_proj"] }

The pkg_config feature enables the use of pkg-config when linking against libproj – note that pkg-config must be available on your system.

The bundled_proj feature allows you to link against a PROJ included with (and built from source by) the proj-sys crate, upon which this crate is built. To do so, enable the bundled_proj Cargo feature. This feature requires sqlite3 to be present on your system.

Examples

Convert from NAD 83 US Survey Feet to NAD 83 Meters Using EPSG Codes

```rust use proj::Proj;

extern crate geotypes; use geotypes::Point;

let from = "EPSG:2230"; let to = "EPSG:26946"; let fttom = Proj::newknowncrs(&from, &to, None).unwrap(); let result = fttom .convert(Point::new(4760096.421921, 3744293.729449)) .unwrap(); assertalmosteq(result.x(), 1450880.29); assertalmosteq(result.y(), 1141263.01); ```

Note that as of v5.0.0, PROJ uses the pipeline operator, which allows an arbitrary number of steps in a conversion. The example below works as follows:

• define the operation as a pipeline operation
• define step 1 as an inverse transform, yielding geodetic coordinates
• define step 2 as a forward transform to projected coordinates, yielding metres.

Convert from NAD 83 US Survey Feet to NAD 83 Meters Using the pipeline Operator

```rust use proj::Proj;

extern crate geotypes; use geotypes::Point;

let fttom = Proj::new(" +proj=pipeline +step +inv +proj=lcc +lat1=33.88333333333333 +lat2=32.78333333333333 +lat0=32.16666666666666 +lon0=-116.25 +x0=2000000.0001016 +y0=500000.0001016001 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +nodefs +step +proj=lcc +lat1=33.88333333333333 +lat2=32.78333333333333 +lat0=32.16666666666666 +lon0=-116.25 +x0=2000000 +y0=500000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +nodefs ").unwrap(); // The Presidio, approximately let result = fttom.convert(Point::new(4760096.421921, 3744293.729449)).unwrap(); asserteq!(result.x(), 1450880.29); asserteq!(result.y(), 1141263.01); ```

Inverse Projection from Stereo70 to Geodetic

```rust use proj::Proj;

extern crate geotypes; use geotypes::Point;

// Carry out an inverse projection from Pulkovo 1942(58) / Stereo70 (EPSG 3844) // into geodetic lon and lat coordinates (in radians) let stereo70 = Proj::new(" +proj=sterea +lat0=46 +lon0=25 +k=0.99975 +x0=500000 +y0=500000 +ellps=krass +towgs84=33.4,-146.6,-76.3,-0.359,-0.053,0.844,-0.84 +units=m +nodefs ").unwrap(); let rp = stereo70.project( Point::new(500119.70352012233, 500027.77896348457), true ).unwrap(); asserteq!(rp, Point::new(0.436332, 0.802851)); ```

Bulk Transformations

The Proj::convert_array() method is available for bulk conversions. It accepts a mutable slice (or anything that can Deref to a mutable slice) of Point<T: Float> or Into<Point<T: Float> elements.

License

Licensed under either of

• Apache License, Version 2.0 (LICENSE-APACHE or http://www.apache.org/licenses/LICENSE-2.0)
• MIT license (LICENSE-MIT or http://opensource.org/licenses/MIT)

at your option.