Re: [gdal-dev] divide a big shape file using ogr2ogr

2012-02-21 Thread ahmet temiz
Thank you Chaitanya It worked. I had tried spat. But it couldn't achive. Could you tell me why spat can not be used to part shp file because of its big size ? regards 2012/2/18 Chaitanya kumar CH chaitanya...@gmail.com: Ahmet, If you want to divide the shapefile to reduce file size or the

[gdal-dev] Hillshade + Topographic Map??

2012-02-21 Thread TJMartin
Hi I have managed to create a hillshade using gdaldem and I have my original topographic rasters. However I do not want the standard grey hillshad or apply a colour ramp. Is there a way to combine the hillshade and topo map together - ie add a transparency or opacity to the hillshade. I have

Re: [gdal-dev] Hillshade + Topographic Map??

2012-02-21 Thread Frank Broniewski
Hi, the standard procedure for creating topo-maps is to overlay your land coverage map with a semi-transparent hillshading. I don't know if there's a tool for combining/adding two raster maps together in gdal, but this is a rather standard procedure in any GIS. In QGIS, for example, you can

Re: [gdal-dev] Hillshade + Topographic Map??

2012-02-21 Thread Volker Wichmann
The workflow with QGIS described by Frank is also working with SAGA GIS. SAGA also allows you to export a permament combination as image file (tif, jpg, png ... including world file and kml) with the 'Export Image' module. Volker On 02/21/2012 02:14 PM, Frank Broniewski wrote: Hi, the

[gdal-dev] Re: Hillshade + Topographic Map??

2012-02-21 Thread Rutger
Hey, TJMartin wrote ... Is there a way to combine the hillshade and topo map together - ie add a transparency or opacity to the hillshade. ... A common way is to convert your (i assume) RGB topo map to HSV color space and replace the intensity band (v) with your hillshade. This is for

Re: [gdal-dev] Re: Hillshade + Topographic Map??

2012-02-21 Thread Joaquim Luis
On 21-02-2012 14:11, Rutger wrote: Hey, TJMartin wrote ... Is there a way to combine the hillshade and topo map together - ie add a transparency or opacity to the hillshade. ... A common way is to convert your (i assume) RGB topo map to HSV color space and replace the intensity band (v)

Re: [gdal-dev] divide a big shape file using ogr2ogr

2012-02-21 Thread Frank Warmerdam
On 12-02-21 12:29 AM, ahmet temiz wrote: Thank you Chaitanya It worked. I had tried spat. But it couldn't achive. Could you tell me why spat can not be used to part shp file because of its big size ? Ahmet, Roughly speaking -spat could be used for this purpose, but any features crossing the

Re: [gdal-dev] nmake.opt compile rule changes

2012-02-21 Thread Jeff McKenna
Thanks for the improvements Frank - no problems here with trunk. -jeff -- Jeff McKenna MapServer Consulting and Training Services http://www.gatewaygeomatics.com/ On 12-02-20 8:39 PM, Frank Warmerdam wrote: Folks, If you don't build on Windows using visual studio NMAKE files you can

Re: [gdal-dev] Re: Hillshade + Topographic Map??

2012-02-21 Thread tim martin
Hi Everyone Many thanks for your help. Lots of recommendations and I have tried each one with limited success 1) QGIS - works but i need a permanent copy. 2) GIMP - problem with this is GIMP does not recognise my rasters. I have expanded to color table to rgb and it still wouldnt open the file.

[gdal-dev] Postgis Raster driver improvements.

2012-02-21 Thread David Zwarg
Hello, I would like to help add write support and more to the Postgis Raster driver. Currently, there is a loader script in the autotest/data directory for the driver -- when write support is complete, will that script no longer be necessary? I don't see any corresponding test loading scripts

Re: [gdal-dev] Postgis Raster driver improvements.

2012-02-21 Thread Frank Warmerdam
On Tue, Feb 21, 2012 at 10:38 AM, David Zwarg dzwarg+g...@azavea.com wrote: Hello, I would like to help add write support and more to the Postgis Raster driver.  Currently, there is a loader script in the autotest/data directory for the driver -- when write support is complete, will that

Re: [gdal-dev] Re: Hillshade + Topographic Map??

2012-02-21 Thread Joaquim Luis
5) Looked at GMT and I cannot find the options to load two rasters on top of each other. Tim, You don't load the two grids on top of one-another. Illumination is NOT like transparency. It uses the gradients of the grid, projects it into the direction of light and than change the color in

[gdal-dev] learning GDAL with C/C++

