<div dir="ltr"><div>Hi Even,</div><div><br></div>The original file is a Landsat-8 image, I can confirm that it has tags set in the imd metadata domain with gdalinfo -mdd imd LC80040242014245LGN00_B3.TIF but I'm creating a brand new raster, not making a copy, so I'm not sure why the metadata is being transferred. The following basic example still has the .IMD file created.<div><div><br></div><div><div>import numpy</div><div>from osgeo import gdal</div><div><br></div><div>band = "C:/Users/deschampsb/Desktop/test/LC80040242014245LGN00/LC80040242014245LGN00_B3.TIF"</div><div>outpath = "C:/Users/deschampsb/Desktop/test/LC80040242014245LGN00/LC80040242014245LGN00_B3_test.TIF"</div><div><br></div><div>ds = gdal.Open(band)</div><div>proj = ds.GetProjection()</div><div>geot = ds.GetGeoTransform()</div><div>dn = ds.ReadAsArray().astype(numpy.float32)</div><div>ds = None</div><div><br></div><div>driver = gdal.GetDriverByName("GTiff")</div><div>ds_out = driver.Create(outpath, dn.shape[1], dn.shape[0], 1, gdal.GDT_Float32)</div><div>ds_out.GetRasterBand(1).WriteArray(dn)</div><div>ds_out.SetProjection(proj)</div><div>ds_out.SetGeoTransform(geot)</div><div>ds_out.GetRasterBand(1).SetNoDataValue(0)</div><div>ds_out = None</div></div><div><br></div><div>However, if I write to a different directory, or more curiously to a name with a different prefix, it does not create the file.</div></div><div><br></div><div><div>outpath = "C:/Users/deschampsb/Desktop/test/LC80040242014245LGN00/newname.TIF"</div></div><div><div>outpath = "C:/Users/deschampsb/Desktop/test/new/LC80040242014245LGN00_B3_test.TIF"</div></div><div><br></div><div>Is it looking up other files in the directory?</div><div><br></div><div>Thanks,</div><div><br></div><div>Benjamin</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jul 14, 2016 at 4:28 PM, Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@spatialys.com" target="_blank">even.rouault@spatialys.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Benjamin,<br>
<span class=""><br>
><br>
> I am seeing .IMD files being created first by driver.Create() with the<br>
> GTIFF driver, and then also by gdalwarp (doing a simple reprojection), with<br>
> both input and output in GTIFF format.<br>
><br>
> I've never seen these files before... Can I prevent them from being<br>
> written?<br>
<br>
</span>If they are created, it means that the IMD metadata domain is set on your<br>
GeoTIFF files. Presumably because it is also set in the source images you use.<br>
There's no option to suppress writing the IMD files when the IMD metadata<br>
domain is set on a GeoTIFF file. So either make sure it is not set, or suppress<br>
the .imd as a postprocessing.<br>
<br>
But just doing a plain driver.Create() without setting metadata should not<br>
result in a .imd file being written.<br>
<span class=""><br>
><br>
> I'm using GDAL 2.0.2 with the Python bindings, and GDAL 2.1.0 for the<br>
> utilities.<br>
<br>
</span>Odd setup. I presume you know what you are doing ;-)<br>
<span class="HOEnZb"><font color="#888888"><br>
Even<br>
<br>
--<br>
Spatialys - Geospatial professional services<br>
<a href="http://www.spatialys.com" rel="noreferrer" target="_blank">http://www.spatialys.com</a><br>
</font></span></blockquote></div><br></div>