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,

[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

[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

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 amanda.m.henn...@boeing.com@lists.osgeo.org Envoyé par :

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

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

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],

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 amanda.m.henn...@boeing.com@lists.osgeo.org Envoyé par : gdal-dev-boun...@lists.osgeo.org 23/09/2009 01:39 PM A gdal-dev@lists.osgeo.org gdal-dev@lists.osgeo.org cc Objet RE: [gdal-dev] Translating between UTM and

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],

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 amanda.m.henn...@boeing.com 23/09/2009 01:44 PM A steve.tout...@inspq.qc.ca steve.tout...@inspq.qc.ca cc gdal-dev@lists.osgeo.org gdal-dev@lists.osgeo.org,

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

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

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] 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]],

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

[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

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,

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

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

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