[gdal-dev] GDAL 1.11.3 ISIS3 Driver

Just built GDAL 1.11.3 using the following: ./configure --prefix=/home/travis/miniconda/envs/_build --with-geos=/home/travis/miniconda/envs/_build/bin/geos-config --with-static-proj4=/home/travis/miniconda/envs/_build --with-hdf5=/home/travis/miniconda/envs/_build --with-xerces=/home/travis/minico

Re: [gdal-dev] Proj4 Transformation Error 6

Super cool. Thanks! On Fri, Jun 12, 2015 at 10:29 AM, Kyle Shannon wrote: > Jay, > > > On Fri, Jun 12, 2015 at 11:21 AM, Kyle Shannon wrote: > > Jay, > > > > > > On Fri, Jun 12, 2015 at 10:59 AM, Jay L. wrote: > >> Using GDAL 1.11.2 (Anaconda Pyth

[gdal-dev] Proj4 Transformation Error 6

Using GDAL 1.11.2 (Anaconda Python osgeo binstar install). I have a WKT projection: 'PROJCS["Mercator",GEOGCS["GCS_Moon_2000",DATUM["D_Moon_2000",SPHEROID["Moon_2000_IAU_IAG",1737400.0,0.0]],PRIMEM["Reference_Meridian",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Mercator"],PARAMETER["False

Re: [gdal-dev] Standard Dev. Environment?

o update the README file in the repo since the > procedure will be easier (since a few things I was using have been directly > integrated as vagrant plugin). > > I plan to finish this during the next week. > > Alan > > > On 14-01-10 01:36 PM, Daniel Morissette wrote

Re: [gdal-dev] Standard Dev. Environment?

ress' > > > /home/vagrant/gdal-1.10.1/libgdal.so: undefined reference to > > > `jpeg_read_scanlines' > > > > > > So lib-jpeg is seen. This doesn't make sense since I have > --without-jpeg > > > specified in the config. I looked in the c

Re: [gdal-dev] Standard Dev. Environment?

ave missed it - does a flag to specify the full PATH to libjpeg exist? Thanks, Jay On Fri, Jan 10, 2014 at 7:24 AM, Even Rouault wrote: > Selon "Jay L." : > > > Frank, > > > > Thanks for the heads up. I am using the --without-jasper flag as part of > &g

Re: [gdal-dev] Standard Dev. Environment?

per? Normally a default build ("./configure") should > work pretty well on a precise system. I'm guessing you have enabled jasper > because it is key to what you want to do? If not --without-jasper should > be straight forward. > > Best regards, > Frank > &g

[gdal-dev] Standard Dev. Environment?

List, I wonder if a standard dev. environment exists? I am working to extend a driver and am having trouble getting an environment working that will compile. Current setup is a VagrantVM, Ubuntu 12.04 32-bit with source downloaded from SVN. Attempting to build with minimized drivers

[gdal-dev] Python Bindings ReadAsArray and ReadRaster Errors

List, I am attempting to open a multi-band image using the PDS driver. I have a working, GDAL readable label that I can run gdalinfo and gdal_translate on without error. I am attempting to access a number of bands and operate on them using NumPy. band.ReadAsArray() is failing silently in proces

Re: [gdal-dev] 32 and 64-bit python bindings concurrently?

Evan, Thanks for the info. I used universal, purely as a descriptor of the functionality I was hoping for. I will have to experiment with installing GDAL in separate locations to see if I can get both installations to get along. Thanks! Jay On Mon, Aug 20, 2012 at 12:24 PM, Even Rouault wrote

[gdal-dev] 32 and 64-bit python bindings concurrently?

List, I have a number of python scripts leveraging gdal, numpy, and scipy. It is necessary, if possible, to support machines which utilize both 32 and 64 bit python installations. The newest ArcGIS10.1 background geoprocessing is a perfect example of this, where processes run within arc utilize

[gdal-dev] Potential thrashing on Band #

List, I am receiving the "Potential thrashing on bans x of message doing some image manipulation using the python bindings. I am manually setting the cache to 1GB and calling band.FlushCache(). I am also reading the image at either the block size or multiples of the block size. What other tech

Re: [gdal-dev] GDAL - Python - ArcGIS10.1

ng of the system > PATH. > Or alternatively, you can run the GDAL console (to set the environment > correctly) and then start python.exe from this command prompt. > > Best regards, > > Tamas > > > 2012/8/6 Jay L. > >> All, >> >> I am trying

[gdal-dev] GDAL - Python - ArcGIS10.1

All, I am trying to get the gdal python bindings working with an ArcGIS 10.1 installation. I grabbed the core from gisinternals.com/sdk (thanks Tamas!) and installed both the 1500 core files (MVS 2008) and the python 2.7 bindings. I appended the path as necessary and added the GDAL_DATA director

Re: [gdal-dev] Capture stderr in interactive Python shell

try this: from osgeo import gdal gdal.UseExceptions() gdal.Open('noexist.tif') This is from: http://trac.osgeo.org/gdal/wiki/PythonGotchas On Sun, Aug 5, 2012 at 5:41 PM, Mike Toews wrote: > I would like to capture or show GDAL error messages in an interactive > Python shell. For instance, if

Re: [gdal-dev] Python Bindings - ReadArray direct to ctypes array?

. Is that test failing on OS X then? [1] http://svn.osgeo.org/gdal/sandbox/rouault/gdal-openjpeg-trunk/autotest/gcore/testnonboundtoswig.py Thanks, Jay On Mon, Jul 30, 2012 at 8:19 AM, Jay L. wrote: > Chris, > > Thanks for the link / info. My issue is not working with the ctypes &g

Re: [gdal-dev] Python Bindings - ReadArray direct to ctypes array?

inaccuracies as I am trying to improve my fluency in GDAL! Jay On Sun, Jul 29, 2012 at 2:20 PM, Chris Barker wrote: > On Fri, Jul 27, 2012 at 8:40 PM, Jay L. wrote: > > > Currently the workflow is to open the dataset and grab the first band, as > > usual. I then read in the arra

Re: [gdal-dev] Python Bindings - ReadArray direct to ctypes array?

Even, Thanks so much! I will give this a try and get it implemented into my code. This is also a wonderful primer on ctypes integration with existing Python scripts. Thanks, Jay On Sat, Jul 28, 2012 at 2:17 AM, Even Rouault wrote: > Le samedi 28 juillet 2012 05:40:57, Jay L. a éc

[gdal-dev] Python Bindings - ReadArray direct to ctypes array?

List, I am working on a python script using the multiprocessing module with shared memory ctype arrays. Currently the workflow is to open the dataset and grab the first band, as usual. I then read in the array, create an empty ctype array, and use memmove to move the numpy array to my ctypes arra

[gdal-dev] Upgraded to 1.9 - Python Bindings Error

List, Just upgraded to 1.9.0 on my windows 7 machine and I can no longer access gdal via the python bindings. Command line (ogr2ogr, gdalinfo, etc. are functioning as expected). Specifically, running 'from osgeo import gdal', returns: File "C:\Python26\ArcGIS10.0\lib\site-packages\osgeo\__init.

Re: [gdal-dev] GDAL Tiling

Akshay, 1) An AOI can be created using the -projwin option in gdal_translate. You can use that option when creating a new raster with a new tiles size. 2) Using the Python bindings it is possible to read a segment of a raster into memory without allocating memory. I am unsure using another othe

