RE: [gdal-dev] gdal_warp

2010-06-21 Thread John Donovan
I've got a slightly older version of the data you're using (OS 10k B&W, but with 2 bits per pixel rather than four), and I've used a two-stage process to achieve what you want. First, I used gdal_translate to expand the image to 8-bit greyscale: gdal_translate -co "TFW=YES" -co "TILED=YES" -expan

RE: [gdal-dev] How to represent multi-dimensional array

2010-06-21 Thread Jason Roberts
> > The trick is your netCDF has to meet a bunch > > of constraints for ArcGIS to recognize it. It has to have square cells. > > bingo -- that was one of our key problems -- wait -- they have to be > "square" -- rectangular won't do? arrgg! Oops, I'm pretty sure you're right, they can be rectangu

Re: [gdal-dev] How to represent multi-dimensional array

2010-06-21 Thread Christopher Barker
Jason Roberts wrote: The trick is your netCDF has to meet a bunch of constraints for ArcGIS to recognize it. It has to have square cells. bingo -- that was one of our key problems -- wait -- they have to be "square" -- rectangular won't do? arrgg! Oops, I'm pretty sure you're right, they can b

[gdal-dev] Use of "LIKE" in OGR layer attribute filters

2010-06-21 Thread Chris Bailey
We are using OGRLayer SetAttributeFilter to implement some search capabilities. One feature that we want to allow is using the SQL "Like" operand in the where clause to allow wildcard searches or searches on partial fields. The documentation says that "Like" is not supported but it appears to

Re: [gdal-dev] Use of "LIKE" in OGR layer attribute filters

2010-06-21 Thread Frank Warmerdam
On Mon, Jun 21, 2010 at 4:52 PM, Chris Bailey wrote: > > We are using OGRLayer SetAttributeFilter to implement some search > capabilities.  One feature that we want to allow is using the SQL "Like" > operand in the where clause to allow wildcard searches or searches on > partial fields. > > The do

[gdal-dev] Re: Use of "LIKE" in OGR layer attribute filters

2010-06-21 Thread Chris Bailey
This is the documentation I was referring to: http://www.gdal.org/ogr/classOGRLayer.html#cb2c6cc5fa3577df5be538284c1b0dde Set a new attribute query. This method sets the attribute query string to be used when fetching features via the GetNextFeature() method. Only features for which the query e

[gdal-dev] Raster to array, then array to Raster (while preserving projection info)

2010-06-21 Thread Discourse Maps
I am using "convertRastertoArray" (modifying the array) and then using "convertArraytoRaster". The problem is that when I convert back to raster, the projection information is lost. The projection of my input is: WGS 84, UTM 17N. (EPSG = 32617) Can someone provide some code as to the cl

Re: [gdal-dev] How to represent multi-dimensional array

2010-06-21 Thread Michael Sumner
I find this categorization helpful: http://en.wikipedia.org/wiki/Regular_grid Any GIS will do regular and cartesian grids (but usually bound to 2D), and with programming constructs you can handle rectlinear or curvilinear grids of any dimension, but it's not well supported in high-level software

Re: [gdal-dev] Raster to array, then array to Raster (while preserving projection info)

2010-06-21 Thread Chaitanya kumar CH
GDALCreateCopy() can help you with creating a raster with the same data including the projection and stuff. You can edit the pixels without bothering about the projection. Use GDALRasterIO() to read and write to that raster. http://www.gdal.org/gdal_tutorial.html http://www.gdal.org/gdal_8h.html