001: //$HeadURL: https://svn.wald.intevation.org/svn/deegree/base/trunk/src/org/deegree/tools/raster/ArcInfo2xyz.java $
002: /*---------------- FILE HEADER ------------------------------------------
003:
004: This file is part of deegree.
005: Copyright (C) 2001-2008 by:
006: EXSE, Department of Geography, University of Bonn
007: http://www.giub.uni-bonn.de/deegree/
008: lat/lon GmbH
009: http://www.lat-lon.de
010:
011: This library is free software; you can redistribute it and/or
012: modify it under the terms of the GNU Lesser General Public
013: License as published by the Free Software Foundation; either
014: version 2.1 of the License, or (at your option) any later version.
015:
016: This library is distributed in the hope that it will be useful,
017: but WITHOUT ANY WARRANTY; without even the implied warranty of
018: MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
019: Lesser General Public License for more details.
020:
021: You should have received a copy of the GNU Lesser General Public
022: License along with this library; if not, write to the Free Software
023: Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
024:
025: Contact:
026:
027: Andreas Poth
028: lat/lon GmbH
029: Aennchenstr. 19
030: 53177 Bonn
031: Germany
032: E-Mail: poth@lat-lon.de
033:
034: Prof. Dr. Klaus Greve
035: Department of Geography
036: University of Bonn
037: Meckenheimer Allee 166
038: 53115 Bonn
039: Germany
040: E-Mail: greve@giub.uni-bonn.de
041:
042: ---------------------------------------------------------------------------*/
043: package org.deegree.tools.raster;
044:
045: import java.io.BufferedReader;
046: import java.io.File;
047: import java.io.FileReader;
048: import java.io.PrintWriter;
049: import java.util.Properties;
050:
051: import org.deegree.framework.util.StringTools;
052: import org.deegree.io.arcinfo_raster.ArcInfoTextRasterReader;
053: import org.deegree.model.coverage.grid.WorldFile;
054: import org.deegree.model.spatialschema.Envelope;
055:
056: /**
057: * converts an ArcInfor raster file (text format):
058: * <pre>
059: * ncols 2404
060: * nrows 2307
061: * xllcorner 2627130
062: * yllcorner 5686612
063: * cellsize 10
064: * NODATA_value -9999
065: * 20636 20593 20569 20573 20571 20564 20564 20558 ...
066: * ...
067: * </pre>
068: * into its XYZ representation:
069: * <pre>
070: * 2627130.0 5709682.0 206.36
071: * 2627140.0 5709682.0 205.93
072: * 2627150.0 5709682.0 205.69
073: * 2627160.0 5709682.0 205.73
074: * 2627170.0 5709682.0 205.71
075: * 2627180.0 5709682.0 205.64
076: * 2627190.0 5709682.0 205.64
077: * 2627200.0 5709682.0 205.58
078: * </pre>
079: *
080: *
081: * @version $Revision: 9346 $
082: * @author <a href="mailto:poth@lat-lon.de">Andreas Poth</a>
083: * @author last edited by: $Author: apoth $
084: *
085: * @version 1.0. $Revision: 9346 $, $Date: 2007-12-27 08:39:07 -0800 (Thu, 27 Dec 2007) $
086: *
087: * @since 2.0
088: */
089: public class ArcInfo2xyz {
090:
091: public static void main(String[] args) throws Exception {
092:
093: Properties map = new Properties();
094: for (int i = 0; i < args.length; i += 2) {
095: map.put(args[i], args[i + 1]);
096: }
097: File arcInfo = new File(map.getProperty("-inFile"));
098: File xyz = new File(map.getProperty("-outFile"));
099:
100: ArcInfoTextRasterReader reader = new ArcInfoTextRasterReader(
101: arcInfo);
102: WorldFile wf = reader.readMetadata();
103: Envelope env = wf.getEnvelope();
104: double res = wf.getResx();
105:
106: BufferedReader in = null;
107: PrintWriter pw = null;
108: try {
109:
110: in = new BufferedReader(new FileReader(arcInfo));
111: for (int i = 0; i < 6; i++) {
112: // skip header
113: in.readLine();
114: }
115:
116: pw = new PrintWriter(xyz);
117:
118: String line = null;
119: double y = env.getMax().getY();
120: while ((line = in.readLine()) != null) {
121: System.out.print(y + "\r");
122: float[] data = StringTools.toArrayFloat(line, " \t");
123: double x = env.getMin().getX();
124: for (int i = 0; i < data.length; i++) {
125: pw.print(x);
126: pw.print(' ');
127: pw.print(y);
128: pw.print(' ');
129: pw.print(data[i]);
130: pw.println();
131: x = x + res;
132: }
133: y = y - res;
134: }
135: } catch (Exception e) {
136: throw e;
137: } finally {
138: pw.close();
139: in.close();
140: }
141:
142: }
143:
144: }
|