PROJ

High-level Rust bindings for the latest stable version of PROJ (7.1.x), compatible with the Georust ecosystem. Includes network grid download functionality.

Requirements

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

Two features are available:

proj = { version = "0.19.0", features = ["pkg_config"] }
proj = = { version = "0.19.0", 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 use a statically-linked libproj included with (and built from source by) the proj-sys crate. Note that this feature requires Sqlite3 and libtiff 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(); assertapproxeq!(result.x() as f64, 1450880.2910605003); assertapproxeq!(result.y() as f64, 1141263.0111604529); ```

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 ", false).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() and Proj::project_array()methods are available for bulk conversions. Both accept mutable slices (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.