Tai,
You can use the spatial filter of the ogr layer to filter the point
features to speed things up. You will still have to test if the point is
within the polygon but it will be faster.
Also, make sure the test is like point->Within(polygon) and not
polygon->Within(point)
OGRLayer::SetSpatialF
Tai,
You can use OGR_G_Within()
OGR_G_WIthin(point,polygon) will return true if the point is within the
polygon.
You can access various other geometry functions from the OGRGeometry class.
http://www.gdal.org/ogr/ogr__api_8h.html#a80ba9e516ccae381ae7240cc9684eb72
On Fri, Jul 20, 2012 at 7:04 AM
Hi everyone,
Do you know how to identify whether a point (OGRPoint) is belong to a
polygon (OGRPolygon) ? (the polygon covers the point)
Thanks and regards,
Tai
--
View this message in context:
http://osgeo-org.1560.n6.nabble.com/How-to-identify-a-point-is-belong-to-a-polygon-tp4989501.html
Se
Hello Everyone
Thanks a lot Jose.
Your Script running okYour solution is very clever.
Congratulations.
--
View this message in context:
http://osgeo-org.1560.n6.nabble.com/Convert-to-temperature-in-degrees-C-in-MODIS-Ocean-Color-tp4989237p4989478.html
Sent from the GDAL
On 19 July 2012 15:25, Jose Gomez-Dans wrote:
> Hi,
>
>
> On 19 July 2012 04:46, sigologo wrote:
>
>> I'm trying import A20030012003008.L3m_8D_NSST_4, with gdal with this
>> script:
>>
>> for file in A*L3m_8D_NSST_4 ; do
>> echo "map: $file"
>>
>> gdal_translate -a_srs "+init=epsg:4326" -a_n
Hi,
Unfortunately the GMT5 CookBook is no longer in html so I have to refer
you the section 4.21 - "Options for grids and images read via GDAL" - of
http://www.soest.hawaii.edu/gmt5/gmt/pdf/GMT_Docs.pdf
where you can see this example
"... Be warned, however, that things are not yet completed
I mean it requires the gdal python bindings and python installed.
You don't need python skills - it's a command-line program. Although
you may need to understand the numpy syntax if you have a complex
formula.
In your case it would be something like
gdal_calc.py -A input1.tif --outfile=result.ti
Hi List,
We've started the MapWindow Summer of Code 2012 again.
It's a bit late, but we don't have an end date ;)
This year we're focusing on the C++ ActiveX control MapWinGIS.
And because a lot of tasks are related to GDAL, I thought I mail this list
as well.
If you are an experience C++ develo
Hi,
On 19 July 2012 04:46, sigologo wrote:
> I'm trying import A20030012003008.L3m_8D_NSST_4, with gdal with this
> script:
>
> for file in A*L3m_8D_NSST_4 ; do
> echo "map: $file"
>
> gdal_translate -a_srs "+init=epsg:4326" -a_nodata 65535 \
> -a_ullr -180 90 180 -90 -co "COMPRESS=PACK
Thanks Etienne
But python is not my best skill in this session
I will try!
--
View this message in context:
http://osgeo-org.1560.n6.nabble.com/Convert-to-temperature-in-degrees-C-in-MODIS-Ocean-Color-tp4989237p4989377.html
Sent from the GDAL - Dev mailing list archive at Nabble.co
You can try gdal_calc.py, although that requires python support
svn.osgeo.org/gdal/trunk/gdal/swig/python/scripts/gdal_calc.py
On Thu, Jul 19, 2012 at 12:46 AM, sigologo wrote:
> Hello Everyone
>
>
> I'm trying import A20030012003008.L3m_8D_NSST_4, with gdal with this script:
>
> for file in A*L
Even Rouault wrote:
>> Windows 7, 64-bit, SATA disk and 2x3 GHz is converting Finland.osm.pbf in
>> about
>> 8 minutes for me. But execution time does not increase at all linearly.
>> Germany.osm.pbf is about 10 times larger in filesize but ogr2ogr had to work
>> about 8 hours with it, thus rough
Even Rouault wrote:
>> Amenity is included in my osmconf.ini and it can be used inside -sql.
>> However, your example gives me this
>>
>> ERROR 1: 'amenity' not recognised as an available field.
>> FAILURE: SetAttributeFilter(amenity='toilets') failed.
> Hum, I think I know what's wrong. I see th
>
> Windows 7, 64-bit, SATA disk and 2x3 GHz is converting Finland.osm.pbf in
> about
> 8 minutes for me. But execution time does not increase at all linearly.
> Germany.osm.pbf is about 10 times larger in filesize but ogr2ogr had to work
> about 8 hours with it, thus roughly 70 times longer. Obvio
Even,
On Tue, Jul 17, 2012 at 10:20 AM, Even Rouault wrote:
>
>
> You could just remove the tweaking done with GDAL_CACHEMAX and -wm and
> leave
> them to their default values. In my experience, a so big value for
> GDAL_CACHEMAX
> is of no use (no need to use more than 500 MB), and when merging
> Amenity is included in my osmconf.ini and it can be used inside -sql.
> However, your example gives me this
>
> ERROR 1: 'amenity' not recognised as an available field.
> FAILURE: SetAttributeFilter(amenity='toilets') failed.
Hum, I think I know what's wrong. I see that, even if a subset of lay
Even Rouault mines-paris.org> writes:
>
> Selon Jukka Rahkonen mmmtike.fi>:
...
>
> > It is
> > much faster to convert all the points than select only a part of those. The
> > error message suggest that ogr2ogr is inspecting the input file more closely
> > than necessary for this use case:
17 matches
Mail list logo
