Re: [gdal-dev] help resampling raster to coarser scale

2011-08-05 Thread Chaitanya kumar CH
Etienne, You can find a lot of sample code in the gdal test scripts[1]. There is a test script for regenerating overviews in gcore/tiff_ovr.py [2]. You can apply a filter using a VRT file[3]. It supports filters of various sizes. [1]: http://trac.osgeo.org/gdal/browser/trunk/autotest/ [2]: http:

Re: [gdal-dev] help resampling raster to coarser scale

2011-08-05 Thread Etienne
Thanks Anton This might seem obvious to many of you, but how could I go about applying a filter using the gdal utilities/functions/python code? Etienne - Original Message - From: Anton Korosov To: gdal-dev@lists.osgeo.org Cc: Sent: Friday, August 5, 2011 3:27:41 AM Subject: Re: [gdal

Re: [gdal-dev] help resampling raster to coarser scale

2011-08-05 Thread Etienne
That's an interesting idea, thanks! I looked around the code and found that the overviews are created in the GDALRegenerateOverviews() function in file gcore/overview.cpp .  It reads: [...] it can also be used to generate downsampled images in one file from another outside the overview architec

Re: [gdal-dev] problem rasterizing categorical data from shapefile - 0 values at some intersections

2011-08-05 Thread Etienne
Chaitanya, I forgot to mention that the entire shapefile (much bigger) has many more polygons with 32 possible values (land cover type), so -with isn't practical in this case. What I meant by  'most representative value' was, to have the value which represents the dominant land cover type in

Re: [gdal-dev] problem rasterizing categorical data from shapefile - 0 values at some intersections

2011-08-05 Thread Chaitanya kumar CH
I suggested the -where option because your shapefile had only two features that are (almost)inverse of each other. You can select only one of them. You might want to use the -at option with this. What do you mean by 'most representative value'? On Fri, Aug 5, 2011 at 5:40 PM, Etienne wrote: > H

Re: [gdal-dev] How to use gdalinfo for HDF

2011-08-05 Thread Anton Korosov
Hi! Looks like a syntax error. There should be a colon after EOS_GRID Try gdalinfo HDF4_EOS:EOS_GRID:"C:\filename.hdf":POSTEL:data1 At least that works at my Linux: gdalinfo HDF4_EOS:EOS_SWATH:"/host/gdaltest/MOD021KM.A2010105.2120.005.2010106075131.hdf":MODIS_SWATH_Type_L1B:EV_1KM_RefSB On

Re: [gdal-dev] problem rasterizing categorical data from shapefile - 0 values at some intersections

2011-08-05 Thread Etienne
Hi and thanks! I set the init value to 255, and the overlapping pixels are set to that value (instead of the default of 0).  Unfortunately, using -where resulted in the same gaps not being filled, and a combination of -where and -at selects too many pixels (the borders all overlap). Looking

[gdal-dev] How to use gdalinfo for HDF

2011-08-05 Thread Luisa Peña
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_GRID"filename.hdf":POSTEL:data1 and I get all the information But, if I have

Re: [gdal-dev] Checking if a label exists in a HDF file

2011-08-05 Thread Anton Korosov
Hi Antonio! In python you'll probably have to open the dataset, loop through all bands, fetch description (or metadata) from each band and then decide if the required band exists. f = gdal.Open('HDF4_EOS:EOS_SWATH:"MOD021KM.A2010105.2120.005.2010106075131.hdf":MODIS_SWATH_Type_L1B:EV_1KM_Ref

Re: [gdal-dev] help resampling raster to coarser scale

2011-08-05 Thread Vincent Schut
On 08/05/2011 03:57 AM, Etienne wrote: Hi all, I am trying to resample a raster with categorical data (land surface vegetation/cover) from a high resolution (500m) to a coarser resolution. I would like to do something similar to what the GRASS operator r.resamp.stats does, using the modal c

Re: [gdal-dev] Checking if a label exists in a HDF file

2011-08-05 Thread António Rocha
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 che

Re: [gdal-dev] problem rasterizing categorical data from shapefile - 0 values at some intersections

2011-08-05 Thread Chaitanya kumar CH
Etienne, Rasters of type Byte cannot hold values below 0. So, the nodata value of -1 is not practical. It's better to choose 0 or 255 depending on your data or choose a different data type. Your problem arises because the two geometries in the shapefile overlap in some places and leave some gaps