GDAL的gdalwarp类源码

/****************************************************************************** * $Id: gdalwarp.cpp 21298 2010-12-20 10:58:34Z rouault $ * * Project: High Performance Image Reprojector * Purpose: Test program for high performance warper API. * Author: Frank Warmerdam <warmerdam@pobox.com> * ****************************************************************************** * Copyright (c) 2002, i3 - information integration and imaging * Fort Collin, CO * * Permission is hereby granted, free of charge, to any person obtaining a * copy of this software and associated documentation files (the "Software"), * to deal in the Software without restriction, including without limitation * the rights to use, copy, modify, merge, publish, distribute, sublicense, * and/or sell copies of the Software, and to permit persons to whom the * Software is furnished to do so, subject to the following conditions: * * The above copyright notice and this permission notice shall be included * in all copies or substantial portions of the Software. * * THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS * OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, * FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL * THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER * LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING * FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER * DEALINGS IN THE SOFTWARE. ****************************************************************************/ #include "gdalwarper.h" #include "cpl_string.h" #include "ogr_spatialref.h" #include "ogr_api.h" CPL_CVSID("$Id: gdalwarp.cpp 21298 2010-12-20 10:58:34Z rouault $"); static void LoadCutline( const char *pszCutlineDSName, const char *pszCLayer, const char *pszCWHERE, const char *pszCSQL, void **phCutlineRet ); static void TransformCutlineToSource( GDALDatasetH hSrcDS, void *hCutline, char ***ppapszWarpOptions, char **papszTO ); static GDALDatasetH GDALWarpCreateOutput( char **papszSrcFiles, const char *pszFilename, const char *pszFormat, char **papszTO, char ***ppapszCreateOptions, GDALDataType eDT ); static double dfMinX=0.0, dfMinY=0.0, dfMaxX=0.0, dfMaxY=0.0; static double dfXRes=0.0, dfYRes=0.0; static int bTargetAlignedPixels = FALSE; static int nForcePixels=0, nForceLines=0, bQuiet = FALSE; static int bEnableDstAlpha = FALSE, bEnableSrcAlpha = FALSE; static int bVRT = FALSE; /******************************************************************************/ /*! \page gdalwarp gdalwarp image reprojection and warping utility \section wsynopsis SYNOPSIS \htmlonly Usage: \endhtmlonly \verbatim gdalwarp [--help-general] [--formats] [-s_srs srs_def] [-t_srs srs_def] [-to "NAME=VALUE"] [-order n] [-tps] [-rpc] [-geoloc] [-et err_threshold] [-te xmin ymin xmax ymax] [-tr xres yres] [-tap] [-ts width height] [-wo "NAME=VALUE"] [-ot Byte/Int16/...] [-wt Byte/Int16] [-srcnodata "value [value...]"] [-dstnodata "value [value...]"] -dstalpha [-r resampling_method] [-wm memory_in_mb] [-multi] [-q] [-cutline datasource] [-cl layer] [-cwhere expression] [-csql statement] [-cblend dist_in_pixels] [-crop_to_cutline] [-of format] [-co "NAME=VALUE"]* [-overwrite] srcfile* dstfile \endverbatim \section wdescription DESCRIPTION <p> The gdalwarp utility is an image mosaicing, reprojection and warping utility. The program can reproject to any supported projection, and can also apply GCPs stored with the image if the image is "raw" with control information. <p> <dl> <dt> <b>-s_srs</b> <em>srs def</em>:</dt><dd> source spatial reference set. The coordinate systems that can be passed are anything supported by the OGRSpatialReference.SetFromUserInput() call, which includes EPSG PCS and GCSes (ie. EPSG:4296), PROJ.4 declarations (as above), or the name of a .prf file containing well known text.</dd> <dt> <b>-t_srs</b> <em>srs_def</em>:</dt><dd> target spatial reference set. The coordinate systems that can be passed are anything supported by the OGRSpatialReference.SetFromUserInput() call, which includes EPSG PCS and GCSes (ie. EPSG:4296), PROJ.4 declarations (as above), or the name of a .prf file containing well known text.</dd> <dt> <b>-to</b> <em>NAME=VALUE</em>:</dt><dd> set a transformer option suitable to pass to GDALCreateGenImgProjTransformer2(). </dd> <dt> <b>-order</b> <em>n</em>:</dt><dd> order of polynomial used for warping (1 to 3). The default is to select a polynomial order based on the number of GCPs.</dd> <dt> <b>-tps</b>:</dt><dd>Force use of thin plate spline transformer based on available GCPs.</dd> <dt> <b>-rpc</b>:</dt> <dd>Force use of RPCs.</dd> <dt> <b>-geoloc</b>:</dt><dd>Force use of Geolocation Arrays.</dd> <dt> <b>-et</b> <em>err_threshold</em>:</dt><dd> error threshold for transformation approximation (in pixel units - defaults to 0.125).</dd> <dt> <b>-te</b> <em>xmin ymin xmax ymax</em>:</dt><dd> set georeferenced extents of output file to be created (in target SRS).</dd> <dt> <b>-tr</b> <em>xres yres</em>:</dt><dd> set output file resolution (in target georeferenced units)</dd> <dt> <b>-tap</b>:</dt><dd> (GDAL >= 1.8.0) (target aligned pixels) align the coordinates of the extent of the output file to the values of the -tr, such that the aligned extent includes the minimum extent.</dd> <dt> <b>-ts</b> <em>width height</em>:</dt><dd> set output file size in pixels and lines. If width or height is set to 0, the other dimension will be guessed from the computed resolution. Note that -ts cannot be used with -tr</dd> <dt> <b>-wo</b> <em>"NAME=VALUE"</em>:</dt><dd> Set a warp options. The GDALWarpOptions::papszWarpOptions docs show all options. Multiple <b>-wo</b> options may be listed.</dd> <dt> <b>-ot</b> <em>type</em>:</dt><dd> For the output bands to be of the indicated data type.</dd> <dt> <b>-wt</b> <em>type</em>:</dt><dd> Working pixel data type. The data type of pixels in the source image and destination image buffers.</dd> <dt> <b>-r</b> <em>resampling_method</em>:</dt><dd> Resampling method to use. Available methods are: <dl> <dt><b>near</b></dt>: <dd>nearest neighbour resampling (default, fastest algorithm, worst interpolation quality).</dd> <dt><b>bilinear</b></dt>: <dd>bilinear resampling.</dd> <dt><b>cubic</b></dt>: <dd>cubic resampling.</dd> <dt><b>cubicspline</b></dt>: <dd>cubic spline resampling.</dd> <dt><b>lanczos</b></dt>: <dd>Lanczos windowed sinc resampling.</dd> </dl> <dt> <b>-srcnodata</b> <em>value [value...]</em>:</dt><dd> Set nodata masking values for input bands (different values can be supplied for each band). If more than one value is supplied all values should be quoted to keep them together as a single operating system argument. Masked values will not be used in interpolation. Use a value of <tt>None</tt> to ignore intrinsic nodata settings on the source dataset.</dd> <dt> <b>-dstnodata</b> <em>value [value...]</em>:</dt><dd> Set nodata values for output bands (different values can be supplied for each band). If more than one value is supplied all values should be quoted to keep them together as a single operating system argument. New files will be initialized to this value and if possible the nodata value will be recorded in the output file.</dd> <dt> <b>-dstalpha</b>:</dt><dd> Create an output alpha band to identify nodata (unset/transparent) pixels. </dd> <dt> <b>-wm</b> <em>memory_in_mb</em>:</dt><dd> Set the amount of memory (in megabytes) that the warp API is allowed to use for caching.</dd> <dt> <b>-multi</b>:</dt><dd> Use multithreaded warping implementation. Multiple threads will be used to process chunks of image and perform input/output operation simultaneously.</dd> <dt>