Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-12 Thread Mateusz Loskot
Jason Roberts wrote: > Mateusz, > > I'm not an expert in this area, but I think that big performance > gains can be obtained by using a spatial index. Yes, likely true. > For example, consider a situation where you want to clip out a study > region from the full resolution GSHHS shoreline data

Re: [gdal-dev] Memory use in GDALDriver::CreateCopy()

2010-01-12 Thread Greg Coats
I find that GDAL version 1.6.3, released 2009/11/19, gdal_translate fully supports reading and writing a 150 GB GeoTiff image 260,000 columns by 195,000 rows by RGB. Greg $ tiffinfo dcua0002_BigTiffYES.tif Image Width: 26 Image Length: 195000 Tile Width: 512 Tile Length: 512 Resolutio

RE: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-12 Thread Jason Roberts
Mateusz, I'm not an expert in this area, but I think that big performance gains can be obtained by using a spatial index. For example, consider a situation where you want to clip out a study region from the full resolution GSHHS shoreline database, a polygon layer. The shoreline polygons have very

Re: [gdal-dev] Memory use in GDALDriver::CreateCopy()

2010-01-12 Thread Even Rouault
Ozy, Did you try with gdal_translate -of NITF src.tif output.tif -co BLOCKSIZE=128 ? Does it give similar results ? I'm a bit surprised that you even managed to read a 40Kx100K large NITF file organized as scanlines. There was a limit until very recently that prevented to read blocks whose o

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-12 Thread Mateusz Loskot
Jason Roberts wrote: > By integrating with GEOS, OGR can perform various spatial operations on > individual geometries, such as buffer, intersection, union, and so on. Is > there a library that efficiently performs these kinds of operations on > entire OGRLayers? For example, this library would hav

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-12 Thread Stephen Woodbridge
Jason Roberts wrote: Doug, Thanks for these suggestions. It looks like PostGIS and SpatialLite both provide a SQL-based approach for accomplishing what I need. Both look promising and I will dig into them in more detail. It might be less than optimal to load data into one of these, exe

RE: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-12 Thread Jason Roberts
Emilio, Thanks for the suggestion of pysal. It does look interesting, but as you speculated, it seems to not aim to include the traditional spatial operators. Instead it looks like a collection of various interesting algorithms, implemented in Python on top of SciPy, NumPy, spatialindex, and Rtree

RE: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-12 Thread Jason Roberts
Doug, Thanks for these suggestions. It looks like PostGIS and SpatialLite both provide a SQL-based approach for accomplishing what I need. Both look promising and I will dig into them in more detail. It might be less than optimal to load data into one of these, execute the desired spatial q

Re: [gdal-dev] Projection issue with a netcdf file

2010-01-12 Thread Frank Warmerdam
Kyle Shannon wrote: After a little digging, it appears that the netCDF driver does not support the tag: Lambert_Conformal:standard_parallel it supports: Lambert_Conformal:standard_parallel_1 and: Lambert_Conformal:standard_parallel_2 I am not sure if the tag is out of specification, or if the

Re: [gdal-dev] Projection issue with a netcdf file

2010-01-12 Thread Kyle Shannon
After a little digging, it appears that the netCDF driver does not support the tag: Lambert_Conformal:standard_parallel it supports: Lambert_Conformal:standard_parallel_1 and: Lambert_Conformal:standard_parallel_2 I am not sure if the tag is out of specification, or if the driver needs to be upd

Re: [gdal-dev] When GDAL 1.7.0 branch?

2010-01-12 Thread Mateusz Loskot
Frank Warmerdam wrote: > Mateusz Loskot wrote: >> Hi, >> >> When the upcoming 1.7.0 will get its own branch in SVN? > > Mateusz, > > My normal practice is to produce a 1.7 branch at the point the first > RC is prepared. Frank, Great. Thanks! By the way, 1.7.0b2 successfully tested with basic s

Re: [gdal-dev] When GDAL 1.7.0 branch?

2010-01-12 Thread Frank Warmerdam
Mateusz Loskot wrote: Hi, When the upcoming 1.7.0 will get its own branch in SVN? Mateusz, My normal practice is to produce a 1.7 branch at the point the first RC is prepared. Best regards, -- ---+-- I set the clouds in