2012-02-21 Thread joseph.nalluri
Hello All, I am very pleased to find this forum. This is my first post. Please bear with me. In my project, I am trying to read in .bil file (hyperspectral image) into my C/C++ code. ( I haven't done this yet, because I am still trying to see, how can I do it). After I read in, I run some

Re: [gdal-dev] Re: Hillshade + Topographic Map??

2012-02-21 Thread Donovan Cameron
Have you reviewed the article from Tim Sutton: A workflow for creating beautiful relief shaded dems using GDAL[1]? It has some valuable pointers. [1] http://linfiniti.com/2010/12/a-workflow-for-creating-beautiful-relief-shaded-dems-using-gdal/ Donovan On Tue, Feb 21, 2012 at 12:29 PM,

Re: [gdal-dev] learning GDAL with C/C++

2012-02-21 Thread Frank Warmerdam
Joseph, I'd start with http://www.gdal.org/gdal_tutorial.html, the general purpose API tutorial. It isn't comprehensive, but I think it is fairly good for the areas covered. Note that GDAL does not directly help with distributed processing. I assume this is what you mean by running on a

[gdal-dev] gdaladdo, overviews and NODATA

2012-02-21 Thread Armin Burger
Hi all I have a question regarding the gdaladdo tool and NODATA pixels: When using resampling like average it seems that along the border of DATA and NODATA pixels there happens an averaging also of some of the original NODATA pixels (looks like 1-2 pixel width). This way in the overviews

Re: [gdal-dev] gdaladdo, overviews and NODATA

2012-02-21 Thread Frank Warmerdam
Armin, Reviewing the overview building code in gdal/gcore/overview.cpp it clearly tries to exclude nodata values from the averaging calculation. What does the gdalinfo report look like for the file before you run gdaladdo? Best regards, On Tue, Feb 21, 2012 at 12:36 PM, Armin Burger

Re: [gdal-dev] gdal post - boundary polygone for rectified images

2012-02-21 Thread Eli Adam
Gianluca, Or if just there is a way to create a vector file showing the boundary of the rectified image. What I am searching is a script in gdal style in order to automate the procedure for several rasters without further user input. OGR has tindex which makes a shapefile outline of

Re: [gdal-dev] gdaladdo, overviews and NODATA

2012-02-21 Thread Armin Burger
Frank the original images showed ... Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=PIXEL ... Band 1 Block=3000x1 Type=Byte, ColorInterp=Red Band 2 Block=3000x1 Type=Byte, ColorInterp=Green Band 3 Block=3000x1 Type=Byte, ColorInterp=Blue Then I tried to set the NODATA

Re: [gdal-dev] gdaladdo, overviews and NODATA

2012-02-21 Thread Frank Warmerdam
On Tue, Feb 21, 2012 at 1:17 PM, Armin Burger armin.bur...@gmx.net wrote: Then I tried to set the NODATA value to 0 using Gdal-Python, and afterwards they had Band 1 Block=3000x1 Type=Byte, ColorInterp=Red  NoData Value=0 Band 2 Block=3000x1 Type=Byte, ColorInterp=Green  NoData Value=0

Re: [gdal-dev] Postgis Raster driver improvements.

2012-02-21 Thread David Zwarg
Hi Frank, Yes, I am modernizing some of the testing scripts for PostGISRaster, and adding write support in the near future. I've been in touch with and am coordinating with Jorge. I think it was just the test script that used the old name, btw. The driver is truly called PostGISRaster, just the

[gdal-dev] gdalwarp extent issues

2012-02-21 Thread Jay L.
List, I am attempting to reproject 8 gtiffs, each in a local equirectangular projection, back to gcs for use in geoserver. I immediately went to gdalwarp to perform this task and am having an issue. The input files are stored in positive longitude, while the output SRS is in negative longitude.

[gdal-dev] comparing two rasters

2012-02-21 Thread jdmorgan
Hello GDAL guru's, I am working on a python script where I read in two rasters of similar extent and resolution.Then I re-assign any values that are greater that zero to a 1.Next, I compare to the rasters and attempt to create a third resulting raster which has 1's everywhere that the two

[gdal-dev] Using SWIG from a Driver

2012-02-21 Thread Michael Speth
Greetings GDAL devs, I am interested in developing a GDAL driver that calls out to a Java library that accesses a database that we are developing using SWIG's Java Bindings. Are there any examples of existing GDAL drivers using the SWIG bindings? If not, how might I use the SWIG library

Re: [gdal-dev] comparing two rasters

2012-02-21 Thread Chaitanya kumar CH
Derek, Can you explain the following lines towards the bottom of the script? data1[data10]=1 ... data2[data20]=1 On Wed, Feb 22, 2012 at 8:46 AM, jdmorgan jdmor...@unca.edu wrote: Hello GDAL guru’s, I am working on a python script where I read in two rasters of similar extent and

[gdal-dev] compiling gdal on Obuntu via Parallels

2012-02-21 Thread Stuart Edwards
Hi This seems to be an old problem (see http://www.osgeo.org/pipermail/gdal-dev/2009-October/022308.html) but when I try to compile GDAL 1.9.0 from source using Obuntu 10.04 through Parallels Desktop 6.0.12106 on OS X 10.6.8 (don't ask why this convoluted approach - long ESRI related story)