<div dir="ltr">Hello,<div><br></div><div> I have created a virtual raster of DEM data covering approximately 72,000 km^2. The DEM has a 1 metre pixel. The source for the virtual raster is about 95 geoTIFF files with irregular shapes. I have a SHP file with a 10 km grid of tiles, with a total of 853 polygons (see image below). I am trying to cut the virtual raster into 853 tiles based on the SHP file. Each tile will be 10 km x 10 km. The tiles in the SHP file are exactly 10,000 metres x 10,000 metres. So I expect the DEM tiles to be exactly 10,000 pixels x 10,000 pixels. However, when I run the QGIS Easy Raster Splitter plugin the output geoTIFF DEM tiles are 9,999 pixels x 10,000 pixels. I have checked all 853 tiles and they are all 9,999 x 10,000 pixels. Here are the first few lines reported by gdalinfo. </div><div><br></div><blockquote style="margin:0 0 0 40px;border:none;padding:0px"><div><i>Driver: GTiff/GeoTIFF</i></div><div><i>Files: nb_2015_16_18_dem_2460000_7480000.tif</i></div><div><i style="background-color:rgb(255,255,0)">Size is 9999, 10000</i></div><div><i>Coordinate System is:</i></div></blockquote><div style="font-style:italic"><i><br></i></div> Here is an image of the 1 metre gap between the 10 km tiles. I set the scale in QGIS to 1:1,000 for this image. <div><br><div><img src="cid:ii_m2zaqkdi2" alt="image.png" width="186" height="193"></div><div><br></div><div> These 1 pixel (1 metre) gaps appear as vertical gaps with a 10 km spacing when several tiles are displayed together.</div><div><br></div><div> Both the SHP file and the source DEM GeoTIFFs all use the same projected coordinate system - EPSG 2953. The
Easy Raster Splitter plugin uses gdalwarp in the background. I suspect gdalwarp is doing some kind of reprojection when it calculates the intersection between the individual tiles in the SHP and the source VRT. But I am unsure. Does anybody have a suggestion?<br><div><div><img src="cid:ii_m2z9cldk1" alt="image.png" width="562" height="451"><br></div><span class="gmail_signature_prefix"><div>Thanks,</div><div><span class="gmail_signature_prefix">Bernie.</span></div>-- </span><br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Bernie Connors, P.Eng<br></div><div>B.Sc Geodesy & Geomatics Engineering</div><div>Fredericton, NB</div><div><a href="mailto:bernie.connors@unb.ca" target="_blank">bernie.connors@unb.ca</a></div></div></div></div></div></div></div>