Re: [gdal-dev] GDAL Tiling

For a geotiff I would use the following to change block size. gdal_translate -of GTiff -co "TILED=YES" -co "BLOCKXSIZE=256" -co "BLOCKYSIZE=256" Replace 256with whatever integer you want. I realize that this is outside of your program, but it looks like you are going to need to modify the inp

Re: [gdal-dev] gdalwarp extent issues

uis wrote: > > > On Wed, Feb 22, 2012 at 12:54 AM, Jay L. wrote: > >> 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 a

[gdal-dev] gdalwarp extent issues

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.

Re: [gdal-dev] Re: 254 into 255

Thanks for posting this! I had to no idea that this functionality existed. Is this documented somewhere? On Thu, Oct 27, 2011 at 4:58 AM, Chaitanya kumar CH wrote: > J, > > A vrt file[1] can implement your requirements. It can create a 'lookup > table' to translate your pixel values. > Start by

Re: [gdal-dev] Re: Gdal_translate, question from an amateur.

I need to convert all the 255 255 255 into 254 254 254. Is this possible? > PS. Do appologize for lack of my language. > > Thanks > J > > On 23 Oct 2011, at 15:26, "Jay L. [via OSGeo.org]" <[hidden > email]<http://user/SendEmail.jtp?type=node&node=6922377&

