Revision: 5017
http://sourceforge.net/p/jump-pilot/code/5017
Author: michaudm
Date: 2016-09-16 18:58:48 +0000 (Fri, 16 Sep 2016)
Log Message:
-----------
Switch to new ProjUtils, rename old ProjUtils to ProjUtilsOld
Modified Paths:
--------------
core/trunk/src/com/vividsolutions/jump/workbench/imagery/ReferencedImageFactoryFileLayerLoader.java
core/trunk/src/org/openjump/core/ccordsys/utils/SRSInfo.java
core/trunk/src/org/openjump/core/rasterimage/TiffTags.java
core/trunk/src/org/openjump/core/ui/io/file/DataSourceFileLayerLoader.java
core/trunk/src/org/openjump/core/ui/plugin/layer/NewLayerPropertiesPlugIn.java
core/trunk/src/org/openjump/core/ui/plugin/raster/RasterImageLayerPropertiesPlugIn.java
Added Paths:
-----------
core/trunk/src/jumptest/ccordsys/ProjUtilsOldTest.java
core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils.java
core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtilsOld.java
Removed Paths:
-------------
core/trunk/src/jumptest/ccordsys/ProjUtilsTest.java
core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils.java
core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils2.java
Modified:
core/trunk/src/com/vividsolutions/jump/workbench/imagery/ReferencedImageFactoryFileLayerLoader.java
===================================================================
---
core/trunk/src/com/vividsolutions/jump/workbench/imagery/ReferencedImageFactoryFileLayerLoader.java
2016-09-11 13:44:55 UTC (rev 5016)
+++
core/trunk/src/com/vividsolutions/jump/workbench/imagery/ReferencedImageFactoryFileLayerLoader.java
2016-09-16 18:58:48 UTC (rev 5017)
@@ -121,7 +121,12 @@
.getStyle(SRIDStyle.class);
int prjSRID = 0;
try {
- prjSRID = ProjUtils.SRID(layer);
+ String prjSRIDString =
ProjUtils.getSRSInfoFromLayerStyleOrSource(layer).getCode();
+ if (prjSRIDString.matches("\\d+")) {
+ prjSRID = Integer.parseInt(prjSRIDString);
+ } else {
+ prjSRID = 0;
+ }
} catch (Exception e) {
prjSRID = 0;
}
Copied: core/trunk/src/jumptest/ccordsys/ProjUtilsOldTest.java (from rev 5014,
core/trunk/src/jumptest/ccordsys/ProjUtilsTest.java)
===================================================================
--- core/trunk/src/jumptest/ccordsys/ProjUtilsOldTest.java
(rev 0)
+++ core/trunk/src/jumptest/ccordsys/ProjUtilsOldTest.java 2016-09-16
18:58:48 UTC (rev 5017)
@@ -0,0 +1,93 @@
+package jumptest.ccordsys;
+
+import junit.framework.TestCase;
+import org.openjump.core.ccordsys.utils.SRSInfo;
+import org.openjump.core.rasterimage.TiffTags;
+
+import java.io.File;
+
+/**
+ * Created by UMichael on 10/09/2016.
+ */
+public class ProjUtilsOldTest extends TestCase {
+ public static void main(String[] args) throws Exception {
+ /*
+ File src1 = new File("D:\\DATA\\FORMATS\\GeoTIFF\\gdal_eg\\cea.tif");
+ File src2 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\GeogToWGS84GeoKey\\GeogToWGS84GeoKey5.tif");
+ File src3 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\intergraph\\albers27.tif");
+ File src4 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\intergraph\\lamb_conf.tif");
+ File src5 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\intergraph\\utm.tif");
+ File src6 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\made_up\\bogota.tif");
+ File src7 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\made_up\\ilatlon_float.tif");
+ File src8 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\made_up\\lcc-datum.tif");
+ File src9 = new File("D:\\DATA\\FORMATS\\GeoTIFF\\misc\\tjpeg.tif");
+ File src10 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\other\\erdas_spnad83.tif");
+ File src11 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\pci_eg\\oblqmerc.tif");
+ File src12 = new File("D:\\DATA\\FORMATS\\GeoTIFF\\usgs\\ilatlon.tif");
+ File src13 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\usgs\\o41078a1.tif");
+ File src14 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\zi_imaging\\image0.tif");
+ File src15 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\zi_imaging\\tp_image0.tif");
+
+ SRSInfo info1 = TiffTags.readMetadata(src1).getSRSInfo();
+ System.out.printf("%s%n%n", info1);
+
+ System.out.println("");
+ SRSInfo info2 = TiffTags.readMetadata(src2).getSRSInfo();
+ System.out.printf("%s%n%n", info2);
+
+ System.out.println("");
+ SRSInfo info3 = TiffTags.readMetadata(src3).getSRSInfo();
+ System.out.printf("%s%n%n", info3);
+
+ System.out.println("");
+ SRSInfo info4 = TiffTags.readMetadata(src4).getSRSInfo();
+ System.out.printf("%s%n%n", info4);
+
+ System.out.println("");
+ SRSInfo info5 = TiffTags.readMetadata(src5).getSRSInfo();
+ System.out.printf("%s%n%n", info5);
+
+ System.out.println("");
+ SRSInfo info6 = TiffTags.readMetadata(src6).getSRSInfo();
+ System.out.printf("%s%n%n", info6);
+
+ SRSInfo info7 = TiffTags.readMetadata(src7).getSRSInfo();
+ System.out.printf("%s%n%n", info7);
+
+ SRSInfo info8 = TiffTags.readMetadata(src8).getSRSInfo();
+ System.out.printf("%s%n%n", info8);
+
+ SRSInfo info9 = TiffTags.readMetadata(src9).getSRSInfo();
+ System.out.printf("%s%n%n", info9);
+
+ SRSInfo info10 = TiffTags.readMetadata(src10).getSRSInfo();
+ System.out.printf("%s%n%n", info10);
+
+ SRSInfo info11 = TiffTags.readMetadata(src11).getSRSInfo();
+ System.out.printf("%s%n%n", info11);
+
+ SRSInfo info12 = TiffTags.readMetadata(src12).getSRSInfo();
+ System.out.printf("%s%n%n", info12);
+
+ SRSInfo info13 = TiffTags.readMetadata(src13).getSRSInfo();
+ System.out.printf("%s%n%n", info13);
+
+ SRSInfo info14 = TiffTags.readMetadata(src14).getSRSInfo();
+ System.out.printf("%s%n%n", info14);
+
+ SRSInfo info15 = TiffTags.readMetadata(src15).getSRSInfo();
+ System.out.printf("%s%n%n", info15);
+ //System.out.printf("%s%n", md.getEnvelope());
+ */
+
+ }
+
+ public void testGetSRSInfoFromGeoTiff() throws Exception {
+ File src1 = new File("D:\\DATA\\FORMATS\\GeoTIFF\\gdal_eg\\cea.tif");
+ SRSInfo info1 = TiffTags.readMetadata(src1).getSRSInfo();
+ System.out.printf("%s%n%n", info1);
+ assertEquals("cea.tiff contains unnamed NAD27", "unnamed NAD27",
info1.getDescription());
+ }
+
+
+}
Deleted: core/trunk/src/jumptest/ccordsys/ProjUtilsTest.java
===================================================================
--- core/trunk/src/jumptest/ccordsys/ProjUtilsTest.java 2016-09-11 13:44:55 UTC
(rev 5016)
+++ core/trunk/src/jumptest/ccordsys/ProjUtilsTest.java 2016-09-16 18:58:48 UTC
(rev 5017)
@@ -1,93 +0,0 @@
-package jumptest.ccordsys;
-
-import junit.framework.TestCase;
-import org.openjump.core.ccordsys.utils.SRSInfo;
-import org.openjump.core.rasterimage.TiffTags;
-
-import java.io.File;
-
-/**
- * Created by UMichael on 10/09/2016.
- */
-public class ProjUtilsTest extends TestCase {
- public static void main(String[] args) throws Exception {
- /*
- File src1 = new File("D:\\DATA\\FORMATS\\GeoTIFF\\gdal_eg\\cea.tif");
- File src2 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\GeogToWGS84GeoKey\\GeogToWGS84GeoKey5.tif");
- File src3 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\intergraph\\albers27.tif");
- File src4 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\intergraph\\lamb_conf.tif");
- File src5 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\intergraph\\utm.tif");
- File src6 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\made_up\\bogota.tif");
- File src7 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\made_up\\ilatlon_float.tif");
- File src8 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\made_up\\lcc-datum.tif");
- File src9 = new File("D:\\DATA\\FORMATS\\GeoTIFF\\misc\\tjpeg.tif");
- File src10 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\other\\erdas_spnad83.tif");
- File src11 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\pci_eg\\oblqmerc.tif");
- File src12 = new File("D:\\DATA\\FORMATS\\GeoTIFF\\usgs\\ilatlon.tif");
- File src13 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\usgs\\o41078a1.tif");
- File src14 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\zi_imaging\\image0.tif");
- File src15 = new
File("D:\\DATA\\FORMATS\\GeoTIFF\\zi_imaging\\tp_image0.tif");
-
- SRSInfo info1 = TiffTags.readMetadata(src1).getSRSInfo();
- System.out.printf("%s%n%n", info1);
-
- System.out.println("");
- SRSInfo info2 = TiffTags.readMetadata(src2).getSRSInfo();
- System.out.printf("%s%n%n", info2);
-
- System.out.println("");
- SRSInfo info3 = TiffTags.readMetadata(src3).getSRSInfo();
- System.out.printf("%s%n%n", info3);
-
- System.out.println("");
- SRSInfo info4 = TiffTags.readMetadata(src4).getSRSInfo();
- System.out.printf("%s%n%n", info4);
-
- System.out.println("");
- SRSInfo info5 = TiffTags.readMetadata(src5).getSRSInfo();
- System.out.printf("%s%n%n", info5);
-
- System.out.println("");
- SRSInfo info6 = TiffTags.readMetadata(src6).getSRSInfo();
- System.out.printf("%s%n%n", info6);
-
- SRSInfo info7 = TiffTags.readMetadata(src7).getSRSInfo();
- System.out.printf("%s%n%n", info7);
-
- SRSInfo info8 = TiffTags.readMetadata(src8).getSRSInfo();
- System.out.printf("%s%n%n", info8);
-
- SRSInfo info9 = TiffTags.readMetadata(src9).getSRSInfo();
- System.out.printf("%s%n%n", info9);
-
- SRSInfo info10 = TiffTags.readMetadata(src10).getSRSInfo();
- System.out.printf("%s%n%n", info10);
-
- SRSInfo info11 = TiffTags.readMetadata(src11).getSRSInfo();
- System.out.printf("%s%n%n", info11);
-
- SRSInfo info12 = TiffTags.readMetadata(src12).getSRSInfo();
- System.out.printf("%s%n%n", info12);
-
- SRSInfo info13 = TiffTags.readMetadata(src13).getSRSInfo();
- System.out.printf("%s%n%n", info13);
-
- SRSInfo info14 = TiffTags.readMetadata(src14).getSRSInfo();
- System.out.printf("%s%n%n", info14);
-
- SRSInfo info15 = TiffTags.readMetadata(src15).getSRSInfo();
- System.out.printf("%s%n%n", info15);
- //System.out.printf("%s%n", md.getEnvelope());
- */
-
- }
-
- public void testGetSRSInfoFromGeoTiff() throws Exception {
- File src1 = new File("D:\\DATA\\FORMATS\\GeoTIFF\\gdal_eg\\cea.tif");
- SRSInfo info1 = TiffTags.readMetadata(src1).getSRSInfo();
- System.out.printf("%s%n%n", info1);
- assertEquals("cea.tiff contains unnamed NAD27", "unnamed NAD27",
info1.getDescription());
- }
-
-
-}
Deleted: core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils.java
===================================================================
--- core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils.java
2016-09-11 13:44:55 UTC (rev 5016)
+++ core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils.java
2016-09-16 18:58:48 UTC (rev 5017)
@@ -1,1151 +0,0 @@
-package org.openjump.core.ccordsys.utils;
-
-import java.awt.Point;
-import java.awt.image.BufferedImage;
-import java.io.File;
-import java.io.IOException;
-import java.io.InputStream;
-import java.io.InputStreamReader;
-import java.net.URISyntaxException;
-import java.util.ArrayList;
-import java.util.HashMap;
-import java.util.Iterator;
-import java.util.List;
-import java.util.Scanner;
-
-import javax.imageio.ImageIO;
-
-import org.apache.commons.imaging.formats.tiff.TiffField;
-import org.apache.commons.imaging.formats.tiff.TiffImageMetadata;
-import org.apache.commons.imaging.formats.tiff.TiffImageParser;
-import org.apache.commons.io.FilenameUtils;
-import org.openjump.core.ccordsys.srid.SRIDStyle;
-import org.openjump.core.rasterimage.GeoTiffConstants;
-
-import com.sun.media.jai.codec.FileSeekableStream;
-import com.sun.media.jai.codec.TIFFDirectory;
-import com.sun.media.jai.codec.TIFFField;
-import com.vividsolutions.jts.geom.Coordinate;
-import com.vividsolutions.jts.geom.Envelope;
-import com.vividsolutions.jump.I18N;
-import com.vividsolutions.jump.feature.Feature;
-import com.vividsolutions.jump.feature.FeatureCollection;
-import com.vividsolutions.jump.io.datasource.DataSourceQuery;
-import com.vividsolutions.jump.util.FileUtil;
-import com.vividsolutions.jump.workbench.Logger;
-import com.vividsolutions.jump.workbench.imagery.ImageryLayerDataset;
-import com.vividsolutions.jump.workbench.imagery.ReferencedImageStyle;
-import com.vividsolutions.jump.workbench.model.Layer;
-import
com.vividsolutions.jump.workbench.ui.plugin.datastore.DataStoreDataSource;
-
-public class ProjUtils {
-
- private static final String PROJECTION_UNSPECIFIED = I18N
-
.get("org.openjump.core.ui.plugin.raster.RasterImageLayerPropertiesPlugIn.unknown_projection");
- private static final String USER_DEFINED = I18N
-
.get("org.openjump.core.ui.plugin.layer.LayerPropertiesPlugIn.User_defined");
- private static final String NOT_RECOGNIZED = I18N
-
.get("org.openjump.core.ui.plugin.layer.LayerPropertiesPlugIn.Not_recognized");
-
- /*
- * Giuseppe Aruta [23_3_2016] This class is used to recognize file
- * projection. There are different methods A) a method to decode projection
- * information from GeoTIFF metadata. B) a method to decode projection info
- * from auxiliary files (.proj and .aux.xml).
- * http://landsathandbook.gsfc.nasa.gov/pdfs/geotiff_spec.pdf
- * http://www.sno.phy.queensu.ca/~phil/exiftool/TagNames/GeoTiff.html
- * http://www.remotesensing.org/geotiff/spec/geotiff6.html
- *
- * the following datasets have been used to test this library: QGIS
- * (http://qgis.org/downloads/data/) OSGEO data samples (
- * http://download.osgeo.org/geotiff/samples/)
- */
-
- /**
- * Test method - to get GeoTiff envelope in case of no srid def.
- *
- * @param fileSourcePath
- * @return envelope
- */
- public static Envelope GeoTiffEnvelope(String fileSourcePath)
- throws Exception {
- BufferedImage image = ImageIO.read(new File(fileSourcePath));
- Point imageDimensions = new Point(image.getWidth(), image.getHeight());
- Envelope env = null;
- Coordinate tiePoint = null, pixelOffset = null, pixelScale = null;
- double[] doubles;
- FileSeekableStream fileSeekableStream = new FileSeekableStream(
- fileSourcePath);
- TIFFDirectory tiffDirectory = new TIFFDirectory(fileSeekableStream, 0);
-
- TIFFField[] availTags = tiffDirectory.getFields();
-
- for (TIFFField availTag : availTags) {
- if (availTag.getTag() == GeoTiffConstants.ModelTiepointTag) {
- doubles = availTag.getAsDoubles();
- if (doubles.length != 6) {
- throw new Exception(
- "unsupported value for ModelTiepointTag ("
- + GeoTiffConstants.ModelTiepointTag + ")");
- }
- if (doubles[0] != 0 || doubles[1] != 0 || doubles[2] != 0) {
- if (doubles[2] == 0)
- pixelOffset = new Coordinate(doubles[0], doubles[1]);
- else
- pixelOffset = new Coordinate(doubles[0], doubles[1],
- doubles[2]);
- }
- if (doubles[5] == 0)
- tiePoint = new Coordinate(doubles[3], doubles[4]);
- else
- tiePoint = new Coordinate(doubles[3], doubles[4],
- doubles[5]);
-
- } else if (availTag.getTag() ==
GeoTiffConstants.ModelPixelScaleTag) {
- // Karteneinheiten pro pixel x bzw. y
- doubles = availTag.getAsDoubles();
- if (doubles[2] == 0)
- pixelScale = new Coordinate(doubles[0], doubles[1]);
- else
- pixelScale = new Coordinate(doubles[0], doubles[1],
- doubles[2]);
- } else {
- }
- }
-
- fileSeekableStream.close();
-
- if (tiePoint != null && pixelScale != null) {
- Coordinate upperLeft;
- Coordinate lowerRight;
-
- if (pixelOffset == null) {
- upperLeft = tiePoint;
- } else {
- upperLeft = new Coordinate(tiePoint.x
- - (pixelOffset.x * pixelScale.x), tiePoint.y
- - (pixelOffset.y * pixelScale.y));
- }
-
- lowerRight = new Coordinate(upperLeft.x
- + (imageDimensions.x * pixelScale.x), upperLeft.y
- - (imageDimensions.y * pixelScale.y));
-
- env = new Envelope(upperLeft, lowerRight);
- }
- return env;
- }
-
- /**
- * Test method to read the entire GeoKeyDirectoryTag, only used to study
- * GeoKeyDirectoryTag structure
- *
- * @param fileSourcePath
- * @return GeoKeyDirectoryTag
- * @throws IOException
- * @throws URISyntaxException
- */
- @SuppressWarnings("static-access")
- public static String readGeoTiffGeoKeyDirectoryTag(String fileSourcePath)
- throws IOException, URISyntaxException {
- String prjname = "";
- File tiffFile = new File(fileSourcePath);
- try {
- TiffImageParser parser = new TiffImageParser();
- TiffImageMetadata metadata = (TiffImageMetadata) parser
- .getMetadata(tiffFile);
- if (metadata != null) {
- List<TiffField> tiffFields = metadata.getAllFields();
- GeoTiffConstants constants = new GeoTiffConstants();
- for (TiffField tiffField : tiffFields) {
- if (tiffField.getTag() == constants.GeoKeyDirectoryTag) {
- String GeoDirTag = tiffField.getValueDescription();
- prjname = GeoDirTag;
- }
- }
- } else {
- prjname = "GeoKeyDirectoryTagis empty";
- }
- } catch (Exception ex) {
- prjname = PROJECTION_UNSPECIFIED;
- }
- return prjname;
- }
-
- /**
- * - Read SRS from GeoTIFF tag - This method gets projection srid code and
- * projection info from a geotiff file. It first scans GeoKeyDirectoryTag
to
- * get either geographic/geocentric (2048 - GeographicTypeGeoKey),
projected
- * (3072 - ProjectedCSTypeGeoKey) or vertical (4096 - VerticalCSTypeGeoKey)
- * info. If no key ID is identified, it scans for GeoAsciiParamsTag
- * projection definition. Last choice, it search for an auxiliary file
- *
- * @param fileSourcePath
- * . eg. "c\documents\folder\image.tif"
- * @return <String> projection srid code as string. eg "32632"
- * @throws IOException
- * @throws URISyntaxException
- */
- @SuppressWarnings("static-access")
- public static String readSRSFromGeoTiffFile(String fileSourcePath)
- throws IOException, URISyntaxException {
- String GeoDirTag = "";
- String GeoDirTag2 = "";
- String prjname = "";
- File tiffFile = new File(fileSourcePath);
- try {
- TiffImageParser parser = new TiffImageParser();
- TiffImageMetadata metadata = (TiffImageMetadata) parser
- .getMetadata(tiffFile);
- if (metadata != null) {
- List<TiffField> tiffFields = metadata.getAllFields();
- GeoTiffConstants constants = new GeoTiffConstants();
-
- String ID = "";
- int start;
- for (TiffField tiffField : tiffFields) {
- if (tiffField.getTag() == constants.GeoKeyDirectoryTag) {
- GeoDirTag = tiffField.getValueDescription();
- }
- if (tiffField.getTag() == constants.GeoAsciiParamsTag) {
- GeoDirTag2 = tiffField.getStringValue().replaceAll(
- "[\\t\\n\\r\\_\\|]", " ");
- }
- if (tiffField.getTag() == constants.ModelTiepointTag) {
- }
- }
- if (GeoDirTag.contains("3072")) {
- start = GeoDirTag.indexOf("3072");
- ID = GeoDirTag.substring(start);
- String[] parts = ID.split(",");
- String part1 = parts[3];
- if (!part1.contains("32767")) {
- prjname = getSRSFromWkt(part1.replaceAll(" ", ""));
- } else {
- prjname = "SRID: " + USER_DEFINED + " - " + GeoDirTag2;
- }
- } else if (!GeoDirTag.contains("3072")
- & GeoDirTag.contains("4096")) {
- start = GeoDirTag.indexOf("4096");
- ID = GeoDirTag.substring(start);
- String[] parts = ID.split(",");
- String part1 = parts[3];
- if (!part1.contains("32767")) {
- prjname = getSRSFromWkt(part1.replaceAll(" ", ""));
- } else {
- prjname = "SRID: " + USER_DEFINED + " - " + GeoDirTag2;
- }
- } else if (!GeoDirTag.contains("3072")
- & !GeoDirTag.contains("4096")
- & GeoDirTag.contains("2048")) {
- start = GeoDirTag.indexOf("2048");
- ID = GeoDirTag.substring(start);
- String[] parts = ID.split(",");
- String part1 = parts[3];
- if (!part1.contains("32767")) {
- prjname = getSRSFromWkt(part1.replaceAll(" ", ""));
- } else {
- prjname = "SRID: " + USER_DEFINED + " - " + GeoDirTag2;
- }
-
- } else if (!GeoDirTag.contains("4096")
- & !GeoDirTag.contains("3072")
- & !GeoDirTag.contains("2048")) {
- if (!GeoDirTag2.isEmpty())
- // It gets "0" which is a non defined projection
- prjname = "SRID: " + NOT_RECOGNIZED + " - "
- + GeoDirTag2;
- else
- prjname = readSRSFromAuxiliaryFile(fileSourcePath);
-
- }
- }
-
- } catch (Exception ex) {
- prjname = PROJECTION_UNSPECIFIED;
- }
- return prjname;
- }
-
- /**
- * - Read SRS from auxiliary file - Method to get a SRS (SRID code + SRID
- * definition) scanning the aux projection file (AUX.XML or PRJ file) for a
- * search string (SRID code or SRID definition). It scans into the registry
- * file (srid.txt) to find a correspondence between the search string and
- * lines of the srid.txt. If the source string corresponds as substring to
a
- * line, it returns the complete line as string. For instance, search
- * strings like "NAD83 UTM zone 10N" or "26910" both return
- * "SRID:26910 - NAD83 UTM zone 10N".
- *
- * @param fileSourcePath auxiliary file path
- * @return <String> SRID and Projection definition
- * @throws URISyntaxException
- * @throws IOException
- */
-
- public static String readSRSFromAuxiliaryFile(String fileSourcePath)
- throws URISyntaxException, IOException {
- InputStream is = ProjUtils.class.getResourceAsStream("srid.txt");
- InputStreamReader isr = new InputStreamReader(is);
- String projectSourceFilePrj = "";
- String projectSourceRFilePrj = "";
- String projectSourceRFileAux = "";
- String textProj = "";
- String prjname = "";
- String SRSDef = PROJECTION_UNSPECIFIED;
- Scanner scanner;
- // --- it reads an auxiliary file and decode a possible proj
- // --- definition to a simple string. Ex. "WGS 84 UTM Zone 32"
- int pos = fileSourcePath.lastIndexOf('.');
- // .shp, .dxf, .asc, .flt files
- projectSourceFilePrj = fileSourcePath.substring(0, pos) + ".prj";
- // image files
- projectSourceRFilePrj = fileSourcePath + ".prj";
- projectSourceRFileAux = fileSourcePath + ".aux.xml";
- List<String> fileList = new ArrayList<String>();
- fileList.add(projectSourceFilePrj);
- fileList.add(projectSourceRFilePrj);
- fileList.add(projectSourceRFileAux);
-
- if (fileList.isEmpty())
- SRSDef = PROJECTION_UNSPECIFIED;
-
- String type = FilenameUtils.getExtension(fileSourcePath).toUpperCase();
-
- if (type.equals("SHP") || type.equals("DXF") || type.equals("ASC")
- || type.equals("FLT") || type.equals("ADF")
- || type.equals("GRD") || type.equals("BIL")) {
- if (new File(projectSourceFilePrj).exists()) {
- scanner = new Scanner(new File(projectSourceFilePrj));
- textProj = scanner.nextLine();
- scanner.close();
- prjname = decodeProjDescription(textProj);
- } else {
- SRSDef = PROJECTION_UNSPECIFIED;
- }
- }
-
- else {
- if ((new File(projectSourceRFileAux).exists())
- & (new File(projectSourceRFilePrj).exists())) {
- scanner = new Scanner(new File(projectSourceRFileAux));
- textProj = scanner.useDelimiter("\\A").next();
- // scanner.close();
- if (textProj.contains("<WKT>") || textProj.contains("<SRS>")) {
- prjname = decodeProjDescription(textProj);
- } else {
- scanner = new Scanner(new File(projectSourceRFilePrj));
- textProj = scanner.nextLine();
- scanner.close();
- prjname = decodeProjDescription(textProj);
- }
- }
- if ((new File(projectSourceRFileAux).exists())
- & !(new File(projectSourceRFilePrj).exists())) {
- scanner = new Scanner(new File(projectSourceRFileAux));
- textProj = scanner.useDelimiter("\\A").next();
- // scanner.close();
- if (textProj.contains("<WKT>") || textProj.contains("<SRS>")) {
- prjname = decodeProjDescription(textProj);
- } else {
- SRSDef = PROJECTION_UNSPECIFIED;
- }
- } else if (!(new File(projectSourceRFileAux).exists())
- & (new File(projectSourceRFilePrj).exists())) {
- scanner = new Scanner(new File(projectSourceRFilePrj));
- textProj = scanner.nextLine();
- // scanner.close();
- prjname = decodeProjDescription(textProj);
- }
-
- else if (!(new File(projectSourceRFileAux).exists())
- & (!new File(projectSourceRFilePrj).exists())) {
- SRSDef = PROJECTION_UNSPECIFIED;
-
- }
- }
- // --- it extracts from proj register file all the info related
- // --- to the previous string (SRSDef). Ex.
- // --- "EPSG:32632 - WGS 84 UTM zone 32"
- if (!prjname.isEmpty()) {
- scanner = new Scanner(isr);
- try {
- while (scanner.hasNextLine()) {
- scanner.useDelimiter("\\n");
- String line = scanner.nextLine();
- String line2 = line.replaceAll("[\\t\\n\\r\\_]", "")
- .replace(";[", " [");
- if (line2.toLowerCase().contains(
- "<" + prjname.toLowerCase() + ">")) {
- int start = line2.indexOf('<');
- int end = line2.indexOf('>', start);
- String def = line2.substring(start + 1, end);
- try {
- int srid = Integer.parseInt(def);
- // 1) WKID <32768 or >5999999
- // will result in an AUTHORITY name of "EPSG".
- // 2) A WKID in range between 33000 and 199999
- // will result in an AUTHORITY name of "ESRI".
- // 3) A WKID in range between 200000 and 209199
- // will result in an AUTHORITY name of "CUSTOM".
- //
(http://help.arcgis.com/en/arcgisserver/10.0/apis/soap/whnjs.htm#SOAP_Geometry_FindSRByWKID.htm)
- if (srid < 32768 || srid > 5999999) {
- // EPSG code between 0 and 32767
- SRSDef = "EPSG:"
- + line2.replaceAll("[<\\>]", " ")
- .replaceAll(";", " - ");
- // ESRI codes range
- } else if (srid > 32999 && srid < 200000) {
- SRSDef = "ESRI:"
- + line2.replaceAll("[<\\>]", " ")
- .replaceAll(";", " - ");
- // Other no EPSG or ESRI codes
- } else {
- SRSDef = "CUSTOM:"
- + line2.replaceAll("[<\\>]", " ")
- .replaceAll(";", " - ");
- }
- // Non numeral SRID like IGNF
- } catch (NumberFormatException e) {
- SRSDef = "SRID:"
- + line2.replaceAll("[<\\>]", " ")
- .replaceAll(";", " - ");
- }
- break;
- } else {
- // --- If no SRSDef is recognized into the register, it
- // --- returns a proj string into a more readable text
- SRSDef =
readableFormatWKTCode(getWktProjDefinition(textProj));
- }
- }
- scanner.close();
- } catch (Exception e) {
- SRSDef = PROJECTION_UNSPECIFIED;
- }
- } else {
- SRSDef = PROJECTION_UNSPECIFIED;
- }
- return SRSDef;
- }
-
- /**
- * Method to show an OGC WKT string in a more readable style
- *
- * @param WKT OGC WKT from auxiliary proj file
- * @return Readable string
- */
- public static String readableFormatWKTCode(String WKT) {
- String HROGC = "";
- // String add_spaces = String.format("%" + count_add++ + "s", "");
- HROGC = WKT.replace(",GEOGCS", ",<br>" + "GEOCS")
- .replace(",DATUM", ",<br>" + "DATUM")
- .replace(",SPHEROID", ",<br>" + "SPEROID")
- .replace("],", "],<br>");
- return HROGC;
-
- }
-
- /**
- * Decode a OGC string to get a unique SRS string definition. This method
is
- * able to understand some WKT common aliases, like OGC WKT and ESRI
WKTode.
- * For instance: "WGS 84 / UTM zone 32", "WGS 1984 UTM zone 32" and
- * "WGS_84_UTM_Zone_32" are converted to the same string
- * "WGS 84 UTM zone 32"
- *
- * @param textProj
- * <String> - OGC/ESRI/other WKT code
- * @return <String> - SRS definition
- */
- private static String decodeProjDescription(String textProj) {
- String prjname = "";
- try {
- // Workaround if aux.xml has been download from web.
- // convert HTML quotes ["] to ["]
- textProj = textProj.replaceAll(""", "\"");
- int start = textProj.indexOf("[\"");
- int end = textProj.indexOf("\",", start);
- prjname = textProj.substring(start + 2, end);
- // The following set of replacements allows to "harmonize" OGC,
ESRI
- // and
- // few other WKT projection definitions
- prjname = prjname.replaceAll("_", " ").replace(" / ", " ")
- .replaceAll("\\bft US\\b", "(ftUS)")
- .replaceAll("\\bftUS\\b", "(ftUS)")
- .replaceAll("\\bft\\b", "(ft)").replaceAll("feet", "ft")
- .replaceAll("WGS 1984", "WGS 84")
- .replaceAll("NAD 1983 UTM", "NAD83 UTM")
- .replaceAll("HARN", "(HARN)")
- .replaceAll("\\bCSRS98\\b", "(CSRS98)")
- .replaceAll("CSRS", "(CSRS)")
- .replaceAll("\\bNSRS2007\\b", "(NSRS2007)")
- .replaceAll("\\bNAD27_76\\b", "NAD27(76)")
- .replaceAll("\\bCGQ77\\b", " (CGQ77)")
- .replaceAll("\\bED77\\b", "(ED77)")
- .replaceAll("\\b1942 83\\b", "1942(83)")
- .replaceAll("\\b1942 58\\b", "1942(58)")
- .replaceAll("\\bSegara Jakarta\\b", "Segara (Jakarta)")
- .replaceAll("\\bRome\\b", "(Rome)")
- .replaceAll("\\bParis\\b", "(Paris)")
- .replaceAll("\\bFerro\\b", "(Ferro)");
-
- } catch (Exception ex) {
- // If there is other info than a WKT definition in the aux file
- prjname = NOT_RECOGNIZED;
- }
- return prjname;
- }
-
- /**
- * returns OGC WKT string located between projection tags (<WKT> or <SRS>)
- * in a projection auxiliary file (AUX.XML)
- *
- * @param textProj
- * string
- * @return OGC WKT string
- */
- private static String getWktProjDefinition(String textProj) {
- String prjname = "";
- try {
- if (textProj.contains("<WKT>")) {
- int start = textProj.indexOf("<WKT>");
- int end = textProj.indexOf("</WKT>", start);
- prjname = textProj.substring(start, end);
- } else if (textProj.contains("<SRS>")) {
- int start = textProj.indexOf("<SRS>");
- int end = textProj.indexOf("</SRS>", start);
- prjname = textProj.substring(start, end);
- } else
- prjname = textProj;
- } catch (Exception ex) {
- prjname = textProj;
- }
- return prjname;
- }
-
- /**
- * It returns the path name of the auxiliary file (AUX.XML or PRJ file)
- * where a projection code is located
- *
- * @param fileSourcePath auxiliary file path
- * @return <String> path name of projection auxiliary file
- * @throws IOException
- */
- public static String getAuxiliaryProjFilePath(String fileSourcePath)
- throws IOException {
- String projectSourceFilePrj = "";
- String projectSourceRFilePrj = "";
- String projectSourceRFileAux = "";
- String textProj = "";
- String filename = "";
- Scanner scanner;
- int pos = fileSourcePath.lastIndexOf('.');
- projectSourceFilePrj = fileSourcePath.substring(0, pos) + ".prj";
- projectSourceRFileAux = fileSourcePath + ".aux.xml";
- projectSourceRFilePrj = fileSourcePath + ".prj";
- String type = FilenameUtils.getExtension(fileSourcePath).toUpperCase();
-
- if (type.equals("SHP") || type.equals("DXF") || type.equals("ASC")
- || type.equals("FLT") || type.equals("ADF")
- || type.equals("GRD") || type.equals("BIL")) {
- if (new File(projectSourceFilePrj).exists()) {
- filename = projectSourceFilePrj;
- } else {
- filename = "";
- }
- } else {
-
- if ((new File(projectSourceRFileAux).exists())
- & (new File(projectSourceRFilePrj).exists())) {
- scanner = new Scanner(new File(projectSourceRFileAux));
- textProj = scanner.useDelimiter("\\A").next();
- scanner.close();
- if (textProj.contains("<WKT>") || textProj.contains("<SRS>")) {
- filename = projectSourceRFileAux;
- } else {
- filename = projectSourceRFilePrj;
- }
- } else if ((new File(projectSourceRFileAux).exists())
- & !(new File(projectSourceRFilePrj).exists())) {
- scanner = new Scanner(new File(projectSourceRFileAux));
- textProj = scanner.useDelimiter("\\A").next();
- scanner.close();
- if (textProj.contains("<WKT>") || textProj.contains("<SRS>")) {
- filename = projectSourceRFileAux;
- } else {
- filename = "";
- }
- } else if (!(new File(projectSourceRFileAux).exists())
- & (new File(projectSourceRFilePrj).exists())) {
- filename = projectSourceRFilePrj;
- } else if (!(new File(projectSourceRFileAux).exists())
- & !(new File(projectSourceRFilePrj).exists())) {
- filename = "";
- }
- }
- return filename;
- }
-
- /**
- * Method to get a SRS (SRID code + SRID definition) using a search string.
- * It scans into the srid list (srid.txt) to find a correspondence between
- * the search string and lines of the srid.txt. If the source string
- * corresponds as substring to a line, it returns the complete line as
- * string.The code is biunivocal: it can use as searchQuerry either SRID
- * code ("26910" or Project definition ("NAD83 UTM zone 10N"). For
instance,
- * searche querries like "NAD83 UTM zone 10N" or "26910" both return
- * "SRID:26910 - NAD83 UTM zone 10N".
- *
- * @param searchQuery search quaery
- * @return <String> SRID and Projection definition
- * @throws URISyntaxException
- * @throws IOException
- */
-
- public static String getSRSFromWkt(String searchQuery)
- throws URISyntaxException, IOException {
-
- InputStream is = ProjUtils.class.getResourceAsStream("srid.txt");
- InputStreamReader isr = new InputStreamReader(is);
- String SRSDef = "";
- Scanner scanner = null;
- try {
- scanner = new Scanner(isr);
- while (scanner.hasNextLine()) {
- scanner.useDelimiter("\\n");
- String line = scanner.nextLine();
- String line2 = line.replaceAll("[\\t\\n\\r\\_]", "").replace(
- ";[", " [");
- if (line2.toLowerCase().contains(
- "<" + searchQuery.toLowerCase() + ">")) {
- int start = line2.indexOf('<');
- int end = line2.indexOf('>', start);
- String def = line2.substring(start + 1, end);
- try {
- int srid = Integer.parseInt(def);// get WKID
- // 1) WKID <32768 or >5999999
- // will result in an AUTHORITY name of "EPSG".
- // 2) A WKID in range between 33000 and 199999
- // will result in an AUTHORITY name of "ESRI".
- // 3) A WKID in range between 200000 and 209199
- // will result in an AUTHORITY name of "CUSTOM".
- //
(http://help.arcgis.com/en/arcgisserver/10.0/apis/soap/whnjs.htm#SOAP_Geometry_FindSRByWKID.htm)
- if (srid < 32768 || srid > 5999999) {
- // EPSG code between 0 and 32767
- SRSDef = "EPSG:"
- + line2.replaceAll("[<\\>]", " ")
- .replaceAll(";", " - ");
- // ESRI codes range
- } else if (srid > 32999 && srid < 200000) {
- SRSDef = "ESRI:"
- + line2.replaceAll("[<\\>]", " ")
- .replaceAll(";", " - ");
- // Other no EPSG or ESRI codes
- } else {
- SRSDef = "CUSTOM:"
- + line2.replaceAll("[<\\>]", " ")
- .replaceAll(";", " - ");
- }
- // Non numeral SRID like IGNF
- } catch (NumberFormatException e) {
- SRSDef = "SRID:"
- + line2.replaceAll("[<\\>]", " ").replaceAll(
- ";", " - ");
- }
- break;
- } else {
-
- SRSDef = searchQuery;
- }
- }
- } finally {
- try {
- if (scanner != null)
- scanner.close();
- } catch (Exception e) {
- System.err.println("Exception while closing scanner "
- + e.toString());
- }
- }
- return SRSDef;
- }
-
- /**
- * Check if selected file is a GeoTIFF. This java code comes from Deegree
- * project org.deegree.tools.raster.MergeRaster
- * (https://github.com/camptocamp
- * /secureOWS/blob/master/owsproxyserver/src/org
- * /deegree/tools/raster/MergeRaster.java)
- *
- * @param fileSourcePath
- * @return true
- * @throws IOException
- */
-
- public static boolean isGeoTIFF(String fileSourcePath) throws IOException {
- FileSeekableStream fileSeekableStream = new FileSeekableStream(
- fileSourcePath);
- TIFFDirectory tifDir = new TIFFDirectory(fileSeekableStream, 0);
- // definition of a geotiff
- if (tifDir.getField(GeoTiffConstants.ModelPixelScaleTag) == null
- && tifDir.getField(GeoTiffConstants.ModelTransformationTag) ==
null
- && tifDir.getField(GeoTiffConstants.ModelTiepointTag) == null
- && tifDir.getField(GeoTiffConstants.GeoKeyDirectoryTag) == null
- && tifDir.getField(GeoTiffConstants.GeoDoubleParamsTag) == null
- && tifDir.getField(GeoTiffConstants.GeoAsciiParamsTag) ==
null) {
- return false;
- } else {
- // is a geotiff and possibly might need to be treated as raw data
- TIFFField bitsPerSample = tifDir.getField(258);
- if (bitsPerSample != null) {
- int samples = bitsPerSample.getAsInt(0);
- if (samples == 16)
- new Integer(16);
- }
- // check the EPSG number
- TIFFField ff =
tifDir.getField(GeoTiffConstants.GeoKeyDirectoryTag);
- if (ff == null) {
- return false;
- }
- char[] ch = ff.getAsChars();
- // resulting HashMap, containing the key and the array of values
- HashMap<Integer, int[]> geoKeyDirectoryTag = new HashMap<Integer,
int[]>(
- ff.getCount() / 4);
- // array of values. size is 4-1.
- int keydirversion, keyrevision, minorrevision, numberofkeys = -99;
- for (int i = 0; i < ch.length; i = i + 4) {
- int[] keys = new int[3];
- keydirversion = ch[i];
- keyrevision = ch[i + 1];
- minorrevision = ch[i + 2];
- numberofkeys = ch[i + 3];
- keys[0] = keyrevision;
- keys[1] = minorrevision;
- keys[2] = numberofkeys;
- geoKeyDirectoryTag.put(new Integer(keydirversion), keys);
- }
- int[] content = new int[3];
- if (geoKeyDirectoryTag.containsKey(new Integer(
- GeoTiffConstants.ModelTiepointTag))) {
- content = (int[]) geoKeyDirectoryTag.get(new Integer(
- GeoTiffConstants.ModelTiepointTag));
- // TIFFTagLocation
- if (content[0] == 0) {
- // return Value_Offset key = content[2];
- } else {
- // TODO other TIFFTagLocation that GeoKeyDirectoryTag
- }
- } else {
- Logger.warn("Can't check EPSG codes, make sure it is ok!");
- }
- return true;
- }
- }
-
- /**
- * - Read SRID from GeoTIFF tag - This method gets projection srid code
from
- * a geotiff file. It first scans GeoKeyDirectoryTag to get either
- * geographic/geocentric (2048 - GeographicTypeGeoKey), projected (3072 -
- * ProjectedCSTypeGeoKey) or vertical (4096 - VerticalCSTypeGeoKey) info.
If
- * no key ID is identified, it scans for GeoAsciiParamsTag projection
- * definition. Last choice, it search for an auxiliary file
- *
- * @param fileSourcePath
- * . eg. "c\documents\folder\image.tif"
- * @return <String> projection srid code as string. eg "32632"
- * @throws IOException
- * @throws URISyntaxException
- */
- @SuppressWarnings("static-access")
- public static String readSRIDFromGeoTiffFile(String fileSourcePath)
- throws IOException, URISyntaxException {
- String GeoDirTag = "";
- String GeoDirTag2 = "";
- String prjname = "";
- File tiffFile = new File(fileSourcePath);
- try {
- TiffImageParser parser = new TiffImageParser();
- TiffImageMetadata metadata = (TiffImageMetadata) parser
- .getMetadata(tiffFile);
- if (metadata != null) {
- List<TiffField> tiffFields = metadata.getAllFields();
- GeoTiffConstants constants = new GeoTiffConstants();
-
- String ID = "";
- int start;
- for (TiffField tiffField : tiffFields) {
- if (tiffField.getTag() == constants.GeoKeyDirectoryTag) {
- GeoDirTag = tiffField.getValueDescription();
- }
- if (tiffField.getTag() == constants.GeoAsciiParamsTag) {
- GeoDirTag2 = tiffField.getStringValue().replaceAll(
- "[\\t\\n\\r\\_\\|]", " ");
- }
- if (tiffField.getTag() == constants.ModelTiepointTag) {
- }
- }
- if (GeoDirTag.contains("3072")) {
- start = GeoDirTag.indexOf("3072");
- ID = GeoDirTag.substring(start);
- String[] parts = ID.split(",");
- String part1 = parts[3];
- if (!part1.contains("32767")) {
- prjname = part1.replaceAll(" ", "");
- } else {
- prjname = "0";
- }
- } else if (!GeoDirTag.contains("3072")
- & GeoDirTag.contains("4096")) {
- start = GeoDirTag.indexOf("4096");
- ID = GeoDirTag.substring(start);
- String[] parts = ID.split(",");
- String part1 = parts[3];
- if (!part1.contains("32767")) {
- prjname = part1.replaceAll(" ", "");
- } else {
- prjname = "0";
- }
- } else if (!GeoDirTag.contains("3072")
- & !GeoDirTag.contains("4096")
- & GeoDirTag.contains("2048")) {
- start = GeoDirTag.indexOf("2048");
- ID = GeoDirTag.substring(start);
- String[] parts = ID.split(",");
- String part1 = parts[3];
- if (!part1.contains("32767")) {
- prjname = part1.replaceAll(" ", "");
- } else {
- prjname = "0";
- }
-
- } else if (!GeoDirTag.contains("4096")
- & !GeoDirTag.contains("3072")
- & !GeoDirTag.contains("2048")) {
- if (!GeoDirTag2.isEmpty())
- // It gets "0" which is a non defined projection
- prjname = "0";
- else
- prjname = readSRIDFromAuxiliaryFile(fileSourcePath);
-
- }
- }
-
- } catch (Exception ex) {
- prjname = PROJECTION_UNSPECIFIED;
- }
- return prjname;
- }
-
- /**
- * - Read SRID from auxiliary file - Method to get SRID code from auxiliary
- * projection file (AUX.XML or PRJ file). It scans into the registry file
- * (srid.txt) to find a correspondence between the search string (auxiliary
- * layer) and lines of the srid.txt. If the source string corresponds as
- * substring to a line, it returns its SRID. For instance, search strings
- * like "NAD83 UTM zone 10N" returns "26910".
- *
- * @param fileSourcePath auxiliary file path
- * @return SRID as String
- * @throws URISyntaxException
- * @throws IOException
- */
-
- public static String readSRIDFromAuxiliaryFile(String fileSourcePath)
- throws URISyntaxException, IOException {
- InputStream is = ProjUtils.class.getResourceAsStream("srid.txt");
- InputStreamReader isr = new InputStreamReader(is);
- String projectSourceFilePrj = "";
- String projectSourceRFilePrj = "";
- String projectSourceRFileAux = "";
- String textProj = "";
- String prjname = "";
- String SRSDef = PROJECTION_UNSPECIFIED;
- Scanner scanner;
- // --- it reads an auxiliary file and decode a possible proj
- // --- definition to a simple string. Ex. "WGS 84 UTM Zone 32"
- int pos = fileSourcePath.lastIndexOf('.');
- // .shp, .dxf, .asc, .flt files
- projectSourceFilePrj = fileSourcePath.substring(0, pos) + ".prj";
- // image files
- projectSourceRFilePrj = fileSourcePath + ".prj";
- projectSourceRFileAux = fileSourcePath + ".aux.xml";
- List<String> fileList = new ArrayList<String>();
- fileList.add(projectSourceFilePrj);
- fileList.add(projectSourceRFilePrj);
- fileList.add(projectSourceRFileAux);
-
- if (fileList.isEmpty())
- SRSDef = PROJECTION_UNSPECIFIED;
-
- String type = FilenameUtils.getExtension(fileSourcePath).toUpperCase();
-
- if (type.equals("SHP") || type.equals("DXF") || type.equals("ASC")
- || type.equals("FLT") || type.equals("ADF")
- || type.equals("GRD") || type.equals("BIL")) {
- if (new File(projectSourceFilePrj).exists()) {
- scanner = new Scanner(new File(projectSourceFilePrj));
- textProj = scanner.nextLine();
- scanner.close();
- prjname = decodeProjDescription(textProj);
- } else {
- SRSDef = PROJECTION_UNSPECIFIED;
- }
- }
-
- else {
- if ((new File(projectSourceRFileAux).exists())
- & (new File(projectSourceRFilePrj).exists())) {
- scanner = new Scanner(new File(projectSourceRFileAux));
- textProj = scanner.useDelimiter("\\A").next();
- // scanner.close();
- if (textProj.contains("<WKT>") || textProj.contains("<SRS>")) {
-
- prjname = decodeProjDescription(textProj);
- } else {
- scanner = new Scanner(new File(projectSourceRFilePrj));
- textProj = scanner.nextLine();
- scanner.close();
- prjname = decodeProjDescription(textProj);
- }
- }
- if ((new File(projectSourceRFileAux).exists())
- & !(new File(projectSourceRFilePrj).exists())) {
- scanner = new Scanner(new File(projectSourceRFileAux));
- textProj = scanner.useDelimiter("\\A").next();
- // scanner.close();
- if (textProj.contains("<WKT>") || textProj.contains("<SRS>")) {
- prjname = decodeProjDescription(textProj);
- } else {
- SRSDef = PROJECTION_UNSPECIFIED;
- }
- } else if (!(new File(projectSourceRFileAux).exists())
- & (new File(projectSourceRFilePrj).exists())) {
- scanner = new Scanner(new File(projectSourceRFilePrj));
- textProj = scanner.nextLine();
- // scanner.close();
- prjname = decodeProjDescription(textProj);
- }
-
- else if (!(new File(projectSourceRFileAux).exists())
- & (!new File(projectSourceRFilePrj).exists())) {
- SRSDef = "0";
-
- }
- }
- // --- it extracts from proj register file all the info related
- // --- to the previous string (SRSDef). Ex.
- // --- "EPSG:32632 - WGS 84 UTM zone 32"
- if (!prjname.isEmpty()) {
- scanner = new Scanner(isr);
- try {
- while (scanner.hasNextLine()) {
- scanner.useDelimiter("\\n");
- String line = scanner.nextLine();
- if (line.toLowerCase().contains(
- "<" + prjname.toLowerCase() + ">")) {
- int start = line.indexOf('<');
- int end = line.indexOf('>', start);
- String def = line.substring(start + 1, end);
- SRSDef = def;
- break;
- } else {
- // --- If no SRSDef is recognized into the register, it
- // --- returns a proj string into a more readable text
- SRSDef = "0";
- }
- }
- scanner.close();
- } catch (Exception e) {
- SRSDef = "0";
- }
- } else {
- SRSDef = "0";
- }
- return SRSDef;
- }
-
- /**
- * Method to get SRID from a layer from Style or from auxiliary file. First
- * scans SRIDStyle, than auxiliary file or GeoTIFF tag. If SRID does not
- * exist, it returns 0.
- *
- * @param layer
- * @return SRID
- * @throws Exception
- */
- public static int SRID(Layer layer) throws Exception {
- String fileSourcePath = "";
- int projection = 0;
- String extension = "";
- SRIDStyle sridStyle = (SRIDStyle) layer.getStyle(SRIDStyle.class);
- final int oldSRID = sridStyle.getSRID();
- // if (layers.length == 1) {
- // First we check if a SRID (Spatial Reference Identifier)
- // code has been recorded by OJ (eg. Spatialite)
- // it excludes 0 from the range of search as it can
- // be consider as "no SRID"
- if (oldSRID > 0) {
- projection = oldSRID;
- // If no SRID has been identified. it checks
- // projection info into external auxiliary files (.prj,
- // aux.xml) or as Geotiff tag
- } else {// Check if selected layer is related to an image file
- if (isImageFileLayer(layer)) {
- FeatureCollection featureCollection = layer
- .getFeatureCollectionWrapper();
- String sourcePathImage = null;
- for (Iterator<?> i = featureCollection.iterator();
i.hasNext();) {
- Feature feature = (Feature) i.next();
- sourcePathImage = (String) feature
- .getString(ImageryLayerDataset.ATTR_URI);
- sourcePathImage = sourcePathImage.substring(5);
- File f = new File(sourcePathImage);
- String filePath = f.getAbsolutePath();
- String filePath1 = filePath.replace("%20", " ");
- fileSourcePath = filePath1;
-
- }
- extension =
FileUtil.getExtension(fileSourcePath).toUpperCase();
- if ((extension.equals("TIF") || extension.equals("TIFF"))) {
- // If TIFF file is a geotiff, it scans into
- // embedded tag
- if (ProjUtils.isGeoTIFF(fileSourcePath)) {
-
- projection = Integer.parseInt(ProjUtils
- .readSRIDFromGeoTiffFile(fileSourcePath));
- // If the TIF file is not a GeiTIFF it looks
- // for a proj code into aux files
- } else {
- projection = Integer.parseInt(ProjUtils
- .readSRIDFromAuxiliaryFile(fileSourcePath));
- }
- } else {
- if (fileSourcePath != null) {
- projection = Integer.parseInt(ProjUtils
- .readSRIDFromAuxiliaryFile(fileSourcePath));
- }
- }// Check if source file is is a file-based vector
- } else { // Only Vector files
- if (!isDataBaseLayer(layer)) {
- DataSourceQuery dsq = layer.getDataSourceQuery();
- String sourceClass = "";
- String sourcePath = "";
- String dsqSourceClass = dsq.getDataSource().getClass()
- .getName();
- if (sourceClass.equals("")) {
- sourceClass = dsqSourceClass;
- }
- Object fnameObj = dsq.getDataSource().getProperties()
- .get("File");
- sourcePath = fnameObj.toString();
- fileSourcePath = sourcePath;
- projection = Integer.parseInt(ProjUtils
- .readSRIDFromAuxiliaryFile(fileSourcePath));
- } else {
-
- projection = 0;
- }
- }
- }
-
- return projection;
-
- }
-
- /**
- * Method to get SRID from a layer file from auxiliary files (.prj or .aux)
- * or GeoTIFFed tag. If the auxiliary file SRID does not exist, it returns
- * 0.
- *
- * @param layer
- * @return SRID
- * @throws Exception
- */
- public static int SRIDFromFile(Layer layer) throws Exception {
- String fileSourcePath = "";
- int SRID = 0;
- String extension = "";
- // First check if selected Layer has an image datasource file
- if (isImageFileLayer(layer)) {
- FeatureCollection featureCollection = layer
- .getFeatureCollectionWrapper();
- String sourcePathImage = null;
- for (Iterator<?> i = featureCollection.iterator(); i.hasNext();) {
- Feature feature = (Feature) i.next();
- sourcePathImage = (String) feature
- .getString(ImageryLayerDataset.ATTR_URI);
- sourcePathImage = sourcePathImage.substring(5);
- File f = new File(sourcePathImage);
- String filePath = f.getAbsolutePath();
- String filePath1 = filePath.replace("%20", " ");
- fileSourcePath = filePath1;
-
- }
- extension = FileUtil.getExtension(fileSourcePath).toUpperCase();
- if ((extension.equals("TIF") || extension.equals("TIFF"))) {
- // If TIFF file is a geotiff, it scans into
- // embedded tag
- if (ProjUtils.isGeoTIFF(fileSourcePath)) {
-
- SRID = Integer.parseInt(ProjUtils
- .readSRIDFromGeoTiffFile(fileSourcePath));
- // If the TIF file is not a GeiTIFF it looks
- // for a proj code into aux files
- } else {
- SRID = Integer.parseInt(ProjUtils
- .readSRIDFromAuxiliaryFile(fileSourcePath));
- }
- // For all other image file types, not TIF
- } else {
- if (fileSourcePath != null) {
- SRID = Integer.parseInt(ProjUtils
- .readSRIDFromAuxiliaryFile(fileSourcePath));
- }
- }
- // Than check if source file is is a file-based vector
- } else {
- // Only Vector files, excluding databases:
- // OpenJUMP has different approches to detect
- // SRIDs for them
- if (!isDataBaseLayer(layer)) {
- DataSourceQuery dsq = layer.getDataSourceQuery();
- String sourceClass = "";
- String sourcePath = "";
- String dsqSourceClass = dsq.getDataSource().getClass()
- .getName();
- if (sourceClass.equals("")) {
- sourceClass = dsqSourceClass;
- }
- Object fnameObj = dsq.getDataSource().getProperties()
- .get("File");
- sourcePath = fnameObj.toString();
- fileSourcePath = sourcePath;
- SRID = Integer.parseInt(ProjUtils
- .readSRIDFromAuxiliaryFile(fileSourcePath));
- } else {
- SRID = 0;
- }
- }
- return SRID;
- }
-
- // Boolean. Selected layer is related to an image file
- private static boolean isImageFileLayer(Layer layer) {
- if (layer.getStyle(ReferencedImageStyle.class) != null
- && (layer.getDescription() != null)) {
- return true;
- } else {
- return false;
- }
- }
-
- // Boolean. Selected layer is related to a database
- private static boolean isDataBaseLayer(Layer layer) {
- DataSourceQuery dsq = layer.getDataSourceQuery();
- if (dsq == null || dsq.getDataSource() instanceof DataStoreDataSource)
{
- return true;
- } else {
- return false;
- }
- }
-
-}
Copied: core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils.java (from
rev 5014, core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils2.java)
===================================================================
--- core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils.java
(rev 0)
+++ core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils.java
2016-09-16 18:58:48 UTC (rev 5017)
@@ -0,0 +1,331 @@
+package org.openjump.core.ccordsys.utils;
+
+import com.vividsolutions.jump.I18N;
+import com.vividsolutions.jump.feature.Feature;
+import com.vividsolutions.jump.feature.FeatureCollection;
+import com.vividsolutions.jump.io.datasource.DataSourceQuery;
+import com.vividsolutions.jump.util.FileUtil;
+import com.vividsolutions.jump.workbench.imagery.ImageryLayerDataset;
+import com.vividsolutions.jump.workbench.imagery.ReferencedImageStyle;
+import com.vividsolutions.jump.workbench.model.Layer;
+import
com.vividsolutions.jump.workbench.ui.plugin.datastore.DataStoreDataSource;
+import org.apache.commons.imaging.ImageReadException;
+import org.apache.commons.io.FilenameUtils;
+import org.openjump.core.ccordsys.srid.SRIDStyle;
+import org.openjump.core.rasterimage.TiffTags;
+
+
+import java.io.File;
+import java.io.IOException;
+import java.net.URISyntaxException;
+import java.util.Iterator;
+import java.util.Scanner;
+
+/**
+ * Giuseppe Aruta [23_3_2016] This class is used to recognize file
+ * projection. There are different methods A) a method to decode projection
+ * information from GeoTIFF metadata. B) a method to decode projection info
+ * from auxiliary files (.proj and .aux.xml).
+ * http://landsathandbook.gsfc.nasa.gov/pdfs/geotiff_spec.pdf
+ * http://www.sno.phy.queensu.ca/~phil/exiftool/TagNames/GeoTiff.html
+ * http://www.remotesensing.org/geotiff/spec/geotiff6.html
+ *
+ * the following datasets have been used to test this library: QGIS
+ * (http://qgis.org/downloads/data/) OSGEO data samples (
+ * http://download.osgeo.org/geotiff/samples/)
+ */
+public class ProjUtils {
+
+ private static final String PROJECTION_UNSPECIFIED = I18N
+
.get("org.openjump.core.ui.plugin.raster.RasterImageLayerPropertiesPlugIn.unknown_projection");
+ private static final String USER_DEFINED = I18N
+
.get("org.openjump.core.ui.plugin.layer.LayerPropertiesPlugIn.User_defined");
+ private static final String NOT_RECOGNIZED = I18N
+
.get("org.openjump.core.ui.plugin.layer.LayerPropertiesPlugIn.Not_recognized");
+
+
+ /**
+ * - Read SRS from GeoTIFF tag - This method gets projection srid code and
+ * projection info from a geotiff file. It first scans GeoKeyDirectoryTag
to
+ * get either geographic/geocentric (2048 - GeographicTypeGeoKey),
projected
+ * (3072 - ProjectedCSTypeGeoKey) or vertical (4096 - VerticalCSTypeGeoKey)
+ * info. If no key ID is identified, it scans for GeoAsciiParamsTag
+ * projection definition. Last choice, it search for an auxiliary file
+ *
+ * @param fileSourcePath
+ * . eg. "c\documents\folder\image.tif"
+ * @return the projection srid as a string. eg "32632"
+ * @throws IOException
+ * @throws URISyntaxException
+ */
+ @SuppressWarnings("static-access")
+ public static SRSInfo getSRSInfoFromGeoTiff(String fileSourcePath)
+ throws IOException, URISyntaxException,
TiffTags.TiffReadingException, ImageReadException {
+ return TiffTags.readMetadata(new File(fileSourcePath)).getSRSInfo();
+ }
+
+
+ /**
+ * - Read SRS from auxiliary file - Method to get a SRS (SRID code + SRID
+ * definition) scanning the aux projection file (AUX.XML or PRJ file) for a
+ * search string (SRID code or SRID definition). It scans into the registry
+ * file (srid.txt) to find a correspondence between the search string and
+ * lines of the srid.txt. If the source string corresponds as substring to
a
+ * line, it returns the complete line as string. For instance, search
+ * strings like "NAD83 UTM zone 10N" or "26910" both return
+ * "SRID:26910 - NAD83 UTM zone 10N".
+ *
+ * @param fileSourcePath auxiliary file path
+ * @return SRSInfo and Projection definition
+ * @throws URISyntaxException
+ * @throws IOException
+ */
+
+ public static SRSInfo getSRSInfoFromAuxiliaryFile(String fileSourcePath)
+ throws URISyntaxException, IOException {
+
+ // --- it reads an auxiliary file and decode a possible proj
+ // --- definition to a simple string. Ex. "WGS 84 UTM Zone 32"
+ int pos = fileSourcePath.lastIndexOf('.');
+ // .shp, .dxf, .asc, .flt files
+ String projectSourceFilePrj = fileSourcePath.substring(0, pos) +
".prj";
+ // image files
+ String projectSourceRFilePrj = fileSourcePath + ".prj";
+ String projectSourceRFileAux = fileSourcePath + ".aux.xml";
+
+ String type = FilenameUtils.getExtension(fileSourcePath).toUpperCase();
+
+ SRSInfo srsInfo = new SRSInfo();
+ String textProj = "";
+ if (type.matches("SHP|DXF|ASC|FLT|ADF|GRD|BIL")) {
+ if (new File(projectSourceFilePrj).exists()) {
+ Scanner scanner = new Scanner(new File(projectSourceFilePrj));
+ textProj = scanner.nextLine();
+ srsInfo = new SRSInfo().setSource(projectSourceFilePrj);
+ scanner.close();
+ }
+ } else if (new File(projectSourceRFileAux).exists()) {
+ Scanner scanner = new Scanner(new File(projectSourceRFileAux));
+ textProj = scanner.useDelimiter("\\A").next();
+ if (!textProj.contains("<WKT>") && !textProj.contains("<SRS>") &&
+ new File(projectSourceRFilePrj).exists()) {
+ Scanner scanner2 = new Scanner(new
File(projectSourceRFilePrj));
+ textProj = getWktProjDefinition(scanner2.nextLine());
+ srsInfo = new SRSInfo().setSource(projectSourceRFileAux);
+ scanner2.close();
+ }
+ } else if (new File(projectSourceRFilePrj).exists()) {
+ Scanner scanner = new Scanner(new File(projectSourceRFilePrj));
+ textProj = scanner.nextLine();
+ srsInfo = new SRSInfo().setSource(projectSourceRFilePrj);
+ scanner.close();
+ }
+
+ String prjName = decodeProjDescription(textProj);
+ srsInfo.setDescription(prjName);
+ srsInfo.complete();
+ return srsInfo;
+ }
+
+
+ /**
+ * Method to show an OGC WKT string in a more readable style
+ *
+ * @param WKT OGC WKT from auxiliary proj file
+ * @return Readable string
+ */
+ public static String readableFormatWKTCode(String WKT) {
+ String HROGC = "";
+ // String add_spaces = String.format("%" + count_add++ + "s", "");
+ HROGC = WKT.replace(",GEOGCS", ",<br>" + "GEOCS")
+ .replace(",DATUM", ",<br>" + "DATUM")
+ .replace(",SPHEROID", ",<br>" + "SPHEROID")
+ .replace("],", "],<br>");
+ return HROGC;
+
+ }
+
+ /**
+ * Decode a OGC string to get a unique SRS string definition. This method
is
+ * able to understand some WKT common aliases, like OGC WKT and ESRI
WKTCode.
+ * For instance: "WGS 84 / UTM zone 32", "WGS 1984 UTM zone 32" and
+ * "WGS_84_UTM_Zone_32" are converted to the same string
+ * "WGS 84 UTM zone 32"
+ *
+ * @param textProj
+ * <String> - OGC/ESRI/other WKT code
+ * @return <String> - SRS definition
+ */
+ private static String decodeProjDescription(String textProj) {
+ String prjname = "";
+ try {
+ // Workaround if aux.xml has been download from web.
+ // convert HTML quotes ["] to ["]
+ textProj = textProj.replaceAll(""", "\"");
+ int start = textProj.indexOf("[\"");
+ int end = textProj.indexOf("\",", start);
+ prjname = textProj.substring(start + 2, end);
+ // The following set of replacements allows to "harmonize" OGC,
ESRI
+ // and
+ // few other WKT projection definitions
+ prjname = prjname.replaceAll("_", " ").replace(" / ", " ")
+ .replaceAll("\\bft US\\b", "(ftUS)")
+ .replaceAll("\\bftUS\\b", "(ftUS)")
+ .replaceAll("\\bft\\b", "(ft)").replaceAll("feet", "ft")
+ .replaceAll("WGS 1984", "WGS 84")
+ .replaceAll("NAD 1983 UTM", "NAD83 UTM")
+ .replaceAll("HARN", "(HARN)")
+ .replaceAll("\\bCSRS98\\b", "(CSRS98)")
+ .replaceAll("CSRS", "(CSRS)")
+ .replaceAll("\\bNSRS2007\\b", "(NSRS2007)")
+ .replaceAll("\\bNAD27_76\\b", "NAD27(76)")
+ .replaceAll("\\bCGQ77\\b", " (CGQ77)")
+ .replaceAll("\\bED77\\b", "(ED77)")
+ .replaceAll("\\b1942 83\\b", "1942(83)")
+ .replaceAll("\\b1942 58\\b", "1942(58)")
+ .replaceAll("\\bSegara Jakarta\\b", "Segara (Jakarta)")
+ .replaceAll("\\bRome\\b", "(Rome)")
+ .replaceAll("\\bParis\\b", "(Paris)")
+ .replaceAll("\\bFerro\\b", "(Ferro)");
+
+ } catch (Exception ex) {
+ // If there is other info than a WKT definition in the aux file
+ prjname = NOT_RECOGNIZED;
+ }
+ return prjname;
+ }
+
+ /**
+ * returns OGC WKT string located between projection tags (<WKT> or <SRS>)
+ * in a projection auxiliary file (AUX.XML)
+ *
+ * @param textProj
+ * string
+ * @return OGC WKT string
+ */
+ private static String getWktProjDefinition(String textProj) {
+ String prjname = "";
+ try {
+ if (textProj.contains("<WKT>")) {
+ int start = textProj.indexOf("<WKT>");
+ int end = textProj.indexOf("</WKT>", start);
+ prjname = textProj.substring(start, end);
+ } else if (textProj.contains("<SRS>")) {
+ int start = textProj.indexOf("<SRS>");
+ int end = textProj.indexOf("</SRS>", start);
+ prjname = textProj.substring(start, end);
+ } else
+ prjname = textProj;
+ } catch (Exception ex) {
+ prjname = textProj;
+ }
+ return prjname;
+ }
+
+
+ /**
+ * Method to get SRID from a layer from Style or from auxiliary file. First
+ * scans SRIDStyle, than auxiliary file or GeoTIFF tag. If SRID does not
+ * exist, it returns 0.
+ *
+ * @param layer
+ * @return SRID
+ * @throws Exception
+ */
+ public static SRSInfo getSRSInfoFromLayerStyleOrSource(Layer layer) throws
Exception {
+ SRSInfo srsInfo = new SRSInfo();
+
+ // First we check if a SRID != 0 has already been recorded for this
Layer
+ SRIDStyle sridStyle = (SRIDStyle) layer.getStyle(SRIDStyle.class);
+ final int oldSRID = sridStyle.getSRID();
+ if (oldSRID > 0) {
+ srsInfo.setCode(Integer.toString(oldSRID));
+ srsInfo.complete();
+ // If no SRID has been identified, search for a SRID in geotiff
tagg
+ // or into an auxiliary file
+ } else {
+ srsInfo = getSRSInfoFromLayerSource(layer);
+ }
+
+ return srsInfo;
+
+ }
+
+ /**
+ * Method to get SRID from a layer file from auxiliary files (.prj or .aux)
+ * or GeoTIFFed tag. If the auxiliary file SRID does not exist, it returns
+ * 0.
+ *
+ * @param layer
+ * @return SRID
+ * @throws Exception
+ */
+ public static SRSInfo getSRSInfoFromLayerSource(Layer layer) throws
Exception {
+ String fileSourcePath = "";
+ SRSInfo srsInfo = new SRSInfo();
+ // Raster layer case
+ if (isImageFileLayer(layer)) {
+ FeatureCollection featureCollection = layer
+ .getFeatureCollectionWrapper();
+ String sourcePathImage = null;
+ for (Iterator<?> i = featureCollection.iterator(); i.hasNext();) {
+ Feature feature = (Feature) i.next();
+ sourcePathImage =
feature.getString(ImageryLayerDataset.ATTR_URI);
+ if (sourcePathImage != null && sourcePathImage.length()>5) {
+ sourcePathImage = sourcePathImage.substring(5);
+ File f = new File(sourcePathImage);
+ if (f.exists()) {
+ fileSourcePath = f.getAbsolutePath().replace("%20", "
");
+ break;
+ }
+ }
+ }
+ String extension =
FileUtil.getExtension(fileSourcePath).toUpperCase();
+ if ((extension.equals("TIF") || extension.equals("TIFF"))) {
+ // If TIFF file is a geotiff, it scans into embedded tag
+ TiffTags.TiffMetadata metadata = TiffTags.readMetadata(new
File(fileSourcePath));
+ if (metadata.isGeoTiff()) {
+ srsInfo = metadata.getSRSInfo();
+ // If the TIF file is not a GeoTIFF it looks
+ // for a proj code into aux files
+ } else {
+ srsInfo = getSRSInfoFromAuxiliaryFile(fileSourcePath);
+ }
+ } else {
+ srsInfo = getSRSInfoFromAuxiliaryFile(fileSourcePath);
+ }
+ }
+ // Vector layer case
+ else {
+ if (!isDataBaseLayer(layer)) {
+ DataSourceQuery dsq = layer.getDataSourceQuery();
+ Object fnameObj =
dsq.getDataSource().getProperties().get("File");
+ fileSourcePath = fnameObj.toString();
+ srsInfo = getSRSInfoFromAuxiliaryFile(fileSourcePath);
+ }
+ }
+ return srsInfo;
+ }
+
+ // Boolean. Selected layer is related to an image file
+ private static boolean isImageFileLayer(Layer layer) {
+ if (layer.getStyle(ReferencedImageStyle.class) != null
+ && (layer.getDescription() != null)) {
+ return true;
+ } else {
+ return false;
+ }
+ }
+
+ // Boolean. Selected layer is related to a database
+ private static boolean isDataBaseLayer(Layer layer) {
+ DataSourceQuery dsq = layer.getDataSourceQuery();
+ if (dsq == null || dsq.getDataSource() instanceof DataStoreDataSource)
{
+ return true;
+ } else {
+ return false;
+ }
+ }
+
+}
+
Deleted: core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils2.java
===================================================================
--- core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils2.java
2016-09-11 13:44:55 UTC (rev 5016)
+++ core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils2.java
2016-09-16 18:58:48 UTC (rev 5017)
@@ -1,330 +0,0 @@
-package org.openjump.core.ccordsys.utils;
-
-import com.vividsolutions.jump.I18N;
-import com.vividsolutions.jump.feature.Feature;
-import com.vividsolutions.jump.feature.FeatureCollection;
-import com.vividsolutions.jump.io.datasource.DataSourceQuery;
-import com.vividsolutions.jump.util.FileUtil;
-import com.vividsolutions.jump.workbench.imagery.ImageryLayerDataset;
-import com.vividsolutions.jump.workbench.imagery.ReferencedImageStyle;
-import com.vividsolutions.jump.workbench.model.Layer;
-import
com.vividsolutions.jump.workbench.ui.plugin.datastore.DataStoreDataSource;
-import org.apache.commons.imaging.ImageReadException;
-import org.apache.commons.io.FilenameUtils;
-import org.openjump.core.ccordsys.srid.SRIDStyle;
-import org.openjump.core.rasterimage.TiffTags;
-
-
-import java.io.File;
-import java.io.IOException;
-import java.net.URISyntaxException;
-import java.util.Iterator;
-import java.util.Scanner;
-
-/**
- * Giuseppe Aruta [23_3_2016] This class is used to recognize file
- * projection. There are different methods A) a method to decode projection
- * information from GeoTIFF metadata. B) a method to decode projection info
- * from auxiliary files (.proj and .aux.xml).
- * http://landsathandbook.gsfc.nasa.gov/pdfs/geotiff_spec.pdf
- * http://www.sno.phy.queensu.ca/~phil/exiftool/TagNames/GeoTiff.html
- * http://www.remotesensing.org/geotiff/spec/geotiff6.html
- *
- * the following datasets have been used to test this library: QGIS
- * (http://qgis.org/downloads/data/) OSGEO data samples (
- * http://download.osgeo.org/geotiff/samples/)
- */
-public class ProjUtils2 {
-
- private static final String PROJECTION_UNSPECIFIED = I18N
-
.get("org.openjump.core.ui.plugin.raster.RasterImageLayerPropertiesPlugIn.unknown_projection");
- private static final String USER_DEFINED = I18N
-
.get("org.openjump.core.ui.plugin.layer.LayerPropertiesPlugIn.User_defined");
- private static final String NOT_RECOGNIZED = I18N
-
.get("org.openjump.core.ui.plugin.layer.LayerPropertiesPlugIn.Not_recognized");
-
-
- /**
- * - Read SRS from GeoTIFF tag - This method gets projection srid code and
- * projection info from a geotiff file. It first scans GeoKeyDirectoryTag
to
- * get either geographic/geocentric (2048 - GeographicTypeGeoKey),
projected
- * (3072 - ProjectedCSTypeGeoKey) or vertical (4096 - VerticalCSTypeGeoKey)
- * info. If no key ID is identified, it scans for GeoAsciiParamsTag
- * projection definition. Last choice, it search for an auxiliary file
- *
- * @param fileSourcePath
- * . eg. "c\documents\folder\image.tif"
- * @return the projection srid as a string. eg "32632"
- * @throws IOException
- * @throws URISyntaxException
- */
- @SuppressWarnings("static-access")
- public static SRSInfo getSRSInfoFromGeoTiff(String fileSourcePath)
- throws IOException, URISyntaxException,
TiffTags.TiffReadingException, ImageReadException {
- return TiffTags.readMetadata(new File(fileSourcePath)).getSRSInfo();
- }
-
-
- /**
- * - Read SRS from auxiliary file - Method to get a SRS (SRID code + SRID
- * definition) scanning the aux projection file (AUX.XML or PRJ file) for a
- * search string (SRID code or SRID definition). It scans into the registry
- * file (srid.txt) to find a correspondence between the search string and
- * lines of the srid.txt. If the source string corresponds as substring to
a
- * line, it returns the complete line as string. For instance, search
- * strings like "NAD83 UTM zone 10N" or "26910" both return
- * "SRID:26910 - NAD83 UTM zone 10N".
- *
- * @param fileSourcePath auxiliary file path
- * @return SRSInfo and Projection definition
- * @throws URISyntaxException
- * @throws IOException
- */
-
- public static SRSInfo getSRSInfoFromAuxiliaryFile(String fileSourcePath)
- throws URISyntaxException, IOException {
-
- // --- it reads an auxiliary file and decode a possible proj
- // --- definition to a simple string. Ex. "WGS 84 UTM Zone 32"
- int pos = fileSourcePath.lastIndexOf('.');
- // .shp, .dxf, .asc, .flt files
- String projectSourceFilePrj = fileSourcePath.substring(0, pos) +
".prj";
- // image files
- String projectSourceRFilePrj = fileSourcePath + ".prj";
- String projectSourceRFileAux = fileSourcePath + ".aux.xml";
-
- String type = FilenameUtils.getExtension(fileSourcePath).toUpperCase();
-
- SRSInfo srsInfo = new SRSInfo();
- String textProj = "";
- if (type.matches("SHP|DXF|ASC|FLT|ADF|GRD|BIL")) {
- if (new File(projectSourceFilePrj).exists()) {
- Scanner scanner = new Scanner(new File(projectSourceFilePrj));
- textProj = scanner.nextLine();
- srsInfo = new SRSInfo().setSource(projectSourceFilePrj);
- scanner.close();
- }
- } else if (new File(projectSourceRFileAux).exists()) {
- Scanner scanner = new Scanner(new File(projectSourceRFileAux));
- textProj = scanner.useDelimiter("\\A").next();
- if (!textProj.contains("<WKT>") && !textProj.contains("<SRS>") &&
- new File(projectSourceRFilePrj).exists()) {
- Scanner scanner2 = new Scanner(new
File(projectSourceRFilePrj));
- textProj = getWktProjDefinition(scanner2.nextLine());
- srsInfo = new SRSInfo().setSource(projectSourceRFileAux);
- scanner2.close();
- }
- } else if (new File(projectSourceRFilePrj).exists()) {
- Scanner scanner = new Scanner(new File(projectSourceRFilePrj));
- textProj = scanner.nextLine();
- srsInfo = new SRSInfo().setSource(projectSourceRFilePrj);
- scanner.close();
- }
-
- String prjName = decodeProjDescription(textProj);
- srsInfo.setDescription(prjName);
- srsInfo.complete();
- return srsInfo;
- }
-
-
- /**
- * Method to show an OGC WKT string in a more readable style
- *
- * @param WKT OGC WKT from auxiliary proj file
- * @return Readable string
- */
- public static String readableFormatWKTCode(String WKT) {
- String HROGC = "";
- // String add_spaces = String.format("%" + count_add++ + "s", "");
- HROGC = WKT.replace(",GEOGCS", ",<br>" + "GEOCS")
- .replace(",DATUM", ",<br>" + "DATUM")
- .replace(",SPHEROID", ",<br>" + "SPHEROID")
- .replace("],", "],<br>");
- return HROGC;
-
- }
-
- /**
- * Decode a OGC string to get a unique SRS string definition. This method
is
- * able to understand some WKT common aliases, like OGC WKT and ESRI
WKTCode.
- * For instance: "WGS 84 / UTM zone 32", "WGS 1984 UTM zone 32" and
- * "WGS_84_UTM_Zone_32" are converted to the same string
- * "WGS 84 UTM zone 32"
- *
- * @param textProj
- * <String> - OGC/ESRI/other WKT code
- * @return <String> - SRS definition
- */
- private static String decodeProjDescription(String textProj) {
- String prjname = "";
- try {
- // Workaround if aux.xml has been download from web.
- // convert HTML quotes ["] to ["]
- textProj = textProj.replaceAll(""", "\"");
- int start = textProj.indexOf("[\"");
- int end = textProj.indexOf("\",", start);
- prjname = textProj.substring(start + 2, end);
- // The following set of replacements allows to "harmonize" OGC,
ESRI
- // and
- // few other WKT projection definitions
- prjname = prjname.replaceAll("_", " ").replace(" / ", " ")
- .replaceAll("\\bft US\\b", "(ftUS)")
- .replaceAll("\\bftUS\\b", "(ftUS)")
- .replaceAll("\\bft\\b", "(ft)").replaceAll("feet", "ft")
- .replaceAll("WGS 1984", "WGS 84")
- .replaceAll("NAD 1983 UTM", "NAD83 UTM")
- .replaceAll("HARN", "(HARN)")
- .replaceAll("\\bCSRS98\\b", "(CSRS98)")
- .replaceAll("CSRS", "(CSRS)")
- .replaceAll("\\bNSRS2007\\b", "(NSRS2007)")
- .replaceAll("\\bNAD27_76\\b", "NAD27(76)")
- .replaceAll("\\bCGQ77\\b", " (CGQ77)")
- .replaceAll("\\bED77\\b", "(ED77)")
- .replaceAll("\\b1942 83\\b", "1942(83)")
- .replaceAll("\\b1942 58\\b", "1942(58)")
- .replaceAll("\\bSegara Jakarta\\b", "Segara (Jakarta)")
- .replaceAll("\\bRome\\b", "(Rome)")
- .replaceAll("\\bParis\\b", "(Paris)")
- .replaceAll("\\bFerro\\b", "(Ferro)");
-
- } catch (Exception ex) {
- // If there is other info than a WKT definition in the aux file
- prjname = NOT_RECOGNIZED;
- }
- return prjname;
- }
-
- /**
- * returns OGC WKT string located between projection tags (<WKT> or <SRS>)
- * in a projection auxiliary file (AUX.XML)
- *
- * @param textProj
- * string
- * @return OGC WKT string
- */
- private static String getWktProjDefinition(String textProj) {
- String prjname = "";
- try {
- if (textProj.contains("<WKT>")) {
- int start = textProj.indexOf("<WKT>");
- int end = textProj.indexOf("</WKT>", start);
- prjname = textProj.substring(start, end);
- } else if (textProj.contains("<SRS>")) {
- int start = textProj.indexOf("<SRS>");
- int end = textProj.indexOf("</SRS>", start);
- prjname = textProj.substring(start, end);
- } else
- prjname = textProj;
- } catch (Exception ex) {
- prjname = textProj;
- }
- return prjname;
- }
-
-
- /**
- * Method to get SRID from a layer from Style or from auxiliary file. First
- * scans SRIDStyle, than auxiliary file or GeoTIFF tag. If SRID does not
- * exist, it returns 0.
- *
- * @param layer
- * @return SRID
- * @throws Exception
- */
- public static SRSInfo getSRSInfoFromLayerStyleOrSource(Layer layer) throws
Exception {
- SRSInfo srsInfo = new SRSInfo();
-
- // First we check if a SRID != 0 has already been recorded for this
Layer
- SRIDStyle sridStyle = (SRIDStyle) layer.getStyle(SRIDStyle.class);
- final int oldSRID = sridStyle.getSRID();
- if (oldSRID > 0) {
- srsInfo.setCode(Integer.toString(oldSRID));
- srsInfo.complete();
- // If no SRID has been identified, search for a SRID in geotiff
tagg
- // or into an auxiliary file
- } else {
- srsInfo = getSRSInfoFromLayerSource(layer);
- }
-
- return srsInfo;
-
- }
-
- /**
- * Method to get SRID from a layer file from auxiliary files (.prj or .aux)
- * or GeoTIFFed tag. If the auxiliary file SRID does not exist, it returns
- * 0.
- *
- * @param layer
- * @return SRID
- * @throws Exception
- */
- public static SRSInfo getSRSInfoFromLayerSource(Layer layer) throws
Exception {
- String fileSourcePath = "";
- SRSInfo srsInfo = new SRSInfo();
- // Raster layer case
- if (isImageFileLayer(layer)) {
- FeatureCollection featureCollection = layer
- .getFeatureCollectionWrapper();
- String sourcePathImage = null;
- for (Iterator<?> i = featureCollection.iterator(); i.hasNext();) {
- Feature feature = (Feature) i.next();
- sourcePathImage =
feature.getString(ImageryLayerDataset.ATTR_URI);
- if (sourcePathImage != null && sourcePathImage.length()>5) {
- sourcePathImage = sourcePathImage.substring(5);
- File f = new File(sourcePathImage);
- if (f.exists()) {
- fileSourcePath = f.getAbsolutePath().replace("%20", "
");
- break;
- }
- }
- }
- String extension =
FileUtil.getExtension(fileSourcePath).toUpperCase();
- if ((extension.equals("TIF") || extension.equals("TIFF"))) {
- // If TIFF file is a geotiff, it scans into embedded tag
- if (ProjUtils.isGeoTIFF(fileSourcePath)) {
- srsInfo = TiffTags.readMetadata(new
File(fileSourcePath)).getSRSInfo();
- // If the TIF file is not a GeoTIFF it looks
- // for a proj code into aux files
- } else {
- srsInfo = getSRSInfoFromAuxiliaryFile(fileSourcePath);
- }
- } else {
- srsInfo = getSRSInfoFromAuxiliaryFile(fileSourcePath);
- }
- }
- // Vector layer case
- else {
- if (!isDataBaseLayer(layer)) {
- DataSourceQuery dsq = layer.getDataSourceQuery();
- Object fnameObj =
dsq.getDataSource().getProperties().get("File");
- fileSourcePath = fnameObj.toString();
- srsInfo = getSRSInfoFromAuxiliaryFile(fileSourcePath);
- }
- }
- return srsInfo;
- }
-
- // Boolean. Selected layer is related to an image file
- private static boolean isImageFileLayer(Layer layer) {
- if (layer.getStyle(ReferencedImageStyle.class) != null
- && (layer.getDescription() != null)) {
- return true;
- } else {
- return false;
- }
- }
-
- // Boolean. Selected layer is related to a database
- private static boolean isDataBaseLayer(Layer layer) {
- DataSourceQuery dsq = layer.getDataSourceQuery();
- if (dsq == null || dsq.getDataSource() instanceof DataStoreDataSource)
{
- return true;
- } else {
- return false;
- }
- }
-
-}
-
Copied: core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtilsOld.java (from
rev 5005, core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtils.java)
===================================================================
--- core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtilsOld.java
(rev 0)
+++ core/trunk/src/org/openjump/core/ccordsys/utils/ProjUtilsOld.java
2016-09-16 18:58:48 UTC (rev 5017)
@@ -0,0 +1,1151 @@
+package org.openjump.core.ccordsys.utils;
+
+import java.awt.Point;
+import java.awt.image.BufferedImage;
+import java.io.File;
+import java.io.IOException;
+import java.io.InputStream;
+import java.io.InputStreamReader;
+import java.net.URISyntaxException;
+import java.util.ArrayList;
+import java.util.HashMap;
+import java.util.Iterator;
+import java.util.List;
+import java.util.Scanner;
+
+import javax.imageio.ImageIO;
+
+import org.apache.commons.imaging.formats.tiff.TiffField;
+import org.apache.commons.imaging.formats.tiff.TiffImageMetadata;
+import org.apache.commons.imaging.formats.tiff.TiffImageParser;
+import org.apache.commons.io.FilenameUtils;
+import org.openjump.core.ccordsys.srid.SRIDStyle;
+import org.openjump.core.rasterimage.GeoTiffConstants;
+
+import com.sun.media.jai.codec.FileSeekableStream;
+import com.sun.media.jai.codec.TIFFDirectory;
+import com.sun.media.jai.codec.TIFFField;
+import com.vividsolutions.jts.geom.Coordinate;
+import com.vividsolutions.jts.geom.Envelope;
+import com.vividsolutions.jump.I18N;
+import com.vividsolutions.jump.feature.Feature;
+import com.vividsolutions.jump.feature.FeatureCollection;
+import com.vividsolutions.jump.io.datasource.DataSourceQuery;
+import com.vividsolutions.jump.util.FileUtil;
+import com.vividsolutions.jump.workbench.Logger;
+import com.vividsolutions.jump.workbench.imagery.ImageryLayerDataset;
+import com.vividsolutions.jump.workbench.imagery.ReferencedImageStyle;
+import com.vividsolutions.jump.workbench.model.Layer;
+import
com.vividsolutions.jump.workbench.ui.plugin.datastore.DataStoreDataSource;
+
+public class ProjUtilsOld {
+
+ private static final String PROJECTION_UNSPECIFIED = I18N
+
.get("org.openjump.core.ui.plugin.raster.RasterImageLayerPropertiesPlugIn.unknown_projection");
+ private static final String USER_DEFINED = I18N
+
.get("org.openjump.core.ui.plugin.layer.LayerPropertiesPlugIn.User_defined");
+ private static final String NOT_RECOGNIZED = I18N
+
.get("org.openjump.core.ui.plugin.layer.LayerPropertiesPlugIn.Not_recognized");
+
+ /*
+ * Giuseppe Aruta [23_3_2016] This class is used to recognize file
+ * projection. There are different methods A) a method to decode projection
+ * information from GeoTIFF metadata. B) a method to decode projection info
+ * from auxiliary files (.proj and .aux.xml).
+ * http://landsathandbook.gsfc.nasa.gov/pdfs/geotiff_spec.pdf
+ * http://www.sno.phy.queensu.ca/~phil/exiftool/TagNames/GeoTiff.html
+ * http://www.remotesensing.org/geotiff/spec/geotiff6.html
+ *
+ * the following datasets have been used to test this library: QGIS
+ * (http://qgis.org/downloads/data/) OSGEO data samples (
+ * http://download.osgeo.org/geotiff/samples/)
+ */
+
+ /**
+ * Test method - to get GeoTiff envelope in case of no srid def.
+ *
+ * @param fileSourcePath
+ * @return envelope
+ */
+ public static Envelope GeoTiffEnvelope(String fileSourcePath)
+ throws Exception {
+ BufferedImage image = ImageIO.read(new File(fileSourcePath));
+ Point imageDimensions = new Point(image.getWidth(), image.getHeight());
+ Envelope env = null;
+ Coordinate tiePoint = null, pixelOffset = null, pixelScale = null;
+ double[] doubles;
+ FileSeekableStream fileSeekableStream = new FileSeekableStream(
+ fileSourcePath);
+ TIFFDirectory tiffDirectory = new TIFFDirectory(fileSeekableStream, 0);
+
+ TIFFField[] availTags = tiffDirectory.getFields();
+
+ for (TIFFField availTag : availTags) {
+ if (availTag.getTag() == GeoTiffConstants.ModelTiepointTag) {
+ doubles = availTag.getAsDoubles();
+ if (doubles.length != 6) {
+ throw new Exception(
+ "unsupported value for ModelTiepointTag ("
+ + GeoTiffConstants.ModelTiepointTag + ")");
+ }
+ if (doubles[0] != 0 || doubles[1] != 0 || doubles[2] != 0) {
+ if (doubles[2] == 0)
+ pixelOffset = new Coordinate(doubles[0], doubles[1]);
+ else
+ pixelOffset = new Coordinate(doubles[0], doubles[1],
+ doubles[2]);
+ }
+ if (doubles[5] == 0)
+ tiePoint = new Coordinate(doubles[3], doubles[4]);
+ else
+ tiePoint = new Coordinate(doubles[3], doubles[4],
+ doubles[5]);
+
+ } else if (availTag.getTag() ==
GeoTiffConstants.ModelPixelScaleTag) {
+ // Karteneinheiten pro pixel x bzw. y
+ doubles = availTag.getAsDoubles();
+ if (doubles[2] == 0)
+ pixelScale = new Coordinate(doubles[0], doubles[1]);
+ else
+ pixelScale = new Coordinate(doubles[0], doubles[1],
+ doubles[2]);
+ } else {
+ }
+ }
+
+ fileSeekableStream.close();
+
+ if (tiePoint != null && pixelScale != null) {
+ Coordinate upperLeft;
+ Coordinate lowerRight;
+
+ if (pixelOffset == null) {
+ upperLeft = tiePoint;
+ } else {
+ upperLeft = new Coordinate(tiePoint.x
+ - (pixelOffset.x * pixelScale.x), tiePoint.y
+ - (pixelOffset.y * pixelScale.y));
+ }
+
+ lowerRight = new Coordinate(upperLeft.x
+ + (imageDimensions.x * pixelScale.x), upperLeft.y
+ - (imageDimensions.y * pixelScale.y));
+
+ env = new Envelope(upperLeft, lowerRight);
+ }
+ return env;
+ }
+
+ /**
+ * Test method to read the entire GeoKeyDirectoryTag, only used to study
+ * GeoKeyDirectoryTag structure
+ *
+ * @param fileSourcePath
+ * @return GeoKeyDirectoryTag
+ * @throws IOException
+ * @throws URISyntaxException
+ */
+ @SuppressWarnings("static-access")
+ public static String readGeoTiffGeoKeyDirectoryTag(String fileSourcePath)
+ throws IOException, URISyntaxException {
+ String prjname = "";
+ File tiffFile = new File(fileSourcePath);
+ try {
+ TiffImageParser parser = new TiffImageParser();
+ TiffImageMetadata metadata = (TiffImageMetadata) parser
+ .getMetadata(tiffFile);
+ if (metadata != null) {
+ List<TiffField> tiffFields = metadata.getAllFields();
+ GeoTiffConstants constants = new GeoTiffConstants();
+ for (TiffField tiffField : tiffFields) {
+ if (tiffField.getTag() == constants.GeoKeyDirectoryTag) {
+ String GeoDirTag = tiffField.getValueDescription();
+ prjname = GeoDirTag;
+ }
+ }
+ } else {
+ prjname = "GeoKeyDirectoryTagis empty";
+ }
+ } catch (Exception ex) {
+ prjname = PROJECTION_UNSPECIFIED;
+ }
+ return prjname;
+ }
+
+ /**
+ * - Read SRS from GeoTIFF tag - This method gets projection srid code and
+ * projection info from a geotiff file. It first scans GeoKeyDirectoryTag
to
+ * get either geographic/geocentric (2048 - GeographicTypeGeoKey),
projected
+ * (3072 - ProjectedCSTypeGeoKey) or vertical (4096 - VerticalCSTypeGeoKey)
+ * info. If no key ID is identified, it scans for GeoAsciiParamsTag
+ * projection definition. Last choice, it search for an auxiliary file
+ *
+ * @param fileSourcePath
+ * . eg. "c\documents\folder\image.tif"
+ * @return <String> projection srid code as string. eg "32632"
+ * @throws IOException
+ * @throws URISyntaxException
+ */
+ @SuppressWarnings("static-access")
+ public static String readSRSFromGeoTiffFile(String fileSourcePath)
+ throws IOException, URISyntaxException {
+ String GeoDirTag = "";
+ String GeoDirTag2 = "";
+ String prjname = "";
+ File tiffFile = new File(fileSourcePath);
+ try {
+ TiffImageParser parser = new TiffImageParser();
+ TiffImageMetadata metadata = (TiffImageMetadata) parser
+ .getMetadata(tiffFile);
+ if (metadata != null) {
+ List<TiffField> tiffFields = metadata.getAllFields();
+ GeoTiffConstants constants = new GeoTiffConstants();
+
+ String ID = "";
+ int start;
+ for (TiffField tiffField : tiffFields) {
+ if (tiffField.getTag() == constants.GeoKeyDirectoryTag) {
+ GeoDirTag = tiffField.getValueDescription();
+ }
+ if (tiffField.getTag() == constants.GeoAsciiParamsTag) {
+ GeoDirTag2 = tiffField.getStringValue().replaceAll(
+ "[\\t\\n\\r\\_\\|]", " ");
+ }
+ if (tiffField.getTag() == constants.ModelTiepointTag) {
+ }
+ }
+ if (GeoDirTag.contains("3072")) {
+ start = GeoDirTag.indexOf("3072");
+ ID = GeoDirTag.substring(start);
+ String[] parts = ID.split(",");
+ String part1 = parts[3];
+ if (!part1.contains("32767")) {
+ prjname = getSRSFromWkt(part1.replaceAll(" ", ""));
+ } else {
+ prjname = "SRID: " + USER_DEFINED + " - " + GeoDirTag2;
+ }
+ } else if (!GeoDirTag.contains("3072")
+ & GeoDirTag.contains("4096")) {
+ start = GeoDirTag.indexOf("4096");
+ ID = GeoDirTag.substring(start);
+ String[] parts = ID.split(",");
+ String part1 = parts[3];
+ if (!part1.contains("32767")) {
+ prjname = getSRSFromWkt(part1.replaceAll(" ", ""));
+ } else {
+ prjname = "SRID: " + USER_DEFINED + " - " + GeoDirTag2;
+ }
+ } else if (!GeoDirTag.contains("3072")
@@ Diff output truncated at 100000 characters. @@
------------------------------------------------------------------------------
_______________________________________________
Jump-pilot-devel mailing list
Jump-pilot-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/jump-pilot-devel
