Hi All,<div><br></div><div>Merging 2 files with a geographical CRS results in the output being shifted a quarter-pixel to the east and south. The datum is WGS84 (altough not explicitly named), I get the same result if I explicitly set it to WGS84.</div>
<div>SImilar data in modis sinusoidal projection gets merged fine with the same command.</div><div><br></div><div>I am using the following command:</div><div><div>gdalwarp -srcnodata 0 -dstnodata 0 -co COMPRESS=LZW MCD45monthly.burndate.Win06.2010-byte.tif MCD45monthly.burndate.Win05.2010-byte.tif MCD45monthly.burndate.sam.2010-Byte.tif</div>
<div><br></div><div>Files can be found at: </div><div><a href="http://www.openclimgeo.org/pub/">http://www.openclimgeo.org/pub/</a></div><div><br></div><div><br></div><div>Any help would be appreciated! Here is the gdalinfo for the 3 files.</div>
<div><br></div><div>Etienne</div><div><br></div>
</div><div><div>tourigny@supernova: /data/research/work/modis/mcd45/gtiff/out/upload $ gdalinfo MCD45monthly.burndate.Win05.2010-byte.tif</div><div>Driver: GTiff/GeoTIFF</div><div>Files: MCD45monthly.burndate.Win05.2010-byte.tif</div>
<div>Size is 10923, 5234</div><div>Coordinate System is:</div><div>GEOGCS[,</div><div> DATUM["unknown",</div><div> SPHEROID["unnamed",6378137,298.257223563]],</div><div> PRIMEM["Greenwich",0],</div>
<div> UNIT["degree",0.0174532925199433]]</div><div>Origin = (-81.997802734374957,12.997802734374995)</div><div>Pixel Size = (0.004394531250000,-0.004394531250000)</div><div>Metadata:</div><div> AREA_OR_POINT=Area</div>
<div>Image Structure Metadata:</div><div> COMPRESSION=LZW</div><div> INTERLEAVE=BAND</div><div>Corner Coordinates:</div><div>Upper Left ( -81.9978027, 12.9978027) ( 81d59'52.09"W, 12d59'52.09"N)</div>
<div>Lower Left ( -81.9978027, -10.0031738) ( 81d59'52.09"W, 10d 0'11.43"S)</div><div>Upper Right ( -33.9963379, 12.9978027) ( 33d59'46.82"W, 12d59'52.09"N)</div><div>Lower Right ( -33.9963379, -10.0031738) ( 33d59'46.82"W, 10d 0'11.43"S)</div>
<div>Center ( -57.9970703, 1.4973145) ( 57d59'49.45"W, 1d29'50.33"N)</div><div>Band 1 Block=10923x1 Type=Byte, ColorInterp=Gray</div><div> NoData Value=0</div><div>tourigny@supernova: /data/research/work/modis/mcd45/gtiff/out/upload $ gdalinfo MCD45monthly.burndate.Win06.2010-byte.tif</div>
<div>Driver: GTiff/GeoTIFF</div><div>Files: MCD45monthly.burndate.Win06.2010-byte.tif</div><div>Size is 10240, 5689</div><div>Coordinate System is:</div><div>GEOGCS[,</div><div> DATUM["unknown",</div><div> SPHEROID["unnamed",6378137,298.257223563]],</div>
<div> PRIMEM["Greenwich",0],</div><div> UNIT["degree",0.0174532925199433]]</div><div>Origin = (-78.997802734374972,-10.002197265624996)</div><div>Pixel Size = (0.004394531250000,-0.004394531250000)</div>
<div>Metadata:</div><div> AREA_OR_POINT=Area</div><div>Image Structure Metadata:</div><div> COMPRESSION=LZW</div><div> INTERLEAVE=BAND</div><div>Corner Coordinates:</div><div>Upper Left ( -78.9978027, -10.0021973) ( 78d59'52.09"W, 10d 0' 7.91"S)</div>
<div>Lower Left ( -78.9978027, -35.0026855) ( 78d59'52.09"W, 35d 0' 9.67"S)</div><div>Upper Right ( -33.9978027, -10.0021973) ( 33d59'52.09"W, 10d 0' 7.91"S)</div><div>Lower Right ( -33.9978027, -35.0026855) ( 33d59'52.09"W, 35d 0' 9.67"S)</div>
<div>Center ( -56.4978027, -22.5024414) ( 56d29'52.09"W, 22d30' 8.79"S)</div><div>Band 1 Block=10240x1 Type=Byte, ColorInterp=Gray</div><div> NoData Value=0</div><div>tourigny@supernova: /data/research/work/modis/mcd45/gtiff/out/upload $ gdalinfo MCD45monthly.burndate.sam.2010-Byte.tif </div>
<div>Driver: GTiff/GeoTIFF</div><div>Files: MCD45monthly.burndate.sam.2010-Byte.tif</div><div>Size is 10923, 10923</div><div>Coordinate System is:</div><div>GEOGCS[,</div><div> DATUM["unknown",</div><div> SPHEROID["unnamed",6378137,298.257223563]],</div>
<div> PRIMEM["Greenwich",0],</div><div> UNIT["degree",0.0174532925199433]]</div><div>Origin = (-81.997802734374957,12.997802734374995)</div><div>Pixel Size = (0.004394531250000,-0.004394531250000)</div>
<div>Metadata:</div><div> AREA_OR_POINT=Area</div><div>Image Structure Metadata:</div><div> COMPRESSION=LZW</div><div> INTERLEAVE=BAND</div><div>Corner Coordinates:</div><div>Upper Left ( -81.9978027, 12.9978027) ( 81d59'52.09"W, 12d59'52.09"N)</div>
<div>Lower Left ( -81.9978027, -35.0036621) ( 81d59'52.09"W, 35d 0'13.18"S)</div><div>Upper Right ( -33.9963379, 12.9978027) ( 33d59'46.82"W, 12d59'52.09"N)</div><div>Lower Right ( -33.9963379, -35.0036621) ( 33d59'46.82"W, 35d 0'13.18"S)</div>
<div>Center ( -57.9970703, -11.0029297) ( 57d59'49.45"W, 11d 0'10.55"S)</div><div>Band 1 Block=10923x1 Type=Byte, ColorInterp=Gray</div><div> NoData Value=0</div></div><div><br></div>