[gdal-dev] gdaladdo errors for overviews on vrt file

2013-03-06 Thread Pietro Rossin
Hi all I'm trying to make overviews on a vrt file. Windows XP x32, dual opteron 4Gb RAM, Osgeo4W. This file is a virtual dataset, a mosaic of 993 tiff images 1bit b/w. gdalinfo for one of these tiff is: Driver: GTiff/GeoTIFF Files: 109021.tif

[gdal-dev] Help me about gdal2tiles

2013-03-06 Thread 市川 雄二
Hi. I want to make data as TMS from NOAA GLOBE, but am in trouble at gdal2tiles.py stage. Here are what I did. 1. installed GDAL library (1.9.2) by HomeBrew on MacOS X(10.8.2) 2. installed gdal by pip on python2.7.3 3. downloaded all data from http://www.ngdc.noaa.gov/mgg/topo/gltiles.html and

Re: [gdal-dev] OGR distance different than in ArcGIS

2013-03-06 Thread Hare, Trent
Uli, Unless your neighbors are very localized (in an undistorted projection), I would recommend using some geodesic methods (in ArcMap or OGR). For ArcMap, you can beat this free add-on (Tools for Graphics and Shapes) which does geodesic lengths, areas, and azimuths no matter the input projecti

Re: [gdal-dev] OGR distance different than in ArcGIS

2013-03-06 Thread Uli Strötz
Thanks for your answers. I tried number two and the files are projected in Albers (meters) (see attachment). The results are still the same as being unprojected. Is there a way to tell GDAL to only use the Projected Coordinate System instead of the CRS? _

Re: [gdal-dev] OGR distance different than in ArcGIS

2013-03-06 Thread Frank Warmerdam
On Wed, Mar 6, 2013 at 3:59 PM, Uli Strötz wrote: > Hi Frank, > > The distance is according to ArcGIS 500m. According to Python(OGR) > 29215.31. I was thinking the same, that GDAL gives the distance according > to the Geographic Coordinate System. Is there a way it uses the units of > the Project

Re: [gdal-dev] Python API

2013-03-06 Thread Frank Warmerdam
On Wed, Mar 6, 2013 at 11:01 AM, Jeremy Palmer wrote: > Hi All, > > I've got a few questions about the support and future of some method calls > within the OGR Python API. > > 1/ What is preferred way to method to access geometries from features? is > it feature.geometry() or feature.GetGeometryR

Re: [gdal-dev] OGR distance different than in ArcGIS

2013-03-06 Thread Uli Strötz
Hi Frank, The distance is according to ArcGIS 500m. According to Python(OGR) 29215.31. I was thinking the same, that GDAL gives the distance according to the Geographic Coordinate System. Is there a way it uses the units of the Projected Coordinate System (which would be meters)? Or do you know of

Re: [gdal-dev] OGR distance different than in ArcGIS

2013-03-06 Thread Frank Warmerdam
Uli, You might want to elaborate on in what ways the results differ. Perhaps for a concrete example pair of geometries. The GDAL/OGR Distance is going to give planar distance in "decimal degrees" if everything is in geographic coordinates without regard to ellipsoidal geometries. I suspect the

[gdal-dev] OGR distance different than in ArcGIS

2013-03-06 Thread Uli Strötz
I want to determine the shortest distance between a point and a line. I use the following Python script to do so. If I compare the result of the script with the result of ArcMap's Near Analysis, I get a completely different results. The shapefiles are in the CRS GCS_North_American_1983 and projecte

Re: [gdal-dev] nodata in .bil format

2013-03-06 Thread Joaquim Luis
Jan, My mail used copy-paste so the spaces are the same as printed by gdalinfo but I don't think that is important as it merely reflects how gdalinfo prints the info. Actually I crossed this problem when my gdalread mex lib showed that the nodata for the .bil case was 1e10 instead of the -999

Re: [gdal-dev] nodata in .bil format

2013-03-06 Thread Jan Heckman
Hi, There should certainly not be a space between Nodata and Value as in your mail below. There are two formats relevant for (at least) arcgis: .bil. In this case the keyword is NoData .flt. In this case the keyword is NoData_Value. Capitalization, afaik, is irrelevant. So you could modify the .hdr

Re: [gdal-dev] WFS Axis order issue with projected coordinate system

2013-03-06 Thread Jeremy Palmer
Thanks Even. So does this ticket account for this problem? http://trac.osgeo.org/gdal/ticket/4329 Or do changes need to be applied in the WFS driver as well? Yes the ISO/OGC people have made my job very hard here in NZ. As NZ's leading geospatial government department LINZ is tasked with pushi

[gdal-dev] Python API

2013-03-06 Thread Jeremy Palmer
Hi All, I've got a few questions about the support and future of some method calls within the OGR Python API. 1/ What is preferred way to method to access geometries from features? is it feature.geometry() or feature.GetGeometryRef()? "geometry()" seems funny as it's not well documented and is

[gdal-dev] nodata in .bil format

2013-03-06 Thread Joaquim Luis
Hi, While playing a bit with Lastools I came across this behavior. I create outputs in both .bil and .asc ESRI formats The weird thing is that beeing the .hdr header file common (I mean exactly equal) to both outputs the .bil format does not recognize the nodata value. Is this a knwon issue?

