Hey Andrew,<br><br>Are your files in the same projection? I know gshhs is projected in wgs84 by default, and I see your raster is in British National Grid. AFAIK, gdal_rasterize assumes the vector & raster data sets are in the same projection. <br>
<br>-Jamie<br><br><div class="gmail_quote">On Mon, Nov 10, 2008 at 11:11 AM, Andrew Brooks <span dir="ltr"><<a href="mailto:arb@sat.dundee.ac.uk">arb@sat.dundee.ac.uk</a>></span> wrote:<br><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hello<br>
<br>
I'm trying to overlay coastlines on an image using gdal_rasterize but<br>
can't get it to do anything at all.<br>
<br>
% gdal_rasterize gshhs_land.shp warped.tif<br>
Usage: gdal_rasterize [-b band] [-i]<br>
[-burn value] | [-a attribute_name] | [-3d]<br>
[-l layername]* [-where expression] [-sql select_statement]<br>
<src_datasource> <dst_filename><br>
<br>
The files are ok:<br>
<br>
% ogrinfo gshhs_land.shp<br>
INFO: Open of `gshhs_land.shp' using driver `ESRI Shapefile' successful.<br>
1: gshhs_land (Polygon)<br>
<br>
% gdalinfo warped.tif<br>
Driver: GTiff/GeoTIFF<br>
Size is 396, 396<br>
Coordinate System is: PROJCS["OSGB 1936 / British National Grid", ...<br>
Band 1 Block=396x10 Type=UInt16, ColorInterp=Gray<br>
<br>
Andrew<br>
<br>
P.S. Unrelated to this problem there's also another problem:<br>
<br>
% gdal_rasterize -l gshhs_land gshhs_land.shp warped.tif<br>
To many of -3d, -burn and -a specified.<br>
<br>
"To" should be "Too", but since I've not specified any of -3d,<br>
-burn or -a then I'm not sure what it's complaining about...<br>
<br>
_______________________________________________<br>
gdal-dev mailing list<br>
<a href="mailto:gdal-dev@lists.osgeo.org">gdal-dev@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/mailman/listinfo/gdal-dev" target="_blank">http://lists.osgeo.org/mailman/listinfo/gdal-dev</a><br>
</blockquote></div><br>