001 //$HeadURL: svn+ssh://rbezema@svn.wald.intevation.org/deegree/base/branches/2.2_testing/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 17:39:07 +0100 (Do, 27 Dez 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( arcInfo );
101 WorldFile wf = reader.readMetadata();
102 Envelope env = wf.getEnvelope();
103 double res = wf.getResx();
104
105 BufferedReader in = null;
106 PrintWriter pw = null;
107 try {
108
109 in = new BufferedReader( new FileReader( arcInfo ) );
110 for ( int i = 0; i < 6; i++ ) {
111 // skip header
112 in.readLine();
113 }
114
115 pw = new PrintWriter( xyz );
116
117 String line = null;
118 double y = env.getMax().getY();
119 while ( ( line = in.readLine() ) != null ) {
120 System.out.print( y + "\r" );
121 float[] data = StringTools.toArrayFloat( line, " \t" );
122 double x = env.getMin().getX();
123 for ( int i = 0; i < data.length; i++ ) {
124 pw.print( x );
125 pw.print( ' ' );
126 pw.print( y );
127 pw.print( ' ' );
128 pw.print( data[i] );
129 pw.println();
130 x = x + res;
131 }
132 y = y - res;
133 }
134 } catch ( Exception e ) {
135 throw e;
136 } finally {
137 pw.close();
138 in.close();
139 }
140
141 }
142
143 }