Re: [gdal-dev] append mode in postgresql

2013-03-06 Thread Even Rouault
Selon "Yves Jacolin (Free)" : > thanks Even. All information is in the pg_driver page. I really need to read > it more carefully :( > > There is just a very small issue, the doc said about SCHEMA: > "Use the -nln option of ogr2ogr instead, or better the active_schema " > > But there is no example

Re: [gdal-dev] append mode in postgresql

2013-03-06 Thread Yves Jacolin (Free)
thanks Even. All information is in the pg_driver page. I really need to read it more carefully :( There is just a very small issue, the doc said about SCHEMA: "Use the -nln option of ogr2ogr instead, or better the active_schema " But there is no example about SCHEMA :/ Y. Le mercredi 6 mars 201

Re: [gdal-dev] append mode in postgresql

2013-03-06 Thread Even Rouault
Selon "Yves Jacolin (Free)" : > Le mercredi 6 mars 2013 15:38:11 Yves Jacolin a écrit : > > Hello, > > > > I am trying to import EDIGEO files in append mode. The doc said to use > > -append flag with ogr2ogr commande. > > > > ogr2ogr -gt 65536 -f "PostgreSQL" "PG:dbname=test_edigeo host='localhos

Re: [gdal-dev] EDIGEO format

2013-03-06 Thread Yves Jacolin (Free)
Le jeudi 28 février 2013 18:51:17 Jürgen E. Fischer a écrit : > Hi Yves, > > On Thu, 28. Feb 2013 at 16:16:02 +0100, Yves Jacolin (Free) wrote: > > I am working on EDIGEO files, is this something I have to expect: > Not sure what you question is. > > > $ ogrinfo *.THF > > Had to open data source

Re: [gdal-dev] append mode in postgresql

2013-03-06 Thread Yves Jacolin (Free)
Le mercredi 6 mars 2013 15:38:11 Yves Jacolin a écrit : > Hello, > > I am trying to import EDIGEO files in append mode. The doc said to use > -append flag with ogr2ogr commande. > > ogr2ogr -gt 65536 -f "PostgreSQL" "PG:dbname=test_edigeo host='localhost'" > - lco "GEOMETRY_NAME=the_geom" -lco "

[gdal-dev] append mode in postgresql

2013-03-06 Thread Yves Jacolin (Free)
Hello, I am trying to import EDIGEO files in append mode. The doc said to use -append flag with ogr2ogr commande. ogr2ogr -gt 65536 -f "PostgreSQL" "PG:dbname=test_edigeo host='localhost'" - lco "GEOMETRY_NAME=the_geom" -lco "SCHEMA=test" -lco "OVERWRITE=YES" -nln EDIGEO_parcelle E0001.THF PA

Re: [gdal-dev] INSPIRE compliant metadata

2013-03-06 Thread laura0
Thanks for the answers, Tom, can you suggest me links (if any) to help me to follow your idea? Thanks, Laura -- View this message in context: http://osgeo-org.1560.n6.nabble.com/INSPIRE-compliant-metadata-tp5038779p5038802.html Sent from the GDAL - Dev mailing list archive at Nabble.com. _

Re: [gdal-dev] INSPIRE compliant metadata

2013-03-06 Thread Kralidis,Tom [Ontario]
One approach could be a script to use the GDAL/OGR libraries to generate the XML document, given that there are metadata elements that you can't derive just from the data. ..Tom -Original Message- From: gdal-dev-boun...@lists.osgeo.org on behalf of laura0 Sent: Wed 06-Mar-13 08:29 To

Re: [gdal-dev] INSPIRE compliant metadata

2013-03-06 Thread Margherita Di Leo
Hi Laura, On Wed, Mar 6, 2013 at 2:29 PM, laura0 wrote: > Hello, > I am using the GDAL Utility Program gdalinfo to report metadata about a > raster file (netCDF, grib, hdf). I would like to transform these > information > collected using the inspire directive (ogc csw). Is there a tool to > auto

[gdal-dev] INSPIRE compliant metadata

2013-03-06 Thread laura0
Hello, I am using the GDAL Utility Program gdalinfo to report metadata about a raster file (netCDF, grib, hdf). I would like to transform these information collected using the inspire directive (ogc csw). Is there a tool to automate this translation? I have found CatMDEdit but it does not support t

[gdal-dev] New Driver: DDS

2013-03-06 Thread Alan Boudreault
Hi, I'd like to get your hints on a new DDS driver. Please, take a look at my brief ticket and let me know if there is anything wrong or if it's ok to commit. http://trac.osgeo.org/gdal/ticket/5017 Regards, Alan -- Alan Boudreault http://www.mapgears.com/ ___

Re: [gdal-dev] GDAL fails on reading rset via Virtual File

2013-03-06 Thread Even Rouault
Selon Jayesh Chaudhary : > Folks, > > I am using GDAL Virtual File mechanism to read data off a stream (something > like what's shown below) > > VSILFILE* fpTemp = VSIFileFromMemBuffer(szTempFileName, (GByte*) > s->file_size(), FALSE); > this->gdal_data_ = (GDALDataset*) GDALOpen(szTempFileName, G