[gdal-dev] Getting data values from a first GeoTIFF and do some map algebra type operations with a second GeoTIFF

2011-11-23 Thread jdmorgan
Hello, I am new to GDAL and raster programming in general and hoping that I can get some guidance. I have GeoTIFF which are projected in WGS84 and a single band.I am trying to start simply with the concept of getting out the data values from the GeoTIFF from one file so that I compare them to

[gdal-dev] gdalinfo on win 7 64-bit

2011-11-27 Thread jdmorgan
Hello, I am running gdal on windows 7 64-bit which I installed via OSGeo4W. When I run gdalinfo on a raster GeoTiff file I get the following runtime error: RS6034 An application has made an attempt to load the C runtime library incorrectly. Any ideas one why this might be occuring, or how I

[gdal-dev] gdalinfo on win 7 64-bit

2011-11-28 Thread jdmorgan
Please disregard my question on this topic. I was able to resolve this issue on my own. For the record what I did was open the OSGeo4W installer and re-install all of the libraries which seems to have cleaned on some dll conflict I am guessing. Thanks and I am digging this forum. -- JD Morga

[gdal-dev] editor that color codes mapfile syntax?

2011-12-13 Thread jdmorgan
Hello, I was wondering if anyone knows of an editor that color codes mapfile syntax? I use notepad++ currently for python, but it would be great to have sytax coloring for .map files. Thanks, -- JDM Morgan @ NEMAC ___ gdal-dev mailing list gdal-de

[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 inpu

Re: [gdal-dev] comparing two rasters

2012-02-22 Thread jdmorgan
ta1>0]=1 ... data2[data2>0]=1 On Wed, Feb 22, 2012 at 8:46 AM, jdmorgan <mailto: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 resolution.Then I re-assign any values that are greater t

Re: [gdal-dev] comparing two rasters

2012-02-23 Thread jdmorgan
: " + str(data2) result_bools = np.logical_and((data1 > 0), (data2 > 0)) result_ints = np.array(result_bools, dtype=int) print "result_ints: " + str(result_ints) dst_ds.GetRasterBand(1).WriteArray(result_ints, 0, r) dst_ds = None On Wed, Feb 22, 2012 at 5:26 PM, jdmorg

[gdal-dev] color-code/class the raster in gdal/python

2012-02-23 Thread jdmorgan
Hello, I am creating a raster file (.img) in python/gdal using.GetRasterBand(1).WriteArray it is seems to be working well.Where I could use some help is in knowing how to color-code/class the resulting raster.For instance my raster is only ones and zeros.I would like to make the ones white an

Re: [gdal-dev] comparing two rasters

2012-02-23 Thread jdmorgan
method SetRasterColorTable(). You can find some sample python code in GDAL's test suite[1][2]. [1]: http://trac.osgeo.org/gdal/browser/trunk/autotest/gcore/colortable.py [2]: http://trac.osgeo.org/gdal/browser/trunk/autotest/gcore/tiff_write.py#L916 On Thu, Feb 23, 2012 at 7:45 PM, jdm

[gdal-dev] possible to take this unprojected/unknown coordinate system data and assign it to the second known ones projected/coordinate system

2012-03-27 Thread jdmorgan
Hello, I was hoping someone in the GDAL community could help me out, as I am stuck on a particular issue of dealing with a source dataset (geotiff) that has no projection information.Here are the gdalinfo results on the coordinate system: Coordinate System is: /LOCAL_CS[" Geocoding informat

Re: [gdal-dev] possible to take this unprojected/unknown coordinate system data and assign it to the second known ones projected/coordinate system

2012-03-31 Thread jdmorgan
in worrying about datum so do something like: gdal_translate -a_srs WGS84 in.tif out.tif After that you should be able to use gdalwarp to reproject normally as long as you know the coordinate system you want to go to. Best regards, Frank On Tue, Mar 27, 2012 at 11:20 AM, jdmorgan wrote

[gdal-dev] georeference and project a jpeg

2012-05-24 Thread jdmorgan
Hello, I am attempting to georeference and project a jpeg image using only gdal tool.I have georeferenced the jpeg using the following command: gdal_translate -a_srs EPSG:102719-of HFA -gcp 0.0 0.0 -82.5586187839508 35.59414007259327 -gcp 0.0 410.0 -82.55858659744263 35.5937998255945 -gcp 52

[gdal-dev] Re: gdal-dev Digest, Vol 96, Issue 55

2012-05-25 Thread jdmorgan
eoreference and project a jpeg To: jdmorgan Cc: gdal-dev@lists.osgeo.org Message-ID: Content-Type: text/plain; charset="utf-8" Derek, Please provide the gdalinfo output of the img file. Also, try this: Create a vrt file instead of a HFA using the same gdal_translate command. See

Re: [gdal-dev] georeference and project a jpeg

2012-05-25 Thread jdmorgan
ee if gdalwarp works with the vrt file. On Fri, May 25, 2012 at 7:17 AM, jdmorgan <mailto:jdmor...@unca.edu>> wrote: Hello, I am attempting to georeference and project a jpeg image using only gdal tool.I have georeferenced the jpeg using the following command: gdal_tr

[gdal-dev] zonal statistics with gdal

2012-06-19 Thread jdmorgan
Hello, Is there a way to do something along the lines of zonal statistics with gdal.Basically, I have a county shape file and a classified raster dataset.I would like to count the number of raster pixels/county (FIPS).Is this possible with GDAL/OGR?If so can someone point me in the right dire

Re: [gdal-dev] zonal statistics with gdal

2012-06-23 Thread jdmorgan
.ca.gov project. Petteri On Tue, 19 Jun 2012, jdmorgan wrote: Hello, Is there a way to do something along the lines of zonal statistics with gdal.Basically, I have a county shape file and a classified raster dataset.I would like to count the number of raster pixels/county (FIPS).Is this possibl

[gdal-dev] add string field to img

2012-07-24 Thread jdmorgan
Hello, I was hoping someone could help me out. I have a raster .img file. with existing pixel values. Those values are integer, but are classed to string values not currently in the same .img file. For example pixel value 2 means ClassTwo. Is it possible to add the string values (a new fiel