Re: [gdal-dev] Re: Gdal_translate, question from an amateur.

Or are you trying to convert only some of the pixels to white? On Sun, Oct 23, 2011 at 10:02 AM, Brian Case wrote: > joolek, > > I noticed you used the word "strip". Are you trying to remove portions > of the image that are all black or white? > > > On Sun, 2011-10-23 at 03:01 -0700, joolek wrot

Re: [gdal-dev] Gdal_translate, question from an amateur.

I'm not sure I fully understand the question. Are you attempting to convert all pixels with RGB value 72, 255, 25 to white? Or are you trying to scale the image to 8-bit? Sounds like the former. In this case -scale is not the way to go. I'm honestly not sure what is, short of reading each band

[gdal-dev] GTiff INTERLEAVE=BAND

I have a series of raster images (GTiff) that are causing ArcMap10 to hang. After trying a number of different things (32bit to 8bit, translate to another format, compression, tiling, etc.) I suspect that the issue could be with the interleave. Gdalinfo shows that INTERLEAVE=BAND. The GTiff driv

Re: [gdal-dev] TRI, TPI, and Ruggedness

Evan, Thanks for the info. Jay On Sat, Oct 8, 2011 at 12:29 PM, Even Rouault wrote: > Le samedi 08 octobre 2011 17:26:58, Jay L. a écrit : > > What is the box size on these indices (TRI, TPI, Roughness)? > > 3x3 window > > > Is it user modifiable? > > No > &g

[gdal-dev] TRI, TPI, and Ruggedness

What is the box size on these indices (TRI, TPI, Roughness)? Is it user modifiable? Finally, in your experience, which index have you selected for a specific application and why? (Cross-posted this question to gis.stackexchange.com) Thanks, Jay ___ gd

Re: [gdal-dev] Fwd: gdal_merge.py - Missing images - sans attachment

e vrt file into a single image. The result > is usually the same. > > > On Fri, Oct 7, 2011 at 9:12 PM, Jay L. wrote: > >> Sent with an example attachment, which is awaiting approval. Here is the >> mail without attachement. >> >> >> I have 8 large

Re: [gdal-dev] Re: gdal_merge.py - Missing images

ge the original IMGs- > > gdalbuildvrt merge.vrt *.img > > -then- > > gdal_translate merge.vrt output.tif > > > > Message: 4 > > Date: Fri, 7 Oct 2011 11:42:52 -0400 > > From: "Jay L." > > Subject: [gdal-dev] Fwd: gdal_merge.py - Mis

[gdal-dev] Fwd: gdal_merge.py - Missing images - sans attachment

Sent with an example attachment, which is awaiting approval. Here is the mail without attachement. I have 8 large Gtiffs which I am trying to mosaic. Tiffs were created from .IMG files downloaded from PDS. I converted to tif using: gdal_translate -of GTiff -scale -co "TILED=YES" -co "

Re: [gdal-dev] Problem with type of data in ESRI SHAPEFILE

What about something like this? The block in the middle copies your fields and field definitions. The def statement handles copying the individual fields for each feature. Your script then goes in an modifies the geometry. This does roughly the same thing that Frank suggested, but is more scri

Re: [gdal-dev] gdal GeoTiff parsing gml

How about turning the GTiff into a point and adding z either into the geometry or into the attribute table? That functionality would be accessible through any of the APIs. Another idea would be to use polygonize to create polygons of your raster. Here is a script

Re: [gdal-dev] tiled vector rendering

I believe that geoserver's geowebcache will perform similar functionality. Here is a URL to geowebcahce documentation: http://docs.geoserver.org/stable/en/user/geowebcache/index.html Best of luck! On Wed, Aug 10, 2011 at 8:33 AM, Jean-François Gigand wrote: > Hi Rashad, > > There is MapServer

Re: [gdal-dev] Raster size and ReadAsArray()

