[gdal-dev] Using gdal_translate to convert to KML overlay
Dear all, I'm using a similar command line like this to convert a geotiff RASTER to KMZ but gdal_translate.exe -of KMLSUPEROVERLAY c:\in.tif c:\out.kmz -co FORMAT=JPEG All I get is a totally black image. Any tips on how to solve this? REgards, Antonio ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Question about using crop_to_cutline
Hi Sorry Joaquim but I'm really trying to find a solution based on GDAL. Thanks Antonio I have a close functionality in Mirone but it was developed mainly for grids, though it works also with images. With a grid/image loaded, load the polygon too. Next right-click on the polygon line and select ROI Crop Tools - Crop Image (or any of the other related options) and it's done. One thing you may not like is that the resulting image is not of the same size (extent) as the original. This is so because I (don't remember anymore why) decided to save only the extent corresponding to the BoundingBox of the polygon. Joaquim ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Question about using crop_to_cutline
Greetings I have a raster image and a shapefile polygon and I wanted to crop the raster for just to have valid pixels inside the polygon. I have used crop_to_cutline but here is a problem: It produces a raster map but pixels does not geographically match with the original raster map. So, i produce a raster map but I'm not able to compare with original image since there is a small displacement. How can I obtain a raster map without this mismatch? Thank you Regards, Antonio GDAL version 1.9.1 ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Renaming a column's name in a Shapefile
Hi I have tried qwhat you said: First: ogrinfo -al -so GPS_Project.shp INFO: Open of `GPS_Project.shp' using driver `ESRI Shapefile' successful. Layer name: GPS_Project Geometry: Point Feature Count: 135 Extent: (643712.710081, 7794078.321332) - (649554.409908, 7796057.897650) Layer SRS WKT: PROJCS[WGS_1984_UTM_Zone_36S, GEOGCS[GCS_WGS_1984, DATUM[WGS_1984, SPHEROID[WGS_1984,6378137.0,298.257223563]], PRIMEM[Greenwich,0.0], UNIT[Degree,0.0174532925199433]], PROJECTION[Transverse_Mercator], PARAMETER[False_Easting,50.0], PARAMETER[False_Northing,1000.0], PARAMETER[Central_Meridian,33.0], PARAMETER[Scale_Factor,0.9996], PARAMETER[Latitude_Of_Origin,0.0], UNIT[Meter,1.0]] ele: Real (19.15) time: String (80.0) magvar: Real (19.15) geoidheigh: Real (19.15) name: String (80.0) cmt: String (80.0) desc_: String (80.0) src: String (80.0) link1_href: String (80.0) link1_text: String (80.0) link1_type: String (80.0) link2_href: String (80.0) link2_text: String (80.0) link2_type: String (80.0) sym: String (80.0) type: String (80.0) fix: String (80.0) sat: Integer (10.0) hdop: Real (19.15) vdop: Real (19.15) pdop: Real (19.15) ageofdgpsd: Real (19.15) dgpsid: Integer (10.0) label: String (80.0) And then: ogr2ogr -sql SELECT [ele AS ele, time AS t, magvar AS magvar, geoidheigh AS geoidheigh, name AS name, cmt AS cmt, desc_ AS desc_, src AS src] FROM GPS_Project.shp, new.shp GPS_Project.shp And I got: ERROR 1: SQL Expression Parsing Error: syntax error What is wrong in my expression? THanks and sorry Antonio Chaitanya kumar CH wrote: António, OGR comes as a part of GDAL. OGR utilities can be used like GDAL utilities. First, run ogrinfo on the shapefile with the options -al and -so to determine the field names in it. Next, use the ogr2ogr utility like this: ogr2ogr -sql SELECT [field1a AS field1b, field2a AS field2b, field 3a AS field3b] FROM shapefileA shapefileB.shp shapefileA.shp This command changes the fields field1a, field2a and field3a in the original shapefile, shapefileA.shp, to field1b, field2b and field3b respectively in the new shapefile, shapefile.shp. 2011/10/26 António Rocha antonio.ro...@deimos.com.pt mailto:antonio.ro...@deimos.com.pt Hi THanks for the feedback. I'm a bit lost on this. Is there any examples and tutorials on using sql commands for OGR? I'm sorry but I'm not use to use Vectorial data and neither OGR (only GDAL) Thanks Antonio Chaitanya kumar CH wrote: António, The easiest way is to use the -sql option in ogr2ogr [1] and set a new name with the AS clause in the SQL statement [2]. [1]: http://www.gdal.org/ogr2ogr.html [2]: http://trac.osgeo.org/gdal/wiki/rfc28_sqlfunc 2011/10/26 António Rocha antonio.ro...@deimos.com.pt mailto:antonio.ro...@deimos.com.pt mailto:antonio.ro...@deimos.com.pt mailto:antonio.ro...@deimos.com.pt Greetings I have to shage a column name in a Shapefile. I'm not much familiar with OGR. How can I do that in OGR? THanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6574 (20111025) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org mailto:gdal-dev@lists.osgeo.org mailto:gdal-dev@lists.osgeo.org mailto:gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev -- Best regards, Chaitanya kumar CH. +91-9494447584 tel:%2B91-9494447584 17.2416N 80.1426E __ Information from ESET NOD32 Antivirus, version of virus signature database 6575 (20111026) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com -- Best regards, Chaitanya kumar CH. +91-9494447584 17.2416N 80.1426E -- António Rocha Aerospace Defence Portugal DEIMOS Engenharia, SA Av. D. Joao II, Lt 1.17.01 - 10º, Ed. Torre Zen 1998-023 Lisboa, Portugal Ph: +351 21 893 3010 Fax: +351 21 896 9099 http://www.deimos.pt __ Information from ESET NOD32 Antivirus, version of virus signature database 6582 (20111028) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Renaming a column's name in a Shapefile
Hi Ralf Thanks for the feedback but it did not work. I tried: ogr2ogr -f esri shapefile -sql SELECT ele, time AS t, magvar,geoidheigh, name, cmt, desc, src FROM GPS_Project new.shp GPS_Project.shp And I got: ERROR 1: SQL Expression Parsing Error: syntax error Even I tried: ogr2ogr -f esri shapefile -sql SELECT ele, magvar,geoidheigh, name, cmt, desc, src FROM GPS_Project new.shp GPS_Project.shp And I got: ERROR 1: SQL Expression Parsing Error: syntax error Any other suggestion? Thanks Antonio Ralf Suhr wrote: ogr2ogr -f esri shapefile -sql SELECT ele, time AS t, magvar, geoidheigh, name, cmt, desc, src FROM GPS_Project new.shp GPS_Project.shp Gr Ralf On Freitag 28 Oktober 2011 13:54:37 António Rocha wrote: ogr2ogr -sql SELECT [ele AS ele, time AS t, magvar AS magvar, geoidheigh AS geoidheigh, name AS name, cmt AS cmt, desc_ AS desc_, src AS src] FROM GPS_Project.shp, new.shp GPS_Project.shp __ Information from ESET NOD32 Antivirus, version of virus signature database 6582 (20111028) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Checking if a label exists in a HDF file
Hi Anton Thanks for the feedbvack. Unfortunely I'm not much familiar with bash Script. I'm running this in Python (without Gdalpy). Is there a way for Python? __ Information from ESET NOD32 Antivirus, version of virus signature database 6351 (20110804) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Checking if a label exists in a HDF file
Greetings I have a Python Script where I call gdalwarp binary (in Windows) to warp an image. But, before I run the gdalwarp I need to check if the file exists (DONE) and if a defined HDF layer exists. How can I check, using gdal, in a Script if a certain layer exists? THanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6348 (20110803) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] usage of subprocess.call
Greetings In a Python script I'm using this to run an external BIN: p=subprocess.call(command) where command is my command to run the bin. p will be equal tot he return values of the bin. (so far so good) But I'm getting an error now while running bin: File C:/EIS/scripts/mycopy, line 577, in module main() File C:/EIS/scripts/mycopy, line 278, in main p=subprocess.call(command) File C:\EIS\Python25\lib\subprocess.py, line 444, in call return Popen(*popenargs, **kwargs).wait() File C:\EIS\Python25\lib\subprocess.py, line 594, in __init__ errread, errwrite) File C:\EIS\Python25\lib\subprocess.py, line 816, in _execute_child startupinfo) WindowsError: [Error 193] My problem is not the error but, is there a way to get the error and stop the process in a more smoother way instead of stopping so drastically? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6286 (20110712) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Using gdalwarp with HDF files- problem reading proejction
Greetings I have a dataset (in HDF) and I need to convert its projection. When I do gdalinfo I get this information: Coordinate System is: GEOGCS[Unknown datum based upon the Clarke 1866 ellipsoid, DATUM[Not specified (based on Clarke 1866 spheroid), SPHEROID[Clarke 1866,6378206.4,294.9786982139006, AUTHORITY[EPSG,7008]]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433]] Origin = (-10.0013888,45.0013890) Pixel Size = (0.0027778,-0.0027778) But this is odd. I was suppose to get Plate-Carrée with geographic Lat/Long and Ref ellipsoid WGS 84 Spherical Radius= Equatorial Radius=Re= 6378,14km) Why am I no obtaining this information? Nevertheless, has this any impact or shall I use *-s_srs and define it do Plate Carée (EPSG:??) ? *Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6272 (20110707) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Using gdal to generate thumbnails
Greetings I have a couple of Gtif Images and I want to create thumbnails. Is it possible to do that with GDAL? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 6119 (20110513) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Re: [GRASS-dev] Error with r.to.vect- Solved?
Hello Markus Thanks. Just one final question: Is there a way to know which is the limit for vector? My idea is to check if raster, to be converted, is too big to be converted before I try to convert? THanks Markus Neteler wrote: Hi António (cc Martin) you can get winGRASS 7 here: http://wingrass.fsv.cvut.cz/grass70/ but I dunno if LFS is enabled in that built. Maybe Martin can tell us. cheers Markus 2011/3/29 António Rocha antonio.ro...@deimos.com.pt: Hi Markus The problem is that I'm running in Windows and in GRASS6.4.0stable release. So, there is no solution for Windows? Thanks Antonio Markus Neteler wrote: 2011/3/29 António Rocha antonio.ro...@deimos.com.pt: Hello Markus Unfortunely, I get this error when I run other datasets. In this case, a landCover with all CELLS and spatial resolution of 30. The problem is that North carolina sample is too small. r.to.vect -v input=FULL@PERMANENT output=FULL feature=area I get this: Default driver / database set to: driver: dbf database: $GISDBASE/$LOCATION_NAME/$MAPSET/dbf/ Extracting areas... Building topology for vector map PT_National_2000_FULL... R4243000ERROR: G_realloc: unable to allocate 16976004 bytes at struct_alloc.c:133 You mean that it is too big :) The file in question is lib/vector/diglib/struct_alloc.c Please try GRASS 7 instead of GRASS 6 which offers large file support for vector data (to some extent). You need to configure it with this options --enable-largefile enable support for large files (LFS) and compile it. Or use the weekly Linux binary snapshot which has LFS enabled. Markus __ Information from ESET NOD32 Antivirus, version of virus signature database 5995 (20110329) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5995 (20110329) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5995 (20110329) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com __ Information from ESET NOD32 Antivirus, version of virus signature database 5995 (20110329) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Re: Warping ECMWF's grib file
Driver: GRIB/GRIdded Binary (.grb) Files: c:\Test_areas\area8.grib Size is 15, 5 Coordinate System is: GEOGCS[Coordinate System imported from GRIB file, DATUM[unknown, SPHEROID[Sphere,6367470,0]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433]] Origin = (169.500,6.000) Pixel Size = (1.500,-1.500) Corner Coordinates: Upper Left ( 169.500, 6.000) (169d30' 0.00E, 6d 0' 0.00N) Lower Left ( 169.500, -1.500) (169d30' 0.00E, 1d30' 0.00S) Upper Right ( 192.000, 6.000) (192d 0' 0.00E, 6d 0' 0.00N) Lower Right ( 192.000, -1.500) (192d 0' 0.00E, 1d30' 0.00S) Center ( 180.750, 2.250) (180d45' 0.00E, 2d15' 0.00N) So it seems some kind of a conflict with the coordinates. Right? Any tip on how to solve it? Thanks António Rocha wrote: Greetings I have a GRIB file with longitude between 170 and 190 (attached to this email) and, when I try to warp it to Geotiff the outcome cannot be displayed. Can anyone give me a tip on how to warp this file? (the problem is not the grib but the longitude coordinates) Thankls Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5996 (20110329) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Error running gdal16 binary
Greetings While I was running GDAL16 binary, in Windows OS, for a grib file I got this error after a few lines: ERROR 2: CPLMalloc(): Out of memory allocating -469696517 bytes. This application has requested the Runtime to terminate it in an unusual way. Please contact the application's support team for more information. What might be happening here? because I have used with similar GRIBS (ECMWF ERA Interim dataset) and I never got this problem Thanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5898 (20110223) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Error running gdal16 binary
Hi Jurgen But in my case I need the binaries in order to integrate it in a package that I'm developing. So I need standalone Windows Binaries. Is there any of it in OSGEO4W? I mean can I pick them up and remove from osgeo4w folder? Thanks Antonio On Wed, 23. Feb 2011 at 09:30:09 +, António Rocha wrote: / Hello Chaitanya // Yesterday I was searching Windows binary gdal18 version and I didn't // find it. Where can I get it? / OSGeo4W[1] has GDAL 1.8. Jürgen [1] http://trac.osgeo.org/osgeo4w/ __ Information from ESET NOD32 Antivirus, version of virus signature database 5898 (20110223) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] using gdalwarp: defining target_spatial resolution
Greetings I'm using a bash script to reproject an image, using gdalwarp, to a set of Coordinates/Projection systems. I'm using -tr flag to define output spatial resolution. Since, i'm using UTM WGS84 I have this information in Meters (e.g. 30) so far, so good. The thing is that i need to reproject to this system +proj=longlat +no_defs +a=6378137 +rf=298.257223563 +towgs84=0.000,0.000,0.000 So, I'm getting an error because this system uses degrees not meters as units. and, I will probably get more problems in the near future with other systems My questions are: 1- from PROJ.4 format, how can I tell the units of the target coordinate/projection system 2- Since it's in degrees, does anyone suggest any (really good) tool to convert meters to degs? In order to parameterize -tr parameter? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5832 (20110130) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Strange printed lines for binary gdalinfo16 in Windows
Hi Thanks for the information. Where can I get a newer verasion of GDAL-BIn for Windows? THanks Chaitanya kumar CH wrote: António, This is not a bug. This is the debugging information printed by the grib driver. The code to print this has been commented out by Even a while ago. Try a later version of GDAL for less verbosity. __ Information from ESET NOD32 Antivirus, version of virus signature database 5772 (20110109) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Strange printed lines for binary gdalinfo16 in Windows
Greetings I have just downloaded from here http://download.osgeo.org/gdal/win32/1.6/gdalwin32exe160.zip, GDAL16 for Windows. When I do gdalinfo.exe for a GRIB file before it prints the header of the file it prints: ID 1 Class 14 Type 9 Stream 1025 Ver 0001, Octet-50 0, Octet-51 0 SectLen 52 repeatedly (nearly 1400). Is this a bug? Because I imagine that this is not an expected print... Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5766 (20110107) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Error: EPSG code not recognized
Greetings I'm using gdal15 and gdal16 and I'm trying to use gdalwarp with the following command: gdalwarp -s_srs EPSG:9823 -t_srs EPSG:3763 file.nc file.tif and I get the error: ERROR 6: EPSG PCS/GCS code 9823 not found in EPSG support files. Is this a valid EPSG coordinate system? ERROR 1: Translating source or target SRS failed: EPSG:9823 The thing is that this EPSG code exists and is this one: http://www.remotesensing.org/geotiff/proj_list/equirectangular.html Is this a bug? or is there other alternative? Thanks Antonio R. __ Information from ESET NOD32 Antivirus, version of virus signature database 5676 (20101205) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Error: EPSG code not recognized
Hi Jean-Claude So, how can I know which EPSG CRS codes are use this projection? Thanks Antonio Message: 2 Date: Mon, 06 Dec 2010 12:40:28 +0100 From: Jean-Claude Repetto jrepe...@free.fr Subject: Re: [gdal-dev] Error: EPSG code not recognized To: gdal-dev@lists.osgeo.org Message-ID: 4cfccbac.5050...@free.fr Content-Type: text/plain; charset=ISO-8859-1; format=flowed Le 06/12/2010 11:21, António Rocha a écrit : Greetings I'm using gdal15 and gdal16 and I'm trying to use gdalwarp with the following command: gdalwarp -s_srs EPSG:9823 -t_srs EPSG:3763 file.nc file.tif and I get the error: ERROR 6: EPSG PCS/GCS code 9823 not found in EPSG support files. Is this a valid EPSG coordinate system? ERROR 1: Translating source or target SRS failed: EPSG:9823 The thing is that this EPSG code exists and is this one: http://www.remotesensing.org/geotiff/proj_list/equirectangular.html Is this a bug? or is there other alternative? Thanks Antonio R. Hi Antonio, EPSG:9823 is not a coordinate reference system, it is only a projection code. Jean-Claude -- __ Information from ESET NOD32 Antivirus, version of virus signature database 5677 (20101206) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Running gdalwarp over ECMWF's NetCDF
Greetings I have runned the same gdalwarp command and I got two different results Command: gdalwarp -t_srs EPSG:3763 2009.nc 2009.tif In Ubuntu (with GDAL 1.7.1, released 2010/02/08) It procuded the raster file In Windows (with OSGEO4w and GDAL1.5.4) I got the eerror: Creating output file that is 3Px4L. Processing input file 2009.nc A target coordinate system was specified, but there is no source coordinate system. Consider using -s_srs option to provide a source coordinate system. Operation terminated Is this a limitation of GDAL1.5.4? Thanks __ Information from ESET NOD32 Antivirus, version of virus signature database 5667 (20101202) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] ECMWF's GRIBs and NetCDFs
Hi I'm using ECMWF's GRIB files and, I realized that there is a small error in geolocation (that can be fixed with gdal_translate). But, since I changed to winGRASS (that uses gdal1.5) I'm not able to use ECMWF's GRIB. Since ECMWF's data is also provided in NetCDF, I would like to know if there is any known issue on using NetCDF data in GDAL1.5? (I didn't find anything). Thanks Antonio R. __ Information from ESET NOD32 Antivirus, version of virus signature database 5661 (20101130) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Create a Google Earth overlay based on Geotiff
Greetings I have a goerreferenced Geotiff image, already in EPSG:4326, and I need to automatically create an KML with an overlay from this geotiff. Is it possible using GDAL? If so, how can I do that? Or can anyone indicate me an Wiki that explain this? I have found this link: http://code.google.com/apis/kml/articles/raster.html but it creates empty doc.kml Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5428 (20100906) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Create a Google Earth overlay based on Geotiff
Hi there It worked (my problem was that I was using gdal_translate instead of gdalwarp to reproject the images to EPSG:4326). But now I have 1 little question: When I do the reprojection of the image from my UTM to EPSG:4326 It creates a black-NO VALUE area around the image due to projection. In order to auxiliate the visualization is it possible to turn this no-data into transparent? Chaitanya kumar CH wrote: António, Refer to http://www.gdal.org/gdal2tiles.html and http://www.gdal.org/gdal_utilities.html 2010/9/6 António Rocha antonio.ro...@deimos.com.pt mailto:antonio.ro...@deimos.com.pt Greetings I have a goerreferenced Geotiff image, already in EPSG:4326, and I need to automatically create an KML with an overlay from this geotiff. Is it possible using GDAL? If so, how can I do that? Or can anyone indicate me an Wiki that explain this? I have found this link: http://code.google.com/apis/kml/articles/raster.html but it creates empty doc.kml Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5428 (20100906) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org mailto:gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev -- Best regards, Chaitanya kumar CH. /tʃaɪθənjə/ /kʊmɑr/ +91-9494447584 17.2416N 80.1426E __ Information from ESET NOD32 Antivirus, version of virus signature database 5428 (20100906) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Retrieve number of layers with gdalpy
Greetings I have a file Geotiff format (and other in GRIB) and I need to know how many layers of data I have in each one of them by using gdalPy. is it possible? (I don't want to use gdalinfo , I just want to get this value directly) Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5356 (20100810) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Using gdalwarp over NCEP - keep pixel resolution
Greetings I'm doing some SRS conversion of GRIB files using gdalwarp (to WGS-UTM) but depending on the zone I'm obtaining different Pixel Size. I know I can use the -tr parameter but in this case I know the source_resolution (1,5 in NCEP) but I don't know to what it corresponds in WGSUTM. Has anyone has a tip on how to figure this correspondence in order to use that value in the Gdalwarp expression? Thanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5228 (20100625) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Gdalwarp Landsat image- Segmentation fault
Hi I'm trying to Warp a Landsat Geotiff band (in UTM coordinates) to a specific set of coordinates for Portugal, using this command: gdalwarp -t_srs EPSG:3763 L71203033_0332601_B10.TIF B10.tif log.log Not only the image size doesn't match, but pixel size was changed (a bit) and I got a segmentation error and TIFFWriteEncodedStrip and TIFFReadDirectory Errors. What am I doing wrong? This is the error log that I obtained Creating output file that is 8229P x 7218L. Processing input file L71203033_0332601_B10.TIF. 0...10...20...30...40...50ERROR 1: TIFFFetchDirectory:B10.tif: Can not read TIFF directory count ERROR 1: TIFFReadDirectory:Failed to read directory at offset 8 ERROR 1: TIFFWriteEncodedStrip:Must set ImageWidth before writing data ERROR 1: TIFFWriteEncodedStrip:Must set ImageWidth before writing data ERROR 1: TIFFWriteEncodedStrip:Must set ImageWidth before writing data ERROR 1: TIFFWriteEncodedStrip:Must set ImageWidth before writing data ERROR 1: TIFFWriteEncodedStrip:Must set ImageWidth before writing data (..) More than 1000 errors or warnings have been reported. No more will be reported from now. ...60...70...80...90... Segmentation fault __ Information from ESET NOD32 Antivirus, version of virus signature database 5143 (20100525) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Gdalwarp Landsat image- Segmentation fault
Greetings Franz This may sound ackward but I was running it at a shared folder (between Windows and my Ubuntu)... But If I run at my ubuntu's user area Desktop it runs perfectly well... Thanks for your help. This is really strange. Sorry for bothering you with a stupid thing like this Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5144 (20100525) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Producing GML files from a Shapefile
Greetings I have a standard ESRI Shapefile and I want to produce a GML file. Is it possible using OGR? (I know that it's possilbe to produce KML) Thanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5134 (20100521) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Obtain more information with Ogrinfo
Greetings I need to obtain a little bit more information regarding a Shapefile. I have done ogrinfo filename but I only obtained: INFO: Open of `file.shp' using driver `ESRI Shapefile' successful. 1: training_data_AA_NAC_N_M (Polygon) I mean it only identified me the sub-dataset and the data-format. I need information regarding projection/coordinates. How can I obtain it? Thanks Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 5100 (20100510) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Converting GEOS to Gtiff lat/lon
Greetings Bas has been helping me in TRANSLATE and WARPING MSG data, in HDF5 format, to a georrectified format that I can use in GRASS. First I do this: gdal_translate -of GTiff -a_srs +proj=geos +lon_0=0.0 +h=35785831 +x_0=0.0 +y_0=0.0 -a_ullr -924124.2 5421728.5 4179561.6 3468466.1 HDF5:HDF5_MSG_200802170100://ABC area_georeferenced.tif Everything runs smoothly. Then I have to do: gdalwarp -s_srs +proj=geos +lon_0=0.0 +h=35785831 +x_0=0.0 +y_0=0.0 -t_srs +proj=latlong area_georeferenced.tif euro_area_latlong.tif For some datasets I get an image with different size (instead of 1701, 651 is 1721, 595) and in another cases I get the following error: Creating output file that is 1721P x 595L. Processing input file euro_area_georeferenced.tif. 0...10...20...30...40...50...60...70...80...90...ERROR 1: TIFFReadDirectory:euro_area_latlong.tif: Can not read TIFF directory count ERROR 1: TIFFWriteEncodedStrip:euro_area_latlong.tif: Must set ImageWidth before writing data ERROR 1: TIFFWriteEncodedStrip:euro_area_latlong.tif: Must set ImageWidth before writing data ERROR 1: TIFFWriteEncodedStrip:euro_area_latlong.tif: Must set ImageWidth before writing data (...) Why am I obtaining a smaller dataset and, why I get this error in some datasets? I even tried with different fileformats and the error is the same. Bas is using GDAL1.6dev in FWTOOLS and I'm using 1.5.4 in Ubuntu OS (to work with GRASS) Thank you Best regards Antonio Rocha __ Information from ESET NOD32 Antivirus, version of virus signature database 4890 (20100223) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Supported Remote Sensing instruments/satellites
Greetings all I'm a GRASS user and, following a few emails from Franz Schiller in GRASS users mailing list, I decided to send an email to GDAL mailing list regarding the following question: Which native Remote Sensing instruments/satellites data formats are supported by GDAL? I mean, some instruments/satellites are distributed in standard formats (e.g. N1, Geotiff) but others not really. So my question is if you have a table, just like this (http://www.gdal.org/formats_list.html), but instead of being organized by formats it's by Satellites/nstruments? Thanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 4812 (20100128) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Supported Remote Sensing instruments/satellites
Greetings all I'm a GRASS user and, following a few emails from Franz Schiller in GRASS users mailing list, I decided to send an email to GDAL mailing list regarding the following question: Which native Remote Sensing instruments/satellites data formats are supported by GDAL? I mean, some instruments/satellites are distributed in standard formats (e.g. N1, Geotiff) but others not really. So my question is if you have a table, just like this (http://www.gdal.org/formats_list.html), but instead of being organized by formats it's by Satellites/nstruments? Thanks Best regards, Antonio __ Information from ESET NOD32 Antivirus, version of virus signature database 4790 (20100120) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev