[gdal-dev] gdal2tiles.py and ProjNatOriginLongGeoKey

[Klokan Petr Přidal]

Zachary,

do you really need a huge big from your files?

In case you need to use your merged data for generating tiles by
gdal2tiles (I guess you do because you posted your message to this
thread) you can use GDAL VRT driver for merging - which is just XML
file describing where is which file located in the huge result raster.

VRT merging can be done simply by writing your own XML, format is
described here:
http://www.gdal.org/gdal_vrttut.html

You can also use one of the existing utilities:
- gdalbuildvrt from GDAL 1.7 dev (GDAL svn trunk)
- gdal_vrtmerge
(http://svn.osgeo.org/gdal/trunk/gdal/swig/python/samples/gdal_vrtmerge.py)
- preprocess.py from MapTiler SVN
(http://code.google.com/p/maptiler/source/browse/trunk/maptiler/preprocess.py)

Best,

Klokan

On Mon, Apr 13, 2009 at 10:36 PM, Zachary Fisk <zacharyfisk at gmail.com> wrote:
>  I'm trying to merge about 345 tifs NED 10M data for the continental U.S. I
> processed the data from ArcInfo Grid to tif (see bottom for the output of
> one of the tifs from gdalinfo) using gdal_translate. The error below is
> generated from the attempted merge, how I do build GDAL with BigTiff
> support? Also I noticed my windows version of the tool supports the
> --optfile parameter, but not the linux python script version. Can this
> command handle all 345 files on one command? Are there other best practices
> methods that I should be aware of, I prefer to process on the linux system
> since they, but am open to processing on windows. Whatever works.
>
>
> Sample command parameters I'm passing to gdal_merge.py
> gdal_merge.py -o big2.tif -of GTIFF -co BIGTIFF=YES -co TILE=YES -co
> COMPRESS=DEFLATE -co PREDICTOR=2 -co BLOCKXSIZE=128 -co BLOCKYSIZE=128
> demgrid_101_1.tif demgrid_101_2.tif demgrid_101_3.tif demgrid_102_1.tif
> demgrid_102_2.tif demgrid_102_3.tif demgrid_103_1.tif demgrid_103_2.tif
> demgrid_103_3.tif demgrid_104_1.tif
> Warning 6: BIGTIFF requested, but GDAL built without BigTIFF
> enabled libtiff, request ignored.
> Segmentation fault
>
>
>
> Output of one of the images passed to the gdal_merge.py from gdalinfo
> gdalinfo demgrid_101_1.tif
> Driver: GTiff/GeoTIFF
> Files: demgrid_101_1.tif
> Size is 32400, 10800
> Coordinate System is:
> GEOGCS["NAD83",
>     DATUM["North_American_Datum_1983",
>         SPHEROID["GRS 1980",6378137,298.2572221010042,
>             AUTHORITY["EPSG","7019"]],
>         AUTHORITY["EPSG","6269"]],
>     PRIMEM["Greenwich",0],
>     UNIT["degree",0.0174532925199433],
>     AUTHORITY["EPSG","4269"]]
> Origin = (-120.000000000014793,47.000000000075204)
> Pixel Size = (0.000092592592600,-0.000092592592600)
> Metadata:
>   AREA_OR_POINT=Area
> Image Structure Metadata:
>   INTERLEAVE=BAND
> Corner Coordinates:
> Upper Left  (-120.0000000,  47.0000000) (120d 0'0.00"W, 47d 0'0.00"N)
> Lower Left  (-120.0000000,  46.0000000) (120d 0'0.00"W, 46d 0'0.00"N)
> Upper Right (-117.0000000,  47.0000000) (117d 0'0.00"W, 47d 0'0.00"N)
> Lower Right (-117.0000000,  46.0000000) (117d 0'0.00"W, 46d 0'0.00"N)
> Center      (-118.5000000,  46.5000000) (118d30'0.00"W, 46d30'0.00"N)
> Band 1 Block=32400x1 Type=Int16, ColorInterp=Gray
>   NoData Value=-3.4028234663852886e+38
>   Metadata:
>     LAYER_TYPE=athematic
>
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>



-- 
http://blog.klokan.cz/
http://www.maptiler.org/
http://www.oldmapsonline.org/