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 <[email protected]> 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 > [email protected] > http://lists.osgeo.org/mailman/listinfo/gdal-dev > -- http://blog.klokan.cz/ http://www.maptiler.org/ http://www.oldmapsonline.org/ _______________________________________________ gdal-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/gdal-dev