[gdal-dev] get projection zone and datum using gdal functions

2010-01-12 Thread Shawn GONG
Hi list, I use gdal2tiles.py to generate KML for Google Earth. gdal2tiles.py requires gdal_translate and gdalwarp to convert images to EPSG:4326 If an optical image is in UTM, then "gdalwarp -s_srs '+proj=utm +zone=## ...' Has anybody extracted projection info using Gdal Python function? I can

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-12 Thread Emilio Mayorga
Hi Jason, This may not be quite what you have in mind, but check out the PySAL (Open Source Python Library for Spatial Analytical Functions) project: http://geodacenter.asu.edu/pysal I've never used it, and have only looked at a recent presentation (http://conference.scipy.org/static/wiki/rey_pys

Re: [gdal-dev] Projection issue with a netcdf file

2010-01-12 Thread Christopher Barker
Kyle Shannon wrote: ncdump -h reports the same correct standard parallels (25) whereas Coordinate System reports 0 for both standard parallels. It appears the projection is a Lambert Conformal Conic with 2 standard parallels, but they are the same effectively making it a lcc with 1 standard pa

Re: [gdal-dev] Projection issue with a netcdf file

2010-01-12 Thread Christopher Barker
Kyle Shannon wrote: I am downloading a netcdf file from a thredds data server. It contains 4 grib datasets of weather data. When I run gdalinfo on one of these subdatasets, the standard parallels in the netCDF metadata (marked by astrisks) differs from the standard parallels in the reported C

Re: [gdal-dev] Projection issue with a netcdf file

2010-01-12 Thread Kyle Shannon
ncdump -h reports the same correct standard parallels (25) whereas Coordinate System reports 0 for both standard parallels. It appears the projection is a Lambert Conformal Conic with 2 standard parallels, but they are the same effectively making it a lcc with 1 standard parallel. # =

Re: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-12 Thread Doug_Newcomb
Jason, If you're working with vector data, why not throw the data into Postgresql/Postgis, http://postgis.refractions.net, and use the spatial operators there to select/buffer/intersect the vector geometries as you describe. http://postgis.refractions.net/documentation/manual-1.4 /ch07.html

RE: [gdal-dev] Open source vector geoprocessing libraries?

2010-01-12 Thread Jason Roberts
Frank, Thanks for your thoughts on this. > I'd like to see something along this line happen. I to do it efficiently > it would be necessary to dig into GEOS past the C interface so that > a spatial index on a collection of features can be maintained over time > rather than created and discarded

[gdal-dev] Projection issue with a netcdf file

2010-01-12 Thread Kyle Shannon
I am downloading a netcdf file from a thredds data server. It contains 4 grib datasets of weather data. When I run gdalinfo on one of these subdatasets, the standard parallels in the netCDF metadata (marked by astrisks) differs from the standard parallels in the reported Coordinate system: k...@

[gdal-dev] Memory use in GDALDriver::CreateCopy()

2010-01-12 Thread ozy sjahputera
I was trying to make a copy of a very large NITF image (about 40Kx100K pixels) using GDALDriver::CreateCopy(). The new file was set to have different block-size (input was a scanline image, output is to have a 128x128 blocksize). The program keeps getting killed by the system (Linux). I monitor the

[gdal-dev] gdal_translate JPEG to GeoTiff

2010-01-12 Thread Marco Afonso
Hi all, I'm new to gdal tools. I'm trying to transform a JPEG world image to GeoTiff. My jpeg world map image is cropped to fit on -180.0, -59.1, 180, 90 (minx, miny, maxx, maxy - EPSG:4326). Here are the steps: % gdal_translate -of GTiff -a_srs EPSG:4326 -co "TILED=YES" -co "JPEG_QUALITY=100"

[gdal-dev] When GDAL 1.7.0 branch?

2010-01-12 Thread Mateusz Loskot
Hi, When the upcoming 1.7.0 will get its own branch in SVN? Best regards, -- Mateusz Loskot, http://mateusz.loskot.net ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev