RE: [gdal-dev] GetStatistics and signed 8-bit IMGs and GeoTIFFs

2009-09-23 Thread Jason Roberts
Ok, here is how a GDAL caller can detect that a file is signed 8-bit: >>> band.GetMetadata('IMAGE_STRUCTURE') {'PIXELTYPE': 'SIGNEDBYTE'} Thanks to Frank W for giving me a clue about this a week ago that I promptly forgot about. Sorry about that. So. does anyone know whether band.GetSt

Re: [gdal-dev] Re: gdal_rasterize with a NetCDF file?

2009-09-23 Thread Roger André
Hey Scott, You could extract that layer from the netCDF, burn it with gdal_rasterize as a geotiff, convert the geoTiff into a GMT grid, and then possibly add that back into the netCDF (somehow). Sounds like a pain though. Roger -- On Mon, Sep 21, 2009 at 11:49 AM, Scott Lewis wrote: > Ah, tha

RE: [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Henneke, Amanda M
Just to document the conclusion, I did figure out if I copied the datum shift grids and pointed my PROJ_LIB environment variable to that directory that seemed to fix the problem for me. Thanks everyone for you assistance! I really appreciate it! Amanda M. Henneke IG/DB Tools - TSS Capability C

Re: [gdal-dev] Output MapInfo File with CoordSys Information

2009-09-23 Thread Reinaldo escada Chohfi
Hi Frank, I tested the program below but it does not out put the .tab or .mif/.mid with the coordinate system.  It is always non-earth.  Could you please check it against your gdal/ogr installation to see if something wrong happened during my compilation of gdal 1.6.2? Thanks a lot, Rei

[gdal-dev] RasterIO antialiased subsampling

2009-09-23 Thread Tamas Szekeres
Folks, I'm not quite sure whether this issue have already been raised or not, but I'd be interested in having on-the fly antialiasing support when the raster is subsampled within a RasterIO operation. I'm just trying to display OpenStreetMap rasters provided by the GDAL TMS mini driver in MapServe

RE: [gdal-dev] GetStatistics and signed 8-bit IMGs and GeoTIFFs

2009-09-23 Thread Jason Roberts
Someone suggested that I try this with 1.6.2 without saying whether it might be a bug or by design. The behavior is the same with 1.6.2, using binaries produced by Tamas Szekeres (thanks Tamas, for providing those on your website). Back to the original question: is it a bug or by design? How do

[gdal-dev] GetStatistics and signed 8-bit IMGs and GeoTIFFs

2009-09-23 Thread Jason Roberts
Is band.GetStatistics() supposed to work with signed 8-bit integer rasters with the HFA or GTiff driver? With GDAL 1.6.0, I can successfully create a signed 8-bit raster with the HFA or GTiff driver using the PIXELTYPE=SIGNEDBYTE option: >>> a = numpy.array([[-128, -1, 0, 1, 127]], dtype='i

RE: [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Craig Miller
I looked at Amanda's files. It looks like it is getting reprojected from UTM to Geographic, but that the datum isn't getting transformed from NAD27 to WGS84. I suggested she verify her PROJ environment variables/NAD data files. Craig > -Original Message- > From: gdal-dev-boun...@lists.

RE: [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Henneke, Amanda M
Even- Thanks! I do get something different from you: -117 40.6527598726051 0 I'm using FWTools 2.3.0, so is there something I can do then to fix this problem since I didn't build my own libraries? Thanks, Amanda Amanda M. Henneke IG/DB Tools - TSS Capability Center When One tugs at a sin

RE: [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Even Rouault
Amanda, You didn't specify why you believe that the coordinates you get are not right. Are you validating against another reference tool ? Could you give an example of input, the output you get and the output you expect. Here's what I've tried with the gdaltransform utility (in latest FWTools 2.4

RE: [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Steve . Toutant
Maybe the problem is in your data. Di you try to reproject zone 11 nad27 with the wkt nad83 below? "Henneke, Amanda M" 23/09/2009 01:44 PM A "steve.tout...@inspq.qc.ca" cc "gdal-dev@lists.osgeo.org" , "gdal-dev-boun...@lists.osgeo.org" Objet RE: [gdal-dev] Translating between UTM and Ge

RE: [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Henneke, Amanda M
PROJCS["NAD_1983_UTM_Zone_12N", GEOGCS["GCS_North_American_1983", DATUM["North_American_Datum_1983", SPHEROID["GRS_1980",6378137,298.257222101]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["cen

RE: [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Steve . Toutant
What is the wkt you are using for zone 12 nad83? "Henneke, Amanda M" @lists.osgeo.org Envoyé par : gdal-dev-boun...@lists.osgeo.org 23/09/2009 01:39 PM A "gdal-dev@lists.osgeo.org" cc Objet RE: [gdal-dev] Translating between UTM and Geographic Here is some more information: So, h

RE: [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Henneke, Amanda M
Here is some more information: So, here is the WKT I'm getting from my code for the "problem" source file: PROJCS["NAD_1927_UTM_Zone_11N", GEOGCS["GCS_North_American_1927", DATUM["North_American_Datum_1927", SPHEROID["Clarke_1866",6378206.4,294.9786982]], PRIMEM["Greenwich",0], UNIT

RE: [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Henneke, Amanda M
Craig- I'm trying to have it in Geograhphic/WGS84, at least that's what I'm specifying in the code...I agree I think it is a datum error, but I don't know what to do to make it translate the datum from NAD27 to WGS84 correctly. I don't know if I can send the raster, but I can send the two shap

RE [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Steve . Toutant
The problem might not be in your code but in the WKT projection information (sSpRefParams). Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques biologiques, environnementaux et occupationnels Institut national de santé publique du Québec 945, avenue Wolfe Qué

RE: RE [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Steve . Toutant
I don't know about the code I'm not a programmer. You said UTM Zone 12 with NAD83 works but this Zone 11 NAD27 doesn't. My guess is you are using Zone 11 NAD83 instead NAD27. I might be wrong though. Steve Steve Toutant, M. Sc. Analyste en géomatique Secteur environnement Direction des risques

RE [gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Steve . Toutant
hi, slight shift in the points (to the east) and an extremely slight shift up (to the north). That seems to be a datum error. Make sure you set the good parameters for datum when using nad83 or nad27. Steve "Henneke, Amanda M" @lists.osgeo.org Envoyé par : gdal-dev-boun...@lists.osgeo.org

[gdal-dev] Translating between UTM and Geographic

2009-09-23 Thread Henneke, Amanda M
Hello all- I'm trying to figure out how to reproject some points from a UTM projection to Geographic (WGS84) projection. This works fine when I'm working in UTM Zone 12 with NAD83. When I try the same code for Zone 11 NAD27, I get a slight shift in the points (to the east) and an extremely sl

[gdal-dev] raster extents and Python API

2009-09-23 Thread Kralidis,Tom [Ontario]
Hi, Using the Python API, I am writing a script to get high level information on datasets to feed metadata and publishing. When fetching extents, osgeo.ogr has GetExtent(). Doesn't look like there is anything like this within osgeo.gdal (thanks crschmidt). One would have to call GetGeoTransfor

Re: [gdal-dev] OGRGeometry Distance

2009-09-23 Thread Mateusz Loskot
Christopher Barker wrote: Yilmaz Arslanoglu wrote: I wanted to compute the distance between two points, you need to be clear by what you mean here: I suspect you mean: "The shortest distance between two points following the surface of the earth" In which case, what you are looking for is r

Re: [Gdal-dev] Programming with GDAL in Windows XP using DEV-C++ or CodeBlocks

2009-09-23 Thread Mateusz Loskot
btaylor wrote: Has anyone used GDAL libraries in a windows compiling environment? Yes, I have. Here is a simple example of Windows application (a tiny converter) developed using GDAL: http://github.com/mloskot/workshop/tree/master/sid2jp2/trunk/ Best regards, -- Mateusz Loskot, http://mate