Hi, I&#39;m trying to stitch several map tiles together into a big image and then warp it. However I&#39;m getting error like this from gdalbuildvrt...<br><br>Warning 6: gdalbuildvrt does not support positive NS resolution. Skipping map1.tiff<br>
<br>Searches for this error suggest its got something to do with north not being up. But its not clear what this means. (My maps are in the northern hemisphere).<br><br>I wondered if I am georeferencing the maps incorrectly. Am I right is saying that top-left is (0,0) and the y-axis points downward for pixel coordinates within the map? (I cant find anywhere where this is documented.) I also tried y-axis pointing upward (giving -ve y-values) as this is what quantumGIS seems to use.<br>
<br><br>COMMANDS I USE....<br><br>gdal_translate.exe -a_srs EPSG:27700 <br>-gcp 0 400 241538.5 134628.66666667 <br>-gcp 600 0 244713.5 136745.33333333 <br>-gcp 0 0 241538.5 136745.33333333 <br>-gcp 600 400 244713.5 134628.66666667 <br>
./map1.png ./map1.tiff<br><br>gdal_translate.exe -a_srs EPSG:27700 <br>-gcp 0 400 241538.5 132512 <br>-gcp 600 0 244713.5 134628.66666667 <br>-gcp 0 0 241538.5 134628.66666667 <br>-gcp 600 400 244713.5 132512 <br>./map2.png ./map2.tiff<br>
<br>gdalbuildvrt -overwrite maps_merged.vrt map1.tiff map2.tiff<br><br><br>OUTPUT OF GDALINFO ON MAP1.TIFF....<br><br>Driver: GTiff/GeoTIFF<br>Files: map1.tiff<br>Size is 600, 400<br>Coordinate System is `&#39;<br>GCP Projection =<br>
PROJCS[&quot;OSGB 1936 / British National Grid&quot;,<br>    GEOGCS[&quot;OSGB 1936&quot;,<br>        DATUM[&quot;OSGB_1936&quot;,<br>            SPHEROID[&quot;Airy 1830&quot;,6377563.396,299.3249646000044,<br>                AUTHORITY[&quot;EPSG&quot;,&quot;7001&quot;]],<br>
            AUTHORITY[&quot;EPSG&quot;,&quot;6277&quot;]],<br>        PRIMEM[&quot;Greenwich&quot;,0],<br>        UNIT[&quot;degree&quot;,0.0174532925199433],<br>        AUTHORITY[&quot;EPSG&quot;,&quot;4277&quot;]],<br>    PROJECTION[&quot;Transverse_Mercator&quot;],<br>
    PARAMETER[&quot;latitude_of_origin&quot;,49],<br>    PARAMETER[&quot;central_meridian&quot;,-2],<br>    PARAMETER[&quot;scale_factor&quot;,0.9996012717],<br>    PARAMETER[&quot;false_easting&quot;,400000],<br>    PARAMETER[&quot;false_northing&quot;,-100000],<br>
    UNIT[&quot;metre&quot;,1,<br>        AUTHORITY[&quot;EPSG&quot;,&quot;9001&quot;]],<br>    AUTHORITY[&quot;EPSG&quot;,&quot;27700&quot;]]<br>GCP[  0]: Id=1, Info=<br>          (0,400) -&gt; (241538.5,134628.66666667,0)<br>
GCP[  1]: Id=2, Info=<br>          (600,0) -&gt; (244713.5,136745.33333333,0)<br>GCP[  2]: Id=3, Info=<br>          (0,0) -&gt; (241538.5,136745.33333333,0)<br>GCP[  3]: Id=4, Info=<br>          (600,400) -&gt; (244713.5,134628.66666667,0)<br>
Metadata:<br>  AREA_OR_POINT=Area<br>Image Structure Metadata:<br>  INTERLEAVE=PIXEL<br>Corner Coordinates:<br>Upper Left  (    0.0,    0.0)<br>Lower Left  (    0.0,  400.0)<br>Upper Right (  600.0,    0.0)<br>Lower Right (  600.0,  400.0)<br>
Center      (  300.0,  200.0)<br>Band 1 Block=600x3 Type=Byte, ColorInterp=Red<br>  Mask Flags: PER_DATASET ALPHA<br>Band 2 Block=600x3 Type=Byte, ColorInterp=Green<br>  Mask Flags: PER_DATASET ALPHA<br>Band 3 Block=600x3 Type=Byte, ColorInterp=Blue<br>
  Mask Flags: PER_DATASET ALPHA<br>Band 4 Block=600x3 Type=Byte, ColorInterp=Alpha<br><br><br><br><br><br><br>