<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
</head>
<body bgcolor="#ffffff" text="#000000">
<div class="moz-text-html" lang="x-western">Hi,<br>
<br>
I am trying to use VRT files and gdalwarp to mosaic and tile
multiresolution images. I use two images at different resolutions
which overlap and compute the mosaic using gdalbuildvrt. <br>
The vrt produced has wrong "y resolution and y size", so I first edit
the y coordinates of the vrt (ticket 3198). Then, I use gdalwarp with
cubic interpolation to cut a portion of the mosaic where the input
images overlap. <br>
The image produced shows a slight offset in position and scaling
between the two input images in the mosaic and no cubic interpolation
of the background image.<br>
<br>
I run the following commands:<br>
<b>gdalbuildvrt mosaic.vrt -resolution highest imageL00.jpg imageL06.jpg<br>
gdalwarp mosaic.vrt -r CUBIC -te 0 0 0.03 0.03 cut_vrt.jpg</b><br>
<br>
The same calculation works perfectly when using directly gdalwarp
without vrt:<br>
<br>
<b>gdalwarp imageL00.jpg imageL06.jpg -r CUBIC -te 0 0 0.03 0.03
cut_warp.jpg</b><br>
<br>
Any idea of what I am doing wrong?<br>
<br>
Many thanks,<br>
<br>
Christophe.<br>
<br>
More details on the images used:<br>
<br>
D:\Geo_Images\multi_Resolution_Data\test\bug_report>gdalinfo
imageL06.jpg<br>
Driver: JPEG/JPEG JFIF<br>
Files: imageL06.jpg<br>
imageL06.jpg.aux.xml<br>
Size is 1024, 1024<br>
Coordinate System is:<br>
GEOGCS["WGS 84",<br>
DATUM["WGS_1984",<br>
SPHEROID["WGS 84",6378137,298.257223563,<br>
AUTHORITY["EPSG","7030"]],<br>
AUTHORITY["EPSG","6326"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4326"]]<br>
Origin = (0.000000000000000,0.015625000000000)<br>
Pixel Size = (0.000015258789063,-0.000015258789063)<br>
Metadata:<br>
AREA_OR_POINT=Area<br>
Image Structure Metadata:<br>
SOURCE_COLOR_SPACE=YCbCr<br>
INTERLEAVE=PIXEL<br>
COMPRESSION=JPEG<br>
Corner Coordinates:<br>
Upper Left ( 0.0000000, 0.0156250) ( 0d 0'0.01"E, 0d 0'56.25"N)<br>
Lower Left ( 0.0000000, 0.0000000) ( 0d 0'0.01"E, 0d 0'0.01"N)<br>
Upper Right ( 0.0156250, 0.0156250) ( 0d 0'56.25"E, 0d 0'56.25"N)<br>
Lower Right ( 0.0156250, 0.0000000) ( 0d 0'56.25"E, 0d 0'0.01"N)<br>
Center ( 0.0078125, 0.0078125) ( 0d 0'28.13"E, 0d 0'28.13"N)<br>
Band 1 Block=1024x1 Type=Byte, ColorInterp=Red<br>
Metadata:<br>
LAYER_TYPE=athematic<br>
Image Structure Metadata:<br>
COMPRESSION=JPEG<br>
Band 2 Block=1024x1 Type=Byte, ColorInterp=Green<br>
Metadata:<br>
LAYER_TYPE=athematic<br>
Image Structure Metadata:<br>
COMPRESSION=JPEG<br>
Band 3 Block=1024x1 Type=Byte, ColorInterp=Blue<br>
Metadata:<br>
LAYER_TYPE=athematic<br>
Image Structure Metadata:<br>
COMPRESSION=JPEG<br>
<br>
D:\Geo_Images\multi_Resolution_Data\test\bug_report>gdalinfo
imageL00.jpg<br>
Driver: JPEG/JPEG JFIF<br>
Files: imageL00.jpg<br>
imageL00.jpg.aux.xml<br>
Size is 1024, 1024<br>
Coordinate System is:<br>
GEOGCS["WGS 84",<br>
DATUM["WGS_1984",<br>
SPHEROID["WGS 84",6378137,298.257223563,<br>
AUTHORITY["EPSG","7030"]],<br>
AUTHORITY["EPSG","6326"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4326"]]<br>
Origin = (0.000000000000000,1.000000000000000)<br>
Pixel Size = (0.000976562500000,-0.000976562500000)<br>
Metadata:<br>
AREA_OR_POINT=Area<br>
Image Structure Metadata:<br>
SOURCE_COLOR_SPACE=YCbCr<br>
INTERLEAVE=PIXEL<br>
COMPRESSION=JPEG<br>
Corner Coordinates:<br>
Upper Left ( 0.0000000, 1.0000000) ( 0d 0'0.01"E, 1d 0'0.00"N)<br>
Lower Left ( 0.0000000, 0.0000000) ( 0d 0'0.01"E, 0d 0'0.01"N)<br>
Upper Right ( 1.0000000, 1.0000000) ( 1d 0'0.00"E, 1d 0'0.00"N)<br>
Lower Right ( 1.0000000, 0.0000000) ( 1d 0'0.00"E, 0d 0'0.01"N)<br>
Center ( 0.5000000, 0.5000000) ( 0d30'0.00"E, 0d30'0.00"N)<br>
Band 1 Block=1024x1 Type=Byte, ColorInterp=Red<br>
Metadata:<br>
LAYER_TYPE=athematic<br>
Image Structure Metadata:<br>
COMPRESSION=JPEG<br>
Band 2 Block=1024x1 Type=Byte, ColorInterp=Green<br>
Metadata:<br>
LAYER_TYPE=athematic<br>
Image Structure Metadata:<br>
COMPRESSION=JPEG<br>
Band 3 Block=1024x1 Type=Byte, ColorInterp=Blue<br>
Metadata:<br>
LAYER_TYPE=athematic<br>
Image Structure Metadata:<br>
COMPRESSION=JPEG<br>
<br>
D:\Geo_Images\multi_Resolution_Data\test\bug_report>gdalinfo
mosaic.vrt<br>
Driver: VRT/Virtual Raster<br>
Files: mosaic.vrt<br>
Size is 65536, 65536<br>
Coordinate System is:<br>
GEOGCS["WGS 84",<br>
DATUM["WGS_1984",<br>
SPHEROID["WGS 84",6378137,298.257223563,<br>
AUTHORITY["EPSG","7030"]],<br>
AUTHORITY["EPSG","6326"]],<br>
PRIMEM["Greenwich",0],<br>
UNIT["degree",0.0174532925199433],<br>
AUTHORITY["EPSG","4326"]]<br>
Origin = (0.000000000000000,1.000000000000000)<br>
Pixel Size = (0.000015258789063,-0.000015258789063)<br>
Corner Coordinates:<br>
Upper Left ( 0.0000000, 1.0000000) ( 0d 0'0.01"E, 1d 0'0.00"N)<br>
Lower Left ( 0.0000000, 0.0000000) ( 0d 0'0.01"E, 0d 0'0.01"N)<br>
Upper Right ( 1.0000000, 1.0000000) ( 1d 0'0.00"E, 1d 0'0.00"N)<br>
Lower Right ( 1.0000000, 0.0000000) ( 1d 0'0.00"E, 0d 0'0.01"N)<br>
Center ( 0.5000000, 0.5000000) ( 0d30'0.00"E, 0d30'0.00"N)<br>
Band 1 Block=128x128 Type=Byte, ColorInterp=Red<br>
Band 2 Block=128x128 Type=Byte, ColorInterp=Green<br>
Band 3 Block=128x128 Type=Byte, ColorInterp=Blue<br>
</div>
<br>
<br>
</body>
</html>