<Moses.Gone <at> t-systems.com> writes: > > > > Hi everyone, > I am getting this weird error when I try to convert a vrt to a tif. > > Some context: We are tiling a set of tifs(in a list file). First we generate a vrt from the list. Then we assign a projection on the generated vrt using gdaltranslate. We then clip the output using gdalwarp applying a clipping at the same time to the desired > BBOX(This is achieved using the –te option). In the warp operation we also specify the target SRS(using the –t_srs option). Lastly we tile using our customized gdal_retile.py. All is fine with the projection and clipping, but the resultant tiles are empty-no > data just white images. > > So I tried to backtrack my steps, to check at what point is the data being lost. Generation of the VRT works fine. However, performing gdaltranslate to assign the source SRS on the VRT result into this error: > > bash-2.05$ gdal_translate -of GTiff ffm_dtk200_wdd.vrt test.tif > Input file size is 2147483647, 2147483647 > 0ERROR 1: GTiffCreateCopy():No space for strip arrays > ERROR 1: MissingRequired:TIFF directory is missing required "StripOffsets" field > ERROR 1: MissingRequired:TIFF directory is missing required "StripOffsets" field > > I tried to also convert to AAIGrid format, thinking that something is wrong converting to Tiff, but also got this error: > > bash-2.05$ gdal_translate -of AAIGrid ffm_dtk200_wdd.vrt test3 > Input file size is 2147483647, 2147483647 > 0ERROR 1: CPLMalloc(-4): Silly size requested. > > ERROR 1: ffm_dtk200_wdd.vrt, band 1: The buffer into which the data should be read is null > Segmentation Fault > > Other tiling jobs following a similar method, with different source data works fine though. Could the problem be with the data? If so how comes it worked before?
Hi, Perhaps adding -co BIGTIFF could help. Much better though would be better to skip the whole step of creating a temporary tiff file because it is slow and unnecessary. Use buildvrt with -a_srs and resulting vrt contains projection. Even if it is missing you should be able to tell it directly for gdalwarp as gdalwarp -s_srs epsg:xxx -t_srs epsg:yyyy -te .... input.vrt output.tif -Jukka Rahkonen- _______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
