23 | * This is the base class for all exceptions
24 | * thrown in the Proj4J API.
25 | *
26 | * @author mbdavis
27 | *
28 | */
29 | public class Proj4jException extends RuntimeException
30 | {
31 | public Proj4jException() {
32 | }
33 |
34 | public Proj4jException(String message) {
35 | super(message);
36 | }
37 | }
38 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/ProjectionException.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j;
18 |
19 | import org.locationtech.proj4j.proj.Projection;
20 |
21 | /**
22 | * Signals that an erroneous situation has
23 | * occured during the computation of
24 | * a projected coordinate system value.
25 | *
26 | * @author mbdavis
27 | */
28 | public class ProjectionException extends Proj4jException {
29 |
30 | public static String ERR_17 = "non-convergent inverse meridinal dist";
31 |
32 | public ProjectionException() {}
33 |
34 | public ProjectionException(String message) {
35 | super(message);
36 | }
37 |
38 | public ProjectionException(Projection proj, String message) {
39 | this(proj.toString() + ": " + message);
40 | }
41 | }
42 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/UnknownAuthorityCodeException.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j;
18 |
19 | /**
20 | * Signals that an authority code is unknown
21 | * and cannot be mapped to a CRS definition.
22 | *
23 | * @author mbdavis
24 | */
25 | public class UnknownAuthorityCodeException extends Proj4jException {
26 | public UnknownAuthorityCodeException(String message) {
27 | super(message);
28 | }
29 | }
30 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/UnsupportedParameterException.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j;
18 |
19 | /**
20 | * Signals that a parameter in a CRS specification
21 | * is not currently supported, or unknown.
22 | *
23 | * @author mbdavis
24 | *
25 | */
26 | public class UnsupportedParameterException extends Proj4jException
27 | {
28 |
29 | public UnsupportedParameterException(String message) {
30 | super(message);
31 | }
32 | }
33 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/geodesic/Constants.java:
--------------------------------------------------------------------------------
1 | /**
2 | * Implementation of the net.sf.geographiclib.Constants class
3 | *
4 | * Copyright (c) Charles Karney (2013)
13 | * Define constants specifying the WGS84 ellipsoid.
14 | ***********************************************************************/
15 | public class Constants {
16 | /**
17 | * The equatorial radius of WGS84 ellipsoid (6378137 m).
18 | **********************************************************************/
19 | public static final double WGS84_a = 6378137;
20 | /**
21 | * The flattening of WGS84 ellipsoid (1/298.257223563).
22 | **********************************************************************/
23 | public static final double WGS84_f = 1/298.257223563;
24 |
25 | private Constants() {}
26 | }
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/geodesic/GeographicErr.java:
--------------------------------------------------------------------------------
1 | /**
2 | * Implementation of the net.sf.geographiclib.GeographicErr class
3 | *
4 | * Copyright (c) Charles Karney (2013)
13 | * A class to handle exceptions. It's derived from RuntimeException so it
14 | * can be caught by the usual catch clauses.
15 | **********************************************************************/
16 | public class GeographicErr extends RuntimeException {
17 | /**
18 | * Constructor
19 | *
20 | * @param msg a string message, which is accessible in the catch
21 | * clause via getMessage().
22 | **********************************************************************/
23 | public GeographicErr(String msg) { super(msg); }
24 | }
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/geodesic/Pair.java:
--------------------------------------------------------------------------------
1 | /**
2 | * Implementation of the net.sf.geographiclib.Pair class
3 | *
4 | * Copyright (c) Charles Karney (2013-2020)
13 | * This duplicates the C++ class {@code std::pair
27 | * @param first the first member of the pair.
28 | * @param second the second member of the pair.
29 | **********************************************************************/
30 | public Pair(double first, double second)
31 | { this.first = first; this.second = second; }
32 | /**
33 | * No-argument Constructor
34 | **********************************************************************/
35 | public Pair() {}
36 | }
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/AitoffProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | *******************************************************************************/
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import java.util.Objects;
23 |
24 | import org.locationtech.proj4j.ProjCoordinate;
25 |
26 |
27 | public class AitoffProjection extends PseudoCylindricalProjection {
28 |
29 | protected final static int AITOFF = 0;
30 | protected final static int WINKEL = 1;
31 |
32 | private boolean winkel = false;
33 | private double cosphi1 = 0;
34 |
35 | public AitoffProjection() {
36 | }
37 |
38 | public AitoffProjection(int type, double projectionLatitude) {
39 | this.projectionLatitude = projectionLatitude;
40 | winkel = type == WINKEL;
41 | }
42 |
43 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
44 | double c = 0.5 * lplam;
45 | double d = Math.acos(Math.cos(lpphi) * Math.cos(c));
46 |
47 | if (d != 0) {
48 | out.x = 2. * d * Math.cos(lpphi) * Math.sin(c) * (out.y = 1. / Math.sin(d));
49 | out.y *= d * Math.sin(lpphi);
50 | } else
51 | out.x = out.y = 0.0;
52 | if (winkel) {
53 | out.x = (out.x + lplam * cosphi1) * 0.5;
54 | out.y = (out.y + lpphi) * 0.5;
55 | }
56 | return out;
57 | }
58 |
59 | public void initialize() {
60 | super.initialize();
61 | if (winkel) {
62 | //FIXME
63 | // if (pj_param(P->params, "tlat_1").i)
64 | // if ((cosphi1 = Math.cos(pj_param(P->params, "rlat_1").f)) == 0.)
65 | // throw new IllegalArgumentException("-22")
66 | // else /* 50d28' or acos(2/pi) */
67 | cosphi1 = 0.636619772367581343;
68 | }
69 | }
70 |
71 | public boolean hasInverse() {
72 | return false;
73 | }
74 |
75 | public String toString() {
76 | return winkel ? "Winkel Tripel" : "Aitoff";
77 | }
78 |
79 | @Override
80 | public boolean equals(Object that) {
81 | if (this == that) {
82 | return true;
83 | }
84 | if (that instanceof AitoffProjection) {
85 | AitoffProjection p = (AitoffProjection) that;
86 | return (winkel == p.winkel && super.equals(that));
87 | }
88 | return false;
89 | }
90 |
91 | @Override
92 | public int hashCode() {
93 | return Objects.hash(winkel, super.hashCode());
94 | }
95 | }
96 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/AugustProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class AugustProjection extends Projection {
25 |
26 | private final static double M = 1.333333333333333;
27 |
28 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
29 | double t, c1, c, x1, x12, y1, y12;
30 |
31 | t = Math.tan(.5 * lpphi);
32 | c1 = Math.sqrt(1. - t * t);
33 | c = 1. + c1 * Math.cos(lplam *= .5);
34 | x1 = Math.sin(lplam) * c1 / c;
35 | y1 = t / c;
36 | out.x = M * x1 * (3. + (x12 = x1 * x1) - 3. * (y12 = y1 * y1));
37 | out.y = M * y1 * (3. + 3. * x12 - y12);
38 | return out;
39 | }
40 |
41 | /**
42 | * Returns true if this projection is conformal
43 | */
44 | public boolean isConformal() {
45 | return true;
46 | }
47 |
48 | public boolean hasInverse() {
49 | return false;
50 | }
51 |
52 | public String toString() {
53 | return "August Epicycloidal";
54 | }
55 |
56 | }
57 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/AzimuthalProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j.proj;
18 |
19 | import java.util.Objects;
20 |
21 | import org.locationtech.proj4j.util.ProjectionMath;
22 |
23 | /**
24 | * The superclass for all azimuthal map projections
25 | */
26 | public abstract class AzimuthalProjection extends Projection {
27 |
28 | public final static int NORTH_POLE = 1;
29 | public final static int SOUTH_POLE = 2;
30 | public final static int EQUATOR = 3;
31 | public final static int OBLIQUE = 4;
32 |
33 | protected int mode;
34 | protected double sinphi0, cosphi0;
35 | private double mapRadius = 90.0;
36 |
37 | public AzimuthalProjection() {
38 | this( Math.toRadians(45.0), Math.toRadians(45.0) );
39 | }
40 |
41 | public AzimuthalProjection(double projectionLatitude, double projectionLongitude) {
42 | this.projectionLatitude = projectionLatitude;
43 | this.projectionLongitude = projectionLongitude;
44 | initialize();
45 | }
46 |
47 | public void initialize() {
48 | super.initialize();
49 | if (Math.abs(Math.abs(projectionLatitude) - ProjectionMath.HALFPI) < EPS10)
50 | mode = projectionLatitude < 0. ? SOUTH_POLE : NORTH_POLE;
51 | else if (Math.abs(projectionLatitude) > EPS10) {
52 | mode = OBLIQUE;
53 | sinphi0 = Math.sin(projectionLatitude);
54 | cosphi0 = Math.cos(projectionLatitude);
55 | } else
56 | mode = EQUATOR;
57 | }
58 |
59 | public boolean inside(double lon, double lat) {
60 | return ProjectionMath.greatCircleDistance( Math.toRadians(lon), Math.toRadians(lat), projectionLongitude, projectionLatitude) < Math.toRadians(mapRadius);
61 | }
62 |
63 | /**
64 | * Set the map radius (in degrees). 180 shows a hemisphere, 360 shows the whole globe.
65 | */
66 | public void setMapRadius(double mapRadius) {
67 | this.mapRadius = mapRadius;
68 | }
69 |
70 | public double getMapRadius() {
71 | return mapRadius;
72 | }
73 |
74 | @Override
75 | public boolean equals(Object that) {
76 | if (this == that) {
77 | return true;
78 | }
79 | if (that instanceof AzimuthalProjection) {
80 | AzimuthalProjection p = (AzimuthalProjection) that;
81 | return (
82 | mode == p.mode &&
83 | sinphi0 == p.sinphi0 &&
84 | cosphi0 == p.cosphi0 &&
85 | mapRadius == p.mapRadius &&
86 | super.equals(that));
87 | }
88 | return false;
89 | }
90 |
91 | @Override
92 | public int hashCode() {
93 | return Objects.hash(mode, sinphi0, cosphi0, mapRadius, super.hashCode());
94 | }
95 | }
96 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/BoggsProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class BoggsProjection extends PseudoCylindricalProjection {
26 |
27 | private final static int NITER = 20;
28 | private final static double EPS = 1e-7;
29 | private final static double ONETOL = 1.000001;
30 | private final static double FXC = 2.00276;
31 | private final static double FXC2 = 1.11072;
32 | private final static double FYC = 0.49931;
33 | private final static double FYC2 = 1.41421356237309504880;
34 |
35 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
36 | double theta, th1, c;
37 | int i;
38 |
39 | theta = lpphi;
40 | if (Math.abs(Math.abs(lpphi) - ProjectionMath.HALFPI) < EPS)
41 | out.x = 0.;
42 | else {
43 | c = Math.sin(theta) * Math.PI;
44 | for (i = NITER; i > 0; --i) {
45 | theta -= th1 = (theta + Math.sin(theta) - c) /
46 | (1. + Math.cos(theta));
47 | if (Math.abs(th1) < EPS) break;
48 | }
49 | theta *= 0.5;
50 | out.x = FXC * lplam / (1. / Math.cos(lpphi) + FXC2 / Math.cos(theta));
51 | }
52 | out.y = FYC * (lpphi + FYC2 * Math.sin(theta));
53 | return out;
54 | }
55 |
56 | /**
57 | * Returns true if this projection is equal area
58 | */
59 | public boolean isEqualArea() {
60 | return true;
61 | }
62 |
63 | public boolean hasInverse() {
64 | return false;
65 | }
66 |
67 | public String toString() {
68 | return "Boggs Eumorphic";
69 | }
70 |
71 | }
72 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/CentralCylindricalProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.ProjectionException;
24 | import org.locationtech.proj4j.util.ProjectionMath;
25 |
26 | public class CentralCylindricalProjection extends CylindricalProjection {
27 |
28 | private final static double EPS10 = 1.e-10;
29 |
30 | public CentralCylindricalProjection() {
31 | minLatitude = Math.toRadians(-80);
32 | maxLatitude = Math.toRadians(80);
33 | }
34 |
35 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
36 | if (Math.abs(Math.abs(lpphi) - ProjectionMath.HALFPI) <= EPS10) throw new ProjectionException("F");
37 | out.x = lplam;
38 | out.y = Math.tan(lpphi);
39 | return out;
40 | }
41 |
42 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
43 | out.y = Math.atan(xyy);
44 | out.x = xyx;
45 | return out;
46 | }
47 |
48 | public boolean hasInverse() {
49 | return true;
50 | }
51 |
52 | public String toString() {
53 | return "Central Cylindrical";
54 | }
55 |
56 | }
57 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/CollignonProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.ProjectionException;
24 | import org.locationtech.proj4j.util.ProjectionMath;
25 |
26 | public class CollignonProjection extends Projection {
27 |
28 | private final static double FXC = 1.12837916709551257390;
29 | private final static double FYC = 1.77245385090551602729;
30 | private final static double ONEEPS = 1.0000001;
31 |
32 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
33 | if ((out.y = 1. - Math.sin(lpphi)) <= 0.)
34 | out.y = 0.;
35 | else
36 | out.y = Math.sqrt(out.y);
37 | out.x = FXC * lplam * out.y;
38 | out.y = FYC * (1. - out.y);
39 | return out;
40 | }
41 |
42 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
43 | double lpphi = xyy / FYC - 1.;
44 | if (Math.abs(out.y = 1. - lpphi * lpphi) < 1.)
45 | out.y = Math.asin(lpphi);
46 | else if (Math.abs(lpphi) > ONEEPS) throw new ProjectionException("I");
47 | else out.y = lpphi < 0. ? -ProjectionMath.HALFPI : ProjectionMath.HALFPI;
48 | if ((out.x = 1. - Math.sin(lpphi)) <= 0.)
49 | out.x = 0.;
50 | else
51 | out.x = xyx / (FXC * Math.sqrt(out.x));
52 | out.y = lpphi;
53 | return out;
54 | }
55 |
56 | /**
57 | * Returns true if this projection is equal area
58 | */
59 | public boolean isEqualArea() {
60 | return true;
61 | }
62 |
63 | public boolean hasInverse() {
64 | return true;
65 | }
66 |
67 | public String toString() {
68 | return "Collignon";
69 | }
70 |
71 | }
72 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/ConicProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j.proj;
18 |
19 |
20 |
21 | /**
22 | * A projection formed by projecting the sphere
23 | * onto a cone tangent, or secant, to the sphere
24 | * along any small circle (usually a mid-latitude parallel).
25 | * In the normal aspect (which is oblique for conic projections),
26 | * parallels are projected as concentric arcs of circles,
27 | * and meridians are projected as straight lines
28 | * radiating at uniform angular intervals from the apex of the flattened cone.
29 | */
30 | public abstract class ConicProjection extends Projection {
31 |
32 | public String toString() {
33 | return "Conic";
34 | }
35 |
36 | }
37 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/CrasterProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class CrasterProjection extends Projection {
25 |
26 | private final static double XM = 0.97720502380583984317;
27 | private final static double RXM = 1.02332670794648848847;
28 | private final static double YM = 3.06998012383946546542;
29 | private final static double RYM = 0.32573500793527994772;
30 | private final static double THIRD = 0.333333333333333333;
31 |
32 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
33 | lpphi *= THIRD;
34 | out.x = XM * lplam * (2. * Math.cos(lpphi + lpphi) - 1.);
35 | out.y = YM * Math.sin(lpphi);
36 | return out;
37 | }
38 |
39 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
40 | out.y = 3. * Math.asin(xyy * RYM);
41 | out.x = xyx * RXM / (2. * Math.cos((out.y + out.y) * THIRD) - 1);
42 | return out;
43 | }
44 |
45 | /**
46 | * Returns true if this projection is equal area
47 | */
48 | public boolean isEqualArea() {
49 | return true;
50 | }
51 |
52 | public boolean hasInverse() {
53 | return true;
54 | }
55 |
56 | public String toString() {
57 | return "Craster Parabolic (Putnins P4)";
58 | }
59 |
60 | }
61 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/CylindricalEqualAreaProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.ProjectionException;
24 | import org.locationtech.proj4j.util.ProjectionMath;
25 |
26 | public class CylindricalEqualAreaProjection extends Projection {
27 |
28 | private double qp;
29 | private double[] apa;
30 |
31 | public CylindricalEqualAreaProjection() {
32 | this(0.0, 0.0, 0.0);
33 | }
34 |
35 | public CylindricalEqualAreaProjection(double projectionLatitude, double projectionLongitude, double trueScaleLatitude) {
36 | this.projectionLatitude = projectionLatitude;
37 | this.projectionLongitude = projectionLongitude;
38 | this.trueScaleLatitude = trueScaleLatitude;
39 | initialize();
40 | }
41 |
42 | public void initialize() {
43 | super.initialize();
44 | double t = trueScaleLatitude;
45 |
46 | scaleFactor = Math.cos(t);
47 | if (es != 0) {
48 | t = Math.sin(t);
49 | scaleFactor /= Math.sqrt(1. - es * t * t);
50 | apa = ProjectionMath.authset(es);
51 | qp = ProjectionMath.qsfn(1., e, one_es);
52 | }
53 | }
54 |
55 | public ProjCoordinate project(double lam, double phi, ProjCoordinate xy) {
56 | if (spherical) {
57 | xy.x = scaleFactor * lam;
58 | xy.y = Math.sin(phi) / scaleFactor;
59 | } else {
60 | xy.x = scaleFactor * lam;
61 | xy.y = .5 * ProjectionMath.qsfn(Math.sin(phi), e, one_es) / scaleFactor;
62 | }
63 | return xy;
64 | }
65 |
66 | public ProjCoordinate projectInverse(double x, double y, ProjCoordinate lp) {
67 | if (spherical) {
68 | double t;
69 |
70 | if ((t = Math.abs(y *= scaleFactor)) - EPS10 <= 1.) {
71 | if (t >= 1.)
72 | lp.y = y < 0. ? -ProjectionMath.HALFPI : ProjectionMath.HALFPI;
73 | else
74 | lp.y = Math.asin(y);
75 | lp.x = x / scaleFactor;
76 | } else throw new ProjectionException();
77 | } else {
78 | lp.y = ProjectionMath.authlat(Math.asin( 2. * y * scaleFactor / qp), apa);
79 | lp.x = x / scaleFactor;
80 | }
81 | return lp;
82 | }
83 |
84 | public boolean hasInverse() {
85 | return true;
86 | }
87 |
88 | public boolean isRectilinear() {
89 | return true;
90 | }
91 |
92 | }
93 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/CylindricalProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j.proj;
18 |
19 |
20 |
21 | /**
22 | * A projection
23 | * in which meridians are mapped to equally spaced vertical lines
24 | * and circles of latitude (parallels) are mapped to horizontal lines.
25 | */
26 | public abstract class CylindricalProjection extends Projection {
27 |
28 | public boolean isRectilinear() {
29 | return true;
30 | }
31 |
32 | public String toString() {
33 | return "Cylindrical";
34 | }
35 |
36 |
37 | }
38 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/DenoyerProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 |
25 | public class DenoyerProjection extends Projection {
26 |
27 | public final static double C0 = 0.95;
28 | public final static double C1 = -.08333333333333333333;
29 | public final static double C3 = 0.00166666666666666666;
30 | public final static double D1 = 0.9;
31 | public final static double D5 = 0.03;
32 |
33 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
34 | out.y = lpphi;
35 | out.x = lplam;
36 | double aphi = Math.abs(lplam);
37 | out.x *= Math.cos((C0 + aphi * (C1 + aphi * aphi * C3)) *
38 | (lpphi * (D1 + D5 * lpphi * lpphi * lpphi * lpphi)));
39 | return out;
40 | }
41 |
42 | public boolean parallelsAreParallel() {
43 | return true;
44 | }
45 |
46 | public boolean hasInverse() {
47 | return false;
48 | }
49 |
50 | public String toString() {
51 | return "Denoyer Semi-elliptical";
52 | }
53 |
54 | }
55 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Eckert1Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class Eckert1Projection extends Projection {
25 |
26 | private final static double FC = .92131773192356127802;
27 | private final static double RP = .31830988618379067154;
28 |
29 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
30 | out.x = FC * lplam * (1. - RP * Math.abs(lpphi));
31 | out.y = FC * lpphi;
32 | return out;
33 | }
34 |
35 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
36 | out.y = xyy / FC;
37 | out.x = xyx / (FC * (1. - RP * Math.abs(out.y)));
38 | return out;
39 | }
40 |
41 | public boolean hasInverse() {
42 | return true;
43 | }
44 |
45 | public String toString() {
46 | return "Eckert I";
47 | }
48 |
49 | }
50 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Eckert2Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.ProjectionException;
24 | import org.locationtech.proj4j.util.ProjectionMath;
25 |
26 | public class Eckert2Projection extends Projection {
27 |
28 | private final static double FXC = 0.46065886596178063902;
29 | private final static double FYC = 1.44720250911653531871;
30 | private final static double C13 = 0.33333333333333333333;
31 | private final static double ONEEPS = 1.0000001;
32 |
33 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
34 | out.x = FXC * lplam * (out.y = Math.sqrt(4. - 3. * Math.sin(Math.abs(lpphi))));
35 | out.y = FYC * (2. - out.y);
36 | if ( lpphi < 0.) out.y = -out.y;
37 | return out;
38 | }
39 |
40 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
41 | out.x = xyx / (FXC * ( out.y = 2. - Math.abs(xyy) / FYC) );
42 | out.y = (4. - out.y * out.y) * C13;
43 | if (Math.abs(out.y) >= 1.) {
44 | if (Math.abs(out.y) > ONEEPS) throw new ProjectionException("I");
45 | else
46 | out.y = out.y < 0. ? -ProjectionMath.HALFPI : ProjectionMath.HALFPI;
47 | } else
48 | out.y = Math.asin(out.y);
49 | if (xyy < 0)
50 | out.y = -out.y;
51 | return out;
52 | }
53 |
54 | public boolean hasInverse() {
55 | return true;
56 | }
57 |
58 | public String toString() {
59 | return "Eckert II";
60 | }
61 |
62 | }
63 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Eckert4Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class Eckert4Projection extends Projection {
26 |
27 | private final static double C_x = .42223820031577120149;
28 | private final static double C_y = 1.32650042817700232218;
29 | private final static double RC_y = .75386330736002178205;
30 | private final static double C_p = 3.57079632679489661922;
31 | private final static double RC_p = .28004957675577868795;
32 | private final static double EPS = 1e-7;
33 | private final int NITER = 6;
34 |
35 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
36 | double p, V, s, c;
37 | int i;
38 |
39 | p = C_p * Math.sin(lpphi);
40 | V = lpphi * lpphi;
41 | lpphi *= 0.895168 + V * ( 0.0218849 + V * 0.00826809 );
42 | for (i = NITER; i > 0; --i) {
43 | c = Math.cos(lpphi);
44 | s = Math.sin(lpphi);
45 | lpphi -= V = (lpphi + s * (c + 2.) - p) /
46 | (1. + c * (c + 2.) - s * s);
47 | if (Math.abs(V) < EPS)
48 | break;
49 | }
50 | if (i == 0) {
51 | out.x = C_x * lplam;
52 | out.y = lpphi < 0. ? -C_y : C_y;
53 | } else {
54 | out.x = C_x * lplam * (1. + Math.cos(lpphi));
55 | out.y = C_y * Math.sin(lpphi);
56 | }
57 | return out;
58 | }
59 |
60 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
61 | double c;
62 |
63 | out.y = ProjectionMath.asin(xyy / C_y);
64 | out.x = xyx / (C_x * (1. + (c = Math.cos(out.y))));
65 | out.y = ProjectionMath.asin((out.y + Math.sin(out.y) * (c + 2.)) / C_p);
66 | return out;
67 | }
68 |
69 | public boolean hasInverse() {
70 | return true;
71 | }
72 |
73 | public boolean isEqualArea() {
74 | return true;
75 | }
76 |
77 | public String toString() {
78 | return "Eckert IV";
79 | }
80 |
81 | }
82 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Eckert5Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class Eckert5Projection extends Projection {
25 |
26 | private final static double XF = 0.44101277172455148219;
27 | private final static double RXF = 2.26750802723822639137;
28 | private final static double YF = 0.88202554344910296438;
29 | private final static double RYF = 1.13375401361911319568;
30 |
31 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
32 | out.x = XF * (1. + Math.cos(lpphi)) * lplam;
33 | out.y = YF * lpphi;
34 | return out;
35 | }
36 |
37 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
38 | out.x = RXF * xyx / (1. + Math.cos( out.y = RYF * xyy));
39 | return out;
40 | }
41 |
42 | public boolean hasInverse() {
43 | return true;
44 | }
45 |
46 | public String toString() {
47 | return "Eckert V";
48 | }
49 |
50 | }
51 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Eckert6Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 | /*
17 | Copyright 2011 Martin Davis
18 |
19 | Licensed under the Apache License, Version 2.0 (the "License");
20 | you may not use this file except in compliance with the License.
21 | You may obtain a copy of the License at
22 |
23 | http://www.apache.org/licenses/LICENSE-2.0
24 |
25 | Unless required by applicable law or agreed to in writing, software
26 | distributed under the License is distributed on an "AS IS" BASIS,
27 | WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
28 | See the License for the specific language governing permissions and
29 | limitations under the License.
30 | */package org.locationtech.proj4j.proj;
31 |
32 | import org.locationtech.proj4j.*;
33 |
34 | public class Eckert6Projection extends PseudoCylindricalProjection {
35 |
36 | private static final double n = 2.570796326794896619231321691;
37 | private static final double C_y = Math.sqrt((2) / n);
38 | private static final double C_x = C_y / 2;
39 | private static final int MAX_ITER = 8;
40 | private static final double LOOP_TOL = 1e-7;
41 |
42 | public ProjCoordinate project(double lam, double phi, ProjCoordinate xy) {
43 |
44 | int i;
45 | double k, V;
46 | k = n * Math.sin(phi);
47 | for (i = MAX_ITER; i > 0;) {
48 | phi -= V = (phi + Math.sin(phi) - k) / (1 + Math.cos(phi));
49 | if (Math.abs(V) < LOOP_TOL) {
50 | break;
51 | }
52 | --i;
53 | }
54 | if (i == 0) {
55 | throw new ProjectionException("F_ERROR");
56 | }
57 |
58 | xy.x = C_x * lam * (1 + Math.cos(phi));
59 | xy.y = C_y * phi;
60 | return xy;
61 | }
62 |
63 | public ProjCoordinate projectInverse(double x, double y, ProjCoordinate lp) {
64 | y /= C_y;
65 | lp.y = Math.asin((y + Math.sin(y)) / n);
66 | lp.x = x / (C_x * (1 + Math.cos(y)));
67 | return lp;
68 | }
69 |
70 | public boolean hasInverse() {
71 | return true;
72 | }
73 |
74 | public boolean isEqualArea() {
75 | return true;
76 | }
77 |
78 | public String toString() {
79 | return "Eckert VI";
80 | }
81 | }
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/EulerProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class EulerProjection extends SimpleConicProjection {
23 |
24 | public EulerProjection() {
25 | super( SimpleConicProjection.EULER );
26 | }
27 |
28 | public String toString() {
29 | return "Euler";
30 | }
31 |
32 | }
33 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/FaheyProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class FaheyProjection extends Projection {
25 |
26 | private final static double TOL = 1e-6;
27 |
28 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
29 | out.y = 1.819152 * ( out.x = Math.tan(0.5 * lpphi) );
30 | out.x = 0.819152 * lplam * asqrt(1 - out.x * out.x);
31 | return out;
32 | }
33 |
34 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
35 | out.y = 2. * Math.atan(out.y /= 1.819152);
36 | out.x = Math.abs(out.y = 1. - xyy * xyy) < TOL ? 0. :
37 | xyx / (0.819152 * Math.sqrt(xyy));
38 | return out;
39 | }
40 |
41 | private double asqrt(double v) {
42 | return (v <= 0) ? 0. : Math.sqrt(v);
43 | }
44 |
45 | public boolean hasInverse() {
46 | return true;
47 | }
48 |
49 | public String toString() {
50 | return "Fahey";
51 | }
52 |
53 | }
54 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/FoucautProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class FoucautProjection extends SineTangentSeriesProjection {
23 |
24 | public FoucautProjection() {
25 | super( 2., 2., true );
26 | }
27 |
28 | public String toString() {
29 | return "Foucaut";
30 | }
31 |
32 | }
33 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/FoucautSinusoidalProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.ProjectionException;
24 | import org.locationtech.proj4j.util.ProjectionMath;
25 |
26 | public class FoucautSinusoidalProjection extends Projection {
27 | private double n, n1;
28 |
29 | private final static int MAX_ITER = 10;
30 | private final static double LOOP_TOL = 1e-7;
31 |
32 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
33 | double t;
34 |
35 | t = Math.cos(lpphi);
36 | out.x = lplam * t / (n + n1 * t);
37 | out.y = n * lpphi + n1 * Math.sin(lpphi);
38 | return out;
39 | }
40 |
41 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
42 | double V;
43 | int i;
44 |
45 | if (n != 0) {
46 | out.y = xyy;
47 | for (i = MAX_ITER; i > 0; --i) {
48 | out.y -= V = (n * out.y + n1 * Math.sin(out.y) - xyy ) /
49 | (n + n1 * Math.cos(out.y));
50 | if (Math.abs(V) < LOOP_TOL)
51 | break;
52 | }
53 | if (i == 0)
54 | out.y = xyy < 0. ? -ProjectionMath.HALFPI : ProjectionMath.HALFPI;
55 | } else
56 | out.y = ProjectionMath.asin(xyy);
57 | V = Math.cos(out.y);
58 | out.x = xyx * (n + n1 * V) / V;
59 | return out;
60 | }
61 |
62 | public void initialize() {
63 | super.initialize();
64 | // n = pj_param(params, "dn").f;
65 | if (n < 0. || n > 1.)
66 | throw new ProjectionException("-99");
67 | n1 = 1. - n;
68 | }
69 |
70 | public boolean hasInverse() {
71 | return true;
72 | }
73 |
74 | public String toString() {
75 | return "Foucaut Sinusoidal";
76 | }
77 |
78 | }
79 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/GallProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class GallProjection extends Projection {
25 |
26 | private final static double YF = 1.70710678118654752440;
27 | private final static double XF = 0.70710678118654752440;
28 | private final static double RYF = 0.58578643762690495119;
29 | private final static double RXF = 1.41421356237309504880;
30 |
31 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
32 | out.x = XF * lplam;
33 | out.y = YF * Math.tan(.5 * lpphi);
34 | return out;
35 | }
36 |
37 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
38 | out.x = RXF * xyx;
39 | out.y = 2. * Math.atan(xyy * RYF);
40 | return out;
41 | }
42 |
43 | public boolean hasInverse() {
44 | return true;
45 | }
46 |
47 | public String toString() {
48 | return "Gall (Gall Stereographic)";
49 | }
50 |
51 | }
52 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/GeocentProjection.java:
--------------------------------------------------------------------------------
1 | package org.locationtech.proj4j.proj;
2 |
3 | import org.locationtech.proj4j.ProjCoordinate;
4 | import org.locationtech.proj4j.datum.GeocentricConverter;
5 |
6 | public class GeocentProjection extends Projection {
7 |
8 | @Override
9 | public ProjCoordinate projectRadians(ProjCoordinate src, ProjCoordinate dst) {
10 | GeocentricConverter geocentricConverter = new GeocentricConverter(this.ellipsoid);
11 | geocentricConverter.convertGeodeticToGeocentric(dst);
12 | return dst;
13 | }
14 |
15 | @Override
16 | public ProjCoordinate inverseProjectRadians(ProjCoordinate src, ProjCoordinate dst) {
17 | GeocentricConverter geocentricConverter = new GeocentricConverter(this.ellipsoid);
18 | geocentricConverter.convertGeocentricToGeodetic(dst);
19 | return dst;
20 | }
21 | }
22 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Ginsburg8Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class Ginsburg8Projection extends Projection {
25 |
26 | private final static double Cl = 0.000952426;
27 | private final static double Cp = 0.162388;
28 | private final static double C12 = 0.08333333333333333;
29 |
30 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
31 | double t = lpphi * lpphi;
32 |
33 | out.y = lpphi * (1. + t * C12);
34 | out.x = lplam * (1. - Cp * t);
35 | t = lplam * lplam;
36 | out.x *= (0.87 - Cl * t * t);
37 | return out;
38 | }
39 |
40 | public String toString() {
41 | return "Ginsburg VIII (TsNIIGAiK)";
42 | }
43 |
44 | }
45 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/GoodeProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class GoodeProjection extends Projection {
25 |
26 | private final static double Y_COR = 0.05280;
27 | private final static double PHI_LIM = .71093078197902358062;
28 |
29 | private SinusoidalProjection sinu = new SinusoidalProjection();
30 | private MolleweideProjection moll = new MolleweideProjection();
31 |
32 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
33 | if (Math.abs(lpphi) <= PHI_LIM)
34 | out = sinu.project(lplam, lpphi, out);
35 | else {
36 | out = moll.project(lplam, lpphi, out);
37 | out.y -= lpphi >= 0.0 ? Y_COR : -Y_COR;
38 | }
39 | return out;
40 | }
41 |
42 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
43 | if (Math.abs(xyy) <= PHI_LIM)
44 | out = sinu.projectInverse(xyx, xyy, out);
45 | else {
46 | xyy += xyy >= 0.0 ? Y_COR : -Y_COR;
47 | out = moll.projectInverse(xyx, xyy, out);
48 | }
49 | return out;
50 | }
51 |
52 | public boolean hasInverse() {
53 | return true;
54 | }
55 |
56 | public String toString() {
57 | return "Goode Homolosine";
58 | }
59 |
60 | }
61 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/HammerProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import java.util.Objects;
23 |
24 | import org.locationtech.proj4j.ProjCoordinate;
25 | import org.locationtech.proj4j.ProjectionException;
26 |
27 | public class HammerProjection extends PseudoCylindricalProjection {
28 |
29 | private double w = 0.5;
30 | private double m = 1;
31 | private double rm;
32 |
33 | public HammerProjection() {
34 | }
35 |
36 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate xy) {
37 | double cosphi, d;
38 |
39 | d = Math.sqrt(2./(1. + (cosphi = Math.cos(lpphi)) * Math.cos(lplam *= w)));
40 | xy.x = m * d * cosphi * Math.sin(lplam);
41 | xy.y = rm * d * Math.sin(lpphi);
42 | return xy;
43 | }
44 |
45 | public void initialize() {
46 | super.initialize();
47 | if ((w = Math.abs(w)) <= 0.)
48 | throw new ProjectionException("-27");
49 | else
50 | w = .5;
51 | if ((m = Math.abs(m)) <= 0.)
52 | throw new ProjectionException("-27");
53 | else
54 | m = 1.;
55 | rm = 1. / m;
56 | m /= w;
57 | es = 0.;
58 | }
59 |
60 | /**
61 | * Returns true if this projection is equal area
62 | */
63 | public boolean isEqualArea() {
64 | return true;
65 | }
66 |
67 | // Properties
68 | public void setW( double w ) {
69 | this.w = w;
70 | }
71 |
72 | public double getW() {
73 | return w;
74 | }
75 |
76 | public void setM( double m ) {
77 | this.m = m;
78 | }
79 |
80 | public double getM() {
81 | return m;
82 | }
83 |
84 | public String toString() {
85 | return "Hammer & Eckert-Greifendorff";
86 | }
87 |
88 | @Override
89 | public boolean equals(Object that) {
90 | if (this == that) {
91 | return true;
92 | }
93 | if (that instanceof HammerProjection) {
94 | HammerProjection p = (HammerProjection) that;
95 | return (
96 | m == p.m &&
97 | w == p.w &&
98 | super.equals(that));
99 | }
100 | return false;
101 | }
102 |
103 | @Override
104 | public int hashCode() {
105 | return Objects.hash(m, w, super.hashCode());
106 | }
107 | }
108 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/HatanoProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.ProjectionException;
24 | import org.locationtech.proj4j.util.ProjectionMath;
25 |
26 | public class HatanoProjection extends Projection {
27 |
28 | private final static int NITER = 20;
29 | private final static double EPS = 1e-7;
30 | private final static double ONETOL = 1.000001;
31 | private final static double CN = 2.67595;
32 | private final static double CS = 2.43763;
33 | private final static double RCN = 0.37369906014686373063;
34 | private final static double RCS = 0.41023453108141924738;
35 | private final static double FYCN = 1.75859;
36 | private final static double FYCS = 1.93052;
37 | private final static double RYCN = 0.56863737426006061674;
38 | private final static double RYCS = 0.51799515156538134803;
39 | private final static double FXC = 0.85;
40 | private final static double RXC = 1.17647058823529411764;
41 |
42 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
43 | double th1, c;
44 | int i;
45 |
46 | c = Math.sin(lpphi) * (lpphi < 0. ? CS : CN);
47 | for (i = NITER; i > 0; --i) {
48 | lpphi -= th1 = (lpphi + Math.sin(lpphi) - c) / (1. + Math.cos(lpphi));
49 | if (Math.abs(th1) < EPS) break;
50 | }
51 | out.x = FXC * lplam * Math.cos(lpphi *= .5);
52 | out.y = Math.sin(lpphi) * (lpphi < 0. ? FYCS : FYCN);
53 | return out;
54 | }
55 |
56 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
57 | double th;
58 |
59 | th = xyy * ( xyy < 0. ? RYCS : RYCN);
60 | if (Math.abs(th) > 1.)
61 | if (Math.abs(th) > ONETOL) throw new ProjectionException("I");
62 | else th = th > 0. ? ProjectionMath.HALFPI : - ProjectionMath.HALFPI;
63 | else
64 | th = Math.asin(th);
65 | out.x = RXC * xyx / Math.cos(th);
66 | th += th;
67 | out.y = (th + Math.sin(th)) * (xyy < 0. ? RCS : RCN);
68 | if (Math.abs(out.y) > 1.)
69 | if (Math.abs(out.y) > ONETOL) throw new ProjectionException("I");
70 | else out.y = out.y > 0. ? ProjectionMath.HALFPI : - ProjectionMath.HALFPI;
71 | else
72 | out.y = Math.asin(out.y);
73 | return out;
74 | }
75 |
76 | public boolean hasInverse() {
77 | return true;
78 | }
79 |
80 | public String toString() {
81 | return "Hatano Asymmetrical Equal Area";
82 | }
83 |
84 | }
85 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/KavraiskyVProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class KavraiskyVProjection extends SineTangentSeriesProjection {
23 |
24 | public KavraiskyVProjection() {
25 | super( 1.50488, 1.35439, false );
26 | }
27 |
28 | public String toString() {
29 | return "Kavraisky V";
30 | }
31 |
32 | }
33 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/LagrangeProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import java.util.Objects;
23 |
24 | import org.locationtech.proj4j.ProjCoordinate;
25 | import org.locationtech.proj4j.ProjectionException;
26 | import org.locationtech.proj4j.util.ProjectionMath;
27 |
28 | public class LagrangeProjection extends Projection {
29 |
30 | // Parameters
31 | private double hrw;
32 | private double rw = 1.4;
33 | private double a1;
34 | private double phi1;
35 |
36 | private final static double TOL = 1e-10;
37 |
38 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate xy) {
39 | double v, c;
40 |
41 | if ( Math.abs(Math.abs(lpphi) - ProjectionMath.HALFPI) < TOL) {
42 | xy.x = 0;
43 | xy.y = lpphi < 0 ? -2. : 2.;
44 | } else {
45 | lpphi = Math.sin(lpphi);
46 | v = a1 * Math.pow((1. + lpphi)/(1. - lpphi), hrw);
47 | if ((c = 0.5 * (v + 1./v) + Math.cos(lplam *= rw)) < TOL)
48 | throw new ProjectionException();
49 | xy.x = 2. * Math.sin(lplam) / c;
50 | xy.y = (v - 1./v) / c;
51 | }
52 | return xy;
53 | }
54 |
55 | public void setW( double w ) {
56 | this.rw = w;
57 | }
58 |
59 | public double getW() {
60 | return rw;
61 | }
62 |
63 | public void initialize() {
64 | super.initialize();
65 | if (rw <= 0)
66 | throw new ProjectionException("-27");
67 | hrw = 0.5 * (rw = 1. / rw);
68 | phi1 = projectionLatitude1;
69 | if (Math.abs(Math.abs(phi1 = Math.sin(phi1)) - 1.) < TOL)
70 | throw new ProjectionException("-22");
71 | a1 = Math.pow((1. - phi1)/(1. + phi1), hrw);
72 | }
73 |
74 | /**
75 | * Returns true if this projection is conformal
76 | */
77 | public boolean isConformal() {
78 | return true;
79 | }
80 |
81 | public boolean hasInverse() {
82 | return false;
83 | }
84 |
85 | public String toString() {
86 | return "Lagrange";
87 | }
88 |
89 | @Override
90 | public boolean equals(Object that) {
91 | if (this == that) {
92 | return true;
93 | }
94 | if (that instanceof LagrangeProjection) {
95 | LagrangeProjection p = (LagrangeProjection) that;
96 | return (rw == p.rw) && super.equals(that);
97 | }
98 | return false;
99 | }
100 |
101 | @Override
102 | public int hashCode() {
103 | return Objects.hash(rw, super.hashCode());
104 | }
105 | }
106 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/LambertEqualAreaConicProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.util.ProjectionMath;
23 |
24 | public class LambertEqualAreaConicProjection extends AlbersProjection {
25 |
26 | public LambertEqualAreaConicProjection() {
27 | this( false );
28 | }
29 |
30 | public LambertEqualAreaConicProjection( boolean south ) {
31 | minLatitude = Math.toRadians(0);
32 | maxLatitude = Math.toRadians(90);
33 | projectionLatitude1 = south ? -ProjectionMath.QUARTERPI : ProjectionMath.QUARTERPI;
34 | projectionLatitude2 = south ? -ProjectionMath.HALFPI : ProjectionMath.HALFPI;
35 | initialize();
36 | }
37 |
38 | public String toString() {
39 | return "Lambert Equal Area Conic";
40 | }
41 |
42 | }
43 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/LarriveeProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class LarriveeProjection extends Projection {
25 |
26 | private final static double SIXTH = .16666666666666666;
27 |
28 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
29 | out.x = 0.5 * lplam * (1. + Math.sqrt(Math.cos(lpphi)));
30 | out.y = lpphi / (Math.cos(0.5 * lpphi) * Math.cos(SIXTH * lplam));
31 | return out;
32 | }
33 |
34 | public String toString() {
35 | return "Larrivee";
36 | }
37 |
38 | }
39 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/LaskowskiProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 |
25 | public class LaskowskiProjection extends Projection {
26 |
27 | private final static double a10 = 0.975534;
28 | private final static double a12 = -0.119161;
29 | private final static double a32 = -0.0143059;
30 | private final static double a14 = -0.0547009;
31 | private final static double b01 = 1.00384;
32 | private final static double b21 = 0.0802894;
33 | private final static double b03 = 0.0998909;
34 | private final static double b41 = 0.000199025;
35 | private final static double b23 = -0.0285500;
36 | private final static double b05 = -0.0491032;
37 |
38 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
39 | double l2, p2;
40 |
41 | l2 = lplam * lplam;
42 | p2 = lpphi * lpphi;
43 | out.x = lplam * (a10 + p2 * (a12 + l2 * a32 + p2 * a14));
44 | out.y = lpphi * (b01 + l2 * (b21 + p2 * b23 + l2 * b41) +
45 | p2 * (b03 + p2 * b05));
46 | return out;
47 | }
48 |
49 | public String toString() {
50 | return "Laskowski";
51 | }
52 |
53 | }
54 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/LinearProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j.proj;
18 |
19 | import org.locationtech.proj4j.ProjCoordinate;
20 |
21 |
22 | public class LinearProjection extends Projection {
23 |
24 | public ProjCoordinate project(ProjCoordinate src, ProjCoordinate dst) {
25 | dst.x = src.x;
26 | dst.y = src.y;
27 | return dst;
28 | }
29 |
30 | public void transform(double[] srcPoints, int srcOffset, double[] dstPoints, int dstOffset, int numPoints) {
31 | for (int i = 0; i < numPoints; i++) {
32 | dstPoints[dstOffset++] = srcPoints[srcOffset++];
33 | dstPoints[dstOffset++] = srcPoints[srcOffset++];
34 | }
35 | }
36 |
37 | public ProjCoordinate inverseProject(ProjCoordinate src, ProjCoordinate dst) {
38 | dst.x = src.x;
39 | dst.y = src.y;
40 | return dst;
41 | }
42 |
43 | public void inverseTransform(double[] srcPoints, int srcOffset, double[] dstPoints, int dstOffset, int numPoints) {
44 | for (int i = 0; i < numPoints; i++) {
45 | dstPoints[dstOffset++] = srcPoints[srcOffset++];
46 | dstPoints[dstOffset++] = srcPoints[srcOffset++];
47 | }
48 | }
49 |
50 | public boolean hasInverse() {
51 | return true;
52 | }
53 |
54 | public boolean isRectilinear() {
55 | return true;
56 | }
57 |
58 | public String toString() {
59 | return "Linear";
60 | }
61 |
62 | }
63 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/LongLatProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j.proj;
18 |
19 | import org.locationtech.proj4j.units.Units;
20 |
21 | /**
22 | * A "projection" for geodetic coordinates in Decimal Degrees.
23 | */
24 | public class LongLatProjection extends Projection
25 | {
26 | // TODO: implement projection methods (which are basically just no-ops)
27 | /*
28 |
29 | public Point2D.Double transformRadians( Point2D.Double src, Point2D.Double dst ) {
30 | dst.x = src.x;
31 | dst.y = src.y;
32 | return dst;
33 | }
34 |
35 | */
36 |
37 | public String toString() {
38 | return "LongLat";
39 | }
40 |
41 | public void initialize()
42 | {
43 | // units are always in Decimal Degrees
44 | unit = Units.DEGREES;
45 | totalScale = 1.0;
46 | }
47 |
48 | @Override public Boolean isGeographic() {
49 | return true;
50 | }
51 | }
52 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/LoximuthalProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class LoximuthalProjection extends PseudoCylindricalProjection {
26 |
27 | private final static double FC = .92131773192356127802;
28 | private final static double RP = .31830988618379067154;
29 | private final static double EPS = 1e-8;
30 |
31 | private double phi1;
32 | private double cosphi1;
33 | private double tanphi1;
34 |
35 | public LoximuthalProjection() {
36 | phi1 = Math.toRadians(40.0);//FIXME - param
37 | cosphi1 = Math.cos(phi1);
38 | tanphi1 = Math.tan(ProjectionMath.QUARTERPI + 0.5 * phi1);
39 | }
40 |
41 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
42 | double x;
43 | double y = lpphi - phi1;
44 | if (y < EPS)
45 | x = lplam * cosphi1;
46 | else {
47 | x = ProjectionMath.QUARTERPI + 0.5 * lpphi;
48 | if (Math.abs(x) < EPS || Math.abs(Math.abs(x) - ProjectionMath.HALFPI) < EPS)
49 | x = 0.;
50 | else
51 | x = lplam * y / Math.log( Math.tan(x) / tanphi1 );
52 | }
53 | out.x = x;
54 | out.y = y;
55 | return out;
56 | }
57 |
58 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
59 | double latitude = xyy + phi1;
60 | double longitude;
61 | if (Math.abs(xyy) < EPS)
62 | longitude = xyx / cosphi1;
63 | else if (Math.abs( longitude = ProjectionMath.QUARTERPI + 0.5 * xyy ) < EPS ||
64 | Math.abs(Math.abs(xyx) -ProjectionMath.HALFPI) < EPS)
65 | longitude = 0.;
66 | else
67 | longitude = xyx * Math.log( Math.tan(longitude) / tanphi1 ) / xyy;
68 |
69 | out.x = longitude;
70 | out.y = latitude;
71 | return out;
72 | }
73 |
74 | public boolean hasInverse() {
75 | return true;
76 | }
77 |
78 | public String toString() {
79 | return "Loximuthal";
80 | }
81 |
82 | }
83 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/McBrydeThomasFlatPolarParabolicProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.ProjectionException;
24 | import org.locationtech.proj4j.util.ProjectionMath;
25 |
26 | public class McBrydeThomasFlatPolarParabolicProjection extends Projection {
27 |
28 | private final static double CS = .95257934441568037152;
29 | private final static double FXC = .92582009977255146156;
30 | private final static double FYC = 3.40168025708304504493;
31 | private final static double C23 = .66666666666666666666;
32 | private final static double C13 = .33333333333333333333;
33 | private final static double ONEEPS = 1.0000001;
34 |
35 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
36 | out.y = Math.asin(CS * Math.sin(lpphi));
37 | out.x = FXC * lplam * (2. * Math.cos(C23 * lpphi) - 1.);
38 | out.y = FYC * Math.sin(C13 * lpphi);
39 | return out;
40 | }
41 |
42 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
43 | out.y = xyy / FYC;
44 | if (Math.abs(out.y) >= 1.) {
45 | if (Math.abs(out.y) > ONEEPS) throw new ProjectionException("I");
46 | else out.y = (out.y < 0.) ? -ProjectionMath.HALFPI : ProjectionMath.HALFPI;
47 | } else
48 | out.y = Math.asin(out.y);
49 | out.x = xyx / ( FXC * (2. * Math.cos(C23 * (out.y *= 3.)) - 1.) );
50 | if (Math.abs(out.y = Math.sin(out.y) / CS) >= 1.) {
51 | if (Math.abs(out.y) > ONEEPS) throw new ProjectionException("I");
52 | else out.y = (out.y < 0.) ? -ProjectionMath.HALFPI : ProjectionMath.HALFPI;
53 | } else
54 | out.y = Math.asin(out.y);
55 | return out;
56 | }
57 |
58 | public boolean hasInverse() {
59 | return true;
60 | }
61 |
62 | public String toString() {
63 | return "McBride-Thomas Flat-Polar Parabolic";
64 | }
65 |
66 | }
67 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/McBrydeThomasFlatPolarSine1Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class McBrydeThomasFlatPolarSine1Projection extends SineTangentSeriesProjection {
23 |
24 | public McBrydeThomasFlatPolarSine1Projection() {
25 | super( 1.48875, 1.36509, false );
26 | }
27 |
28 | public String toString() {
29 | return "McBryde-Thomas Flat-Polar Sine (No. 1)";
30 | }
31 |
32 | }
33 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/McBrydeThomasFlatPolarSine2Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class McBrydeThomasFlatPolarSine2Projection extends Projection {
26 |
27 | private final static int MAX_ITER = 10;
28 | private final static double LOOP_TOL = 1e-7;
29 | private final static double C1 = 0.45503;
30 | private final static double C2 = 1.36509;
31 | private final static double C3 = 1.41546;
32 | private final static double C_x = 0.22248;
33 | private final static double C_y = 1.44492;
34 | private final static double C1_2 = 0.33333333333333333333333333;
35 |
36 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
37 | double k, V, t;
38 | int i;
39 |
40 | k = C3 * Math.sin(lpphi);
41 | for (i = MAX_ITER; i > 0; i--) {
42 | t = lpphi / C2;
43 | out.y -= V = (C1 * Math.sin(t) + Math.sin(lpphi) - k) /
44 | (C1_2 * Math.cos(t) + Math.cos(lpphi));
45 | if (Math.abs(V) < LOOP_TOL)
46 | break;
47 | }
48 | t = lpphi / C2;
49 | out.x = C_x * lplam * (1. + 3. * Math.cos(lpphi)/Math.cos(t) );
50 | out.y = C_y * Math.sin(t);
51 | return out;
52 | }
53 |
54 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
55 | double t, s;
56 |
57 | out.y = C2 * (t = ProjectionMath.asin(xyy / C_y));
58 | out.x = xyx / (C_x * (1. + 3. * Math.cos(out.y)/Math.cos(t)));
59 | out.y = ProjectionMath.asin((C1 * Math.sin(t) + Math.sin(out.y)) / C3);
60 | return out;
61 | }
62 |
63 | public boolean hasInverse() {
64 | return true;
65 | }
66 |
67 | public String toString() {
68 | return "McBryde-Thomas Flat-Pole Sine (No. 2)";
69 | }
70 |
71 | }
72 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/MercatorProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class MercatorProjection extends CylindricalProjection {
26 |
27 | public MercatorProjection() {
28 | minLatitude = ProjectionMath.degToRad(-85);
29 | maxLatitude = ProjectionMath.degToRad(85);
30 | }
31 |
32 | public ProjCoordinate project(double lam, double phi, ProjCoordinate out) {
33 | if (spherical) {
34 | out.x = scaleFactor * lam;
35 | out.y = scaleFactor * Math.log(Math.tan(ProjectionMath.QUARTERPI + 0.5 * phi));
36 | } else {
37 | out.x = scaleFactor * lam;
38 | out.y = -scaleFactor * Math.log(ProjectionMath.tsfn(phi, Math.sin(phi), e));
39 | }
40 | return out;
41 | }
42 |
43 | public ProjCoordinate projectInverse(double x, double y, ProjCoordinate out) {
44 | if (spherical) {
45 | out.y = ProjectionMath.HALFPI - 2. * Math.atan(Math.exp(-y / scaleFactor));
46 | out.x = x / scaleFactor;
47 | } else {
48 | out.y = ProjectionMath.phi2(Math.exp(-y / scaleFactor), e);
49 | out.x = x / scaleFactor;
50 | }
51 | return out;
52 | }
53 |
54 | public boolean hasInverse() {
55 | return true;
56 | }
57 |
58 | public boolean isRectilinear() {
59 | return true;
60 | }
61 |
62 | /**
63 | * Returns the ESPG code for this projection, or 0 if unknown.
64 | */
65 | public int getEPSGCode() {
66 | return 9804;
67 | }
68 |
69 | public String toString() {
70 | return "Mercator";
71 | }
72 |
73 | }
74 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/MillerProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 |
26 | public class MillerProjection extends CylindricalProjection {
27 |
28 | public MillerProjection() {
29 | }
30 |
31 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
32 | out.x = lplam;
33 | out.y = Math.log(Math.tan(ProjectionMath.QUARTERPI + lpphi * .4)) * 1.25;
34 | return out;
35 | }
36 |
37 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
38 | out.x = xyx;
39 | out.y = 2.5 * (Math.atan(Math.exp(.8 * xyy)) - ProjectionMath.QUARTERPI);
40 | return out;
41 | }
42 |
43 | public boolean hasInverse() {
44 | return true;
45 | }
46 |
47 | public String toString() {
48 | return "Miller Cylindrical";
49 | }
50 |
51 | }
52 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Murdoch1Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class Murdoch1Projection extends SimpleConicProjection {
23 |
24 | public Murdoch1Projection() {
25 | super( SimpleConicProjection.MURD1 );
26 | }
27 |
28 | public String toString() {
29 | return "Murdoch I";
30 | }
31 |
32 | }
33 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Murdoch2Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class Murdoch2Projection extends SimpleConicProjection {
23 |
24 | public Murdoch2Projection() {
25 | super( SimpleConicProjection.MURD2 );
26 | }
27 |
28 | public String toString() {
29 | return "Murdoch II";
30 | }
31 |
32 | }
33 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Murdoch3Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class Murdoch3Projection extends SimpleConicProjection {
23 |
24 | public Murdoch3Projection() {
25 | super( SimpleConicProjection.MURD3 );
26 | }
27 |
28 | public String toString() {
29 | return "Murdoch III";
30 | }
31 |
32 | }
33 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/NellHProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class NellHProjection extends Projection {
26 |
27 | private final static int NITER = 9;
28 | private final static double EPS = 1e-7;
29 |
30 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
31 | out.x = 0.5 * lplam * (1. + Math.cos(lpphi));
32 | out.y = 2.0 * (lpphi - Math.tan(0.5 *lpphi));
33 | return out;
34 | }
35 |
36 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
37 | double V, c, p;
38 | int i;
39 |
40 | p = 0.5 * xyy;
41 | for (i = NITER; i > 0 ; --i) {
42 | c = Math.cos(0.5 * xyy);
43 | out.y -= V = (xyy - Math.tan(xyy/2) - p)/(1. - 0.5/(c*c));
44 | if (Math.abs(V) < EPS)
45 | break;
46 | }
47 | if (i == 0) {
48 | out.y = p < 0. ? -ProjectionMath.HALFPI : ProjectionMath.HALFPI;
49 | out.x = 2. * xyx;
50 | } else
51 | out.x = 2. * xyx / (1. + Math.cos(xyy));
52 | return out;
53 | }
54 |
55 | public boolean hasInverse() {
56 | return true;
57 | }
58 |
59 | public String toString() {
60 | return "Nell-Hammer";
61 | }
62 |
63 | }
64 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/NellProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class NellProjection extends Projection {
26 |
27 | private final static int MAX_ITER = 10;
28 | private final static double LOOP_TOL = 1e-7;
29 |
30 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
31 | double k, V;
32 | int i;
33 |
34 | k = 2. * Math.sin(lpphi);
35 | V = lpphi * lpphi;
36 | out.y *= 1.00371 + V * (-0.0935382 + V * -0.011412);
37 | for (i = MAX_ITER; i > 0 ; --i) {
38 | out.y -= V = (lpphi + Math.sin(lpphi) - k) /
39 | (1. + Math.cos(lpphi));
40 | if (Math.abs(V) < LOOP_TOL)
41 | break;
42 | }
43 | out.x = 0.5 * lplam * (1. + Math.cos(lpphi));
44 | out.y = lpphi;
45 | return out;
46 | }
47 |
48 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
49 | double th, s;
50 |
51 | out.x = 2. * xyx / (1. + Math.cos(xyy));
52 | out.y = ProjectionMath.asin(0.5 * (xyy + Math.sin(xyy)));
53 | return out;
54 | }
55 |
56 | public boolean hasInverse() {
57 | return true;
58 | }
59 |
60 | public String toString() {
61 | return "Nell";
62 | }
63 |
64 | }
65 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/NicolosiProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class NicolosiProjection extends Projection {
26 |
27 | private final static double EPS = 1e-10;
28 |
29 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
30 | if (Math.abs(lplam) < EPS) {
31 | out.x = 0;
32 | out.y = lpphi;
33 | } else if (Math.abs(lpphi) < EPS) {
34 | out.x = lplam;
35 | out.y = 0.;
36 | } else if (Math.abs(Math.abs(lplam) - ProjectionMath.HALFPI) < EPS) {
37 | out.x = lplam * Math.cos(lpphi);
38 | out.y = ProjectionMath.HALFPI * Math.sin(lpphi);
39 | } else if (Math.abs(Math.abs(lpphi) - ProjectionMath.HALFPI) < EPS) {
40 | out.x = 0;
41 | out.y = lpphi;
42 | } else {
43 | double tb, c, d, m, n, r2, sp;
44 |
45 | tb = ProjectionMath.HALFPI / lplam - lplam / ProjectionMath.HALFPI;
46 | c = lpphi / ProjectionMath.HALFPI;
47 | d = (1 - c * c)/((sp = Math.sin(lpphi)) - c);
48 | r2 = tb / d;
49 | r2 *= r2;
50 | m = (tb * sp / d - 0.5 * tb)/(1. + r2);
51 | n = (sp / r2 + 0.5 * d)/(1. + 1./r2);
52 | double x = Math.cos(lpphi);
53 | x = Math.sqrt(m * m + x * x / (1. + r2));
54 | out.x = ProjectionMath.HALFPI * ( m + (lplam < 0. ? -x : x));
55 | double y = Math.sqrt(n * n - (sp * sp / r2 + d * sp - 1.) /
56 | (1. + 1./r2));
57 | out.y = ProjectionMath.HALFPI * ( n + (lpphi < 0. ? y : -y ));
58 | }
59 | return out;
60 | }
61 |
62 | public String toString() {
63 | return "Nicolosi Globular";
64 | }
65 |
66 | }
67 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/NullProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j.proj;
18 |
19 | import org.locationtech.proj4j.ProjCoordinate;
20 |
21 | /**
22 | * A projection which does nothing. Use this for drawing non-geographical overlays.
23 | */
24 | public class NullProjection extends Projection {
25 |
26 | public NullProjection()
27 | {
28 | initialize();
29 | }
30 |
31 | public ProjCoordinate project( ProjCoordinate src, ProjCoordinate dst ) {
32 | dst.x = src.x;
33 | dst.y = src.y;
34 | return dst;
35 | }
36 |
37 | public ProjCoordinate projectInverse( ProjCoordinate src, ProjCoordinate dst ) {
38 | dst.x = src.x;
39 | dst.y = src.y;
40 | return dst;
41 | }
42 |
43 | public boolean isRectilinear() {
44 | return true;
45 | }
46 |
47 | public String toString() {
48 | return "Null";
49 | }
50 |
51 | }
52 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/PerspectiveConicProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class PerspectiveConicProjection extends SimpleConicProjection {
23 |
24 | public PerspectiveConicProjection() {
25 | super( SimpleConicProjection.PCONIC );
26 | }
27 |
28 | public String toString() {
29 | return "Perspective Conic";
30 | }
31 |
32 | }
33 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/PlateCarreeProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j.proj;
18 |
19 |
20 | public class PlateCarreeProjection extends CylindricalProjection {
21 |
22 | public boolean hasInverse() {
23 | return true;
24 | }
25 |
26 | public boolean isRectilinear() {
27 | return true;
28 | }
29 |
30 | public String toString() {
31 | return "Plate Carr\u00e9e";
32 | }
33 |
34 | }
35 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/PseudoCylindricalProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j.proj;
18 |
19 |
20 | /**
21 | * The superclass for all pseudo-cylindrical projections - eg. sinusoidal
22 | * These are projections where parallels are straight, but meridians aren't
23 | */
24 | public class PseudoCylindricalProjection extends CylindricalProjection {
25 |
26 | public boolean isRectilinear() {
27 | return false;
28 | }
29 |
30 | public String toString() {
31 | return "Pseudo Cylindrical";
32 | }
33 |
34 | }
35 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/PutninsP2Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class PutninsP2Projection extends Projection {
26 |
27 | private final static double C_x = 1.89490;
28 | private final static double C_y = 1.71848;
29 | private final static double C_p = 0.6141848493043784;
30 | private final static double EPS = 1e-10;
31 | private final static int NITER = 10;
32 | private final static double PI_DIV_3 = 1.0471975511965977;
33 |
34 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
35 | double p, c, s, V;
36 | int i;
37 |
38 | p = C_p * Math.sin(lpphi);
39 | s = lpphi * lpphi;
40 | out.y *= 0.615709 + s * ( 0.00909953 + s * 0.0046292 );
41 | for (i = NITER; i > 0; --i) {
42 | c = Math.cos(lpphi);
43 | s = Math.sin(lpphi);
44 | out.y -= V = (lpphi + s * (c - 1.) - p) /
45 | (1. + c * (c - 1.) - s * s);
46 | if (Math.abs(V) < EPS)
47 | break;
48 | }
49 | if (i == 0)
50 | out.y = lpphi < 0 ? - PI_DIV_3 : PI_DIV_3;
51 | out.x = C_x * lplam * (Math.cos(lpphi) - 0.5);
52 | out.y = C_y * Math.sin(lpphi);
53 | return out;
54 | }
55 |
56 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
57 | double c;
58 |
59 | out.y = ProjectionMath.asin(xyy / C_y);
60 | out.x = xyx / (C_x * ((c = Math.cos(out.y)) - 0.5));
61 | out.y = ProjectionMath.asin((out.y + Math.sin(out.y) * (c - 1.)) / C_p);
62 | return out;
63 | }
64 |
65 | public boolean hasInverse() {
66 | return true;
67 | }
68 |
69 | public String toString() {
70 | return "Putnins P2";
71 | }
72 |
73 | }
74 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/PutninsP4Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class PutninsP4Projection extends Projection {
26 |
27 | protected double C_x;
28 | protected double C_y;
29 |
30 | public PutninsP4Projection() {
31 | C_x = 0.874038744;
32 | C_y = 3.883251825;
33 | }
34 |
35 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate xy) {
36 | lpphi = ProjectionMath.asin(0.883883476 * Math.sin(lpphi));
37 | xy.x = C_x * lplam * Math.cos(lpphi);
38 | xy.x /= Math.cos(lpphi *= 0.333333333333333);
39 | xy.y = C_y * Math.sin(lpphi);
40 | return xy;
41 | }
42 |
43 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate lp) {
44 | lp.y = ProjectionMath.asin(xyy / C_y);
45 | lp.x = xyx * Math.cos(lp.y) / C_x;
46 | lp.y *= 3.;
47 | lp.x /= Math.cos(lp.y);
48 | lp.y = ProjectionMath.asin(1.13137085 * Math.sin(lp.y));
49 | return lp;
50 | }
51 |
52 | /**
53 | * Returns true if this projection is equal area
54 | */
55 | public boolean isEqualArea() {
56 | return true;
57 | }
58 |
59 | public boolean hasInverse() {
60 | return true;
61 | }
62 |
63 | public String toString() {
64 | return "Putnins P4";
65 | }
66 |
67 | }
68 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/PutninsP5PProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 |
23 | public class PutninsP5PProjection extends PutninsP5Projection {
24 |
25 | public PutninsP5PProjection() {
26 | A = 1.5;
27 | B = 0.5;
28 | }
29 |
30 | public String toString() {
31 | return "Putnins P5P";
32 | }
33 |
34 | }
35 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/PutninsP5Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class PutninsP5Projection extends Projection {
25 |
26 | protected double A;
27 | protected double B;
28 |
29 | private final static double C = 1.01346;
30 | private final static double D = 1.2158542;
31 |
32 | public PutninsP5Projection() {
33 | A = 2;
34 | B = 1;
35 | }
36 |
37 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate xy) {
38 | xy.x = C * lplam * (A - B * Math.sqrt(1. + D * lpphi * lpphi));
39 | xy.y = C * lpphi;
40 | return xy;
41 | }
42 |
43 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate lp) {
44 | lp.y = xyy / C;
45 | lp.x = xyx / (C * (A - B * Math.sqrt(1. + D * lp.y * lp.y)));
46 | return lp;
47 | }
48 |
49 | public boolean hasInverse() {
50 | return true;
51 | }
52 |
53 | public String toString() {
54 | return "Putnins P5";
55 | }
56 |
57 | }
58 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/QuarticAuthalicProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class QuarticAuthalicProjection extends SineTangentSeriesProjection {
23 |
24 | public QuarticAuthalicProjection() {
25 | super( 2., 2., false );
26 | }
27 |
28 | public String toString() {
29 | return "Quartic Authalic";
30 | }
31 |
32 | }
33 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/RectangularPolyconicProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class RectangularPolyconicProjection extends Projection {
25 |
26 | private double phi0;
27 | private double phi1;
28 | private double fxa;
29 | private double fxb;
30 | private boolean mode;
31 |
32 | private final static double EPS = 1e-9;
33 |
34 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
35 | double fa;
36 |
37 | if (mode)
38 | fa = Math.tan(lplam * fxb) * fxa;
39 | else
40 | fa = 0.5 * lplam;
41 | if (Math.abs(lpphi) < EPS) {
42 | out.x = fa + fa;
43 | out.y = - phi0;
44 | } else {
45 | out.y = 1. / Math.tan(lpphi);
46 | out.x = Math.sin(fa = 2. * Math.atan(fa * Math.sin(lpphi))) * out.y;
47 | out.y = lpphi - phi0 + (1. - Math.cos(fa)) * out.y;
48 | }
49 | return out;
50 | }
51 |
52 | public void initialize() { // rpoly
53 | super.initialize();
54 | /*FIXME
55 | if ((mode = (phi1 = Math.abs(pj_param(params, "rlat_ts").f)) > EPS)) {
56 | fxb = 0.5 * Math.sin(phi1);
57 | fxa = 0.5 / fxb;
58 | }
59 | */
60 | }
61 |
62 | public String toString() {
63 | return "Rectangular Polyconic";
64 | }
65 |
66 | }
67 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/SineTangentSeriesProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import java.util.Objects;
23 |
24 | import org.locationtech.proj4j.ProjCoordinate;
25 | import org.locationtech.proj4j.util.ProjectionMath;
26 |
27 | class SineTangentSeriesProjection extends ConicProjection {
28 | private double C_x;
29 | private double C_y;
30 | private double C_p;
31 | private boolean tan_mode;
32 |
33 | protected SineTangentSeriesProjection( double p, double q, boolean mode ) {
34 | es = 0.;
35 | C_x = q / p;
36 | C_y = p;
37 | C_p = 1/ q;
38 | tan_mode = mode;
39 | initialize();
40 | }
41 |
42 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate xy) {
43 | double c;
44 |
45 | xy.x = C_x * lplam * Math.cos(lpphi);
46 | xy.y = C_y;
47 | lpphi *= C_p;
48 | c = Math.cos(lpphi);
49 | if (tan_mode) {
50 | xy.x *= c * c;
51 | xy.y *= Math.tan(lpphi);
52 | } else {
53 | xy.x /= c;
54 | xy.y *= Math.sin(lpphi);
55 | }
56 | return xy;
57 | }
58 |
59 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate lp) {
60 | double c;
61 |
62 | xyy /= C_y;
63 | c = Math.cos(lp.y = tan_mode ? Math.atan(xyy) : ProjectionMath.asin(xyy));
64 | lp.y /= C_p;
65 | lp.x = xyx / (C_x * Math.cos(lp.y));
66 | if (tan_mode)
67 | lp.x /= c * c;
68 | else
69 | lp.x *= c;
70 | return lp;
71 | }
72 |
73 | public boolean hasInverse() {
74 | return true;
75 | }
76 |
77 | @Override
78 | public boolean equals(Object that) {
79 | if (this == that) {
80 | return true;
81 | }
82 | if (that instanceof SineTangentSeriesProjection) {
83 | SineTangentSeriesProjection p = (SineTangentSeriesProjection) that;
84 | return (
85 | C_x == p.C_x &&
86 | C_y == p.C_y &&
87 | C_p == p.C_p &&
88 | tan_mode == p.tan_mode &&
89 | super.equals(that));
90 | }
91 | return false;
92 | }
93 |
94 | @Override
95 | public int hashCode() {
96 | return Objects.hash(C_x, C_y, C_p, tan_mode, super.hashCode());
97 | }
98 | }
99 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/SinusoidalProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class SinusoidalProjection extends PseudoCylindricalProjection {
26 |
27 | public ProjCoordinate project(double lam, double phi, ProjCoordinate xy) {
28 | xy.x = lam * Math.cos(phi);
29 | xy.y = phi;
30 | return xy;
31 | }
32 |
33 | public ProjCoordinate projectInverse(double x, double y, ProjCoordinate lp) {
34 | lp.x = x / Math.cos(y);
35 | lp.y = y;
36 | return lp;
37 | }
38 |
39 | public double getWidth(double y) {
40 | return ProjectionMath.normalizeLongitude(Math.PI) * Math.cos(y);
41 | }
42 |
43 | public boolean hasInverse() {
44 | return true;
45 | }
46 |
47 | public boolean isEqualArea() {
48 | return true;
49 | }
50 |
51 | public String toString() {
52 | return "Sinusoidal";
53 | }
54 |
55 | }
56 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/TissotProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class TissotProjection extends SimpleConicProjection {
23 |
24 | public TissotProjection() {
25 | super( SimpleConicProjection.TISSOT );
26 | }
27 |
28 | public String toString() {
29 | return "Tissot";
30 | }
31 |
32 | }
33 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/TransverseCylindricalEqualArea.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class TransverseCylindricalEqualArea extends Projection {
25 |
26 | private double rk0;
27 |
28 | public TransverseCylindricalEqualArea() {
29 | initialize();
30 | }
31 |
32 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
33 | out.x = rk0 * Math.cos(lpphi) * Math.sin(lplam);
34 | out.y = scaleFactor * (Math.atan2(Math.tan(lpphi), Math.cos(lplam)) - projectionLatitude);
35 | return out;
36 | }
37 |
38 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
39 | double t;
40 |
41 | out.y = xyy * rk0 + projectionLatitude;
42 | out.x *= scaleFactor;
43 | t = Math.sqrt(1. - xyx * xyx);
44 | out.y = Math.asin(t * Math.sin(xyy));
45 | out.x = Math.atan2(xyx, t * Math.cos(xyy));
46 | return out;
47 | }
48 |
49 | public void initialize() { // tcea
50 | super.initialize();
51 | rk0 = 1 / scaleFactor;
52 | }
53 |
54 | public boolean isRectilinear() {
55 | return false;
56 | }
57 |
58 | public boolean hasInverse() {
59 | return true;
60 | }
61 |
62 | public String toString() {
63 | return "Transverse Cylindrical Equal Area";
64 | }
65 |
66 | }
67 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/TranverseCentralCylindricalProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.ProjectionException;
24 | import org.locationtech.proj4j.util.ProjectionMath;
25 |
26 | public class TranverseCentralCylindricalProjection extends CylindricalProjection {
27 |
28 | public TranverseCentralCylindricalProjection() {
29 | minLongitude = ProjectionMath.degToRad(-60);
30 | maxLongitude = ProjectionMath.degToRad(60);
31 | }
32 |
33 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
34 | double b, bt;
35 |
36 | b = Math.cos(lpphi) * Math.sin(lplam);
37 | if ((bt = 1. - b * b) < EPS10)
38 | throw new ProjectionException("F");
39 | out.x = b / Math.sqrt(bt);
40 | out.y = Math.atan2(Math.tan(lpphi), Math.cos(lplam));
41 | return out;
42 | }
43 |
44 | public boolean isRectilinear() {
45 | return false;
46 | }
47 |
48 | public String toString() {
49 | return "Transverse Central Cylindrical";
50 | }
51 |
52 | }
53 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/UrmaevFlatPolarSinusoidalProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import java.util.Objects;
23 |
24 | import org.locationtech.proj4j.ProjCoordinate;
25 | import org.locationtech.proj4j.ProjectionException;
26 | import org.locationtech.proj4j.util.ProjectionMath;
27 |
28 | public class UrmaevFlatPolarSinusoidalProjection extends Projection {
29 |
30 | private final static double C_x = 0.8773826753;
31 | private final static double Cy = 1.139753528477;
32 |
33 | private double n = 0.8660254037844386467637231707;// wag1
34 | private double C_y;
35 |
36 | public UrmaevFlatPolarSinusoidalProjection() {
37 | }
38 |
39 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
40 | out.y = ProjectionMath.asin(n * Math.sin(lpphi));
41 | out.x = C_x * lplam * Math.cos(lpphi);
42 | out.y = C_y * lpphi;
43 | return out;
44 | }
45 |
46 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
47 | xyy /= C_y;
48 | out.y = ProjectionMath.asin(Math.sin(xyy) / n);
49 | out.x = xyx / (C_x * Math.cos(xyy));
50 | return out;
51 | }
52 |
53 | public boolean hasInverse() {
54 | return true;
55 | }
56 |
57 | public void initialize() { // urmfps
58 | super.initialize();
59 | if (n <= 0. || n > 1.)
60 | throw new ProjectionException("-40");
61 | C_y = Cy / n;
62 | }
63 |
64 | // Properties
65 | public void setN( double n ) {
66 | this.n = n;
67 | }
68 |
69 | public double getN() {
70 | return n;
71 | }
72 |
73 | public String toString() {
74 | return "Urmaev Flat-Polar Sinusoidal";
75 | }
76 |
77 | @Override
78 | public boolean equals(Object that) {
79 | if (this == that) {
80 | return true;
81 | }
82 | if (that instanceof UrmaevFlatPolarSinusoidalProjection) {
83 | UrmaevFlatPolarSinusoidalProjection p = (UrmaevFlatPolarSinusoidalProjection) that;
84 | return (n == p.n) && super.equals(that);
85 | }
86 | return false;
87 | }
88 |
89 | @Override
90 | public int hashCode() {
91 | return Objects.hash(n, super.hashCode());
92 | }
93 | }
94 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/VitkovskyProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class VitkovskyProjection extends SimpleConicProjection {
23 |
24 | public VitkovskyProjection() {
25 | super( SimpleConicProjection.VITK1 );
26 | }
27 |
28 | }
29 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Wagner1Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class Wagner1Projection extends UrmaevFlatPolarSinusoidalProjection {
23 |
24 | public String toString() {
25 | return "Wagner I";
26 | }
27 | }
28 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Wagner2Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 | import org.locationtech.proj4j.util.ProjectionMath;
24 |
25 | public class Wagner2Projection extends Projection {
26 |
27 | private final static double C_x = 0.92483;
28 | private final static double C_y = 1.38725;
29 | private final static double C_p1 = 0.88022;
30 | private final static double C_p2 = 0.88550;
31 |
32 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
33 | out.y = ProjectionMath.asin(C_p1 * Math.sin(C_p2 * lpphi));
34 | out.x = C_x * lplam * Math.cos(lpphi);
35 | out.y = C_y * lpphi;
36 | return out;
37 | }
38 |
39 | public ProjCoordinate projectInverse(double xyx, double xyy, ProjCoordinate out) {
40 | out.y = xyy / C_y;
41 | out.x = xyx / (C_x * Math.cos(out.y));
42 | out.y = ProjectionMath.asin(Math.sin(out.y) / C_p1) / C_p2;
43 | return out;
44 | }
45 |
46 | public boolean hasInverse() {
47 | return true;
48 | }
49 |
50 | public String toString() {
51 | return "Wagner II";
52 | }
53 | }
54 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Wagner3Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class Wagner3Projection extends PseudoCylindricalProjection {
25 |
26 | private final static double TWOTHIRD = 0.6666666666666666666667;
27 |
28 | private double C_x;
29 |
30 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate xy) {
31 | xy.x = C_x * lplam * Math.cos(TWOTHIRD * lpphi);
32 | xy.y = lpphi;
33 | return xy;
34 | }
35 |
36 | public ProjCoordinate projectInverse(double x, double y, ProjCoordinate lp) {
37 | lp.y = y;
38 | lp.x = x / (C_x * Math.cos(TWOTHIRD * lp.y));
39 | return lp;
40 | }
41 |
42 | public void initialize() {
43 | super.initialize();
44 | C_x = Math.cos(trueScaleLatitude) / Math.cos(2.*trueScaleLatitude/3.);
45 | es = 0.;
46 | }
47 |
48 | public boolean hasInverse() {
49 | return true;
50 | }
51 |
52 | public String toString() {
53 | return "Wagner III";
54 | }
55 |
56 | }
57 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Wagner4Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class Wagner4Projection extends MolleweideProjection {
23 |
24 | public Wagner4Projection() {
25 | super( MolleweideProjection.WAGNER4 );
26 | }
27 |
28 | }
29 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Wagner5Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | public class Wagner5Projection extends MolleweideProjection {
23 |
24 | public Wagner5Projection() {
25 | super( MolleweideProjection.WAGNER5 );
26 | }
27 |
28 | }
29 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/Wagner7Projection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 | import org.locationtech.proj4j.ProjCoordinate;
23 |
24 | public class Wagner7Projection extends Projection {
25 |
26 | public ProjCoordinate project(double lplam, double lpphi, ProjCoordinate out) {
27 | double theta, ct, D;
28 |
29 | theta = Math.asin(out.y = 0.90630778703664996 * Math.sin(lpphi));
30 | out.x = 2.66723 * (ct = Math.cos(theta)) * Math.sin(lplam /= 3.);
31 | out.y *= 1.24104 * (D = 1/(Math.sqrt(0.5 * (1 + ct * Math.cos(lplam)))));
32 | out.x *= D;
33 | return out;
34 | }
35 |
36 | /**
37 | * Returns true if this projection is equal area
38 | */
39 | public boolean isEqualArea() {
40 | return true;
41 | }
42 |
43 | public String toString() {
44 | return "Wagner VII";
45 | }
46 |
47 | }
48 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/WerenskioldProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 |
23 | public class WerenskioldProjection extends PutninsP4Projection {
24 |
25 | public WerenskioldProjection() {
26 | C_x = 1;
27 | C_y = 4.442882938;
28 | }
29 |
30 | public String toString() {
31 | return "Werenskiold I";
32 | }
33 |
34 | }
35 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/proj/WinkelTripelProjection.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /*
18 | * This file was semi-automatically converted from the public-domain USGS PROJ source.
19 | */
20 | package org.locationtech.proj4j.proj;
21 |
22 |
23 | public class WinkelTripelProjection extends AitoffProjection {
24 |
25 | public WinkelTripelProjection() {
26 | super( WINKEL, 0.636619772367581343 );
27 | }
28 |
29 | public String toString() {
30 | return "Winkel Tripel";
31 | }
32 |
33 | }
34 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/units/DegreeUnit.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2006, 2017 Jerry Huxtable, Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | package org.locationtech.proj4j.units;
18 |
19 |
20 | public class DegreeUnit extends Unit {
21 |
22 | static final long serialVersionUID = -3212757578604686538L;
23 |
24 | private static AngleFormat format = new AngleFormat(AngleFormat.ddmmssPattern, true);
25 |
26 | public DegreeUnit() {
27 | super("degree", "degrees", "deg", 1);
28 | }
29 |
30 | public double parse(String s) throws NumberFormatException {
31 | try {
32 | return format.parse(s).doubleValue();
33 | }
34 | catch (java.text.ParseException e) {
35 | throw new NumberFormatException(e.getMessage());
36 | }
37 | }
38 |
39 | public String format(double n) {
40 | return format.format(n)+" "+abbreviation;
41 | }
42 |
43 | public String format(double n, boolean abbrev) {
44 | if (abbrev)
45 | return format.format(n)+" "+abbreviation;
46 | return format.format(n);
47 | }
48 |
49 | public String format(double x, double y, boolean abbrev) {
50 | if (abbrev)
51 | return format.format(x)+"/"+format.format(y)+" "+abbreviation;
52 | return format.format(x)+"/"+format.format(y);
53 | }
54 | }
55 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/util/Complex.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 | package org.locationtech.proj4j.util;
17 |
18 | public final class Complex {
19 |
20 | public double r;
21 |
22 | public double i;
23 |
24 | public Complex(double r, double i) {
25 | this.r = r;
26 | this.i = i;
27 | }
28 |
29 | public Complex(Complex that) {
30 | this(that.r, that.i);
31 | }
32 |
33 | public int hashCode() {
34 | return Double.valueOf(r).hashCode() | 37 * Double.valueOf(i).hashCode();
35 | }
36 |
37 | public boolean equals(Object that) {
38 | if (that instanceof Complex) {
39 | return ((Complex) that).r == r && ((Complex)that).i == i;
40 | } else {
41 | return false;
42 | }
43 | }
44 | }
45 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/util/FloatPolarCoordinate.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 | package org.locationtech.proj4j.util;
17 |
18 | import java.io.Serializable;
19 |
20 | public final class FloatPolarCoordinate implements Serializable {
21 |
22 | public float lam, phi;
23 |
24 | public FloatPolarCoordinate(FloatPolarCoordinate that) {
25 | this(that.lam, that.phi);
26 | }
27 |
28 | public FloatPolarCoordinate(float lam, float phi) {
29 | this.lam = lam;
30 | this.phi = phi;
31 | }
32 |
33 | @Override
34 | public int hashCode() {
35 | return new Float(lam).hashCode() | (17 * new Float(phi).hashCode());
36 | }
37 |
38 | @Override
39 | public boolean equals(Object that) {
40 | if (that instanceof FloatPolarCoordinate) {
41 | FloatPolarCoordinate c = (FloatPolarCoordinate) that;
42 | return lam == c.lam && phi == c.phi;
43 | } else {
44 | return false;
45 | }
46 | }
47 | }
48 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/util/IntPolarCoordinate.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 | package org.locationtech.proj4j.util;
17 |
18 | import java.io.Serializable;
19 |
20 | public final class IntPolarCoordinate implements Serializable {
21 |
22 | public int lam, phi;
23 |
24 | public IntPolarCoordinate(IntPolarCoordinate that) {
25 | this(that.lam, that.phi);
26 | }
27 |
28 | public IntPolarCoordinate(int lam, int phi) {
29 | this.lam = lam;
30 | this.phi = phi;
31 | }
32 |
33 | @Override
34 | public String toString() {
35 | return String.format("ILP %x:%x", lam, phi);
36 | }
37 |
38 | @Override
39 | public int hashCode() {
40 | return lam | (17 * phi);
41 | }
42 |
43 | @Override
44 | public boolean equals(Object that) {
45 | if (that instanceof IntPolarCoordinate) {
46 | IntPolarCoordinate c = (IntPolarCoordinate) that;
47 | return lam == c.lam && phi == c.phi;
48 | } else {
49 | return false;
50 | }
51 | }
52 | }
53 |
--------------------------------------------------------------------------------
/core/src/main/java/org/locationtech/proj4j/util/Pair.java:
--------------------------------------------------------------------------------
1 | /*******************************************************************************
2 | * Copyright 2009, 2017 Martin Davis
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 | package org.locationtech.proj4j.util;
17 |
18 | import java.util.Objects;
19 |
20 | public class Pair There is one exception to the above paragraph: whether an object is a geographic or projected CRS.
39 | * Because the type of a Java object cannot change dynamically, whether a CRS is geographic or projected
40 | * is determined at {@code geoapi(CoordinateReferenceSystem)} invocation time.Dependency to a Unit of Measurement library
24 | * This module requires a JSR-385 (Units of Measurement) implementation
25 | * to be present on the class-path or module-path.
26 | * The choice of an implementation is left to the user. Some implementations are
27 | * Indriya,
28 | * Seshat and
29 | * Apache SIS.
30 | * The two latter support EPSG codes for units of measurement.
31 | *
32 | * Mutability
33 | * No information is copied. All methods of the views delegate their work to the PROJ4J implementation.
34 | * Consequently, since PROJ4J objects are mutable, changes to the wrapped PROJ4J object are immediately
35 | * reflected in the view. However, it is not recommended to change a wrapped PROJ4J object as CRS should
36 | * be immutable.
37 | *
38 | * Serialization
43 | * The serialization details are not committed API.
44 | * Serialization is okay for exchanging objects between JVM running the same version of PROJ4J,
45 | * but is not guaranteed to be compatible between different versions of PROJ4J. This module does not define
46 | * {@code serialVersionUID} because the backing PROJ4J objects do not define those UID anyway.
47 | *
48 | * @author Martin Desruisseaux (Geomatys)
49 | */
50 | package org.locationtech.proj4j.geoapi;
51 |
--------------------------------------------------------------------------------
/geoapi/src/main/java/org/locationtech/proj4j/geoapi/spi/package-info.java:
--------------------------------------------------------------------------------
1 | /*
2 | * Copyright 2025, PROJ4J contributors
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 |
17 | /**
18 | * Registration of PROJ4J wrappers as referencing services.
19 | * Developers should not use this package directly, as it may change in any future version.
20 | * In particular, it may be simplified if PROJ4J migrates to Java 9+ module system.
21 | *
22 | * @author Martin Desruisseaux (Geomatys)
23 | */
24 | package org.locationtech.proj4j.geoapi.spi;
25 |
--------------------------------------------------------------------------------
/geoapi/src/main/resources/META-INF/services/org.opengis.referencing.crs.CRSAuthorityFactory:
--------------------------------------------------------------------------------
1 | org.locationtech.proj4j.geoapi.spi.AuthorityFactory
2 |
--------------------------------------------------------------------------------
/geoapi/src/main/resources/META-INF/services/org.opengis.referencing.operation.CoordinateOperationFactory:
--------------------------------------------------------------------------------
1 | org.locationtech.proj4j.geoapi.spi.OperationFactory
2 |
--------------------------------------------------------------------------------
/geoapi/src/test/java/org/locationtech/proj4j/geoapi/ServicesTest.java:
--------------------------------------------------------------------------------
1 | /*
2 | * Copyright 2025, PROJ4J contributors
3 | *
4 | * Licensed under the Apache License, Version 2.0 (the "License");
5 | * you may not use this file except in compliance with the License.
6 | * You may obtain a copy of the License at
7 | *
8 | * http://www.apache.org/licenses/LICENSE-2.0
9 | *
10 | * Unless required by applicable law or agreed to in writing, software
11 | * distributed under the License is distributed on an "AS IS" BASIS,
12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
13 | * See the License for the specific language governing permissions and
14 | * limitations under the License.
15 | */
16 | package org.locationtech.proj4j.geoapi;
17 |
18 | import java.util.Iterator;
19 | import java.util.ServiceLoader;
20 | import org.junit.Test;
21 | import org.opengis.referencing.crs.CRSAuthorityFactory;
22 | import org.opengis.referencing.crs.CoordinateReferenceSystem;
23 | import org.opengis.referencing.operation.CoordinateOperationFactory;
24 | import org.opengis.util.FactoryException;
25 |
26 | import static org.locationtech.proj4j.CoordinateReferenceSystem.CS_GEO;
27 | import static org.junit.Assert.*;
28 |
29 |
30 | /**
31 | * Tests fetching factory instances as services.
32 | *
33 | * @author Martin Desruisseaux (Geomatys)
34 | */
35 | public class ServicesTest {
36 | /**
37 | * Creates a new test case.
38 | */
39 | public ServicesTest() {
40 | }
41 |
42 | /**
43 | * Returns the factory of the given type, making sure that there is exactly one instance.
44 | */
45 | private static