Sorry, I forget to say that I built Gdal with GEOS<br><br><div class="gmail_quote">2012/9/4 Gabriel Fusca <span dir="ltr"><<a href="mailto:gabrielfusca@suremptec.com.ar" target="_blank">gabrielfusca@suremptec.com.ar</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi,<br><br>I'm using Gdal 1.9.<br><br>I have a function that has to perform a intersection of a geometry based on the click that is made on the screen. <br>
<br>The raster file (Geotiff) that has a projected reference system (WGS 84 / UTM zone 21S) and the geometry a georeferenced system (EPSG - 4326 WGS84). <br>
<br>The click that comes from the screen (spatial reference of the image) create a geometry from calling point.Buffer (radius).<br><br>When i want to verify that the geometries intersect the intersects method from OGRGeometry fails, even if the querypoint is transformed into the spatial reference system of the geometry.<br>
<br>Envelopes from the 2 geometries. <br>QueryPoint <br>env. <br>if i generates the querypoint from the spatial reference of the geometry that i want to intersect)<br>minx -61.66, maxx -53.66;<br>miny -41.97, maxy -33.97<br>
else<br>minx 443322.50, maxx 443330.50<br>miny 5798153.59, maxy 5798161.59<br><br>
Geometry env.<br>minx -57.79, maxx -57.59;<br>miny -38.08, maxy-37.92<br><br><br>The code is something like this:<br><br>// spatialfilterwindow image subset transformed into <br>// the spatial reference of the vector to intersect<br>
Extent spatialfilter(spatialfilterwindow);<br>player->SetSpatialFilterRect(spatialfilter.min_.x_,<br> spatialfilter.min_.y_,<br> spatialfilter.max_.x_,<br> spatialfilter.max_.y_);<br>
OGRLayer* player = ...;<br>OGRFeature *pfeature = player->GetNextFeature(); <br>OGRPoint querypoint = OGRPoint(Point.x_, Point.y_);<br>querypoint.assignSpatialReference(new OGRSpatialReference(SrWkt.c_str()));<br>OGRGeometry* pquerypoly = querypoint.Buffer(radius);<br>
<b>while </b>(pfeature) {<br> OGRGeometry *pgeometry = pfeature->GetGeometryRef();<br> <b>if </b>((pgeometry && pgeometry->Intersects(pquerypoly) == TRUE)<br> || (pquerypoly && pquerypoly->Intersects(pgeometry) == TRUE)) {<br>
// do something...<br> } <b>else </b>{<br> // handle the intersection fail<br> }<br> pfeature = player->GetNextFeature();<br>}<br><br>Any idea why it fails?<br><br clear="all">Thanks,<br>Regards,<span class="HOEnZb"><font color="#888888"><br>
<br>-- <br><span><font color="#888888">Gabriel Fusca<br>
SUR Emprendimientos Tecnológicos<br>
<br>Perú 345 Piso 5to Oficina "B" (C1067AAG)<br>Ciudad de Buenos Aires, Argentina<br>Tel. +54 (11) 4342-2976/84<br><a href="mailto:gabrielfusca@suremptec.com.ar" target="_blank">gabrielfusca@suremptec.com.ar</a><br>
<a href="http://www.suremptec.com/" target="_blank">http://www.suremptec.com/</a><br></font></span><br>
</font></span></blockquote></div><br><br clear="all"><br>-- <br><span><font color="#888888">Gabriel Fusca<br>SUR Emprendimientos Tecnológicos<br><br>Perú 345 Piso 5to Oficina "B" (C1067AAG)<br>Ciudad de Buenos Aires, Argentina<br>
Tel. +54 (11) 4342-2976/84<br><a href="mailto:gabrielfusca@suremptec.com.ar" target="_blank">gabrielfusca@suremptec.com.ar</a><br><a href="http://www.suremptec.com/" target="_blank">http://www.suremptec.com/</a><br></font></span><br>