// License: GPL. For details, see LICENSE file.
package org.openstreetmap.josm.data.projection;
import java.io.BufferedReader;
import java.io.FileInputStream;
import java.io.FileNotFoundException;
import java.io.IOException;
import java.io.InputStreamReader;
import java.nio.charset.StandardCharsets;
import java.util.Collection;
import java.util.HashMap;
import java.util.Map;
import java.util.Map.Entry;
import org.junit.Test;
import org.openstreetmap.josm.data.coor.EastNorth;
import org.openstreetmap.josm.data.coor.LatLon;
import org.openstreetmap.josm.gui.preferences.projection.ProjectionChoice;
import org.openstreetmap.josm.gui.preferences.projection.ProjectionPreference;
/**
* Test projections using reference data. (Currently provided by proj.4)
*
* The data file data_nodist/projection-reference-data.csv can be created like this:
* Fist run this file's main method to collect epsg codes and bounds data.
* Then pipe the result into test/generate-proj-data.pl.
*/
public class ProjectionRefTest {
/**
* create a list of epsg codes and bounds to be used by the perl script
* @param args
*/
public static void main(String[] args) {
Map<String, Projection> allCodes = new HashMap<>();
for (ProjectionChoice pc : ProjectionPreference.getProjectionChoices()) {
for (String code : pc.allCodes()) {
Collection<String> pref = pc.getPreferencesFromCode(code);
pc.setPreferences(pref);
Projection p = pc.getProjection();
allCodes.put(code, p);
}
}
for (Entry<String, Projection> e : allCodes.entrySet()) {
System.out.println(String.format("%s %s", e.getKey(), e.getValue().getWorldBoundsLatLon()));
}
}
@Test
public void test() throws IOException, FileNotFoundException {
try (BufferedReader in = new BufferedReader(new InputStreamReader(
new FileInputStream("data_nodist/projection-reference-data.csv"), StandardCharsets.UTF_8))) {
StringBuilder fail = new StringBuilder();
String line;
while ((line = in.readLine()) != null) {
if (line.startsWith("#")) {
continue;
}
String[] f = line.split(",");
String code = f[0];
double lat = Double.parseDouble(f[1]);
double lon = Double.parseDouble(f[2]);
double east = Double.parseDouble(f[3]);
double north = Double.parseDouble(f[4]);
Projection p = Projections.getProjectionByCode(code);
EastNorth en = p.latlon2eastNorth(new LatLon(lat, lon));
String errorEN = String.format("%s (%s): Projecting latlon(%s,%s):%n" +
" expected: eastnorth(%s,%s),%n" +
" but got: eastnorth(%s,%s)!%n",
p.toString(), code, lat, lon, east, north, en.east(), en.north());
final double EPSILON_EN = SwissGridTest.SWISS_EPSG_CODE.equals(code)
? SwissGridTest.EPSILON_APPROX
: 1e-3; // 1 mm accuracy
if (Math.abs(east - en.east()) > EPSILON_EN || Math.abs(north - en.north()) > EPSILON_EN) {
fail.append(errorEN);
}
LatLon ll = p.eastNorth2latlon(new EastNorth(east, north));
String errorLL = String.format("%s (%s): Inverse projecting eastnorth(%s,%s):%n" +
" expected: latlon(%s,%s),%n" +
" but got: latlon(%s,%s)!%n",
p.toString(), code, east, north, lat, lon, ll.lat(), ll.lon());
final double EPSILON_LL = Math.toDegrees(EPSILON_EN / 6378137); // 1 mm accuracy (or better)
if (Math.abs(lat - ll.lat()) > EPSILON_LL || Math.abs(lon - ll.lon()) > EPSILON_LL) {
if (!("yes".equals(System.getProperty("suppressPermanentFailure")) && code.equals("EPSG:21781"))) {
fail.append(errorLL);
}
}
}
if (fail.length() > 0) {
System.err.println(fail.toString());
throw new AssertionError(fail.toString());
}
}
}
}