[gdal-dev] Generating a Google Earth image overlay from a Geotiff image
Greetings I have a raster RGB image (Geotiff) and I want to create an image overlay KML file to be seen on Google Earth. Is it possible to do this in GDAL? If yes, how ? Thanks Best regards, Luisa ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Re: How to use gdalinfo for HDF- Windows
Hey Brent I'm using ESA's GLOBCOVER data.In this case I'm trying to set a GRASS script to run gdal directly. I do: gdalinfo C:\data\GL.hdf and I retrieve all basic information plus HDF's layers Driver: HDF4/Hierarchical Data Format Release 4 Files: c:\data\GL.hdf Size is 512, 512 Coordinate System is `' Metadata: HDFEOSVersion=HDFEOS_V2.12 Subdatasets: SUBDATASET_1_NAME=HDF4_EOS:EOS_GRID:c:\data\GL.hdf:POSTEL:MEAN_1 SUBDATASET_1_DESC=[1800x1800] MEAN_1 POSTEL (16-bit integer) SUBDATASET_2_NAME=HDF4_EOS:EOS_GRID:c:\data\GL.hdf:POSTEL:MEAN_2 SUBDATASET_2_DESC=[1800x1800] MEAN_2 POSTEL (16-bit integer) SUBDATASET_3_NAME=HDF4_EOS:EOS_GRID:c:\data\GL.hdf:POSTEL:MEAN_3 SUBDATASET_3_DESC=[1800x1800] MEAN_3 POSTEL (16-bit integer) SUBDATASET_4_NAME=HDF4_EOS:EOS_GRID:c:\data\GL.hdf:POSTEL:MEAN_4 SUBDATASET_4_DESC=[1800x1800] MEAN_4 POSTEL (16-bit integer) SUBDATASET_5_NAME=HDF4_EOS:EOS_GRID:c:\data\GL.hdf:POSTEL:MEAN_5 (...) then: gdalinfo HDF4_EOS:EOS_GRID:c:\data\GL.hdf:POSTEL:MEAN_1 I get: ERROR 4: 'HDF4_EOS:EOS_GRID:c:\data\GL.hdf:POSTEL:MEAN_1' does not exist in fyle system, and is not recognised as a supported dataset name. gdalinfo failed- unable to open 'HDF4_EOS:EOS_GRID:c:\data\GL.hdf:POSTEL:MEAN_1'' open GDAL Datasets: 1 N DriverIsNULL 512x512x0 But if I do: gdalinfo HDF4_EOS:EOS_GRID::\data\GL.hdf:POSTEL:MEAN_1 it works (without the C from the driver) The problem is that if I'm trying to to use data in D it does not work. Any suggestions? 2011/8/9 Brent Fraser bfra...@geoanalytic.com Luisa, I use GDAL to access MODIS data in HDF4 on Windows like this: gdalinfo HDF4_EOS:EOS_SWATH:G:\Projects\temp\image3\MOD02HKM.A2011186.1640.005.NRT.hdf:MODIS_SWATH_Type_L1B:EV_500_RefSB So I'd suggest using double quotes around your path. Best Regards, Brent Fraser On 8/9/2011 8:24 AM, Luisa Peña wrote: Greetings This seems odd but it does not work with HDF4_EOS:EOS_GRID:C:/filepath only HDF4_EOS:EOS_GRID:filepath And if my MERIS data is not placed in my main hard disk (C:) it does not work. Anyone with some experience on this in WINDOWS? Thanks Luisa ___ gdal-dev mailing listgdal-dev@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Re: How to use gdalinfo for HDF- Windows
Greetings This seems odd but it does not work with HDF4_EOS:EOS_GRID:C:/filepath only HDF4_EOS:EOS_GRID:filepath And if my MERIS data is not placed in my main hard disk (C:) it does not work. Anyone with some experience on this in WINDOWS? Thanks Luisa ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] How to use gdalinfo for HDF
Greetings I have read This webpage: http://www.gdal.org/frmt_hdf4.html But I'm having one difficulty If I have my HDF in my current folder I'm able to do gdalinfo over my file and then over my layer gdalinfo HDF4_EOS:EOS_GRIDfilename.hdf:POSTEL:data1 and I get all the information But, if I have my data in other folder (like C:\data) I do: gdalinfo C:\data\filename.hdf and I get a list of layers and then: gdalinfo HDF4_EOS:EOS_GRIDC:\filename.hdf:POSTEL:data1 I get: Open GDAL datasets: 1 N DriverIsNull 512x512x0 How can I apply gdalinfo (or gdalwarp) for cases like this? THanks ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Reproject Landsat USGS image
Hello I got a reply to my message stating that GDAL ignores a flag in the GEOTIF specification that says whether the coordinates are for the center or the lower-left of each pixel. Is this true? It's exacly what I am obtaining and what Markus Neteler explained me (about the useage of NN method) Can anyone confirm me? Thanks Luisa 2011/3/25 Luisa Peña luisapena1...@gmail.com Back again on this I have talked with Markus Neteler and he told me that if I use gdalwarp with NN method, I will get a pixel shift and he suggested me to use bilinear method. The problem with bilinear is that my pixel values are changed (based on bilinear method) So is there any method to keep the pixels values and without any pixel position shift? Thanks Best regards, Luisa 2011/3/23 Luisa Peña luisapena1...@gmail.com For now GDAL15 (the one that comes with winGRASS GIS) and GDAL16 binary. Why? (I have sent this to the wrong mailing list) Why? 2011/3/22 Kyle Shannon ksshan...@gmail.com What version of gdal are you using? # Kyle Shannon Physical Science Technician RMRS Fire Sciences Lab Fire, Fuels Smoke - RWU 4405 5775 Highway 10 W. Missoula, MT 59808 (406)829-6954 kshan...@fs.fed.us # On Tue, Mar 22, 2011 at 15:36, Luisa Peña luisapena1...@gmail.comwrote: Greetings I have sent an message a couple of days ago ( http://lists.osgeo.org/pipermail/gdal-dev/2011-March/028124.html) about this I used gdalwarp with this formulation gdalpwarp -t_srs EPSG -tr 30 30 input.tif output.tif Well it works fine but It seems to be with a deviation of half pixel on Latitude and half pixel on Longitude when compared with a similar operation in ARCGis Maybe I'm missing some parameter that might be relevant. Any suggestion? Thanks Luisa ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Reproject Landsat USGS image
Back again on this I have talked with Markus Neteler and he told me that if I use gdalwarp with NN method, I will get a pixel shift and he suggested me to use bilinear method. The problem with bilinear is that my pixel values are changed (based on bilinear method) So is there any method to keep the pixels values and without any pixel position shift? Thanks Best regards, Luisa 2011/3/23 Luisa Peña luisapena1...@gmail.com For now GDAL15 (the one that comes with winGRASS GIS) and GDAL16 binary. Why? (I have sent this to the wrong mailing list) Why? 2011/3/22 Kyle Shannon ksshan...@gmail.com What version of gdal are you using? # Kyle Shannon Physical Science Technician RMRS Fire Sciences Lab Fire, Fuels Smoke - RWU 4405 5775 Highway 10 W. Missoula, MT 59808 (406)829-6954 kshan...@fs.fed.us # On Tue, Mar 22, 2011 at 15:36, Luisa Peña luisapena1...@gmail.comwrote: Greetings I have sent an message a couple of days ago ( http://lists.osgeo.org/pipermail/gdal-dev/2011-March/028124.html) about this I used gdalwarp with this formulation gdalpwarp -t_srs EPSG -tr 30 30 input.tif output.tif Well it works fine but It seems to be with a deviation of half pixel on Latitude and half pixel on Longitude when compared with a similar operation in ARCGis Maybe I'm missing some parameter that might be relevant. Any suggestion? Thanks Luisa ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Reproject Landsat USGS image
For now GDAL15 (the one that comes with winGRASS GIS) and GDAL16 binary. Why? (I have sent this to the wrong mailing list) Why? 2011/3/22 Kyle Shannon ksshan...@gmail.com What version of gdal are you using? # Kyle Shannon Physical Science Technician RMRS Fire Sciences Lab Fire, Fuels Smoke - RWU 4405 5775 Highway 10 W. Missoula, MT 59808 (406)829-6954 kshan...@fs.fed.us # On Tue, Mar 22, 2011 at 15:36, Luisa Peña luisapena1...@gmail.comwrote: Greetings I have sent an message a couple of days ago ( http://lists.osgeo.org/pipermail/gdal-dev/2011-March/028124.html) about this I used gdalwarp with this formulation gdalpwarp -t_srs EPSG -tr 30 30 input.tif output.tif Well it works fine but It seems to be with a deviation of half pixel on Latitude and half pixel on Longitude when compared with a similar operation in ARCGis Maybe I'm missing some parameter that might be relevant. Any suggestion? Thanks Luisa ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Reproject Landsat USGS image
Greetings I have sent an message a couple of days ago ( http://lists.osgeo.org/pipermail/gdal-dev/2011-March/028124.html) about this I used gdalwarp with this formulation gdalpwarp -t_srs EPSG -tr 30 30 input.tif output.tif Well it works fine but It seems to be with a deviation of half pixel on Latitude and half pixel on Longitude when compared with a similar operation in ARCGis Maybe I'm missing some parameter that might be relevant. Any suggestion? Thanks Luisa ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Reproject USGS Landsat image
Hi, I'm planning to reproject a Landsat image using gdalwarp but I have a question: - Besiodes the 7 bands, I also have an AUX file (for each band), an RRD for each file and a GCP. Is it possibvle to use gdalwarp to read this GCP file or is it useless since all the required information is already available in each band? Thanks Best regards, Luisa Pena ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Re: Updating GDAL on OSGEO4w (windows)
Hello Jurgen But you're running gdalinfo from the same shell, that you ran gdal16/gdal17/gdaldev and your current directory is _not_ %OSGEO4W_ROOT%\bin? Uhm Maybe I wasn't clear. I run gdal16 in C:/osgeo4w/bin and then I did gdalinfo --version I attached a screenshot that exemplifies it Thanks Luisa attachment: gdal16_01.PNG___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Processing GRIB data to ASCII or text file
Hi there I have several GRIB files and I wanted to do some statistical processing. To do so, I wanted to convert them to ASCII or text file. Is it possible using GDAL? Thanks Best regards, Luisa ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Gdalwarp of GLOBCOVER to UTMSpecific Zone
Hi there I'm a bit of a newbie in GDAL but I tried to warp a GLOBCOVER image with this information: Coordinate System is: GEOGCS[WGS 84, DATUM[WGS_1984, SPHEROID[WGS 84,6378137,298.257223563, AUTHORITY[EPSG,7030]], AUTHORITY[EPSG,6326]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433], AUTHORITY[EPSG,4326]] Into an UTM zone (EPSG:21036) I runned this command: gdalwarp -t_srs EPSG:21036 INPUT.tif OUTPUT.tif And I get this Warning/ERROR: ERROR 1: TIFFFetchDirectory:MZ.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 Segmentation fault And I obtained an OUTPUT.tif with a different size. - What might not be ok in this? - What can I do to obtain an image with the same size? Thanks best regards Luisa ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Spatial Reference syntax
Hello I'm reading a few examples of GDALWARP usage and I have a few questions regarding spatial reference syntax. Is there any Wiki explaining explaining coordinates systems definition to be passed to GDALWARP/other GDAL tool? Thanks ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Warp a file from two separate Lat and Long
Greetings I've used GDAL a few times but this is the first time that I'm lost on what I can do with GDAL to solve my problem. I have 3 separate files, one with my data, other with Lattitude and another with Longitude. They have the same size and, as an example, it means that Pixel (1,1) of data is located in latitude retrieved from Pixel (1,1) from latitude and longitude (pixel (1,1) from longitude). All of them are in WGS84 but without GCP's points. I need to create a georreferenced/rectified Geotiff with this data and, to do that, I need lat and Long. How can I do this? Thanks Best regards, Luisa Peña S.Paulo, Brazil ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Warp a file from two separate Lat and Long
Hi Andrew But How can I create those VRT, I never used those, from my original binaries datasets (HDF)? It seems that VRT is the key for my problem. But I never used it before Thanks for your help Luisa 2010/2/18 Andrew Clegg a...@pml.ac.uk Hi Luisa, I've tried something similar in the past with satellite data. I had 3 files, containing latitude, longitude and value stored as floats. I had some very minor success with the following solution (I got out an image but it seemed to drop to a very low resolution). I'm posting it more in the hope that someone else may be able to spot what I did wrong, because I would still really like to be able to do this. My approach was to create a virtual dataset for each file: sst.vrt (the values): VRTDataset rasterXSize=1296 rasterYSize=2030 Metadata domain=GEOLOCATION MDI key=X_DATASETlon.vrt/MDI MDI key=X_BAND1/MDI MDI key=Y_DATASETlat.vrt/MDI MDI key=Y_BAND1/MDI MDI key=PIXEL_OFFSET0/MDI MDI key=LINE_OFFSET0/MDI MDI key=PIXEL_STEP1/MDI MDI key=LINE_STEP1/MDI /Metadata VRTRasterBand dataType=Float32 band=1 subClass=VRTRawRasterBand SourceFilename relativetoVRT=1sst.img/SourceFilename ImageOffset0/ImageOffset PixelOffset4/PixelOffset LineOffset5184/LineOffset ByteOrderLSB/ByteOrder /VRTRasterBand /VRTDataset lon.vrt (longitude): (latitude is basically the same) VRTDataset rasterXSize=1296 rasterYSize=2030 VRTRasterBand dataType=Float32 band=1 subClass=VRTRawRasterBand SourceFilename relativetoVRT=1lon.img/SourceFilename ImageOffset0/ImageOffset PixelOffset4/PixelOffset LineOffset5184/LineOffset ByteOrderLSB/ByteOrder /VRTRasterBand /VRTDataset I then tried to run gdalwarp on sst.vrt, trying either '-geoloc' or '-to METHOD=GEOLOC_ARRAY'. Regards, Andrew Clegg Plymouth Marine Laboratory Luisa Peña wrote: Greetings I've used GDAL a few times but this is the first time that I'm lost on what I can do with GDAL to solve my problem. I have 3 separate files, one with my data, other with Lattitude and another with Longitude. They have the same size and, as an example, it means that Pixel (1,1) of data is located in latitude retrieved from Pixel (1,1) from latitude and longitude (pixel (1,1) from longitude). All of them are in WGS84 but without GCP's points. I need to create a georreferenced/rectified Geotiff with this data and, to do that, I need lat and Long. How can I do this? Thanks Best regards, Luisa Peña S.Paulo, Brazil ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev Plymouth Marine Laboratory Registered Office: Prospect Place The Hoe Plymouth PL1 3DH Website: www.pml.ac.uk Registered Charity No. 1091222 PML is a company limited by guarantee registered in England Wales company number 4178503 PML is a member of the Plymouth Marine Sciences Partnership Website: www.pmsp.org.uk This e-mail, its content and any file attachments are confidential. If you have received this e-mail in error please do not copy, disclose it to any third party or use the contents or attachments in any way. Please notify the sender by replying to this e-mail or e-mail fori...@pml.ac.ukand then delete the email without making any copies or using it in any other way. The content of this message may contain personal views which are not the views of Plymouth Marine Laboratory unless specifically stated. You are reminded that e-mail communications are not secure and may contain viruses. Plymouth Marine Laboratory accepts no liability for any loss or damage which may be caused by viruses. ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev