001    //$HeadURL: svn+ssh://rbezema@svn.wald.intevation.org/deegree/base/branches/2.2_testing/src/org/deegree/model/coverage/grid/TiffDEMParser.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.model.coverage.grid;
044    
045    import java.awt.image.BufferedImage;
046    import java.awt.image.DataBuffer;
047    import java.io.File;
048    
049    import org.apache.batik.ext.awt.image.codec.FileCacheSeekableStream;
050    import org.apache.batik.ext.awt.image.codec.tiff.TIFFDecodeParam;
051    import org.apache.batik.ext.awt.image.codec.tiff.TIFFImage;
052    
053    /**
054     * Parses 4 channel (32Bit) tiff images as DEM and returns a float matrix containing the DEM heights
055     * 
056     * 
057     * @version $Revision: 9343 $
058     * @author <a href="mailto:poth@lat-lon.de">Andreas Poth</a>
059     * @author last edited by: $Author: apoth $
060     * 
061     * @version 1.0. $Revision: 9343 $, $Date: 2007-12-27 14:30:32 +0100 (Do, 27 Dez 2007) $
062     * 
063     * @since 2.0
064     */
065    public class TiffDEMParser {
066    
067        private DataBuffer db;
068    
069        private float scale = 0;
070    
071        private float offset = 0;
072    
073        private int width = 0;
074    
075        private int height = 0;
076    
077        /**
078         * 
079         * @param dataFile
080         *            image containing DEM data instead color information
081         */
082        public TiffDEMParser( File dataFile ) {
083            this( dataFile, 1, 0 );
084        }
085    
086        /**
087         * 
088         * @param dataFile
089         *            image containing DEM data instead color information
090         * @param scale
091         *            scale factor; newHeight[i][j] = height[i][j] * scale
092         */
093        public TiffDEMParser( File dataFile, float scale ) {
094            this( dataFile, scale, 0 );
095        }
096    
097        /**
098         * 
099         * @param dataFile
100         *            image containing DEM data instead color information
101         * @param scale
102         *            scale factor; newHeight[i][j] = height[i][j] * scale
103         * @param offset
104         *            height offset; newHeight[i][j] = height[i][j] + offset
105         */
106        public TiffDEMParser( File dataFile, float scale, float offset ) {
107            try {
108                FileCacheSeekableStream fss = new FileCacheSeekableStream( dataFile.toURL().openStream() );
109                TIFFImage tiff = new TIFFImage( fss, new TIFFDecodeParam(), 0 );
110                db = tiff.getData().getDataBuffer();
111                width = tiff.getWidth();
112                height = tiff.getHeight();
113            } catch ( Exception e ) {
114                e.printStackTrace();
115            }
116            this.scale = scale;
117            this.offset = offset;
118        }
119    
120        /**
121         * 
122         * @param data
123         *            image containing DEM data instead color information
124         */
125        public TiffDEMParser( BufferedImage data ) {
126            this( data, 1, 0 );
127        }
128    
129        /**
130         * 
131         * @param data
132         *            image containing DEM data instead color information
133         * @param scale
134         *            scale factor; newHeight[i][j] = height[i][j] * scale
135         */
136        public TiffDEMParser( BufferedImage data, float scale ) {
137            this( data, scale, 0 );
138        }
139    
140        /**
141         * 
142         * @param data
143         * @param scale
144         *            scale factor; newHeight[i][j] = height[i][j] * scale
145         * @param offset
146         *            height offset; newHeight[i][j] = height[i][j] + offset
147         */
148        public TiffDEMParser( BufferedImage data, float scale, float offset ) {
149            this.db = data.getRaster().getDataBuffer();
150            this.scale = scale;
151            this.offset = offset;
152            width = data.getWidth();
153            height = data.getHeight();
154        }
155    
156        /**
157         * returns the DEM heights as float matrix
158         * 
159         * @return the DEM heights as float matrix
160         */
161        public float[][] parse() {
162    
163            float[][] terrain = new float[height][width];
164    
165            for ( int j = 0; j < height; j++ ) {
166                for ( int i = 0; i < width; i++ ) {
167                    terrain[j][i] = db.getElemFloat( width * j + i ) * scale + offset;
168                }
169            }
170    
171            return terrain;
172    
173        }
174    
175    }