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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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
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.
# =
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
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
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...@
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
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"
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
23 matches
Mail list logo