[gdal-dev] Real and "raster" coordinates

Hi, I have a georeferenced raster and coordinates of some point (lat, lon) inside raster. Is it possible to convert this lat/lon coordinates to raster coordinates (column and row)? I use GDAL 1.6.3 with Python 2.5 (installed with OSGeo4W installer). -- Alexander Bruy ___

[gdal-dev] Re: How to divide a raster layer into a set of tiles

A short question Why did you not use 1) gdal_merge 2) gdal_translate 3) gdal_retile Btw, if there is no need for reprojection, gdal_retile alone does the job. Did you have problems using these utilities ? Eli Adam writes: Gilles, I do a similar process (build large vrt mosaic, re

[gdal-dev] Maintenance report for Feb 24, 2010 to Mar 23, 2010

Hi, I updated my GDAL maintenance report at http://trac.osgeo.org/gdal/wiki/MaintenanceReportsByChaitany I am also attaching a copy. -- Best regards, Chaitanya kumar CH. /tʃaɪθənjə/ /kʊmɑr/ +91-9848167848 17.2416N 80.1426E GDAL Maint

Re: [gdal-dev] VRT creation with Python

Michael Aye wrote: Dear all, I have some, as of yet unpublished, PDS data I am reading via a well working simple vrt definition looking like this: FC1_2005-06-23T21.55.54.000Z_ID10_019407_F1.img 8704 2 2184 LSB which is saved as a FC1.vrt file. How can I programmatica

[gdal-dev] Re: Gdal / Python / Win7 crash with ReadAsArray()

I am running the 32bit Versions of python and gdal. bye Nop -- View this message in context: http://n2.nabble.com/Gdal-Python-Win7-crash-with-ReadAsArray-tp4398401p4788484.html Sent from the GDAL - Dev mailing list archive at Nabble.com. _

[gdal-dev] Re: How to divide a raster layer into a set of tiles

Gilles, I do a similar process (build large vrt mosaic, reproject the mosaic vrt to another vrt, tile the reprojected mosaic vrt in new projection to tif) although my end tiles are in tif and not compressed. My method is the brute force way and not elegant at all. You may also want to know t

[gdal-dev] gdal2tiles.py "direct warping" issue?

I'm looking through gdal2tiles.py to get a feel for how it works and any potential issues I might encounter. Around line 801 of the script in GDAL 1.7.1 (right after the call to gdal.AutoCreateWarpedVRT) there is a comment that says: TODO: HIGH PRIORITY: Correction of AutoCreateWarpedVRT accordi

[gdal-dev] gdal2tiles.py "direct warping" issue?

I'm looking through gdal2tiles.py to get a feel for what it's doing and any potential issues I might encounter. Around line 801 of the script in GDAL 1.7.1 (right after the call to gdal.AutoCreateWarpedVRT) there is a comment that says: TODO: HIGH PRIORITY: Correction of AutoCreateWarpedVRT acco

[gdal-dev] VRT creation with Python

Dear all, I have some, as of yet unpublished, PDS data I am reading via a well working simple vrt definition looking like this: FC1_2005-06-23T21.55.54.000Z_ID10_019407_F1.img 8704 2 2184 LSB which is saved as a FC1.vrt file. How can I programmatically opening many mor

Re: [gdal-dev] Re: Gdal / Python / Win7 crash with ReadAsArray()

Hi, You didn't mention whether you are running against the x64 version of python.exe or the x86 version? I also see a problem when running the x64 version of python with numpy installed could crash gdal (x64), however without installing numpy it should work properly. > It seems that this is a Wi

[gdal-dev] Re: Gdal / Python / Win7 crash with ReadAsArray()

Hi! Eventually I had some time to look into this again. First thing I tried is updating python to 2.6.5 - no change, still crashes. Then I tried running exactly the same files in a 32bit XP virtual machine. It works there. This is a workaround, but not a solution as I cannot run the whole work

Re: [gdal-dev] Re: GDAL custom build (ECW required)

Frank Warmerdam wrote: > Gilles Bassière wrote: >>> with COMPRESS=JPEG: >>> >>> Creating output file that is 1P x 1L. >>> Processing input file ortho-dept04-2004-la2e.vrt. >>> 0...10...20...30..gdalwarp: tif_jpeg.c:691: JPEGPreDecode: Assertion >>> `sp->cinfo.comm.is_decompressor' failed

[gdal-dev] Re: How can I display the extent of a shapefile when using ogrinfo or from another command line tool?

John Mitchell gmail.com> writes: > > Jukka Rahkonen mmmtike.fi> writes: > > Thanks, > > That worked great!!! Anyone know why I did not receive the below response to my > question to my email account? I had to search via google to find this response. Hi, One reason may be that I did not

Re: [gdal-dev] Re: Bug in coordinate system library? unknown elliptical parameter for BNG coordinate system

Just to remind you: I reported a similar problem when 1.7 just came out. Maybe it helps in identifying the problem. http://www.mail-archive.com/gdal-dev@lists.osgeo.org/msg05605.html Kind regards, Alex On Mar 23 2010, Jean-Claude REPETTO wrote: Andreas Joisten a écrit : As you can see I did

[gdal-dev] Re: Bug in coordinate system library? unknown elliptical parameter for BNG coordinate system

Andreas Joisten a écrit : As you can see I did not specifiy the ellipsoid parameter: gdalwarp -s_srs EPSG:27700 -t_srs EPSG:4326 input_geotiff_in_BNG.tiff output_geotiff.tiff Sorry, but you didn't give this command line in your first message. the gdal error message was: ERROR 6: Failed

Re: [gdal-dev] Bug in coordinate system library? unknown elliptical parameter for BNG coordinate system

Andreas Joisten wrote: Hi Jean-Claude, many thanks for your answer. As you can see I did not specifiy the ellipsoid parameter: gdalwarp -s_srs EPSG:27700 -t_srs EPSG:4326 input_geotiff_in_BNG.tiff output_geotiff.tiff the gdal error message was: ERROR 6: Failed to initialize PROJ.4 with `+p

Re: [gdal-dev] Re: GDAL custom build (ECW required)

Gilles Bassière wrote: with COMPRESS=JPEG: Creating output file that is 1P x 1L. Processing input file ortho-dept04-2004-la2e.vrt. 0...10...20...30..gdalwarp: tif_jpeg.c:691: JPEGPreDecode: Assertion `sp->cinfo.comm.is_decompressor' failed. Aborted Gilles, The likely solution to t

[gdal-dev] Bug in coordinate system library? unknown elliptical parameter for BNG coordinate system

Hi Jean-Claude, many thanks for your answer. As you can see I did not specifiy the ellipsoid parameter: gdalwarp -s_srs EPSG:27700 -t_srs EPSG:4326 input_geotiff_in_BNG.tiff output_geotiff.tiff the gdal error message was: ERROR 6: Failed to initialize PROJ.4 with `+proj=tmerc +lat_0=49 +lon_

[gdal-dev] Re: GDAL custom build (ECW required)

Gilles Bassière wrote: > Hi (again), > > I'm facing problems with GeoTiff compression on a custom build of GDAL. > > My source dataset is a virtual raster composed of ECW tiles built with: > gdalbuildvrt ortho-dept04-2004-la2e.vrt ortho-dept04-2004-la2e/*.ecw > > The error comes when I try to e

Re: [gdal-dev] Problem with +ellps=airy which is used for BNG

Andreas Joisten a écrit : Hi all, The BNG (British Natial Grid) based on the airy Ellipsoid but GDAL did not know this parameter. +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=40 +y_0=-10 +ellps=airy +datum=OSGB36 +units=m +no_defs Probably, gdal does not know this parameter

[gdal-dev] Problem with +ellps=airy which is used for BNG

Hi all, The BNG (British Natial Grid) based on the airy Ellipsoid but GDAL did not know this parameter. +proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=40 +y_0=-10 +ellps=airy +datum=OSGB36 +units=m +no_defs Probably, gdal does not know this parameter: +ellps=airy But gdal kno

[gdal-dev] Re: How to divide a raster layer into a set of tiles?

Gilles Bassière gmail.com> writes: > > Hi, > > I'm trying to reproject an aerial imagery dataset with GDAL. > > The initial dataset is composed of ECW tiles in EPSG:27572 and my goal > is to compute a new set set of ECW tiles projected in EPSG:2154. Hi, May I ask if you are ready with your n

[gdal-dev] Re: How can I display the extent of a shapefile when using ogrinfo or from another command line tool?

Jukka Rahkonen mmmtike.fi> writes: Thanks, That worked great!!! Anyone know why I did not receive the below response to my question to my email account? I had to search via google to find this response. John ___ gdal-dev mailing list gdal-dev@

Re: [gdal-dev] Re: How to divide a raster layer into a set of tiles?

Jukka Rahkonen wrote: > Gilles Bassière gmail.com> writes: > >> Hi, >> >> I'm trying to reproject an aerial imagery dataset with GDAL. >> >> The initial dataset is composed of ECW tiles in EPSG:27572 and my goal >> is to compute a new set set of ECW tiles projected in EPSG:2154. > ... > >> My cu

[gdal-dev] GDAL custom build (ECW required)

Hi (again), I'm facing problems with GeoTiff compression on a custom build of GDAL. My source dataset is a virtual raster composed of ECW tiles built with: gdalbuildvrt ortho-dept04-2004-la2e.vrt ortho-dept04-2004-la2e/*.ecw The error comes when I try to extract and reproject some portions of t

[gdal-dev] Re: How to divide a raster layer into a set of tiles?

Christian Müller wrote: > That is a surprise. > In all my tests I had no memory problems with gdal_retile. Sometimes I > had some memory problems with gdal_merge. The biggest input file for > gdal_retile was a 80 GB Erdas image (not compressed) and it works fine. > I have never testet with an ECW F

[gdal-dev] ERROR 6: Failed to initialize PROJ.4 -- convert from BNG to WGS84; using gdalwarp

Hi all, I want to use gdalwarp to convert a GeoTiff file in BNG (British National Grid) to WGS84, but I got the following error: ERROR 6: Failed to initialize PROJ.4 with `+proj=tmerc +lat_0=49 +lon_0=-2 +k=0. 9996012717 +x_0=40 +y_0=-10 +ellps=airy +datum=OSGB36 +units=m +no_defs '. Th

[gdal-dev] how: dealing with osm data and creating custom shape files

Hi all, thanks to the hints from last week I was able to successfully create a postgresql database with osm data and from there I am able to create shape files containing all the attributes from osm. Thanks to Ralf and Frederik! But now I want to group those data together. Well, I want to create

Re: [gdal-dev] Re: How to divide a raster layer into a set of tiles?

Christian Müller wrote: > gdal_retile.py definitely creates georeferenced tiles, I developed this > utility. > I did some experiments reprojecting the single tiles while retiling, but > it simply does not work. As a consequence I canceled the reprojecting > support. > The solution for your problem

Re: [gdal-dev] How to divide a raster layer into a set of tiles?

Even Rouault wrote: > Le Monday 22 March 2010 19:35:47 Gilles Bassière, vous avez écrit : >> Hi, >> >> The problem with this workflow is that it creates non-georeferenced >> tiles (I don't know if this is the intended behaviour of this command). > > gdal_retile.py *does* create georeferenced tiles

[gdal-dev] Re: How to divide a raster layer into a set of tiles?

That is a surprise. In all my tests I had no memory problems with gdal_retile. Sometimes I had some memory problems with gdal_merge. The biggest input file for gdal_retile was a 80 GB Erdas image (not compressed) and it works fine. I have never testet with an ECW File and my tests run on a Li