| 1 | // License: GPL. For details, see LICENSE file.
|
|---|
| 2 | package org.openstreetmap.josm.data.projection;
|
|---|
| 3 |
|
|---|
| 4 | import java.io.BufferedReader;
|
|---|
| 5 | import java.io.IOException;
|
|---|
| 6 | import java.io.InputStream;
|
|---|
| 7 | import java.io.InputStreamReader;
|
|---|
| 8 | import java.nio.charset.StandardCharsets;
|
|---|
| 9 | import java.util.ArrayList;
|
|---|
| 10 | import java.util.Collection;
|
|---|
| 11 | import java.util.Collections;
|
|---|
| 12 | import java.util.HashMap;
|
|---|
| 13 | import java.util.HashSet;
|
|---|
| 14 | import java.util.List;
|
|---|
| 15 | import java.util.Locale;
|
|---|
| 16 | import java.util.Map;
|
|---|
| 17 | import java.util.Set;
|
|---|
| 18 | import java.util.regex.Matcher;
|
|---|
| 19 | import java.util.regex.Pattern;
|
|---|
| 20 |
|
|---|
| 21 | import org.openstreetmap.josm.Main;
|
|---|
| 22 | import org.openstreetmap.josm.data.coor.EastNorth;
|
|---|
| 23 | import org.openstreetmap.josm.data.coor.LatLon;
|
|---|
| 24 | import org.openstreetmap.josm.data.projection.datum.Datum;
|
|---|
| 25 | import org.openstreetmap.josm.data.projection.datum.GRS80Datum;
|
|---|
| 26 | import org.openstreetmap.josm.data.projection.datum.NTV2GridShiftFileWrapper;
|
|---|
| 27 | import org.openstreetmap.josm.data.projection.datum.WGS84Datum;
|
|---|
| 28 | import org.openstreetmap.josm.data.projection.proj.ClassProjFactory;
|
|---|
| 29 | import org.openstreetmap.josm.data.projection.proj.LambertConformalConic;
|
|---|
| 30 | import org.openstreetmap.josm.data.projection.proj.LonLat;
|
|---|
| 31 | import org.openstreetmap.josm.data.projection.proj.Mercator;
|
|---|
| 32 | import org.openstreetmap.josm.data.projection.proj.Proj;
|
|---|
| 33 | import org.openstreetmap.josm.data.projection.proj.ProjFactory;
|
|---|
| 34 | import org.openstreetmap.josm.data.projection.proj.SwissObliqueMercator;
|
|---|
| 35 | import org.openstreetmap.josm.data.projection.proj.TransverseMercator;
|
|---|
| 36 | import org.openstreetmap.josm.gui.preferences.projection.ProjectionChoice;
|
|---|
| 37 | import org.openstreetmap.josm.gui.preferences.projection.ProjectionPreference;
|
|---|
| 38 | import org.openstreetmap.josm.io.CachedFile;
|
|---|
| 39 | import org.openstreetmap.josm.tools.Pair;
|
|---|
| 40 | import org.openstreetmap.josm.tools.Utils;
|
|---|
| 41 |
|
|---|
| 42 | /**
|
|---|
| 43 | * Class to handle projections
|
|---|
| 44 | *
|
|---|
| 45 | */
|
|---|
| 46 | public final class Projections {
|
|---|
| 47 |
|
|---|
| 48 | private Projections() {
|
|---|
| 49 | // Hide default constructor for utils classes
|
|---|
| 50 | }
|
|---|
| 51 |
|
|---|
| 52 | public static EastNorth project(LatLon ll) {
|
|---|
| 53 | if (ll == null) return null;
|
|---|
| 54 | return Main.getProjection().latlon2eastNorth(ll);
|
|---|
| 55 | }
|
|---|
| 56 |
|
|---|
| 57 | public static LatLon inverseProject(EastNorth en) {
|
|---|
| 58 | if (en == null) return null;
|
|---|
| 59 | return Main.getProjection().eastNorth2latlon(en);
|
|---|
| 60 | }
|
|---|
| 61 |
|
|---|
| 62 | /*********************************
|
|---|
| 63 | * Registry for custom projection
|
|---|
| 64 | *
|
|---|
| 65 | * should be compatible to PROJ.4
|
|---|
| 66 | */
|
|---|
| 67 | static final Map<String, ProjFactory> projs = new HashMap<>();
|
|---|
| 68 | static final Map<String, Ellipsoid> ellipsoids = new HashMap<>();
|
|---|
| 69 | static final Map<String, Datum> datums = new HashMap<>();
|
|---|
| 70 | static final Map<String, NTV2GridShiftFileWrapper> nadgrids = new HashMap<>();
|
|---|
| 71 | static final Map<String, Pair<String, String>> inits = new HashMap<>();
|
|---|
| 72 |
|
|---|
| 73 | static {
|
|---|
| 74 | registerBaseProjection("lonlat", LonLat.class, "core");
|
|---|
| 75 | registerBaseProjection("josm:smerc", Mercator.class, "core");
|
|---|
| 76 | registerBaseProjection("lcc", LambertConformalConic.class, "core");
|
|---|
| 77 | registerBaseProjection("somerc", SwissObliqueMercator.class, "core");
|
|---|
| 78 | registerBaseProjection("tmerc", TransverseMercator.class, "core");
|
|---|
| 79 |
|
|---|
| 80 | ellipsoids.put("airy", Ellipsoid.Airy);
|
|---|
| 81 | ellipsoids.put("mod_airy", Ellipsoid.AiryMod);
|
|---|
| 82 | ellipsoids.put("aust_SA", Ellipsoid.AustSA);
|
|---|
| 83 | ellipsoids.put("bessel", Ellipsoid.Bessel1841);
|
|---|
| 84 | ellipsoids.put("clrk66", Ellipsoid.Clarke1866);
|
|---|
| 85 | ellipsoids.put("clarkeIGN", Ellipsoid.ClarkeIGN);
|
|---|
| 86 | ellipsoids.put("intl", Ellipsoid.Hayford);
|
|---|
| 87 | ellipsoids.put("helmert", Ellipsoid.Helmert);
|
|---|
| 88 | ellipsoids.put("krass", Ellipsoid.Krassowsky);
|
|---|
| 89 | ellipsoids.put("GRS67", Ellipsoid.GRS67);
|
|---|
| 90 | ellipsoids.put("GRS80", Ellipsoid.GRS80);
|
|---|
| 91 | ellipsoids.put("WGS72", Ellipsoid.WGS72);
|
|---|
| 92 | ellipsoids.put("WGS84", Ellipsoid.WGS84);
|
|---|
| 93 |
|
|---|
| 94 | datums.put("WGS84", WGS84Datum.INSTANCE);
|
|---|
| 95 | datums.put("GRS80", GRS80Datum.INSTANCE);
|
|---|
| 96 |
|
|---|
| 97 | nadgrids.put("BETA2007.gsb", NTV2GridShiftFileWrapper.BETA2007);
|
|---|
| 98 | nadgrids.put("ntf_r93_b.gsb", NTV2GridShiftFileWrapper.ntf_rgf93);
|
|---|
| 99 |
|
|---|
| 100 | loadInits();
|
|---|
| 101 | }
|
|---|
| 102 |
|
|---|
| 103 | /**
|
|---|
| 104 | * Plugins can register additional base projections.
|
|---|
| 105 | *
|
|---|
| 106 | * @param id The "official" PROJ.4 id. In case the projection is not supported
|
|---|
| 107 | * by PROJ.4, use some prefix, e.g. josm:myproj or gdal:otherproj.
|
|---|
| 108 | * @param fac The base projection factory.
|
|---|
| 109 | * @param origin Multiple plugins may implement the same base projection.
|
|---|
| 110 | * Provide plugin name or similar string, so it be differentiated.
|
|---|
| 111 | */
|
|---|
| 112 | public static void registerBaseProjection(String id, ProjFactory fac, String origin) {
|
|---|
| 113 | projs.put(id, fac);
|
|---|
| 114 | }
|
|---|
| 115 |
|
|---|
| 116 | public static void registerBaseProjection(String id, Class<? extends Proj> projClass, String origin) {
|
|---|
| 117 | registerBaseProjection(id, new ClassProjFactory(projClass), origin);
|
|---|
| 118 | }
|
|---|
| 119 |
|
|---|
| 120 | public static Proj getBaseProjection(String id) {
|
|---|
| 121 | ProjFactory fac = projs.get(id);
|
|---|
| 122 | if (fac == null) return null;
|
|---|
| 123 | return fac.createInstance();
|
|---|
| 124 | }
|
|---|
| 125 |
|
|---|
| 126 | public static Ellipsoid getEllipsoid(String id) {
|
|---|
| 127 | return ellipsoids.get(id);
|
|---|
| 128 | }
|
|---|
| 129 |
|
|---|
| 130 | public static Datum getDatum(String id) {
|
|---|
| 131 | return datums.get(id);
|
|---|
| 132 | }
|
|---|
| 133 |
|
|---|
| 134 | public static NTV2GridShiftFileWrapper getNTV2Grid(String id) {
|
|---|
| 135 | return nadgrids.get(id);
|
|---|
| 136 | }
|
|---|
| 137 |
|
|---|
| 138 | /**
|
|---|
| 139 | * Get the projection definition string for the given id.
|
|---|
| 140 | * @param id the id
|
|---|
| 141 | * @return the string that can be processed by #{link CustomProjection}.
|
|---|
| 142 | * Null, if the id isn't supported.
|
|---|
| 143 | */
|
|---|
| 144 | public static String getInit(String id) {
|
|---|
| 145 | Pair<String, String> r = inits.get(id.toUpperCase(Locale.ENGLISH));
|
|---|
| 146 | if (r == null) return null;
|
|---|
| 147 | return r.b;
|
|---|
| 148 | }
|
|---|
| 149 |
|
|---|
| 150 | /**
|
|---|
| 151 | * Load +init "presets" from file
|
|---|
| 152 | */
|
|---|
| 153 | private static void loadInits() {
|
|---|
| 154 | Pattern epsgPattern = Pattern.compile("<(\\d+)>(.*)<>");
|
|---|
| 155 | try (
|
|---|
| 156 | InputStream in = new CachedFile("resource://data/projection/epsg").getInputStream();
|
|---|
| 157 | BufferedReader r = new BufferedReader(new InputStreamReader(in, StandardCharsets.UTF_8));
|
|---|
| 158 | ) {
|
|---|
| 159 | String line, lastline = "";
|
|---|
| 160 | while ((line = r.readLine()) != null) {
|
|---|
| 161 | line = line.trim();
|
|---|
| 162 | if (!line.startsWith("#") && !line.isEmpty()) {
|
|---|
| 163 | if (!lastline.startsWith("#")) throw new AssertionError("EPSG file seems corrupted");
|
|---|
| 164 | String name = lastline.substring(1).trim();
|
|---|
| 165 | Matcher m = epsgPattern.matcher(line);
|
|---|
| 166 | if (m.matches()) {
|
|---|
| 167 | inits.put("EPSG:" + m.group(1), Pair.create(name, m.group(2).trim()));
|
|---|
| 168 | } else {
|
|---|
| 169 | Main.warn("Failed to parse line from the EPSG projection definition: "+line);
|
|---|
| 170 | }
|
|---|
| 171 | }
|
|---|
| 172 | lastline = line;
|
|---|
| 173 | }
|
|---|
| 174 | } catch (IOException ex) {
|
|---|
| 175 | throw new RuntimeException(ex);
|
|---|
| 176 | }
|
|---|
| 177 | }
|
|---|
| 178 |
|
|---|
| 179 | private static final Set<String> allCodes = new HashSet<>();
|
|---|
| 180 | private static final Map<String, ProjectionChoice> allProjectionChoicesByCode = new HashMap<>();
|
|---|
| 181 | private static final Map<String, Projection> projectionsByCode_cache = new HashMap<>();
|
|---|
| 182 |
|
|---|
| 183 | static {
|
|---|
| 184 | for (ProjectionChoice pc : ProjectionPreference.getProjectionChoices()) {
|
|---|
| 185 | for (String code : pc.allCodes()) {
|
|---|
| 186 | allProjectionChoicesByCode.put(code, pc);
|
|---|
| 187 | }
|
|---|
| 188 | }
|
|---|
| 189 | allCodes.addAll(inits.keySet());
|
|---|
| 190 | allCodes.addAll(allProjectionChoicesByCode.keySet());
|
|---|
| 191 | }
|
|---|
| 192 |
|
|---|
| 193 | public static Projection getProjectionByCode(String code) {
|
|---|
| 194 | Projection proj = projectionsByCode_cache.get(code);
|
|---|
| 195 | if (proj != null) return proj;
|
|---|
| 196 | ProjectionChoice pc = allProjectionChoicesByCode.get(code);
|
|---|
| 197 | if (pc != null) {
|
|---|
| 198 | Collection<String> pref = pc.getPreferencesFromCode(code);
|
|---|
| 199 | pc.setPreferences(pref);
|
|---|
| 200 | try {
|
|---|
| 201 | proj = pc.getProjection();
|
|---|
| 202 | } catch (Exception e) {
|
|---|
| 203 | String cause = e.getMessage();
|
|---|
| 204 | Main.warn("Unable to get projection "+code+" with "+pc + (cause != null ? ". "+cause : ""));
|
|---|
| 205 | }
|
|---|
| 206 | }
|
|---|
| 207 | if (proj == null) {
|
|---|
| 208 | Pair<String, String> pair = inits.get(code);
|
|---|
| 209 | if (pair == null) return null;
|
|---|
| 210 | String name = pair.a;
|
|---|
| 211 | String init = pair.b;
|
|---|
| 212 | proj = new CustomProjection(name, code, init, null);
|
|---|
| 213 | }
|
|---|
| 214 | projectionsByCode_cache.put(code, proj);
|
|---|
| 215 | return proj;
|
|---|
| 216 | }
|
|---|
| 217 |
|
|---|
| 218 | public static Collection<String> getAllProjectionCodes() {
|
|---|
| 219 | return Collections.unmodifiableCollection(allCodes);
|
|---|
| 220 | }
|
|---|
| 221 |
|
|---|
| 222 | private static String listKeys(Map<String, ?> map) {
|
|---|
| 223 | List<String> keys = new ArrayList<>(map.keySet());
|
|---|
| 224 | Collections.sort(keys);
|
|---|
| 225 | return Utils.join(", ", keys);
|
|---|
| 226 | }
|
|---|
| 227 |
|
|---|
| 228 | /**
|
|---|
| 229 | * Replies the list of projections as string (comma separated).
|
|---|
| 230 | * @return the list of projections as string (comma separated)
|
|---|
| 231 | * @since 8533
|
|---|
| 232 | */
|
|---|
| 233 | public static String listProjs() {
|
|---|
| 234 | return listKeys(projs);
|
|---|
| 235 | }
|
|---|
| 236 |
|
|---|
| 237 | /**
|
|---|
| 238 | * Replies the list of ellipsoids as string (comma separated).
|
|---|
| 239 | * @return the list of ellipsoids as string (comma separated)
|
|---|
| 240 | * @since 8533
|
|---|
| 241 | */
|
|---|
| 242 | public static String listEllipsoids() {
|
|---|
| 243 | return listKeys(ellipsoids);
|
|---|
| 244 | }
|
|---|
| 245 |
|
|---|
| 246 | /**
|
|---|
| 247 | * Replies the list of datums as string (comma separated).
|
|---|
| 248 | * @return the list of datums as string (comma separated)
|
|---|
| 249 | * @since 8533
|
|---|
| 250 | */
|
|---|
| 251 | public static String listDatums() {
|
|---|
| 252 | return listKeys(datums);
|
|---|
| 253 | }
|
|---|
| 254 |
|
|---|
| 255 | /**
|
|---|
| 256 | * Replies the list of nadgrids as string (comma separated).
|
|---|
| 257 | * @return the list of nadgrids as string (comma separated)
|
|---|
| 258 | * @since 8533
|
|---|
| 259 | */
|
|---|
| 260 | public static String listNadgrids() {
|
|---|
| 261 | return listKeys(nadgrids);
|
|---|
| 262 | }
|
|---|
| 263 | }
|
|---|