<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-family:"Calibri","sans-serif";
mso-fareast-language:EN-US;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:3.0cm 2.0cm 3.0cm 2.0cm;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="DA" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">I have the following polygon:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">POLYGON ((539500.00004634191 6314429.5954758488,539499.999 6314429.5964747742,539500.0 6314500.001,539563.49086194986 6314500.001,539574.70556674246 6314491.5008581309,539591.23494657199 6314478.9664945751,539600.001
6314472.1168122729,539600.00014411798 6314460.5601839637,539543.07124903298 6314452.2691510962,539546.30058158806 6314438.619163854,539507.3702219734 6314429.939004397,539505.68004823104 6314429.8589810915,539500.00004634191 6314429.5954758488))<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I have created a memory vector layer and added this polygon as the only feature.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> projection = …<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> geometry = …<o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> ogr_mem_drv = ogr.GetDriverByName(</span><i><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:#00AA00">'Memory'</span></i><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">)</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas"><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas"><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> mem_vector_ds = ogr_mem_drv.CreateDataSource(</span><i><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:#00AA00">''</span></i><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">)</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas"><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> mem_vector_layer = mem_vector_ds.CreateLayer(</span><i><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:#00AA00">''</span></i><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">,
srs=projection)</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas"><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> feature = ogr.Feature(mem_vector_layer.GetLayerDefn())</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas"><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> feature.SetGeometry(geometry.Clone())</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> mem_vector_layer.CreateFeature(feature)<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"><o:p> </o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black;background:yellow;mso-highlight:yellow"> gdal_mem_drv</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">
= gdal.GetDriverByName(</span><i><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:#00AA00">'MEM'</span></i><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">)</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas"><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> mem_raster_ds =
<span style="background:yellow;mso-highlight:yellow">gdal_mem_drv</span>.Create(</span><i><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:#00AA00">''</span></i><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">,5,5,</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:maroon">1</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">,gdal.GDT_Byte)</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas"><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> mem_raster_ds.SetGeoTransform([…])</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas"><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> mem_raster_ds.SetProjection(mem_vector_layer.GetSpatialRef().ExportToWkt())</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas"><o:p></o:p></span></p>
<p class="MsoNormal" style="text-autospace:none"><u><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">
</span></u><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"> err = gdal.RasterizeLayer(mem_raster_ds, [</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:maroon">1</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">],
mem_vector_layer, </span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:blue">None</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">,
</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:blue">None</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">, [</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:maroon">1</span><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">],
[</span><i><u><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:#00AA00">'ALL_TOUCHED=TRUE'</span></u></i><u><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black">])<o:p></o:p></span></u></p>
<p class="MsoNormal"><u><span lang="EN-US" style="font-size:10.0pt;font-family:Consolas;color:black"><o:p><span style="text-decoration:none"> </span></o:p></span></u></p>
<p class="MsoNormal"><span lang="EN-US">When running the last line with gdal.RasterizeLayer(…) my program waits a couple of seconds and then just terminates without any exception or error.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">If I remove vertex number two from the above WKT it works. Vertex 1 and 2 are very close (located in the lower left corner of the polygon).<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">I am running GDAL 1.10.1 and Python 3.2, builds are stable releases from
<a href="http://www.gisinternals.com/sdk/">http://www.gisinternals.com/sdk/</a>.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">What could be the reason for this strange error?<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US">Kind regards, Casper<o:p></o:p></span></p>
</div>
</body>
</html>