[gdal-dev] Re: satellite image processing in Python

2009-02-10 Thread Yann Chemin
Sorry, this was not the version I intended to send initially. Here it is... ___ # Landsat 7 ETM + image downlaoded from GLCF #ftp.glcf.umiacs.umd.edu/glcf/Landsat/WRS2/p127/r049/p127r049_7x20001104.ETM-EarthSat-Orthorectified/ #p127r049 2000/11/04 around Ubon Ra

[gdal-dev] satellite image processing in Python

2009-02-10 Thread Yann Chemin
Hello, Would need a small help to troubleshoot this small code: 1 - When the for loop is over, the last file is not completely closed (~480Mb instead of 509Mb) Is there a way to force it to close? (opening another file forces the previous file to be closed with a proper"?" 509Mb) 2 - A segmenta

Re: Fwd: Re[4]: [gdal-dev] gdalwarp: mercator to lambert azimutal equeal area

2009-02-10 Thread Joaquim Luis
Hi, A couple of things, 1. "Planar world" is not 44.000.000 meters wide. It follows from the meter definitions that it is 40 000 000 meters 2. Mercator maps cannot span to +/- 90 degrees latitude 3. The images in www.oera.net/How2/TextureMaps2.htm are not Mecator 4. "and changing something

Re: Fwd: Re[4]: [gdal-dev] gdalwarp: mercator to lambert azimutal equeal area

2009-02-10 Thread Even Rouault
I think you don't define enough the projections and should specify more parameters. PROJ.4 documentation and the following links may be helpful : http://www.remotesensing.org/geotiff/proj_list/mercator_1sp.html http://www.remotesensing.org/geotiff/proj_list/mercator_2sp.html http://www.remotesen

Fwd: Re[4]: [gdal-dev] gdalwarp: mercator to lambert azimutal equeal area

2009-02-10 Thread jumpjack
I tried georeferencing one of these images ( http://www.oera.net/How2/TextureMaps2.htm ) using MapWindow GIS; the pogram generated a .jgw file, which then I fixed by replacing "," by "." and changing something like 179.12345 to 180.0 . Then I tried again this command: >gdalwarp -s_srs "+proj=me

Re[4]: [gdal-dev] gdalwarp: mercator to lambert azimutal equeal area

2009-02-10 Thread jumpjack
Thanks, but I can't yet understand... I have a map of the whole world. It is 640 pixel wide. "Planar world" is 44.000.000 meters wide so which number should I specify as "A"?? The same applies for "E". About C and F: which coordinates should I specify? Lat and Lon? This would be 90,-180, I g

Re: [gdal-dev] Generating GML from Database - Problem: No decimal values

2009-02-10 Thread Gregor Fikoczek
Hi list, i got some new By generating GML from oracle, i first have to set the NLS_LANG-Parameter: export NLS_LANG=AMERICAN_AMERICA.WE8ISO8859P1 ogrinfo OCI:* -sql "SELECT * from test_decimal_values" Layer name: SELECT * from test_decimal_values Geometry: Unknown (any) Featu

[gdal-dev] Generating GML from Database - Problem: No decimal values

2009-02-10 Thread Gregor Fikoczek
Hi list, i got a problem by generating a GML-File from the Oracle- and MySQL-DB using ogr2ogr, so I alway getting in MySQL/oracle instead of a float-value an rounded integer-value (mostly 0) without any decimal numbers. Some Info: # Structure MySQL: