Opened 10 years ago

Closed 10 years ago

#3432 closed defect (fixed)

gdalbuildvrt does not support positive NS resolution when based off pixel coordinates

Reported by: abarton Owned by: warmerdam
Priority: normal Milestone: 1.8.0
Component: default Version: 1.6.3
Severity: normal Keywords:
Cc: ysiddiqui@…

Description

When coordinates are based on pixels which increase positively as you go down, this causes an error in gdalbuildvrt because it can't handle positive NS resolutions.

gdalbuildvrt does not support positive NS resolution.

Attached is a sample file that demonstrates the issue:

e.g. gdalbuildvrt blank_image.vrt blank_image.tif

Attachments (3)

blank_image.tif (29.5 KB) - added by abarton 10 years ago.
SV20091120_b0.hdf.txt (376 bytes) - added by abarton 10 years ago.
SV20091120_b2.hdf.txt (376 bytes) - added by abarton 10 years ago.

Download all attachments as: .zip

Change History (8)

Changed 10 years ago by abarton

Attachment: blank_image.tif added

comment:1 Changed 10 years ago by Even Rouault

What is the use case of this ? Do you have georeferenced images that can be mosaiced and have positive NS resolution ? Your sample is just a non-georeferenced image, which makes no sense to try mosaicing

comment:2 Changed 10 years ago by abarton

I have a number of VRT bands that compose one imagery product. I'm trying to combine the bands and then geo-reference them. Currently I must geo-reference each band before combining them which introduces complexity and overhead.

comment:3 Changed 10 years ago by Even Rouault

I'm trying to understand your last message, but I'm not sure to have gotten it. My interpretation is that you'll plan to use the -separate option of gdalbuildvrt to combine several ungeoreferenced single-band images into a VRT and then georeference the VRT ? Or is it something else ?

Maybe you could give example of the command lines you would use if gdalbuildvrt did what you want, and give the output of gdalinfo on a few of your input images.

comment:4 Changed 10 years ago by abarton

You are correct, I am trying to combine several non-georeferenced single-band images into a VRT and then georeference the VRT.

I'm attaching a couple of text files with output from gdalinfo on the input images.

My basic use case would be :

gdalbuildvrt -separate <output_vrt> <input_vrt_files>

Changed 10 years ago by abarton

Attachment: SV20091120_b0.hdf.txt added

Changed 10 years ago by abarton

Attachment: SV20091120_b2.hdf.txt added

comment:5 Changed 10 years ago by Even Rouault

Milestone: 1.8.0
Resolution: fixed
Status: newclosed

r18888 /trunk/gdal/apps/gdalbuildvrt.cpp: gdalbuildvrt : support stacking ungeoreferenced images when using -separate, provided they have the same size (#3432)

r18889 /trunk/autotest/utilities/test_gdalbuildvrt.py: Test gdalbuildvrt -separate with ungeoreferenced images (#3432Ã)

Note: See TracTickets for help on using tickets.