    001    //$HeadURL: https://svn.wald.intevation.org/svn/deegree/base/branches/2.4_testing/src/org/deegree/tools/srs/TransformRasterFile.java $
    002    /*----------------------------------------------------------------------------
    003     This file is part of deegree, http://deegree.org/
    004     Copyright (C) 2001-2009 by:
    005       Department of Geography, University of Bonn
    006     and
    007       lat/lon GmbH
    008
    009     This library is free software; you can redistribute it and/or modify it under
    010     the terms of the GNU Lesser General Public License as published by the Free
    011     Software Foundation; either version 2.1 of the License, or (at your option)
    012     any later version.
    013     This library is distributed in the hope that it will be useful, but WITHOUT
    014     ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
    015     FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
    016     details.
    017     You should have received a copy of the GNU Lesser General Public License
    018     along with this library; if not, write to the Free Software Foundation, Inc.,
    019     59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
    020
    021     Contact information:
    022
    023     lat/lon GmbH
    024     Aennchenstr. 19, 53177 Bonn
    025     Germany
    026     http://lat-lon.de/
    027
    028     Department of Geography, University of Bonn
    029     Prof. Dr. Klaus Greve
    030     Postfach 1147, 53001 Bonn
    031     Germany
    032     http://www.geographie.uni-bonn.de/deegree/
    033
    034     e-mail: info@deegree.org
    035    ----------------------------------------------------------------------------*/
    036    package org.deegree.tools.srs;
    037
    038    import java.awt.image.BufferedImage;
    039    import java.util.Properties;
    040
    041    import javax.media.jai.InterpolationNearest;
    042
    043    import org.deegree.datatypes.CodeList;
    044    import org.deegree.framework.util.ImageUtils;
    045    import org.deegree.model.coverage.grid.ImageGridCoverage;
    046    import org.deegree.model.coverage.grid.WorldFile;
    047    import org.deegree.model.crs.CRSFactory;
    048    import org.deegree.model.crs.CoordinateSystem;
    049    import org.deegree.model.crs.GeoTransformer;
    050    import org.deegree.model.spatialschema.Envelope;
    051    import org.deegree.ogcwebservices.SupportedFormats;
    052    import org.deegree.ogcwebservices.SupportedSRSs;
    053    import org.deegree.ogcwebservices.wcs.describecoverage.CoverageOffering;
    054    import org.deegree.ogcwebservices.wcs.describecoverage.DomainSet;
    055    import org.deegree.ogcwebservices.wcs.describecoverage.RangeSet;
    056    import org.deegree.ogcwebservices.wcs.describecoverage.SpatialDomain;
    057
    058    /**
    059     * Tool to transform raster files from one CRS to another.
    060     *
    061     * @author <a href="mailto:tonnhofer@lat-lon.de">Oliver Tonnhofer</a>
    062     * @author last edited by: $Author: mschneider $
    063     *
    064     * @version $Revision: 18195 $, $Date: 2009-06-18 17:55:39 +0200 (Do, 18. Jun 2009) $
    065     */
    066    public class TransformRasterFile {
    067
    068        private final static float DEFAULT_IMAGE_QUALITY = 0.9f;
    069
    070        private final static int DEFAULT_PP_GRID_SIZE = 5;
    071
    072        private final static int DEFAULT_POLY_ORDER = 3;
    073
    074        private static void transformRasterFile( String inFile, String inCRS, String outFile, String outCRS,
    075                                                 Float imageQuality, int ppgridsize, int polynomOrder )
    076                                throws Exception {
    077
    078            System.out.println( "Loading raster " + inFile );
    079            BufferedImage image = ImageUtils.loadImage( inFile );
    080
    081            CoordinateSystem sourceCRS = CRSFactory.create( inCRS );
    082            CoordinateSystem targetCRS = CRSFactory.create( outCRS );
    083
    084            WorldFile worldFile = WorldFile.readWorldFile( inFile, WorldFile.TYPE.CENTER, image );
    085            Envelope inEnvelope = worldFile.getEnvelope();
    086
    087            // create minimal CoverageOffering for ImageGridCoverage
    088            // most parts are not used
    089            DomainSet ds = new DomainSet( new SpatialDomain( new Envelope[] { inEnvelope } ) );
    090            RangeSet rs = new RangeSet( "", "" );
    091            CodeList[] dummyCodeList = new CodeList[] { new CodeList( "", new String[] {} ) };
    092            CodeList[] nativeSRSCodeList = new CodeList[] { new CodeList( "", new String[] { inCRS } ) };
    093
    094            SupportedSRSs supSRSs = new SupportedSRSs( dummyCodeList, dummyCodeList, dummyCodeList, nativeSRSCodeList );
    095
    096            SupportedFormats supFormats = new SupportedFormats( dummyCodeList );
    097
    098            CoverageOffering coverageOffering = new CoverageOffering( "", "", "", null, null, null, ds, rs, supSRSs,
    099                                                                      supFormats, null, null );
    100
    101            ImageGridCoverage igc = new ImageGridCoverage( coverageOffering, inEnvelope, image );
    102
    103            GeoTransformer gt = new GeoTransformer( targetCRS );
    104
    105            Envelope outEnvelope = gt.transform( inEnvelope, sourceCRS, true );
    106
    107            // calculate new output size
    108            // use square pixels for output, ie. the aspect ratio changes
    109            double deltaX = outEnvelope.getWidth();
    110            double deltaY = outEnvelope.getHeight();
    111            double diagSize = Math.sqrt( deltaX * deltaX + deltaY * deltaY );
    112            // pixelSize for calculation of the new image size
    113            double pixelSize = diagSize / Math.sqrt( Math.pow( image.getWidth(), 2 ) + Math.pow( image.getHeight(), 2 ) );
    114            int height = (int) ( deltaY / pixelSize + 0.5 );
    115            int width = (int) ( deltaX / pixelSize + 0.5 );
    116            // realPixelSize for center type world files, etc.
    117            double realPixelSize = diagSize
    118                                   / Math.sqrt( Math.pow( image.getWidth() - 1, 2 ) + Math.pow( image.getHeight() - 1, 2 ) );
    119
    120            System.out.println( "Transforming raster from " + inCRS + " to " + outCRS );
    121            igc = (ImageGridCoverage) gt.transform( igc, outEnvelope, width, height, ppgridsize, polynomOrder,
    122                                                    new InterpolationNearest() );
    123
    124            image = igc.getAsImage( -1, -1 );
    125
    126            System.out.println( "Saving raster " + outFile );
    127            ImageUtils.saveImage( image, outFile, imageQuality );
    128
    129            // save new WorldFile
    130            WorldFile outWorldFile = new WorldFile( realPixelSize, realPixelSize, 0.0f, 0.0f, outEnvelope );
    131            String basename = outFile.substring( 0, outFile.lastIndexOf( "." ) );
    132            WorldFile.writeWorldFile( outWorldFile, basename );
    133
    134        }
    135
    136        private static void printHelpAndExit() {
    137            System.out.println( "Usage: java [...] org.deegree.tools.srs.TransformRasterFile " );
    138            System.out.println( "                  -inFile filename -inCRS crs " );
    139            System.out.println( "                  [-outFile filename] -outCRS crs" );
    140            System.out.println( "                  [-imageQuality 0.X]" );
    141            System.out.println( "                  -passpointGridSize 5 -polynomOrder 3 " );
    142            System.exit( 1 );
    143        }
    144
    145        /**
    146         * @param args
    147         */
    148        public static void main( String[] args ) {
    149
    150            if ( args.length % 2 != 0 )
    151                printHelpAndExit();
    152
    153            Properties map = new Properties();
    154            for ( int i = 0; i < args.length; i += 2 ) {
    155                map.put( args[i], args[i + 1] );
    156            }
    157
    158            String outCRS = (String) map.get( "-outCRS" );
    159            if ( outCRS == null )
    160                printHelpAndExit();
    161
    162            String inCRS = (String) map.get( "-inCRS" );
    163            if ( inCRS == null )
    164                printHelpAndExit();
    165
    166            String inFilename = (String) map.get( "-inFile" );
    167            if ( inFilename == null )
    168                printHelpAndExit();
    169
    170            String outFilename = (String) map.get( "-outFile" );
    171            if ( outFilename == null ) {
    172                String ext = inFilename.substring( inFilename.lastIndexOf( "." ) );
    173                outFilename = inFilename.substring( 0, inFilename.lastIndexOf( "." ) );
    174                outFilename = outFilename + "." + outCRS + ext;
    175            }
    176
    177            String imageQualityString = (String) map.get( "-imageQuality" );
    178            float imageQuality = DEFAULT_IMAGE_QUALITY;
    179            if ( imageQualityString != null ) {
    180                imageQuality = Float.valueOf( imageQualityString );
    181            }
    182            String ppgridsizeString = (String) map.get( "-passpointGridSize" );
    183            int ppgridsize = DEFAULT_PP_GRID_SIZE;
    184            if ( ppgridsizeString != null ) {
    185                ppgridsize = Integer.parseInt( ppgridsizeString );
    186            }
    187
    188            String polynomOrderString = (String) map.get( "-polynomOrder" );
    189            int polynomOrder = DEFAULT_POLY_ORDER;
    190            if ( polynomOrderString != null ) {
    191                polynomOrder = Integer.parseInt( polynomOrderString );
    192            }
    193
    194            try {
    195                transformRasterFile( inFilename, inCRS, outFilename, outCRS, imageQuality, ppgridsize, polynomOrder );
    196            } catch ( Exception e ) {
    197                e.printStackTrace();
    198            }
    199        }
    200
    201    }