Opened 15 years ago
Closed 10 years ago
#3328 closed enhancement (invalid)
gdalwarp doesn't support producing VRTs containing multiple source files
Reported by: | EliL | Owned by: | warmerdam |
---|---|---|---|
Priority: | lowest | Milestone: | |
Component: | Utilities | Version: | unspecified |
Severity: | trivial | Keywords: | |
Cc: |
Description
This is mostly copied from the email list (http://lists.osgeo.org/pipermail/gdal-dev/2010-January/023194.html)
Trying to mosaic multiple source files with gdalwarp to vrt format doesn't work properly.
gdalwarp -t_srs "EPSG:2913" -s_srs "EPSG:2994" -srcnodata "0 0 0" -dstnodata "0 0 0" -of VRT be44123g8b be44123g8c be44123g8a test2.vrt
Creating output file that is 11039P x 15191L. Processing input file be44123g8b.
test2.vrt looks fine with gdalinfo but not in OpenEv. Also, when I look in the .vrt file itself, I only see one reference to the 1 grid and not to the other 2 grids:
<SourceDataset relativeToVRT="1">be44123g8b</SourceDataset>
The same command using tif format works as expected:
gdalwarp -t_srs "EPSG:2913" -s_srs "EPSG:2994" -srcnodata "0 0 0" -dstnodata "0 0 0" -of GTIFF be44123g8b be44123g8c be44123g8a test2.tif
Creating output file that is 11039P x 15191L. Processing input file be44123g8b. 0...10...20...30...40...50...60...70...80...90...100 - done. Processing input file be44123g8c. 0...10...20...30...40...50...60...70...80...90...100 - done. Processing input file be44123g8a. 0...10...20...30...40...50...60...70...80...90...100 - done.
test2.tif looks fine with gdalinfo and in OpenEV.
I also notice that when processing to tif each of the 3 grids are listed. When processing to .vrt only 1 grid is listed.
This doesn't seem very important since Frank's Advice ("Or if the source files are all in the same coordinate system, perhaps you could build a mosaic vrt, and then use gdalwarp to produce a corresponding warped vrt.") worked quite well:
gdalbuildvrt gdalbuildvrt.vrt be44123g8a be44123g8b be44123g8c be44123h7 be44123h8a be44123h8b be44124g1 be44124h1 be45123a7 be45123a8 be45124a1
gdalwarp -t_srs "EPSG:2913" -s_srs "EPSG:2994" -of VRT gdalbuildvrt.vrt gdalbuildvrt2994.vrt
Change History (5)
comment:1 by , 15 years ago
comment:2 by , 13 years ago
Looking back on this, perhaps this should be closed as invalid or wontfix since the functionality is easily achieved with a slightly different approach. What do you think? Is there a good reason that gdalwarp should support -of VRT when the input is not VRT?
comment:3 by , 13 years ago
Yes we can close it. However your last sentence is not correct: gdalwarp does support warping something that is not a VRT into a VRT. What is not support is wrapping multiple non-VRT into VRT
comment:4 by , 10 years ago
Conclusion 3 years ago was that the ticket can be closed. Suggestion: close it now.
comment:5 by , 10 years ago
Resolution: | → invalid |
---|---|
Status: | new → closed |
Doesn't really improve the situation, but at least the user is warned :
r18564 /trunk/gdal/apps/gdalwarp.cpp: gdalwarp: print a warning message when trying to specify multiple source datasets as input of gdalwarp -of VRT