[Gdal-dev] Re: Mosaic in GDAL Warp

[David W. Lowther]

Jason,

Thanks for the reply. Since I'm running on Windows I'd have to set
python up, which I'm sure is doable. However, before I do that I'd like
to know if it is possible to mosaic with gdalwarp.

Anyone?

Dave

-----Original Message-----
From: gdal-dev-admin at remotesensing.org
[mailto:gdal-dev-admin at remotesensing.org] On Behalf Of jason
Sent: Thursday, June 03, 2004 5:20 PM
To: gdal-dev at remotesensing.org
Subject: [Gdal-dev] Re: Mosaic in GDAL Warp

On Thu, 03 Jun 2004 16:33:18 -0500, David W. Lowther wrote:

>  
> 
> Hi List,
> 
>  
> 
> I'm trying to use the mosaic with gdal warp. The docs say "Mosaicing
> into an existing output file is supported if the output file already
> exists."
> 
>  
> 
> I have four files (quarter quads). I'd like them mosaiced and
> reprojected - either in one operation or more. I can make reproject
> happen no prob for one file.
> 
>  
> 
> Here are some failed attempts to do this (note 109999 is custom
defined
> OK Albers Projection):
> 
>  
> 
> E:\GDAL\OpenEV_FW\data>..\bin\gdalwarp -s_srs EPSG:26715 -t_srs
> EPSG:109999 -of GTIFF -co 'TFW=YES' e:\project\corps\38095a21.tif
> e:\project\corps\out\a38095.tif
> 
> Creating output file is that 6440P x 7932L.
> 
> :0...10...20...30...40...50...60...70...80...90...100 - done.
> 
>  
> 
> E:\GDAL\OpenEV_FW\data>..\bin\gdalwarp -s_srs EPSG:26715 -t_srs
> EPSG:109999 -of GTIFF -co 'TFW=YES' e:\project\corps\38095a22.tif
> e:\project\corps\out\a38095.tif
> 
> Output dataset e:\project\corps\out\a38095.tif exists, but some
> commandline options were provided indicating a new dataset should be
> created.  Please delete existing dataset and run again.
> 
>  
> 
> E:\GDAL\OpenEV_FW\data>..\bin\gdalwarp -s_srs EPSG:26715 -t_srs
> EPSG:109999 -co 'TFW=YES' e:\project\corps\38095a22.tif
> e:\project\corps\out\a38095.tif
> 
> Output dataset e:\project\corps\out\a38095.tif exists, but some
> commandline options were provided indicating a new dataset should be
> created.  Please delete existing dataset and run again.
> 
>  
> 
> E:\GDAL\OpenEV_FW\data>..\bin\gdalwarp -s_srs EPSG:26715 -co 'TFW=YES'
> e:\project\corps\38095a22.tif e:\project\corps\out\a38095.tif
> 
> Output dataset e:\project\corps\out\a38095.tif exists, but some
> commandline options were provided indicating a new dataset should be
> created.  Please delete existing dataset and run again.
> 
>  
> 
> E:\GDAL\OpenEV_FW\data>..\bin\gdalwarp -s_srs EPSG:26715
> e:\project\corps\38095a
> 
> 22.tif e:\project\corps\out\a38095.tif
> 
> :0...10...20...30...40...50...60...70...80...90...100 - done.
> 
>  
> 
> E:\GDAL\OpenEV_FW\data>..\bin\gdalwarp -s_srs EPSG:26715
> e:\project\corps\38095a
> 
> 23.tif e:\project\corps\out\a38095.tif
> 
> :0...10...20...30...40...50...60...70...80...90...100 - done.
> 
>  
> 
> E:\GDAL\OpenEV_FW\data>..\bin\gdalwarp -s_srs EPSG:26715
> e:\project\corps\38095a
> 
> 24.tif e:\project\corps\out\a38095.tif
> 
> :0...10...20...30...40...50...60...70...80...90...100 - done.
> 
>  
> 
>  
> 
> Thanks in advance for any help anyone can provide.
> 
>  
> 
>  
> 
> David Lowther
> 
> Senior Software Engineer
> 
> Strategic Consulting International
> 
> (405) 227-9051
> 
> dlowther at scigis.com

Have you checked out the script gdal_merge?

[root at localhost ogr]# gdal_merge.py 
Usage: gdal_merge.py [-o out_filename] [-f out_format] [-co NAME=VALUE]*
                     [-ps pixelsize_x pixelsize_y] [-separate] [-v]
[-pct]
                     [-ul_lr ulx uly lrx lry] [-n nodata_value]
                     input_files

I regularly use it to splice together and chop out a combined section of
georeferenced TIFFs.

_______________________________________________
Gdal-dev mailing list
Gdal-dev at remotesensing.org
http://remotesensing.org/mailman/listinfo/gdal-dev