[gdal-dev] GeoTransform confusion

2015-11-18 Thread Simon Lyngby Kokkendorff
Hi List, I have a geoid in AAIGRID-format, with the following header: ncols 1001 nrows 731 xllcenter 40.0 yllcenter 6045000.0 cellsize 500.0 NODATA_value GDAL will read this just fine, and report: Driver: AAIGrid/Arc/Info ASCII Grid Files: dkgeoid13b.utm32 Size is 1001, 731 Coordinat

Re: [gdal-dev] Limit on the number of subdatasets

2015-11-18 Thread Ivan Lucena
The georaster driver uses GDAL subdataset to browse the database as you can see here: gdal.org/frmt_georaster.html section "Browsing the database for GeoRasters" That is not a big deal because there are other ways go around it. But it can get a little confusing for users accustomed to use that

Re: [gdal-dev] GeoPDF offsets from actual location after gdalwarp

2015-11-18 Thread Even Rouault
Le mercredi 18 novembre 2015 19:42:30, Gane R a écrit : > Here is link to a GeoPDF ( > https://drive.google.com/file/d/0B4shTM1bPS5TMmk2OWxuTGJpY1k/view?usp=shari > ng) this has a mile offset from the current location, > on GDAL gdal-2.0.1 when gdalwarp warps the PDF (poppler driver) to GeoTIF > >

Re: [gdal-dev] Limit on the number of subdatasets

2015-11-18 Thread Even Rouault
Le mercredi 18 novembre 2015 20:22:38, Ivan Lucena a écrit : > Hi there. > > > Is there a limit on the number of subdatasets that GDAL can handle? There's no hard-coded limit AFAIK. Should be memory bound. But performance of various stuff handing lists might not be great with a big number. >

[gdal-dev] GeoPDF offsets from actual location after gdalwarp

2015-11-18 Thread Gane R
Here is link to a GeoPDF ( https://drive.google.com/file/d/0B4shTM1bPS5TMmk2OWxuTGJpY1k/view?usp=sharing) this has a mile offset from the current location, on GDAL gdal-2.0.1 when gdalwarp warps the PDF (poppler driver) to GeoTIF Here is gdalinfo for the file Driver: PDF/Geospatial PDF Files: Ass

Re: [gdal-dev] Standard Parallel 1 - being changed but not displayed?

2015-11-18 Thread Jonathan Moules
Hi Even,Ok, that clears things up a little. I guess the real question I have at this point is - which is correct in this circumstance? The original GeoTIFF? ArcGIS? GDAL's conversion? None of the above? I get the feeling it's all a bit subjective. >From your comments, I'm surprised this does

Re: [gdal-dev] GCPsToGeoTransform returns wrong values

2015-11-18 Thread Bar Ori
I found in did that the GDALGCPsToGeoTransform using the pixels values for calculation: At gdal_misc.cpp : padfGeoTransform[1] = (pasGCPs[1].dfGCPX - pasGCPs[0].dfGCPX) / (pasGCPs[1].dfGCPPixel - pasGCPs[0].dfGCPPixel); padfGeoTransform[0] = pasGCPs[0].dfGCPX - pasGCPs[0].dfGCPPixel * padfGeoTr

Re: [gdal-dev] GCPsToGeoTransform returns wrong values

2015-11-18 Thread Ari Jolma
On 18.11.2015 11:59, Bar Ori wrote: Thanks Ari for your reply. In my situation, sometimes I am working with rectangles and sometimes not. Depends the source of information. I thought to always use the GCPs as follows: In case my shape has rectangular form, so the dataSet.GetGCPCount() wi

Re: [gdal-dev] GCPsToGeoTransform returns wrong values

2015-11-18 Thread Bar Ori
Thanks Ari for your reply. In my situation, sometimes I am working with rectangles and sometimes not. Depends the source of information. I thought to always use the GCPs as follows: In case my shape has rectangular form, so the dataSet.GetGCPCount() will return me “0” and I will use the dataSet.

Re: [gdal-dev] GCPsToGeoTransform returns wrong values

2015-11-18 Thread Ari Jolma
Bar, I think you are setting the GCPs up wrongly if you want them to be at the corners. GCP(LeftDown.X, LeftDown.Y, LeftDown.Z, Height - 0.5, null, "4") does not have 0 before Height-0.5. GCPsToGeoTransform checks the special case of GCPs in rectangular shape and in top-left, top-right, bo