time)? J On Wed, Aug 3, 2011 at 11:19 AM, Even Rouault wrote: > Le mercredi 03 août 2011 17:32:53, Antonio Valentino a écrit : > > Hi Jay, > > > > Il 03/08/2011 16:53, Jay L. ha scritto: > > > I have been working on this problem as well. Initially, the attempt &g

Re: [gdal-dev] Raster size and ReadAsArray()

I have been working on this problem as well. Initially, the attempt was to ReadAsArray small chunks. Unfortunately this is quite inefficient. Someone more knowledgeable will know why, but I suspect it has to do with either thrashing or the fact that full blocks are not being read in (as is the c

Re: [gdal-dev] JP2k, ReadError, Python

! Jay On Sat, Jul 30, 2011 at 10:44 AM, Even Rouault wrote: > Le samedi 30 juillet 2011 19:14:17, Jay L. a écrit : > > Here is the link. As I said, it is 28MB, including the .jp2, .aux.xml, > and > > .prj file. > > > > http://www.personal.psu.edu/jzl5325/CTX.zip >

Re: [gdal-dev] JP2k, ReadError, Python

Here is the link. As I said, it is 28MB, including the .jp2, .aux.xml, and .prj file. http://www.personal.psu.edu/jzl5325/CTX.zip I would have thought that unloading the JPEG2000 driver would have forced the ECW driver to be used. It looks to me as if that driver is not being used although it i

Re: [gdal-dev] JP2k, ReadError, Python

x(False) ? I had found that band.GetMinimum() and band.GetMaximum(), which processed quite quickly, often returned "None" is statistics had not already been calculated for the raster. Using band.ComputeRasterMinMax(), while slower, always allowed me to get a value. Is there a bet

[gdal-dev] JP2k, ReadError, Python

Morning all, I have a rather large JP2 file (created using the Kakadu driver I believe) that I am trying to process with a python script. I am able to open the dataset and grab the first band. To calculate band statistics I use the following code which returns 1000 lines of error before exiting

Re: [gdal-dev] preserving alpha channel data with gdal_transform.py

What format do you want the alphas channel in after running gdal_polygonize.py? I ask because after processing you could open the input and read just the alpha band (band 4 I assume). Here is a short script framework that might get you started. Something like: from osgeo import gdal gdal.open(i

[gdal-dev] Python WriteArray()

I am working on processing some raster images and have been really struggling with processing times. I have tracked the issue down to GDAL's WriteArray() function. I am cutting the image into blocks and processing in small chunks. This is facilitate working with really large images where the ent

[gdal-dev] OGR equivalent to ARCpy MakeFeatureLayer

I have a script, written to use ESRI's ArcPy, which I am attempting to re-write using the ogr/gdal python bindings. One of the parameters that I can pass to ESRI's MakeFeatureLayer function is an SQL Query. This creates a featurelayer (in memory) which contains only those rows which meet the crit

[gdal-dev] OSR Projection Error 6

Morning all, I am struggling to get a coordinate transformation to work in a python script. Below is the code snippet that is causing the following error: "ERROR 6: No translation for 'Azimuthal_Equidistant' to PROJ.4 format is known." I am pulling the input spatial reference from the input sha

Re: [gdal-dev] Enabling GEOS with GDAL on Windows

Chris, Shapely is my backup plan. I have, essentially, wrapped shapely in gdal using the load and dump wkt options before, but would love to be able to access geos without any work arounds. Thanks! Jay On Wed, Sep 22, 2010 at 2:02 PM, Christopher Barker wrote: > Jay L. wrote: > >

[gdal-dev] Enabling GEOS with GDAL on Windows

Afternoon all, I have both FWTools and GDAL 1.6 (with the python bindings) installed (windows7). I have a pair of scripts that I would like to use the GEOS functionality with (intersection and buffer). I believe that FWTools installed GEOS. Is there a way to use the GEOS installed FWTools from

[gdal-dev] Shapefile - Multipolygon with interior rings

Good afternoon, I am working on a program to convert a shapefile from ocentric to ographic latitude. The program reads an ESRI shapefile, creates a new shapefile, copying the attribute data, and performs a latitude "shift" for each point in the geometry. OGR is working wonderfully on all geometry