[GRASS-user] Fwd: PhD vacancy ‘Bayesian geostatistics for mapping soil and agronomic variables’
Forwarded Message From: Heuvelink, Gerard Subject: PhD vacancy ‘Bayesian geostatistics for mapping soil and agronomic variables’ Date: Wed, 8 Jan 2014 15:04:09 + Dear colleagues, We have a PhD vacancy on ‘Bayesian geostatistics for mapping soil and agronomic variables’. Perhaps you know suitable candidates, in which case we hope that you will forward them this email. Please also feel free to forward this message to anyone who might know of suitable candidates. The application deadline is January 20. More information and application instructions at https://www.academictransfer.com/employer/WUR/vacancy/21307/lang/en/. Best wishes, Gerard Gerard B.M. Heuvelink Wageningen University and ISRIC World Soil Information P.O. Box 47 6700 AA Wageningen The Netherlands tel +31 317 483723 Mon Tue Wed tel +31 317 482716 Wed Thu Fri emailgerard.heuvel...@wur.nl http://wageningenur.nl/sgl http://www.isric.org/staff SSbanner2014 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.mapcalc(nulls into 0)
Hi, I am trying to do something really simple but not finding the right way to do it. I need to convert nulls into zeros. I have done this: *r.mapcalc raster2 = if(raster1 = null(), 0, raster1) * *r.mapcalc raster2 = if(raster1 == null(), 0, raster1) ** * But this convert all cells into nulls. I have tried another approach with r.reclass, rules: null/null()/no data/nodata...=0 1=1 2=2 3=2 4=2 5=2 No fun, all I get is an ERROR: null/null()/no data/nodata...=0 What I'm doing wrong? Eloi Ribeiro GIS Analyst 39,45º -0,40º flavors.me/eloiribeiro ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.mapcalc(nulls into 0)
Thank you all! Solved. Eloi Ribeiro GIS Analyst 39,45º -0,40º flavors.me/eloiribeiro 2012/7/11 Alberto Pettazzi alberto.petta...@meteogalicia.es Hi, I think you should do raster2=if(isnull(raster1),0,raster1) El 11/07/12 14:02, Eloi Ribeiro escribió: Hi, I am trying to do something really simple but not finding the right way to do it. I need to convert nulls into zeros. I have done this: *r.mapcalc raster2 = if(raster1 = null(), 0, raster1) * *r.mapcalc raster2 = if(raster1 == null(), 0, raster1) ** * But this convert all cells into nulls. I have tried another approach with r.reclass, rules: null/null()/no data/nodata...=0 1=1 2=2 3=2 4=2 5=2 No fun, all I get is an ERROR: null/null()/no data/nodata...=0 What I'm doing wrong? Eloi Ribeiro GIS Analyst 39,45º -0,40º flavors.me/eloiribeiro ___ grass-user mailing listgrass-user@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user -- Alberto Pettazzi MeteoGalicia - Departamento de Climatología y Observación Consellería de Medio Ambiente, Territorio e Infraestruturas Rúa de Roma, 6 15707 Santiago de Compostela. A Coruña Teléfono: +34-881-999646 e-mail: alberto.petta...@meteogalicia.es ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Load PostGIS layer with tow geometry columns
Dear all, Trying to load a PostGIS layer with tow geometry columns it output's the following error: 'ERROR: Layer layer_name not available'. If one of the tow geometry's is dropped it works fine. Is there a way to define witch geometry I wont to load, something like 'geom_column=the_geom'? This is what I'm trying to do: *v.in.ogr dsn=PG:host=localhost dbname=bd_ceamsig user=postgres layer=sch_ambiente.escv_incendios_todos where=anyo=2009 output=temp* Thanks in advance! Regards, Eloi Ribeiro GIS Analyst 39,45º -4,40º http://eloiribeiro.wordpress.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Load PostGIS layer with tow geometry columns
On Tue, Apr 5, 2011 at 13:17, Martin Landa landa.mar...@gmail.com wrote: Hi, 2011/4/5 Eloi Ribeiro eloi.ribe...@gmail.com: This is what I'm trying to do: v.in.ogr dsn=PG:host=localhost dbname=bd_ceamsig user=postgres layer=sch_ambiente.escv_incendios_todos where=anyo=2009 output=temp missing `format=PostgreSQL` ? Martin Thanks for your quick reply Martin, Actually it's not a format parameter that I miss but a parameter that would define which geometry I want to use in case there is more than one geometry per layer/postgis table. Example. In my table '*escv_incendios_todos'* I have tow filds/columns with geometry, one is called geom_23030 (coord. syst. epsg 23030) and the other geom_4258 (coord. syst. epsg 4258) and a trigger maintaining this tow geometry's synchronized. In the geometry_columns table there are tow entry's for the '*escv_incendios_todos' *table: 3970433;'';sch_ambiente;escv_incendios_todos;geom_23030;2;23030;MULTIPOLYGON 3970434;'';sch_ambiente;escv_incendios_todos;geom_4258;2;4258;MULTIPOLYGON So when I try to load '*escv_incendios_todos'* in GRASS there is no way to define which geometry I want to use. And the error 'ERROR: Layer layer_name not available' jumps out. By now I worked around by deleting the entry with the geometry I'm not interested on the geometry_columns table and then insert it again. Thanks. Regards, ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] OBDC connection with postgresql in ubuntu 10.04 64 bits
Maybe this help you, although is in portuguese. http://eloiribeiro.wordpress.com/2009/01/27/criar-um-conector-odbc-para-postgresql-em-ubuntu/ By, Eloi Ribeiro GIS Analyst 39,45º -4,40º http://eloiribeiro.wordpress.com On Fri, Oct 1, 2010 at 20:13, Ricardo Rodríguez ricardorodo...@gmail.comwrote: hi all, how I can create ODBC connections with postgresql in Ubuntu 10.04 of 64 bits, since the ODBC plugin for postgresql is very different from what is documented to someone already done, I could explain it. thank you all for your time and help Ricardo Rodríguez Univalle ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Convert FCELL to CELL
Hi, How can I convert a raster type FCELL into a raster type CELL? Thanks, -- Eloi Ribeiro GIS Analyst 39,45º -4,40º http://eloiribeiro.wordpress.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Convert FCELL to CELL
2010/6/28 Micha Silver mi...@arava.co.il On 28/06/2010 14:40, Eloi Ribeiro wrote: Hi, How can I convert a raster type FCELL into a raster type CELL? Wouldn't: r.mapcalc new_map=int(original_map) do that? Thanks, -- Eloi Ribeiro GIS Analyst 39,45º -4,40º http://eloiribeiro.wordpress.com This mail was received via Mail-SeCure System. ___ grass-user mailing listgrass-u...@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silverhttp://www.surfaces.co.il/ Arava Development Co. +972-52-3665918 That's it! Thanks Micha! -- Eloi Ribeiro GIS Analyst 39,45º -4,40º http://eloiribeiro.wordpress.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Retrieve pixel resolution using Python Script
Olá António Rocha, This works for me, give it a try: ### import os try: from osgeo import gdal from osgeo import osr except: import gdal import osr def gdalinfo(pasta): '''Le metadados com GDAL''' ficheiros = os.listdir(pasta) for ficheiro in ficheiros: if ficheiro[-3:].lower() == 'tif': # abrir imagem filename = pasta+ficheiro dataset = gdal.Open(filename, gdal.GA_ReadOnly) if dataset is None: print 'Nao foi possivel abrir a imagem %s!' % ficheiro RasterXSize = dataset.RasterXSize RasterYSize = dataset.RasterYSize print RasterXSize,RasterYSize print '\nfim' return gdalinfo('/home/eloi/Geo/raster/20790/') ### 2010/6/28 Mark Seibel msei...@gmail.com If you use gdalinfo on an image, it will output the Pixel Size. Hope that helps. Mark 2010/6/28 António Rocha antonio.ro...@deimos.com.pt: Greetings I'm trying to retrieve an image (not yet imported to GRASS) pixel resolution through a Python Script but so far i was unable to figure out a method to do this. Is this possible? I have used g.proj to retrieve its projection and datum but not the Pixel resolution. THanks for your help Best regards, Antonio R. __ Information from ESET NOD32 Antivirus, version of virus signature database 5234 (20100628) __ The message was checked by ESET NOD32 Antivirus. http://www.eset.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- Eloi Ribeiro GIS Analyst 39,45º -4,40º http://eloiribeiro.wordpress.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user