RE: [GRASS-user] Surface spot interpolates heights for pointis in a feature class from surface

[Moskovitz, Bob]

You should take a look at v.what.rast.  According to the manual:
v.what.rast reads raster value for each point in the vector and updates 
col column in vector attribute table by this value. 
The column should be type number (integer, float, double, ... ). If 
more points have the same category, attribute value is set 
to NULL. If raster values is NULL, attribute value is set to NULL.

-Bob


Bob Moskovitz
Seismic Hazards Mapping Program
California Geological Survey

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.


-Original Message-
From: grass-user-boun...@lists.osgeo.org 
[mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of Susana Arciniegas
Sent: Friday, July 23, 2010 3:26 PM
To: grass-user@lists.osgeo.org
Subject: [GRASS-user] Surface spot interpolates heights for pointis in a 
feature class from surface
Importance: High

I found tool "Tools for point layer (Herramientas para capas de puntos)/ 
Interpolate Raster (Muestrear capa raster)" in SEXTANTE, it is similar to 
"Surface Spot" of ArcGis, I need to found it, in Grass GIS 6.40ORC6, but I can 
not.
The tool interpolates heights for points in a feature class from a surface (ex: 
DTM). The point heights (interpolate heights) are stored in a field that's 
added to the input feature class.
Thank you for your answer
___
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

RE: [GRASS-user] Re: How to merge rasters using feathering/blending ?

[Moskovitz, Bob]

I haven't tried it yet, but what about this:
http://grass.osgeo.org/wiki/GRASS_AddOns#r.smoothpatch

There are other useful addons, so you might want to browse around that
page.

Bob

Bob Moskovitz
Seismic Hazards Mapping Program
California Geological Survey
801 K Street   MS-1231
Sacramento, CA  95814
ph. (916) 323-0043

CONFIDENTIALITY NOTICE: This communication is intended only for the use
of the individual or entity to which it is addressed. This message
contains information from the State of California, California Geological
Survey, which may be privileged, confidential and exempt from disclosure
under applicable law, including the Electronic Communications Privacy
Act. If the reader of this communication is not the intended recipient,
you are hereby notified that any dissemination, distribution, or copying
of this communication is strictly prohibited.



> -Original Message-
> From: grass-user-boun...@lists.osgeo.org [mailto:grass-user-
> boun...@lists.osgeo.org] On Behalf Of kaipi
> Sent: Monday, February 28, 2011 6:27 AM
> To: grass-user@lists.osgeo.org
> Subject: [GRASS-user] Re: How to merge rasters using
> feathering/blending ?
> 
> Gimp is fine if you have a smal number of scenes but I need to process
> large
> amounts of scenes where I need an automated way to feather the seam
> lines.
> 
> --
> View this message in context: http://osgeo-
> org.1803224.n2.nabble.com/How-to-merge-rasters-using-feathering-
> blending-tp6061489p6073453.html
> Sent from the Grass - Users mailing list archive at Nabble.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

[GRASS-user] Interpolate scatter data with faults?

[Moskovitz, Bob]

Hello Grass Users,

 

Googling around, I see that I can interpolate scattered data by using
v.surf.rst and use masks.  I also see that Surfit (see
http://surfit.sourceforge.net/ ) can do it too.  Anyone know of other
methods?

 

Thanks!

 

Bob 

 

Bob Moskovitz

California Geological Survey

Seismic Hazards Zonation Program

 

 

CONFIDENTIALITY NOTICE: This communication is intended only for the use
of the individual or entity to which it is addressed. This message
contains information from the State of California, California Geological
Survey, which may be privileged, confidential and exempt from disclosure
under applicable law, including the Electronic Communications Privacy
Act. If the reader of this communication is not the intended recipient,
you are hereby notified that any dissemination, distribution, or copying
of this communication is strictly prohibited.

 

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Interpolate scatter data with faults?

[Moskovitz, Bob]

Tom,

 

Somehow I did not mention in the body of my message the word "fault".
It must be Friday J  Anyways, I am looking for a way to interpolate data
that takes faults into consideration.  Something like this:
http://surfit.sourceforge.net/surfit/fault__aniso_8tcl-example.html.
I'd use surfit, but it appears to be abandoned (haven't been updated
since 2006).

 

Bob

 

From: Thomas Adams [mailto:thomas.ad...@noaa.gov] 
Sent: Friday, March 23, 2012 2:19 PM
To: Moskovitz, Bob
Cc: grass-user@lists.osgeo.org
Subject: Re: [GRASS-user] Interpolate scatter data with faults?

 

Bob,

 

You may try R and GRASS together, using the R gstat contributed package;
or use GRASS and standalone gstat (http://www.gstat.org/) - both work
very nicely with GRASS.

 

Tom

On Fri, Mar 23, 2012 at 4:26 PM, Moskovitz, Bob
 wrote:

Hello Grass Users,

 

Googling around, I see that I can interpolate scattered data by using
v.surf.rst and use masks.  I also see that Surfit (see
http://surfit.sourceforge.net/ ) can do it too.  Anyone know of other
methods?

 

Thanks!

 

Bob 

 

Bob Moskovitz

California Geological Survey

Seismic Hazards Zonation Program

 

 

CONFIDENTIALITY NOTICE: This communication is intended only for the use
of the individual or entity to which it is addressed. This message
contains information from the State of California, California Geological
Survey, which may be privileged, confidential and exempt from disclosure
under applicable law, including the Electronic Communications Privacy
Act. If the reader of this communication is not the intended recipient,
you are hereby notified that any dissemination, distribution, or copying
of this communication is strictly prohibited.

 


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user





 

-- 

Thomas E Adams
National Weather Service
Ohio River Forecast Center
1901 South State Route 134
Wilmington, OH 45177
 
EMAIL:  thomas.ad...@noaa.gov
 
VOICE:  937-383-0528
FAX:937-383-0033
 

 

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Creating dem using a contour map.

[Moskovitz, Bob]

Hello list,  
 
I have a 20 foot contour map shapefile that has over 47000 line strings.  Many 
of these lines have gaps (see   
 
http://picasaweb.google.com/bobdebm/Problems/photo#5133539789447176242)   in 
many of the contour lines.  I want to use v.surf.rst to make a dem.  Will I 
have problems using this file? 
 
Thanks, 
Bob 

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Intergraph MGGA grid to Grass raster

[Moskovitz, Bob]

Hello list,
 
I am using a GIS called MGE which has a module (called MGGA) that can export 
MGGA grids into Grass raster format.  It produces the directories 
cats,cell,cellhd,cellmisc, and  colr.  Anyone know what I do with this?
 
Thanks in advance,
Bob


Bob Moskovitz 
Seismic Hazard Zonation Project
California Geological Survey 
http://www.conservation.ca.gov/cgs/shzp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.

 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Intergraph MGGA grid to Grass raster

[Moskovitz, Bob]

I put those files in a mapset, but when I try to display it, I get this:
   WARNING: can't get history information from [EMAIL PROTECTED]
   WARNING: can't read range file from [EMAIL PROTECTED]
   ERROR: Unable to read range file.

What should I do now?

Btw, it would be great if text was selectable in dialog boxes...then I won't 
have to type (just cut and paste) the error messages in this email :)

> -Original Message-
> From: [EMAIL PROTECTED] 
> [mailto:[EMAIL PROTECTED] Behalf
> Of Markus Neteler
> Sent: Friday, January 11, 2008 11:41 AM
> To: Moskovitz, Bob
> Cc: Grass-User (E-mail)
> Subject: Re: [GRASS-user] Intergraph MGGA grid to Grass raster
> 
> 
> On Jan 11, 2008 8:06 PM, Moskovitz, Bob
> <[EMAIL PROTECTED]> wrote:
> >
> >
> > Hello list,
> >
> > I am using a GIS called MGE which has a module (called 
> MGGA) that can export
> > MGGA grids into Grass raster format.  It produces the directories
> > cats,cell,cellhd,cellmisc, and  colr.  Anyone know what I 
> do with this?
> 
> Put them into a
> location/mapset/
> directory structure. However, the typical files from 
> PERMANENT are also needed.
> 
> See page bottom of
>  http://download.osgeo.org/grass/grass6_progman/
> for a diagram (rebuilding now, should be ok in 1h from now).
> 
> Markus
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] please help: v.rast.stats

[Moskovitz, Bob]

Hello Jonathan,

I just tried to see if I could grab the latest starspan via cvs.  The cvs 
commands do not seem to be right.  It says:

cvs -d :pserver:[EMAIL PROTECTED]:/cvsroot/starspan login
cvs -d :pserver:[EMAIL PROTECTED]:/cvsroot/starspan checkout modulename 

Shouldn't cvs1 be a fully qualified domain name?

Bob

Bob Moskovitz
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.



> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] Behalf Of Jonathan
> Greenberg
> Sent: Tuesday, March 04, 2008 8:18 PM
> To: Nikos Alexandris
> Cc: grass-user@lists.osgeo.org
> Subject: Re: [GRASS-user] please help: v.rast.stats
> 
> 
> Hey guys, wanted to echo this -- we are continuing development on 
> starspan again -- it works with gdal 1.5, among other 
> improvements, and 
> if you have any feature requests (or if anyone can help us make this 
> into a grass package :), please email me or Carlos Rueda!
> 
> --j
> 
> Nikos Alexandris wrote:
> > On Tue, 2008-03-04 at 17:08 -0800, Dylan Beaudette wrote:
> >> On Tuesday 04 March 2008, Nikos Alexandris wrote:
> >>> On Tue, 2008-03-04 at 12:57 +, christian Brandt wrote:
>  Dear list,
> >> Also-- see the application 'starspan'. works with GRASS, 
> based on GDAL and 
> >> GEOS.
> > 
> > Although Jose Luis Gomez replied to a question of mine 
> about this in the
> > past, I never really gave it a try. Now I realise it's fast and does
> > what I want ;-)
> > 
> > I am still satisfied doing my attempt to script because I learnt 2-3
> > things more.
> > 
> >> Cheers,
> >>
> >> Dylan
> > 
> > 
> > ___
> > grass-user mailing list
> > grass-user@lists.osgeo.org
> > http://lists.osgeo.org/mailman/listinfo/grass-user
> > 
> 
> -- 
> 
> Jonathan A. Greenberg, PhD
> Postdoctoral Scholar
> Center for Spatial Technologies and Remote Sensing (CSTARS)
> University of California, Davis
> One Shields Avenue
> The Barn, Room 250N
> Davis, CA 95616
> Cell: 415-794-5043
> AIM: jgrn307, MSN: [EMAIL PROTECTED], Gchat: jgrn307
> ___
> 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

RE: [GRASS-user] Error with wxpython GUI in RC6

[Moskovitz, Bob]

I think you might need python-wxgtk2.8-dev ?

Bob Moskovitz
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.



> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] Behalf Of 
> Gustavo Martinez
> Sent: Thursday, April 10, 2008 3:02 PM
> To: grass-user@lists.osgeo.org
> Subject: Re: [GRASS-user] Error with wxpython GUI in RC6
> 
> 
> El Thursday 10 April 2008 18:48:15 escribió:
> > dpkg -L python-wxgtk2.8 | grep _core_
> 
> 
> /usr/lib/python2.4/site-packages/wx-2.8-gtk2-unicode/wx/_core_.so
> 
> -- 
> Lic. Gustavo Martínez
> Gestión de Pesquerías
> Sub Secretaría de Pesca de la Nación
> Paseo Colón 982, Anexo Jardín Pesca
> Tel: (54) 11 43492165
> Buenos Aires, República Argentina
> 
> “Que la tierra se vaya haciendo camino ante tus pasos,
> que el viento sople siempre a tus espaldas,
> que el sol brille cálido sobre tu cara,
> que la lluvia caiga suavemente sobre tus campos y,
> hasta tanto volvamos a encontrarnos,
> que Dios te lleve en la palma de su mano.”
> ___
> 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

RE: [GRASS-user] GRASS workshop submissions for FOSS4G2008

[Moskovitz, Bob]

It would be neat to see our borehole database used in a GRASS demo.  This data 
is used in our liquefaction zone maps.  If you want access to data from tens of 
thousand of boreholes, I can give you instructions on how to get it.

-Bob

Bob Moskovitz
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.



> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] Behalf Of Massimiliano
> Cannata
> Sent: Thursday, May 08, 2008 10:29 AM
> To: GRASS user list
> Subject: Re: [GRASS-user] GRASS workshop submissions for FOSS4G2008
> 
> 
> Dear all,
> glad to here that there are willing people.
> Maria Brovelli and Me have submitted the workshop due to the alredy 
> closed deadline.
> The title is:
> Supporting sustainable development: the GRASS opportunity.
> 
> First and foremost thank you for your willingness.
> We will try to take advantages of all interested people, 
> involving you 
> in preparing exercises in the topic of solar energy and risk 
> assessment. 
> Available free dataset are also welcome.
> 
> Massimiliano
> 
> 
> Massimiliano Cannata wrote:
> > Hi,
> > I will go to FOSS4G2008.
> > If anyone is available for sharing the efforts and or for 
> support in 
> > class room,
> > I'm fine to prepare a base workshop on GRASS.
> > Optionally something about DTM and/or risk assessment for 
> > flooding/tsunami
> >
> > Maxi
> >
> >
> >
> > Markus Neteler wrote:
> >> Hi,
> >>
> >> I just got notified that there aren't yet GRASS workshop 
> submissions
> >> for FOSS4G2008. Please consider to submit a proposal (if you go
> >> there):
> >>
> >> http://conference.osgeo.org/index.php/foss4g/2008
> >>
> >> Markus
> >> 
> --
> --
> >>
> >> ___
> >> 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
> >
> 
> 
> -- 
> -
> Dr. Massimiliano Cannata
> Environmental & Geomatic Engineer
> 
> Institute of Earth Scences - SUPSI
> Trevano, C.P. 72, CH-6952 Canobbio, SWITZERLAND
> 
> Tel:+41 (0)58 / 666 62 14  
> Fax:+41 (0)58 / 666 62 09
> E-mail: [EMAIL PROTECTED]
> 
> Web: 
> http://www.ist.supsi.ch
> http://istgis.ist.supsi.ch:8001/geomatica/
> ---
> 
> 
> ___
> 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

RE: [GRASS-user] Plateau / Flatbed detection

[Moskovitz, Bob]

It would be cool to see something like the Topographic Position Index(TPI) in 
GRASS.  

Details here:  http://www.jennessent.com/downloads/TPI-poster-TNC_18x22.pdf

> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] Behalf Of 
> Daniel Victoria
> Sent: Friday, May 16, 2008 12:48 PM
> To: grass
> Subject: [GRASS-user] Plateau / Flatbed detection
> 
> 
> Hi all,
> 
> Has anybody done or heard off a plateau (flatbed) detection? I'm
> trying to automatically detect these features (Brazilian chapadas -
> http://en.wikipedia.org/wiki/Chapada ) using SRTM. Anyone knows if
> it's possible?
> 
> Many thanks
> Daniel
> ___
> 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

RE: [GRASS-user] Plateau / Flatbed detection

[Moskovitz, Bob]

According to the r.terraflow man page, tci is "topographic convergence index, 
defined as the logarithm of the ratio of flow accumulation and local slope".  
Based on what I read about TPI, it does not appear to be the same.

> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] Behalf Of Jonathan
> Greenberg
> Sent: Friday, May 16, 2008 3:51 PM
> Cc: grass
> Subject: Re: [GRASS-user] Plateau / Flatbed detection
> 
> 
> Actually, I think r.terraflow has that!  look for the TCI output!
> 
> --j
> 
> Moskovitz, Bob wrote:
> > It would be cool to see something like the Topographic 
> Position Index(TPI) in GRASS.  
> >
> > Details here:  
> http://www.jennessent.com/downloads/TPI-poster-TNC_18x22.pdf
> >
> >   
> >> -Original Message-
> >> From: [EMAIL PROTECTED]
> >> [mailto:[EMAIL PROTECTED] Behalf Of 
> >> Daniel Victoria
> >> Sent: Friday, May 16, 2008 12:48 PM
> >> To: grass
> >> Subject: [GRASS-user] Plateau / Flatbed detection
> >>
> >>
> >> Hi all,
> >>
> >> Has anybody done or heard off a plateau (flatbed) detection? I'm
> >> trying to automatically detect these features (Brazilian chapadas -
> >> http://en.wikipedia.org/wiki/Chapada ) using SRTM. Anyone knows if
> >> it's possible?
> >>
> >> Many thanks
> >> Daniel
> >> ___
> >> 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
> >
> >   
> 
> -- 
> Jonathan A. Greenberg, PhD
> Postdoctoral Scholar
> Center for Spatial Technologies and Remote Sensing (CSTARS)
> University of California, Davis
> One Shields Avenue
> The Barn, Room 250N
> Davis, CA 95616
> Cell: 415-794-5043
> AIM: jgrn307, MSN: [EMAIL PROTECTED], Gchat: jgrn307
> 
> ___
> 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

RE: [GRASS-user] Plateau / Flatbed detection

[Moskovitz, Bob]

Once you digest all this info, can you summarize it to the list?

Thanks!

> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] Behalf Of 
> Daniel Victoria
> Sent: Friday, May 16, 2008 3:52 PM
> To: grass
> Subject: Re: [GRASS-user] Plateau / Flatbed detection
> 
> 
> Thanks for all the tips guys!
> 
> My first gut reaction was to get the slope and the look for the edges
> of the flatbeds and select only those edges that are close to a very
> flat place.
> 
> The problem is that if there is a big mountain with steep slope
> besides a plain, it will also be characterized as a flatbed.
> 
> I'll take a look at the r.terraflow TCI that Jonathan mentioned and
> all the other references sent
> 
> Cheers
> Daniel
> 
> On Fri, May 16, 2008 at 6:50 PM, Moskovitz, Bob
> <[EMAIL PROTECTED]> wrote:
> > It would be cool to see something like the Topographic 
> Position Index(TPI) in GRASS.
> >
> > Details here:  
> http://www.jennessent.com/downloads/TPI-poster-TNC_18x22.pdf
> >
> >> -Original Message-
> >> From: [EMAIL PROTECTED]
> >> [mailto:[EMAIL PROTECTED] Behalf Of
> >> Daniel Victoria
> >> Sent: Friday, May 16, 2008 12:48 PM
> >> To: grass
> >> Subject: [GRASS-user] Plateau / Flatbed detection
> >>
> >>
> >> Hi all,
> >>
> >> Has anybody done or heard off a plateau (flatbed) detection? I'm
> >> trying to automatically detect these features (Brazilian chapadas -
> >> http://en.wikipedia.org/wiki/Chapada ) using SRTM. Anyone knows if
> >> it's possible?
> >>
> >> Many thanks
> >> Daniel
> >> ___
> >> 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 mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Change NaN to 0 with r.mapcalc(or whatever)?

[Moskovitz, Bob]

Hello All,
 
Is there a way to r.mapcalc or some other command to change NaN values in 
rasters to another value?
 
Bob 

Bob Moskovitz 
Seismic Hazard Zonation Project
California Geological Survey 
http://www.conservation.ca.gov/cgs/shzp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Change NaN to 0 with r.mapcalc(or whatever)?

[Moskovitz, Bob]

Based on what I read from that link I tried:
r.mapcalc "q50_eh_nonan = if(q50_eh != q50_eh,null(),q50_eh)"

It assigned nan to 0.  It would be great if r.null can use such a simple trick.

> -Original Message-
> From: Hamish [mailto:[EMAIL PROTECTED]
> Sent: Tuesday, May 20, 2008 8:34 PM
> To: Grass-User (E-mail); Moskovitz, Bob
> Subject: Re: [GRASS-user] Change NaN to 0 with r.mapcalc(or whatever)?
> 
> 
> Bob wrote:
> > Is there a way to r.mapcalc or some other command to change
> > NaN values in rasters to another value?
> 
> 
> If you mean real "nan" not grass NULL values, see this thread 
> from the grass-dev list from a few days ago:
>   http://thread.gmane.org/gmane.comp.gis.grass.devel/27160
> and the question of where they came from.
> 
> 
> For grass NULL values, r.null is your friend, or 
> isnull(mapname) in r.mapcalc.
> 
> 
> Hamish
> 
> 
> 
> 
>   
> 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Change NaN to 0 with r.mapcalc(or whatever)?

[Moskovitz, Bob]

NaN is an abbreviation for Not a Number.  Here is more info: 
http://en.wikipedia.org/wiki/NaN

I was able to solve my problem by doing this:  r.mapcalc "q50_eh_nonan = 
if(q50_eh != q50_eh,null(),q50_eh)"

> -Original Message-
> From: Daniel Victoria [mailto:[EMAIL PROTECTED]
> Sent: Tuesday, May 20, 2008 5:26 PM
> To: Moskovitz, Bob
> Cc: Grass-User (E-mail)
> Subject: Re: [GRASS-user] Change NaN to 0 with r.mapcalc(or whatever)?
> 
> 
> If by NaN you mean null values, you are looking for the r.null command
> 
> Daniel
> 
> On Tue, May 20, 2008 at 9:11 PM, Moskovitz, Bob
> <[EMAIL PROTECTED]> wrote:
> > Hello All,
> >
> > Is there a way to r.mapcalc or some other command to change 
> NaN values in
> > rasters to another value?
> >
> > Bob
> >
> > Bob Moskovitz
> > Seismic Hazard Zonation Project
> > California Geological Survey
> > http://www.conservation.ca.gov/cgs/shzp
> >
> > CONFIDENTIALITY NOTICE: This communication is intended only 
> for the use of
> > the individual or entity to which it is addressed. This 
> message contains
> > information from the State of California, California 
> Geological Survey,
> > which may be privileged, confidential and exempt from 
> disclosure under
> > applicable law, including the Electronic Communications 
> Privacy Act. If the
> > reader of this communication is not the intended recipient, 
> you are hereby
> > notified that any dissemination, distribution, or copying of this
> > communication is strictly prohibited.
> >
> > ___
> > 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

RE: [GRASS-user] Change NaN to 0 with r.mapcalc(or whatever)?

[Moskovitz, Bob]

Thanks for making this patch available.  I'll update to r31470 and tell you how 
it goes.

> -Original Message-
> From: Glynn Clements [mailto:[EMAIL PROTECTED]
> Sent: Thursday, May 22, 2008 12:27 AM
> To: Moskovitz, Bob
> Cc: [EMAIL PROTECTED]; Grass-User (E-mail)
> Subject: RE: [GRASS-user] Change NaN to 0 with r.mapcalc(or whatever)?
> 
> 
> 
> Moskovitz, Bob wrote:
> 
> > Based on what I read from that link I tried:
> > r.mapcalc "q50_eh_nonan = if(q50_eh != q50_eh,null(),q50_eh)"
> > 
> > It assigned nan to 0.  It would be great if r.null can use 
> such a simple trick.
> 
> Can you try the following patch (committed to SVN trunk as r31470):
> 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] (no subject)

[Moskovitz, Bob]

Hello list,
 
A coworker was asking me if I know of a slope curvature algorithm that lets you 
adjust the kernel size.  Does GRASS have such a program?
 
Bob

Bob Moskovitz 
Seismic Hazard Zonation Project
California Geological Survey 
http://www.conservation.ca.gov/cgs/shzp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.

 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] (no subject)

[Moskovitz, Bob]

Hi Dylan,

It looks like r.param.scale is the module he's looking for...thanks!

Bob

> -Original Message-
> From: Dylan Beaudette [mailto:[EMAIL PROTECTED]
> Sent: Tuesday, June 03, 2008 2:32 PM
> To: grass-user@lists.osgeo.org
> Cc: Moskovitz, Bob
> Subject: Re: [GRASS-user] (no subject)
> 
> 
> Hi Bob,
> 
> Have a look at r.param.scale -- this module allows you to 
> specify the window 
> size.
> 
> Cheers,
> 
> Dylan
> 
> On Tuesday 03 June 2008, Moskovitz, Bob wrote:
> > Hello list,
> >
> > A coworker was asking me if I know of a slope curvature 
> algorithm that lets
> > you adjust the kernel size.  Does GRASS have such a program?
> >
> > Bob
> >
> > Bob Moskovitz
> > Seismic Hazard Zonation Project
> > California Geological Survey
> > http://www.conservation.ca.gov/cgs/shzp prefix = o ns =
> > "urn:schemas-microsoft-com:office:office" />
> >
> > CONFIDENTIALITY NOTICE: This communication is intended only 
> for the use of
> > the individual or entity to which it is addressed. This 
> message contains
> > information from the State of California, California 
> Geological Survey,
> > which may be privileged, confidential and exempt from 
> disclosure under
> > applicable law, including the Electronic Communications 
> Privacy Act. If the
> > reader of this communication is not the intended recipient, 
> you are hereby
> > notified that any dissemination, distribution, or copying of this
> > communication is strictly prohibited.
> 
> 
> 
> -- 
> Dylan Beaudette
> Soil Resource Laboratory
> http://casoilresource.lawr.ucdavis.edu/
> University of California at Davis
> 530.754.7341
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] slope curvature algorithm - adjust the kernel size

[Moskovitz, Bob]

Hi Markus,

I do have the 3rd edition of your book, but it seems like r.param.scale is not 
in your book (at least I don't see a reference in the index).  Since we do have 
lidar, I'd like to take a look at Helena's IEEE paper.  Is it (or a preprint) 
available on the internet?  Also I wanted to know what the illustration is 
showing.  It looks like curvature is on the left and slope on the right which 
are animated to based on some parameter?

Bob

Bob Moskovitz
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.



> -Original Message-
> From: [EMAIL PROTECTED] 
> [mailto:[EMAIL PROTECTED] Behalf
> Of Markus Neteler
> Sent: Wednesday, June 04, 2008 12:28 PM
> To: Moskovitz, Bob
> Cc: Grass-User (E-mail)
> Subject: Re: [GRASS-user] slope curvature algorithm - adjust 
> the kernel
> size
> 
> 
> Hi Bob,
> 
> below answers from Helena and Jaro
> 
> Markus
> 
> On Wed, Jun 4, 2008 at 2:27 AM,  <[EMAIL PROTECTED]> wrote:
> > Jaro Hofierka wrote:
> >> I think r.param.scale is the module that this guy is looking for.
> >
> > that is right. In case he is working with the latest lidar data for
> > California, I think we have a nice example in the book for 
> r.param.scale
> > for the high resolution data showing what neighborhood you 
> need to get
> > meaningful curvature maps from lidar-based DEM.
> >
> >  If he has the bare earth point cloud available he can get 
> even better
> > curvature maps at an appropriate level of detail while keeping the
> > resolution high from v.surf.rst (I think we have that 
> described in the
> > IEEE paper 2005 and here is an illustration)
> > http://skagit.meas.ncsu.edu/~helena/measwork/jockeys/a57tenslpcn.gif
> >
> > Helena
> >
> 
> 
> On Tue, Jun 3, 2008 at 10:23 PM, Moskovitz, Bob
> <[EMAIL PROTECTED]> wrote:
> > Hello list,
> >
> > A coworker was asking me if I know of a slope curvature 
> algorithm that lets
> > you adjust the kernel size.  Does GRASS have such a program?
> >
> > Bob
> >
> > Bob Moskovitz
> > Seismic Hazard Zonation Project
> > California Geological Survey
> > http://www.conservation.ca.gov/cgs/shzp
> >
> 
> -- 
> Open Source Geospatial Foundation
> http://www.osgeo.org/
> http://www.grassbook.org/
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] v.dissolve bugfeature / workaround needed

[Moskovitz, Bob]

Hello Moritz and Peter,

I have the v.dissolve problem too.  I'm running GRASS 6.4.svn (is there a way 
to tell what svn revision I'm using?) on Hardy Heron X86_64.  Hope this info 
tells you something:

 > v.dissolve [EMAIL PROTECTED] output=efz layer=1 column=FaultName
Segmentation fault
dbmi: Protocol error
ERROR: Cannot open select cursor: SELECT cat, FaultName FROM
   alquist_priolo_eq_zones_mrg_ll83 ORDER BY FaultName
ERROR: Unable to open vector map <[EMAIL PROTECTED]> on topology level 2
WARNING: Table  linked to vector map  does not exist
ERROR: Vector map  not found in the current mapset

 > v.info alquist_priolo_eq_zones_mrg_ll83
 ++
 | Layer:   alquist_priolo_eq_zones_mrg_ll83  |
 | Mapset:  PERMANENT |
 | Location:CGS_SW_FaultLines |
 | Database:/home/bobm/Data/Grass |
 | Title: |
 | Map scale:   1:1   |
 | Map format:  native|
 | Name of creator: bobm  |
 | Organization:  |
 | Source date: Fri Nov 17 11:47:54 2006  |
 ||
 |   Type of Map:  vector (level: 2)  |
 ||
 |   Number of points:   0   Number of areas:  2242   |
 |   Number of lines:0   Number of islands:788|
 |   Number of boundaries:   4630Number of faces:  0  |
 |   Number of centroids:2224Number of kernels:0  |
 ||
 |   Map is 3D:  0|
 |   Number of dblinks:  1|
 ||
 | Projection: Lambert Conformal Conic|
 |   N:  4959352.74991602S:  3929916.63677083 |
 |   E:   341967.39773093W:  -438465.59887675 |
 ||
 |   Digitization threshold: 0|
 |   Comments:|
 ||
 ++

 > v.info -c alquist_priolo_eq_zones_mrg_ll83
Displaying column types/names for database connection of layer 1:
INTEGER|cat
DOUBLE PRECISION|ID
CHARACTER|COUNTY
CHARACTER|QUADNAME
CHARACTER|FaultName
CHARACTER|HistoricSR

Bob Moskovitz
Research Analyst I
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.


> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] Behalf Of Moritz Lennert
> Sent: Thursday, August 07, 2008 8:45 AM
> To: Peter Löwe
> Cc: grass-user@lists.osgeo.org
> Subject: Re: [GRASS-user] v.dissolve bugfeature / workaround needed
> 
> 
> On 07/08/08 17:20, Peter Löwe wrote:
> > Hello Milton,
> > 
> > thanks for that beautiful R-hack !
> > Since I am working on a minimal "bread and butter" machine 
> I went for Moritz' solution but it's a neat idea to use R for 
> the task.
> > 
> > Eventually, GIS-hacking should become part of the Olympics 
> and you all should try to get medals :-)
> 
> 
> Actually, I stand corrected. Daniel was right: v.reclass does allow 
> character columns and v.dissolve, therefore, does as well. I 
> just tested 
> with a map as it works with current svn trunk and and old 
> pre-release 6.3:
> 
>  > v.info -c NUTS2_2006
> Displaying column types/names for database connection of layer 1:
> INTEGER|cat
> INTEGER|objectid
> CHARACTER|cntr_code
> INTEGER|stat_levl_
> 

[GRASS-user] r.to.vect and then v.clean using tool=rmarea creates polygons instead of holes.

[Moskovitz, Bob]

Hello Grass-Users,
 
I'm having problems with this sequence of commands:
r.to.vect -s input=lszone output=lszone10 feature=area
v.clean input=lszone10 output=lszone10_clean tool=rmarea thresh=2800
 
Below is the output.
 
Overlaying lszone and lszone10_clean, I see that there are large holes that 
have been replaced with polygons.  Is this a bug?
 
Btw, I've tried this with QGIS on Windows XP as well as GRASS 6.4.svn with the 
same results.
 
Bob
 

==
r.to.vect -s in=lszone out=lszone10 feature=area
Extracting areas...
 100%
 100%
Building topology for vector map ...
Registering primitives...
50074 primitives registered
629002 vertices registered
Building areas...
 100%
27613 areas built
17107 isles built
Attaching islands...
 100%
Attaching centroids...
 100%
Number of nodes: 39568
Number of primitives: 50074
Number of points: 0
Number of lines: 0
Number of boundaries: 33886
Number of centroids: 16188
Number of areas: 27613
Number of isles: 17107
Number of areas without centroid: 11425
r.to.vect complete.
v.clean in=lszone10 out=lszone10_clean tool=rmarea thresh=2800
--
Tool: Threshold
Remove small areas: 2.80e+03
--
Copying vector lines...
Rebuilding parts of topology...
Building topology for vector map ...
Registering primitives...
50074 primitives registered
629002 vertices registered
Building areas...
 100%
27613 areas built
17107 isles built
Attaching islands...
 100%
Attaching centroids...
 100%
Number of nodes: 39568
Number of primitives: 50074
Number of points: 0
Number of lines: 0
Number of boundaries: 33886
Number of centroids: 16188
Number of areas: 27613
Number of isles: 17107
Number of areas without centroid: 11425
--
Tool: Remove small areas
26251 areas of total size 9.34119e+06 removed
--
Rebuilding topology for output vector map...
Building topology for vector map ...
Registering primitives...
4913 primitives registered
222687 vertices registered
Building areas...
 100%
1362 areas built
1127 isles built
Attaching islands...
 100%
Attaching centroids...
 100%
Number of nodes: 4678
Number of primitives: 4913
Number of points: 0
Number of lines: 0
Number of boundaries: 4185
Number of centroids: 728
Number of areas: 1362
Number of isles: 1127
Number of duplicate centroids: 3
Number of areas without centroid: 637
 

Bob Moskovitz 
Research Analyst I 
Seismic Hazard Zonation Project
California Geological Survey 
http://www.conservation.ca.gov/cgs/shzp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] r.to.vect and then v.clean using tool=rmarea createspolygons instead of holes.

[Moskovitz, Bob]

I followed your suggestion, but I'm still getting the same problem.  It seems 
strange that v.to.rast would create a vector file that has flawed topology.  Do 
you know of any alternative to v.clean with rmarea?

Btw, I also tried to use v.generalize with a very small threshold and then 
v.clean.  This method works for me.

> -Original Message-
> From: Alex Mandel [mailto:tech_...@wildintellect.com]
> Sent: Thursday, December 11, 2008 3:14 PM
> To: Moskovitz, Bob
> Cc: Grass-User (E-mail)
> Subject: Re: [GRASS-user] r.to.vect and then v.clean using tool=rmarea
> createspolygons instead of holes.
> 
> 
> Moskovitz, Bob wrote:
> > Hello Grass-Users,
> >  
> > I'm having problems with this sequence of commands:
> > r.to.vect -s input=lszone output=lszone10 feature=area
> > v.clean input=lszone10 output=lszone10_clean 
> tool=rmarea thresh=2800
> >  
> > Below is the output.
> >  
> > Overlaying lszone and lszone10_clean, I see that there are 
> large holes
> > that have been replaced with polygons.  Is this a bug?
> >  
> > Btw, I've tried this with QGIS on Windows XP as well as 
> GRASS 6.4.svn
> > with the same results.
> >  
> > Bob
> >  
> 
> Hmm, I'm wondering if there's a flaw in the topology built by 
> r.to.vect
> 
> What happens if you run v.build before v.clean?
> 
> At second glance I see how if a small polygon gets removed 
> and it shared
> a border with the big polygon, that that border might not exist now.
> maybe add the bpol to tool=rmarea,bpol
> 
> Alex
> 
> 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] r.to.vect and then v.clean using tool=rmarea createspolygons instead of holes.

[Moskovitz, Bob]

Thank you Hamish...your set of commands worked great!

It still bothers me that r.to.vect would produce invalid topology.  Btw, I also 
tried to use v.generalize but my vectyor seem to lose its attribute table.

> -Original Message-
> From: Hamish [mailto:hamis...@yahoo.com]
> Sent: Friday, December 12, 2008 1:17 AM
> To: Alex Mandel; Moskovitz, Bob
> Cc: Grass-User (E-mail)
> Subject: RE: [GRASS-user] r.to.vect and then v.clean using tool=rmarea
> createspolygons instead of holes.
> 
> 
> Bob Moskovitz wrote:
> > I followed your suggestion, but I'm still getting the
> > same problem.  It seems strange that v.to.rast would create
> > a vector file that has flawed topology.  Do you know of any
> > alternative to v.clean with rmarea?
> 
> v.db.addcol area
> v.to.db option=area
> v.extract where="area > 2800"
> ?
> 
> 
> > Btw, I also tried to use v.generalize with a very small
> > threshold and then v.clean.  This method works for me.
> 
> 
> Hamish
> 
> 
> 
>   
> 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] r.to.vect and then v.clean using tool=rmarea createspolygonsinstead of holes.

[Moskovitz, Bob]

Upon closer inspection, I see what you are talking about.  I do want to remove 
holes <= 2800 m, so this method won't work for me.  So, what would cause 
r.to.vect to create areas without centroids?  And why haven't v.clean rmarea 
filled all of my holes since they all don't have centroids?

> -Original Message-
> From: Markus Metz [mailto:markus.metz.gisw...@googlemail.com]
> Sent: Friday, December 12, 2008 5:31 AM
> To: grass-user@lists.osgeo.org; Moskovitz, Bob
> Subject: Re: [GRASS-user] r.to.vect and then v.clean using tool=rmarea
> createspolygonsinstead of holes.
> 
> 
> 
> 
> Hamish wrote:
> > Bob Moskovitz wrote:
> >   
> >> I followed your suggestion, but I'm still getting the
> >> same problem.  It seems strange that v.to.rast would create
> >> a vector file that has flawed topology.  Do you know of any
> >> alternative to v.clean with rmarea?
> >> 
> >
> > v.db.addcol area
> > v.to.db option=area
> > v.extract where="area > 2800"
> > ?
> >   
> That only deletes the centroids of these small areas, the 
> boundaries remain.
> 
>  From your original post:
> vector  has
> 27613 areas
> 17107 isles
> Number of areas without centroid: 11425
> 
> vector  has
> Number of areas: 1362
> Number of isles: 1127
> Number of duplicate centroids: 3
> Number of areas without centroid: 637
> 
> Further on, 26251 areas have been removed.
> v.clean tool=rmarea also removes areas without centroids, that's why
> holes have been replaced by polygons, boundaries of small areas are
> deleted. If you want to merge small areas with adjacent large areas,
> v.clean would do the job, but for all areas, also areas 
> without centroid.
> If you want to keep only larger areas and these unchanged, 
> the method of
> Hamish with v.extract would do that.
> Or you could do
> v.db.addcol map=lszone10 columns="area DOUBLE"
> v.to.db map=lszone10 option=area columns=area
> and then use where="area > 2800" in further operations.
> 
> Markus
> 
> 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] r.to.vect and then v.clean using tool=rmarea createspolygonsinstead of holes.

[Moskovitz, Bob]

How do I pack up my raster so I can let you see what is happening?  I'd like to 
send it to you or any other dev if possible.

> -Original Message-
> From: neteler.os...@gmail.com 
> [mailto:neteler.os...@gmail.com]on Behalf
> Of Markus Neteler
> Sent: Friday, December 12, 2008 1:17 PM
> To: Moskovitz, Bob
> Cc: Markus Metz; grass-user@lists.osgeo.org
> Subject: Re: [GRASS-user] r.to.vect and then v.clean using tool=rmarea
> createspolygonsinstead of holes.
> 
> 
> On Fri, Dec 12, 2008 at 9:01 PM, Moskovitz, Bob
>  wrote:
> > So, what would cause r.to.vect to create areas without centroids?
> 
> I have tried r.to.vect with a Spearfish map and it looks ok:
> 
> g.region rast=trn.sites
> d.mon x0
> d.rast trn.sites
> 
> r.to.vect trn.sites out=trn feature=area
> Extracting areas...
>  100%
>  100%
> Building topology for vector map ...
> Registering primitives...
> 74 primitives registered
> 903 vertices registered
> Building areas...
>  100%
> 20 areas built
> 2 isles built
> Attaching islands...
>  100%
> Attaching centroids...
>  100%
> Number of nodes: 56
> Number of primitives: 74
> Number of points: 0
> Number of lines: 0
> Number of boundaries: 54
> Number of centroids: 20
> Number of areas: 20
> Number of isles: 2
> r.to.vect complete.
> 
> 
> d.vect trn
> 
> d.vect trn disp=topo
> 
> v.out.ascii trn format=standard
> -> boundaries and centroids
> 
> May a potential problem in r.to.vect depend on the shape or the
> boundary lengths? Obviously, a reproducible example is desired.
> 
> Markus
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Aspect direction in r.slope.aspect

[Moskovitz, Bob]

Hello Grass Users,
 
I was just wondering why aspect is measured from the east and increased 
counterclockwise instead of the north and clockwise?  It seems to run counter 
to what is done in other GISs.
 
Bob
 

Bob Moskovitz 
Research Analyst I 
Seismic Hazard Zonation Project
California Geological Survey 
http://www.conservation.ca.gov/cgs/shzp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Aspect direction in r.slope.aspect

[Moskovitz, Bob]

Thank you for the tip.  Btw, I was made aware of this "problem" when a coworker 
took our landslide hazard workflow and tried to apply it to grass.  Everything 
look great, but when we took another look at the various layers, we were 
shocked to see how aspect was defined.  I know...we should have read the 
manuals, but everybody knows how aspect is defined.  Anyways, we can live with 
another step in our process, but it would be great if there was a document that 
shows the differences between GIS systems (such as ESRI, MapInfo, etc), so that 
we can keep surprises at a minimum. :-)

Bob Moskovitz
Research Analyst I
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.



> -Original Message-
> From: Hamish [mailto:hamis...@yahoo.com]
> Sent: Tuesday, January 20, 2009 12:50 PM
> To: Grass-User (E-mail); Moskovitz, Bob
> Subject: Re: [GRASS-user] Aspect direction in r.slope.aspect
> 
> 
> Bob Moskovitz wrote:
> > I was just wondering why aspect is measured from the east
> > and increased counterclockwise instead of the north and
> > clockwise?  It seems to run counter to what is done in other
> > GISs.
> 
> 
> It is just measuring theta positive from the x-axis, in the 
> mathematical
> sense.
> 
> It stays like that for historical reasons, probably been that 
> way for 20+
> years. (aka we can only guess at the original logic)
> 
> Convert between that and the compass (CW from north) convention with
> "90 - angle". Most modules which deal with aspect maps should already
> have an option to read the different type.
> 
> 
> Hamish
> 
> 
> 
>   
> 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Aspect direction in r.slope.aspect

[Moskovitz, Bob]

Sounds like a great idea!  If sometrhing like this can be done, we would have 
one less reason to shy away from FOSS.

> 
> For my 2 cents worth, it seems to make a lot more sense for a  
> *geographic* information system for all output to follow the same  
> spatial organizational standards. I understand the history of 
> the east  
> is 0 convention for parts of GRASS. And I also appreciate the  
> importance of not breaking code-dependent features within versions.  
> However, that does not mean that we should keep a 
> non-standard way of  
> measuring direction for select modules (like r.slope.aspect), while  
> measuring from north for others, simply because the program started  
> out that way. Version changes are a time when we can rethink and  
> standardize different modules that have evolved 
> independently. Scripts  
> are likely to break across the 6/7 transition for other 
> reasons anyway.
> 
> That said, a functionally simple approach that would not be quite so  
> disruptive would be to add a flag to switch to count from east mode  
> (with count from north as the default) or even a flag to switch to  
> count from north mode (with count from east being the default if  
> backward compatibility is indeed very important in this case).
> 
> It is important to get an idea of how many things actually 
> would need  
> to be changed.
> 
> Michael
> 
> ___
> 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] Tips on how to delete areas with no categories?

[Moskovitz, Bob]

Hello Grass-Users,
 
After using r.to.vect, I have areas with no categories.  These areas appear to 
be holes.  This is what the map look like with info from v.what:  
http://bobdebm.googlepages.com/example_area.png/example_area-full;init:.png
 
How do I remove these areas?
 
Thanks in advance.
 
Bob
 

Bob Moskovitz 
Research Analyst I 
Seismic Hazard Zonation Project
California Geological Survey 
http://www.conservation.ca.gov/cgs/shzp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] DAFF/CSIRO soils map

[Moskovitz, Bob]

Richard,  I believe you need to use v.in.e00.  -Bob

> -Original Message-
> From: grass-user-boun...@lists.osgeo.org
> [mailto:grass-user-boun...@lists.osgeo.org]on Behalf Of 
> Richard Chirgwin
> Sent: Tuesday, February 10, 2009 12:36 PM
> To: grass-user@lists.osgeo.org
> Subject: [GRASS-user] DAFF/CSIRO soils map
> 
> 
> Hamish,
> 
> You mentioned this dataset in a previous discussion ...
> > Jarl wrote:
> >   
> >> > The data is an Australian soils atlas including hydraulic
> >> > properties in the dbf. 
> >> 
> >
> > this?  http://www.daff.gov.au/brs/data-tools/daas-download
> >   
> Being curious I took a look ... the data uncompresses as 
> csiro2mgvfc.e00
> but r.in.arc fails on this:
> > Illegal line in header
> > EXP 1 /REFERENCE/NATIONAL/ARCINFO/SOIL/COVERS/CSIRO2MG.E00
> > "yllcorner" field missing from header
> > "xllcorner" field missing from header
> > "nrows" field missing from header
> > "ncols" field missing from header
> > "cellsize" field missing from header
> > "nodata_value" field missing from header
> > Can't get cell header
> >
> What would cause this?
> 
> Richard C
> ___
> 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

RE: [GRASS-user] Tips on how to delete areas with no categories?

[Moskovitz, Bob]

Richard,

I really do need more education on the grass vector side of things.  Thanks for 
the clue...I got it figured out :-)

Bob


> -Original Message-
> From: grass-user-boun...@lists.osgeo.org
> [mailto:grass-user-boun...@lists.osgeo.org]on Behalf Of 
> Richard Chirgwin
> Sent: Tuesday, February 10, 2009 2:41 PM
> To: grass-user@lists.osgeo.org
> Subject: [GRASS-user] Tips on how to delete areas with no categories?
> 
> 
> 
> > Message: 2
> > Date: Tue, 10 Feb 2009 11:40:15 -0800
> > From: "Moskovitz, Bob" 
> > Subject: [GRASS-user] Tips on how to delete areas with no 
> categories?
> > To: "Grass-User (E-mail)" 
> > Message-ID:
> > 
> 
> > Content-Type: text/plain; charset="iso-8859-1"
> >
> > Hello Grass-Users,
> >  
> > After using r.to.vect, I have areas with no categories.  
> These areas appear to be holes.  This is what the map look 
> like with info from v.what:  
> http://bobdebm.googlepages.com/example_area.png/example_area-f
ull;init:.png
>  
> How do I remove these areas?
>  
> Thanks in advance.
>  
> Bob
>  
>
> Bob Moskovitz 
> Research Analyst I 
> Seismic Hazard Zonation Project
> California Geological Survey 
> http://www.conservation.ca.gov/cgs/shzp "urn:schemas-microsoft-com:office:office" />
>   
Bob,

What if you use v.category to attach cats to the empty areas, and then
use v.edit to delete the appropriate category range?

Cheers,
RC
___
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

RE: [GRASS-user] Open source GIS cataloging software?

[Moskovitz, Bob]

Recently, I came across this:  
http://spatialguru.com/ideas/data_cataloguing_background


Bob Moskovitz
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.



> -Original Message-
> From: grass-user-boun...@lists.osgeo.org
> [mailto:grass-user-boun...@lists.osgeo.org]on Behalf Of Jonathan
> Greenberg
> Sent: Wednesday, February 11, 2009 1:16 PM
> To: grass-user@lists.osgeo.org; r-sig-...@stat.math.ethz.ch
> Subject: [GRASS-user] Open source GIS cataloging software?
> 
> 
> Sorry for the cross-posting.  I was wondering if anyone has a unix 
> (preferably) or windows program that can spider a directory, 
> recursively 
> searching for raster and vector data, and create bounding box 
> polygons 
> for each raster/vector it finds, with attributes indicating the 
> path-to-file.  Thanks!
> 
> --j
> 
> -- 
> 
> Jonathan A. Greenberg, PhD
> Postdoctoral Scholar
> Center for Spatial Technologies and Remote Sensing (CSTARS)
> University of California, Davis
> One Shields Avenue
> The Barn, Room 250N
> Davis, CA 95616
> Cell: 415-794-5043
> AIM: jgrn307, MSN: jgrn...@hotmail.com, Gchat: jgrn307 
> 
> ___
> 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

RE: [GRASS-user] Is it possible to work on .mdb geomedia file?

[Moskovitz, Bob]

I'm also interested in using Geomedia .mdb files with grass.  Btw, do you use 
geomedia Grid?

Bob Moskovitz
Research Analyst I
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.




> -Original Message-
> From: grass-user-boun...@lists.osgeo.org
> [mailto:grass-user-boun...@lists.osgeo.org]on Behalf Of Enrico g
> Sent: Thursday, April 09, 2009 8:41 AM
> To: grass-user@lists.osgeo.org
> Subject: [GRASS-user] Is it possible to work on .mdb geomedia file?
> 
> 
> I've searched Internet info for this but I've found very few 
> and dated infos..
> 
> Is it possible to work (import/export) with .mdb geomedia file?
> 
> The only solution I've founded is the one from Safe Software.. anybody
> knows how much it costs?
> 
> Thank you
> ___
> 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

RE: [GRASS-user] exporting areas to shapefiles without filling holes

[Moskovitz, Bob]

Hello Hamish,

I had a possibly related problem with a vector that went through v.clean.  I 
filed this ticket:  http://trac.osgeo.org/grass/ticket/397.

I noticed that layer 0 had the filled holes while layer 1 was the good 
polygons.  v.clean seemed to merge those two layers.  With Markus M's fix, I 
solved my problem.  Which reminds me...should I close this ticket since it 
works for me?

Bob

Bob Moskovitz
Research Analyst I
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.



> -Original Message-
> From: grass-user-boun...@lists.osgeo.org
> [mailto:grass-user-boun...@lists.osgeo.org]on Behalf Of Hamish
> Sent: Thursday, April 09, 2009 12:02 AM
> To: grass list
> Subject: [GRASS-user] exporting areas to shapefiles without filling
> holes
> 
> 
> Hi,
> 
> I've got a grass vector map containing some areas with holes. It's a
> polygon of a river containing some islands; it looks nice in grass,
> centroids with cat in the wet parts, no centroids in the islands.
> There is no DB table.
> 
> after exporting with 'v.out.ogr type=area' to a shapefile, when I load
> the shp file into QGIS 0.8.1 I see that all the holes have 
> been filled in.
> The boundaries are still intact. using qgis's editing tools I can cut
> them out, but it takes a bit of manual labour to fix them all.
> 
> what am I doing wrong?
> 
> 
> thanks,
> Hamish
> ___
> 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] Write python grass program with osgeo4w

[Moskovitz, Bob]

Hello list,
 
A coworker and I have been trying to use the first example in  
 
http://download.osgeo.org/grass/grass6_progman/pythonlib.html.  We discovered 
that "import grass.script as grass" has to be changed to "import grass".  Below 
is the example and the error we are getting is this:
 
C:\>m.test.py
Traceback (most recent call last):
  File "C:\OSGeo4W\apps\grass\GRASS-~1.0SV/scripts\m.test.py", line 31, in 

options, flags = grass.parser()
  File "C:\OSGeo4W\apps\grass\GRASS-~1.0SV\etc\python\grass.py", line 208, in 
parser
os.execvp("g.parser", [name] + argv)
  File "C:\OSGeo4W\apps\Python25\lib\os.py", line 353, in execvp
_execvpe(file, args)
  File "C:\OSGeo4W\apps\Python25\lib\os.py", line 389, in _execvpe
func(fullname, *argrest)
OSError: [Errno 22] Invalid argument
 
 
Any ideas what is going on?
 
Bob
 
 
 
#
#!/usr/bin/env python
 
#%Module
#% description: Checks if vector map is 3D
#% keywords: vector
#%End
#%option
#% key: map
#% type: string
#% gisprompt: old,vector,vector
#% key_desc: name
#% description: Name of vector map 
#% required: yes
#%end
 
import sys
import grass
 
def main():
info = grass.parse_command('v.info',
   flags = 't',
   map = options['map'])
if info['map3d'] == '1':
print 'Vector map is 3D'
else:
print 'Vector map is 2D'
 
return 0
 
if __name__ == "__main__":
options, flags = grass.parser()
sys.exit(main())

 

Bob Moskovitz 
Research Analyst I 
Seismic Hazard Zonation Project
California Geological Survey 
http://www.conservation.ca.gov/cgs/shzp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Write python grass program with osgeo4w

[Moskovitz, Bob]

I see that GRASS in OSGEO4W is only up to r36506.  Hopefully there will be an 
update to OSGEO4W soon...maybe even an update to GRASS 6.4 RC5 :-)


-Original Message-
From: Glynn Clements [mailto:gl...@gclements.plus.com]
Sent: Sat 6/13/2009 7:56 AM
To: Moskovitz, Bob
Cc: Grass-User (E-mail)
Subject: Re: [GRASS-user] Write python grass program with osgeo4w
 

Moskovitz, Bob wrote:

> A coworker and I have been trying to use the first example in
> <http://download.osgeo.org/grass/grass6_progman/pythonlib.html>
> http://download.osgeo.org/grass/grass6_progman/pythonlib.html.  We
> discovered that "import grass.script as grass" has to be changed to
> "import grass".  Below is the example and the error we are getting is
> this:
>  
> C:\>m.test.py
> Traceback (most recent call last):
>   File "C:\OSGeo4W\apps\grass\GRASS-~1.0SV/scripts\m.test.py", line 31, in 
> 
> options, flags = grass.parser()
>   File "C:\OSGeo4W\apps\grass\GRASS-~1.0SV\etc\python\grass.py", line 208, in 
> parser
> os.execvp("g.parser", [name] + argv)
>   File "C:\OSGeo4W\apps\Python25\lib\os.py", line 353, in execvp
> _execvpe(file, args)
>   File "C:\OSGeo4W\apps\Python25\lib\os.py", line 389, in _execvpe
> func(fullname, *argrest)
> OSError: [Errno 22] Invalid argument
>  
>  
> Any ideas what is going on?

This is a bug in grass.py; Windows requires "g.parser.exe" rather than
"g.parser". This was fixed with:

http://trac.osgeo.org/grass/changeset/37075

-- 
Glynn Clements 

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Write python grass program with osgeo4w

[Moskovitz, Bob]

I changed grass.py with the changes in changeset 37075.  Although I'm no longer 
getting any error messages, I am not able to see my script help by simply 
typing "m.test.py --help".  

When I try "python C:\OSGeo4W\apps\grass\grass-6.4.0svn\scripts\m.test.py 
--help", I do see the help.  Without "--help", m.test.py doesn't start up in a 
gui.

What must be done to make python scripts to behave like regular grass commands?

Bob

> -Original Message-
> From: Glynn Clements [mailto:gl...@gclements.plus.com]
> Sent: Saturday, June 13, 2009 7:56 AM
> To: Moskovitz, Bob
> Cc: Grass-User (E-mail)
> Subject: Re: [GRASS-user] Write python grass program with osgeo4w
> 
> 
> 
> Moskovitz, Bob wrote:
> 
> > A coworker and I have been trying to use the first example in
> > <http://download.osgeo.org/grass/grass6_progman/pythonlib.html>
> > http://download.osgeo.org/grass/grass6_progman/pythonlib.html.  We
> > discovered that "import grass.script as grass" has to be changed to
> > "import grass".  Below is the example and the error we are 
> getting is
> > this:
> >  
> > C:\>m.test.py
> > Traceback (most recent call last):
> >   File 
> "C:\OSGeo4W\apps\grass\GRASS-~1.0SV/scripts\m.test.py", line 
> 31, in 
> > options, flags = grass.parser()
> >   File 
> "C:\OSGeo4W\apps\grass\GRASS-~1.0SV\etc\python\grass.py", 
> line 208, in parser
> > os.execvp("g.parser", [name] + argv)
> >   File "C:\OSGeo4W\apps\Python25\lib\os.py", line 353, in execvp
> > _execvpe(file, args)
> >   File "C:\OSGeo4W\apps\Python25\lib\os.py", line 389, in _execvpe
> > func(fullname, *argrest)
> > OSError: [Errno 22] Invalid argument
> >  
> >  
> > Any ideas what is going on?
> 
> This is a bug in grass.py; Windows requires "g.parser.exe" rather than
> "g.parser". This was fixed with:
> 
>   http://trac.osgeo.org/grass/changeset/37075
> 
> -- 
> Glynn Clements 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Problems with GRASS python code in osgeo4w

[Moskovitz, Bob]

Hello List,
 
I'm trying to figure out how to write python scripts using grass.py under 
osgeo4w.  I was hoping that python script (and the binary commands) would 
behave similarly to the grass environment on my Ubuntu box.  You know...just 
type the command and the gui for that command comes up.  What I have to do now 
is type "python c:\osgeo4w\apps\grass\grass-6.4.0svn\scripts\m.dipslope.py" 
instead of just type "m.dipslope.py".  I've also notice that the command gui or 
even the command help does not come up when you just enter the command.  
Ultimately I would like to put my commands in the QGIS GRASS module tree.
 

Bob Moskovitz 
Research Analyst I 
Seismic Hazard Zonation Project
California Geological Survey 
http://www.conservation.ca.gov/cgs/shzp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.

 
 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Problems with GRASS python code in osgeo4w

[Moskovitz, Bob]

Thanks for the info.  I didn't know about the assoc and ftype.  Looking at the 
output of these two commands, I see that my system is set up correctly.  

So, I still don't know why the grass.parser() is not working for me.  I even 
modified my program to use pdb to see what is happening.  The debugger is 
aborted right after I reach os.execvp("g.parser.exe", [name] + argv) in 
grass.py.  But if I use "python 
c:\osgeo4w\apps\grass\grass-6.4.0svn\scripts\m.dipslope.py --help", I get the 
expected results.

Bob Moskovitz
Research Analyst I
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.



> -Original Message-
> From: Glynn Clements [mailto:gl...@gclements.plus.com]
> Sent: Friday, June 19, 2009 7:03 PM
> To: Moskovitz, Bob
> Cc: Grass-User (E-mail)
> Subject: Re: [GRASS-user] Problems with GRASS python code in osgeo4w
> 
> 
> 
> Moskovitz, Bob wrote:
> 
> > I'm trying to figure out how to write python scripts using grass.py
> > under osgeo4w.  I was hoping that python script (and the binary
> > commands) would behave similarly to the grass environment 
> on my Ubuntu
> > box.  You know...just type the command and the gui for that command
> > comes up.  What I have to do now is type "python
> > c:\osgeo4w\apps\grass\grass-6.4.0svn\scripts\m.dipslope.py" 
> instead of
> > just type "m.dipslope.py".  I've also notice that the command gui or
> > even the command help does not come up when you just enter 
> the command.
> > Ultimately I would like to put my commands in the QGIS GRASS module
> > tree.
> 
> 1. Ensure that .py files are configured to be run with the Python
> interpreter. You should be able to double-click on a .py file in
> Explorer and have it run. This should be done by the Python installer,
> but I don't know if the OSGeo4W installer does this.
> 
> If it isn't, you can use e.g.:
> 
>   assoc .py=python.file
>   ftype python.file="C:\Program Files\Python25\python.exe" "%1"
> 
> [Change the pathname to wherever Python is installed.]
> 
> 2. The PATHEXT environment variable needs to contain ".PY" if you want
> to use "m.dispslope" rather than "m.dispslope.py". You can do this
> temporarily with:
> 
>   set PATHEXT=%PATHEXT%;.PY
> 
> or you can make it persistent either using the Control Panel or
> through the registry.
> 
> For the Control Panel in XP, it's:
> 
> My Computer
>  Control Panel
>   System
>Advanced (Tab)
> Environment Variables
> 
> For the registry, the system-wide setting is taken from the PATHEXT
> value of the key:
> 
> HKLM\SYSTEM\CurrentControlSet\Control\Session Manager\Environment
> 
> Or you can override it for the current user by adding a PATHEXT value
> to the key HKCU\Environment
> 
> -- 
> Glynn Clements 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Problems with GRASS python code in osgeo4w

[Moskovitz, Bob]

I see that c:\osgeo4w\apps\grass\grass-6.4.0svn\scripts is already in the PATH 
and "v.report --help" gives me the command usage.  

I do notice that there is no prompting of options or gui forms when I execute 
any grass command without arguments.  I've noticed this behavior on several 
machines.

> -Original Message-
> From: Hamish [mailto:hamis...@yahoo.com]
> Sent: Monday, June 22, 2009 10:42 AM
> To: Glynn Clements; Moskovitz, Bob
> Cc: Grass-User (E-mail)
> Subject: RE: [GRASS-user] Problems with GRASS python code in osgeo4w
> 
> 
> 
> Bob Moskovitz wrote:
> > So, I still don't know why the grass.parser() is not
> > working for me.  I even modified my program to use pdb
> > to see what is happening.  The debugger is aborted
> > right after I reach os.execvp("g.parser.exe", [name] + argv)
> > in grass.py.  But if I use "python
> > c:\osgeo4w\apps\grass\grass-6.4.0svn\scripts\m.dipslope.py --help",
> > I get the expected results.
> 
> 
> try putting 'c:\osgeo4w\apps\grass\grass-6.4.0svn\scripts' in 
> your %PATH%.
> when g.parser is complete it reruns the script but only by 
> its filename.
> If the script isn't in the PATH* it fails at that point because it
> can't find it. The above seems like it should be in the 
> path.. do other
> "official" scripts in that directory work?
> 
> [*] (the current dir is typically also checked in MS Windows 
> although it
> typically is not in UNIX)
> 
> see also GRASS_ADDONS_PATH enviro variable.
> 
> for grass 7 Glynn may have just changed this a couple of days ago??
> 
> 
> 
> Hamish
> 
> 
> 
>   
> 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Problems with GRASS python code in osgeo4w

[Moskovitz, Bob]

There is no etc/python/grass/script/core.py.  Just etc\python\grass.py.  Here 
is a part of a pdb debug session:

> c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(208)parser()
-> if sys.platform == "win32":
(Pdb)
> c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(209)parser()
-> try:
(Pdb)
> c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(210)parser()
-> os.execvp("g.parser.exe", [name] + argv)
(Pdb) p [name] + argv
['c:\\osgeo4w\\apps\\grass\\grass-6.4.0svn\\scripts\\m.dipslope.py', 'c:\\osgeo
4w\\apps\\grass\\grass-6.4.0svn\\scripts\\m.dipslope.py']
(Pdb)

> -Original Message-
> From: Glynn Clements [mailto:gl...@gclements.plus.com]
> Sent: Monday, June 22, 2009 3:30 PM
> To: Moskovitz, Bob
> Cc: Grass-User (E-mail)
> Subject: RE: [GRASS-user] Problems with GRASS python code in osgeo4w
> 
> 
> 
> Moskovitz, Bob wrote:
> 
> > Thanks for the info.  I didn't know about the assoc and 
> ftype.  Looking
> > at the output of these two commands, I see that my system is set up
> > correctly.  
> > 
> > So, I still don't know why the grass.parser() is not 
> working for me.  I
> > even modified my program to use pdb to see what is happening.  The
> > debugger is aborted right after I reach 
> os.execvp("g.parser.exe", [name]
> > + argv) in grass.py.  But if I use "python
> > c:\osgeo4w\apps\grass\grass-6.4.0svn\scripts\m.dipslope.py 
> --help", I
> > get the expected results.
> 
> Does it work if you run the script using its full path, i.e.:
> 
>  c:\osgeo4w\apps\grass\grass-6.4.0svn\scripts\m.dipslope.py --help
> 
> ?
> 
> You could try modifying etc/python/grass/script/core.py to print the
> arguments which are passed to g.parser.
> 
> The code attempts to determine the script's full pathname, as g.parser
> needs this so that it can open the script to read the #% comments
> (PATH won't help here).
> 
> It appears that this part is where the problem lies. If the script is
> run with --help or --ui, g.parser doesn't get around to re-invoking
> the script, so that isn't an issue.
> 
> -- 
> Glynn Clements 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Problems with GRASS python code in osgeo4w

[Moskovitz, Bob]

I also wanted to add that I did run the program like this:  
c:\osgeo4w\apps\grass\grass-6.4.0svn\scripts\m.dipslope.py --help

> -Original Message-
> From: grass-user-boun...@lists.osgeo.org
> [mailto:grass-user-boun...@lists.osgeo.org]on Behalf Of Moskovitz, Bob
> Sent: Monday, June 22, 2009 4:07 PM
> To: Glynn Clements
> Cc: Grass-User (E-mail)
> Subject: RE: [GRASS-user] Problems with GRASS python code in osgeo4w
> 
> 
> There is no etc/python/grass/script/core.py.  Just 
> etc\python\grass.py.  Here is a part of a pdb debug session:
> 
> > c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(208)parser()
> -> if sys.platform == "win32":
> (Pdb)
> > c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(209)parser()
> -> try:
> (Pdb)
> > c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(210)parser()
> -> os.execvp("g.parser.exe", [name] + argv)
> (Pdb) p [name] + argv
> ['c:\\osgeo4w\\apps\\grass\\grass-6.4.0svn\\scripts\\m.dipslop
> e.py', 'c:\\osgeo
> 4w\\apps\\grass\\grass-6.4.0svn\\scripts\\m.dipslope.py']
> (Pdb)
> 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Problems with GRASS python code in osgeo4w

[Moskovitz, Bob]

I debugged a bit deeper with pdb and found this:
(Pdb)
> c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(210)parser()
-> os.execvp("g.parser.exe", [name] + argv)
(Pdb) s
 snip ===
> c:\warmerda\release\apps\python25\lib\os.py(389)_execvpe()
-> func(fullname, *argrest)
(Pdb) p fullname
'c:\\osgeo4w\\bin\\g.parser.exe'
(Pdb) PATH
['c:\\osgeo4w\\bin', 'C:\\OSGeo4W/apps/grass/grass-6.4.0svn\\bin', 'C:\\OSGeo4W/
apps/grass/grass-6.4.0svn\\lib', 'C:\\OSGeo4W\\bin', 'C:\\Perl\\bin\\', 'C:\\WIN
DOWS\\system32', 'C:\\WINDOWS', 'C:\\WINDOWS\\System32\\Wbem', 'C:\\Perl\\bin\\'
,'C:\\WINDOWS\\system32', 'C:\\WINDOWS', 'C:\\WINDOWS\\System32\\Wbem']

There is no c:\osgeo4w\bin\g.parser.exe but there is a 
C:\OSGeo4W\apps\grass\grass-6.4.0svn\bin\g.parser.exe.

> -Original Message-
> From: Moskovitz, Bob 
> Sent: Monday, June 22, 2009 4:17 PM
> To: Moskovitz, Bob; Glynn Clements
> Cc: Grass-User (E-mail)
> Subject: RE: [GRASS-user] Problems with GRASS python code in osgeo4w
> 
> 
> I also wanted to add that I did run the program like this:  
> c:\osgeo4w\apps\grass\grass-6.4.0svn\scripts\m.dipslope.py --help
> 
> > -Original Message-
> > From: grass-user-boun...@lists.osgeo.org
> > [mailto:grass-user-boun...@lists.osgeo.org]on Behalf Of 
> Moskovitz, Bob
> > Sent: Monday, June 22, 2009 4:07 PM
> > To: Glynn Clements
> > Cc: Grass-User (E-mail)
> > Subject: RE: [GRASS-user] Problems with GRASS python code in osgeo4w
> > 
> > 
> > There is no etc/python/grass/script/core.py.  Just 
> > etc\python\grass.py.  Here is a part of a pdb debug session:
> > 
> > > 
> c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(208)parser()
> > -> if sys.platform == "win32":
> > (Pdb)
> > > 
> c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(209)parser()
> > -> try:
> > (Pdb)
> > > 
> c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(210)parser()
> > -> os.execvp("g.parser.exe", [name] + argv)
> > (Pdb) p [name] + argv
> > ['c:\\osgeo4w\\apps\\grass\\grass-6.4.0svn\\scripts\\m.dipslop
> > e.py', 'c:\\osgeo
> > 4w\\apps\\grass\\grass-6.4.0svn\\scripts\\m.dipslope.py']
> > (Pdb)
> > 
> > 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Problems with GRASS python code in osgeo4w

[Moskovitz, Bob]

I noticed that most of the *.bat files have forward slashes.  For example:

C:\OSGeo4W\apps\grass\bin\grass64.bat:  set 
WINGISBASE=/c/OSGeo4W/apps/grass/grass-6.4.0RC4
C:\OSGeo4W\apps\grass\grass-6.4.0svn\bin\d.correlate.bat:   @"%GRASS_SH%" 
-c '"%GISBASE%/scripts/d.correlate" %*'
C:\OSGeo4W\apps\grass\grass-6.4.0svn\etc\Init.bat:  if 
"%GRASS_PROJSHARE%"=="" set GRASS_PROJSHARE=/c/OSGeo4W/share/proj
C:\OSGeo4W\bin\grass64.bat  set 
WINGISBASE=%OSGEO4W_ROOT%/apps/grass/grass-6.4.0svn

Why is there 2 grass64.bat?


Bob Moskovitz
Research Analyst I
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.



> -Original Message-
> From: Glynn Clements [mailto:gl...@gclements.plus.com]
> Sent: Monday, June 22, 2009 7:22 PM
> To: Moskovitz, Bob
> Cc: Grass-User (E-mail)
> Subject: RE: [GRASS-user] Problems with GRASS python code in osgeo4w
> 
> 
> 
> Moskovitz, Bob wrote:
> 
> > There is no etc/python/grass/script/core.py. Just
> > etc\python\grass.py.
> 
> That changed in the last month or so, but I don't think that it has
> any bearing on this.
> 
> > Here is a part of a pdb debug session:
> > 
> > > 
> c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(208)parser()
> 
> I don't like the look of the "grass-~1.0sv" bit.
> 
> > -> if sys.platform == "win32":
> > (Pdb)
> > > 
> c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(209)parser()
> > -> try:
> > (Pdb)
> > > 
> c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(210)parser()
> > -> os.execvp("g.parser.exe", [name] + argv)
> > (Pdb) p [name] + argv
> > 
> ['c:\\osgeo4w\\apps\\grass\\grass-6.4.0svn\\scripts\\m.dipslop
> e.py', 'c:\\osgeo
> > 4w\\apps\\grass\\grass-6.4.0svn\\scripts\\m.dipslope.py']
> > (Pdb)
> 
> That much looks okay.
> 
> > I debugged a bit deeper with pdb and found this:
> > (Pdb)
> > > 
> c:\osgeo4w\apps\grass\grass-~1.0sv\etc\python\grass.py(210)parser()
> > -> os.execvp("g.parser.exe", [name] + argv)
> > (Pdb) s
> >  snip 
> ===
> > > c:\warmerda\release\apps\python25\lib\os.py(389)_execvpe()
> > -> func(fullname, *argrest)
> > (Pdb) p fullname
> > 'c:\\osgeo4w\\bin\\g.parser.exe'
> > (Pdb) PATH
> > ['c:\\osgeo4w\\bin', 
> 'C:\\OSGeo4W/apps/grass/grass-6.4.0svn\\bin', 'C:\\OSGeo4W/
> > apps/grass/grass-6.4.0svn\\lib', 'C:\\OSGeo4W\\bin', 
> 'C:\\Perl\\bin\\', 'C:\\WIN
> > DOWS\\system32', 'C:\\WINDOWS', 
> 'C:\\WINDOWS\\System32\\Wbem', 'C:\\Perl\\bin\\'
> > ,'C:\\WINDOWS\\system32', 'C:\\WINDOWS', 
> 'C:\\WINDOWS\\System32\\Wbem']
> > 
> > There is no c:\osgeo4w\bin\g.parser.exe but there is a 
> C:\OSGeo4W\apps\grass\grass-6.4.0svn\bin\g.parser.exe.
> 
> This is a red herring. os.execvp() iterates through PATH, attempting
> each directory in turn until it succeeds.
> 
> However, the fact that the GRASS bin and lib directories have forward
> slashes might be a problem. Try correcting PATH from the command line
> first.
> 
> -- 
> Glynn Clements 
> 
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Problems with GRASS python code in osgeo4w

[Moskovitz, Bob]

Glynn Clements wrote:

> 1. Ensure that .py files are configured to be run with the Python
> interpreter. You should be able to double-click on a .py file in
> Explorer and have it run. This should be done by the Python installer,
> but I don't know if the OSGeo4W installer does this.
> 
> If it isn't, you can use e.g.:
> 
>   assoc .py=python.file
>   ftype python.file="C:\Program Files\Python25\python.exe" "%1"
> 
> [Change the pathname to wherever Python is installed.]

I just wanted to update you on my progress with this problem.  I still can't 
get my script to work, but I believe I know the cause of the problem.  As you 
can see below I've included an email you sent related to setting up python file 
association.  

BTW, I did file an osgeo4w ticket:  http://trac.osgeo.org/osgeo4w/ticket/95

After looking at the ftype help, it looks like I should have issued the 
command: 
   ftype python.file="C:\Program Files\Python25\python.exe" "%1" %*


I've done that, but I still have the same problem.  I created a simple test 
program that just prints out the scripts arguments:

# printargs.py 
import sys

for arg in sys.argv:
   print arg
###

And here is grass session that illustration the problem:

GRASS 6.4.0svn (Morgan_Hill)> assoc .py
.py=py_auto_file

GRASS 6.4.0svn (Morgan_Hill)> ftype py_auto_file
py_auto_file="C:\OSGeo4W\bin\python.exe" "%1" %*

GRASS 6.4.0svn (Morgan_Hill)> C:\Temp\printargs.py this is a test
C:\Temp\printargs.py

GRASS 6.4.0svn (Morgan_Hill)> python C:\Temp\printargs.py this is a test
C:\Temp\printargs.py
this
is
a
test

Bob Moskovitz
Research Analyst I
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] Problems with GRASS python code in osgeo4w

[Moskovitz, Bob]

See comments below.

Glynn Clements wrote:

> > GRASS 6.4.0svn (Morgan_Hill)> assoc .py
> > .py=py_auto_file
> > 
> > GRASS 6.4.0svn (Morgan_Hill)> ftype py_auto_file
> > py_auto_file="C:\OSGeo4W\bin\python.exe" "%1" %*
> > 
> > GRASS 6.4.0svn (Morgan_Hill)> C:\Temp\printargs.py this is a test
> > C:\Temp\printargs.py
> > 
> > GRASS 6.4.0svn (Morgan_Hill)> python C:\Temp\printargs.py 
> this is a test
> > C:\Temp\printargs.py
> > this
> > is
> > a
> > test
> 
> Okay, so Python is working fine.
> 

If python is working, can you explain why my first example ( 
C:\Temp\printargs.py this is a test ) only prints out sys.argv[0] ?  
I believe that my problem has nothing to do with g.parser.  I believe that the 
arguments to a python script are being stripped 
when you rely on ftype.

> There have been some bugs fixed recently regarding g.parser on
> Windows.
> 

Based on http://trac.osgeo.org/osgeo4w/wiki/pkg-grass, I see that OSGEO4W's 
GRASS (which is about 4 months old) probably don't 
have many of these fixes.  So, I would like to make a request to the GRASS Devs 
to make RC5 available to OSGEO4W.  It's a shame
that the OSGEO4W community can't be involved in squashing the last few bugs :-)

> The first thing to check with a GRASS script is that "script --help"
> generates the help message. If it doesn't, that indicates that either
> g.parser isn't being called or it isn't being passed the script's
> filename correctly.
> 
> If --help works but invoking the script with other arguments doesn't
> work, that indicates that g.parser is failing to re-invoke the script.
> 

script --help   
No output
c:\full\path\to\script --help   
No output
python c:\full\path\to\script --helpHelp 
message

script --ui 
No output
c:\full\path\to\script --ui 
No output
python c:\full\path\to\script --ui  
No output

script --interface-description  No 
output
c:\full\path\to\script --interface-description  No 
output
python c:\full\path\to\script --interface-description   XML 
output

> If the script requires arguments, running the script without arguments
> should bring up the GUI. Some bugs with this were fixed quite
> recently, but I don't know if others remain.

Btw, most (if not all) of the compiled grass commands do not bring up the GUI 
when you run the program without arguments.

> 
> -- 
> Glynn Clements 
> 

Bob Moskovitz
Research Analyst I
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.

___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

RE: [GRASS-user] v.in.e00 vs. v.in.ogr

[Moskovitz, Bob]

I've noticed that using avcimport (which converts .e00 files to binary 
coverages) greatly speeds up import of .e00 files when you use ogr2ogr.  So 
v.in.ogr probably uses avcimport for this reason.

Bob Moskovitz
Research Analyst I
Seismic Hazard Evaluation Project
California Geological Survey
http://gmw.consrv.ca.gov/shmp

CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.



> -Original Message-
> From: grass-user-boun...@lists.osgeo.org
> [mailto:grass-user-boun...@lists.osgeo.org]on Behalf Of Rich Shepard
> Sent: Tuesday, November 17, 2009 4:51 PM
> To: grass-us...@lists.osgeo.org
> Subject: [GRASS-user] v.in.e00 vs. v.in.ogr
> 
> 
>Looking for the v.in.ogr man page I was reminded that 
> there is a v.in.e00.
> Should I use the latter rather than the former?
> 
>Tried running v.in.e00. After entering the .e00 file name 
> and specifying
> 'area' for the type I and saw this fatal error message:
> 
> ERROR: 'avcimport' program not found. Install it first.
> 
>The source file is an ASCII .e00, not compressed and not needing
> comversion from ARC/Info binary to exchange format.
> 
>What have I missed here?
> 
> Rich
> ___
> 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

Re: [GRASS-user] Geomorphons and extension manager

[Moskovitz, Bob@DOC]

When I tried to compile r.geomorphon, I got the errors you see below.  I just 
tried g.extensions with similar results.


$ make MODULE_TOPDIR=/home/bob/dev/gis/grass70

test -d OBJ.x86_64-unknown-linux-gnu || mkdir -p OBJ.x86_64-unknown-linux-gnu

gcc  -g -Wall -Werror-implicit-function-declaration -fno-common -Wextra 
-Wunused   -I/home/bob/dev/gis/grass70/dist.x86_64-unknown-linux-gnu/include 
-I/home/bob/dev/gis/grass70/dist.x86_64-unknown-linux-gnu/include  
-D_FILE_OFFSET_BITS=64  -DPACKAGE=\""grassmods"\"   
-I/home/bob/dev/gis/grass70/dist.x86_64-unknown-linux-gnu/include 
-I/home/bob/dev/gis/grass70/dist.x86_64-unknown-linux-gnu/include -o 
OBJ.x86_64-unknown-linux-gnu/geom.o -c geom.c

geom.c: In function ‘ternary_rotate’:

geom.c:21:5: warning: unused variable ‘res’ [-Wunused-variable]

geom.c: In function ‘radial2cartesian’:

geom.c:188:21: warning: ordered comparison of pointer with integer zero 
[-Wextra]

geom.c: In function ‘extends’:

geom.c:198:37: warning: unused parameter ‘pattern_size’ [-Wunused-parameter]

geom.c: At top level:

geom.c:2:15: warning: ‘dirs’ defined but not used [-Wunused-variable]

gcc  -g -Wall -Werror-implicit-function-declaration -fno-common -Wextra 
-Wunused   -I/home/bob/dev/gis/grass70/dist.x86_64-unknown-linux-gnu/include 
-I/home/bob/dev/gis/grass70/dist.x86_64-unknown-linux-gnu/include  
-D_FILE_OFFSET_BITS=64  -DPACKAGE=\""grassmods"\"   
-I/home/bob/dev/gis/grass70/dist.x86_64-unknown-linux-gnu/include 
-I/home/bob/dev/gis/grass70/dist.x86_64-unknown-linux-gnu/include -o 
OBJ.x86_64-unknown-linux-gnu/main.o -c main.c

main.c: In function ‘main’:

main.c:177:2: warning: format not a string literal and no format arguments 
[-Wformat-security]

main.c:198:2: warning: format not a string literal and no format arguments 
[-Wformat-security]

main.c:207:2: warning: format not a string literal and no format arguments 
[-Wformat-security]

main.c:213:2: warning: format not a string literal and no format arguments 
[-Wformat-security]

main.c:220:3: warning: format not a string literal and no format arguments 
[-Wformat-security]

main.c:234:3: warning: format not a string literal and no format arguments 
[-Wformat-security]

main.c:310:8: warning: format not a string literal and no format arguments 
[-Wformat-security]

main.c:366:5: error: implicit declaration of function ‘radial2cartesian’ 
[-Werror=implicit-function-declaration]

main.c:401:2: error: implicit declaration of function ‘write_contrast_colors’ 
[-Werror=implicit-function-declaration]

main.c:268:20: warning: unused variable ‘formC’ [-Wunused-variable]

main.c:268:13: warning: unused variable ‘formB’ [-Wunused-variable]

main.c:268:6: warning: unused variable ‘formA’ [-Wunused-variable]

main.c:76:22: warning: unused variable ‘radius’ [-Wunused-variable]

main.c:74:11: warning: unused variable ‘n’ [-Wunused-variable]

main.c:74:8: warning: unused variable ‘j’ [-Wunused-variable]

main.c:463:1: warning: control reaches end of non-void function [-Wreturn-type]

cc1: some warnings being treated as errors

make: *** [OBJ.x86_64-unknown-linux-gnu/main.o] Error 1



From: grass-user-boun...@lists.osgeo.org 
[mailto:grass-user-boun...@lists.osgeo.org] On Behalf Of Robert Andrew Brown
Sent: Tuesday, February 26, 2013 9:21 AM
To: grass-user@lists.osgeo.org
Cc: Phillip Owens
Subject: [GRASS-user] Geomorphons and extension manager

To whom it may concern,

I am trying to use extension manager in grass to install geomrphons. However, I 
am having trouble with the extension manager(error below). I emailed the 
geomorphon developer Jarek Jasiewicz(See Below) and he suggest I contact the 
GRASS-Users forum.

Please advise,

-Robert Brown



From: JASIEWICZ Jarosław mailto:jar...@amu.edu.pl>>
Subject: Re: Geomorphons
Date: February 26, 2013 11:56:20 AM EST
To: brown...@purdue.edu

On 02/26/2013 05:21 PM, brown...@purdue.edu wrote:

Dear Sir,

I use grass on my Mac OS X 10.8.2. I have been trying to compile grass to 
install the add-on Geomorphons. However, it is overwhelming to compile grass 
from source code as I have very little experience compiling programs. When I 
saw that your Geomorphons algorithm was part of the grass 7.0 add-ons I was 
pleased and tried to use the extension manager to add it to my grass add-on 
repository, however, I get an several errors which I have posted below. Any 
Suggestions?
#
g.extension extension=r.geomorphon/ operation=add
Fetching  from GRASS-Addons SVN (be patient)...
Compiling...
In file included from geom.c:1:
local_proto.h:1:20: error: stdlib.h: No such file or
directory
local_proto.h:2:20: error: string.h: No such file or
directory
local_proto.h:3:18: error: math.h: No such file or directory
local_proto.h:4:27: error: grass/glocale.h: No such file or
directory
local_proto.h:5:23: error: grass/gis.h: No such file or
directory
local_proto.h:6:26: error: grass/raster.h: No such file or
directory
In fil

[GRASS-user] Interpolate grid(s) based on vector.

[Moskovitz, Bob@DOC]

Hello Grass List,

I have a set of rasters (ground acceleration, modal mag., and distance) that I 
want to interpolate based on the location of my point features (borehole 
location) and update my point feature with the values from this interpolation.  
I see that v.what.rast might help, but I'm not sure what kind of interpolation 
it uses (probably nearest neighbor?).  I also see that v.sample can be used, 
but it does not work with geographic coordinate system locations.  So, I was 
hoping that the list might have other ideas in solving my problem.

Thanks in advance!

Bob


Robert Moskovitz
California Geological Survey
Seismic Hazards Zonation Program


CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

[GRASS-user] Set rasters with inf values to nodata

[Moskovitz, Bob@DOC]

Hello all,

I have a set of geotiffs that have inf (infinity) values that I would like to 
set to nodata.  How can I do that in GRASS or gdal?

Thanks!

Bob


Robert Moskovitz
California Geological Survey
Seismic Hazards Zonation Program


CONFIDENTIALITY NOTICE: This communication is intended only for the use of the 
individual or entity to which it is addressed. This message contains 
information from the State of California, California Geological Survey, which 
may be privileged, confidential and exempt from disclosure under applicable 
law, including the Electronic Communications Privacy Act. If the reader of this 
communication is not the intended recipient, you are hereby notified that any 
dissemination, distribution, or copying of this communication is strictly 
prohibited.


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Set rasters with inf values to nodata

[Moskovitz, Bob@DOC]

Thanks for looking into this.  

Here it is:
R:\SpecialProjects\tsunami work>gdalinfo --version
GDAL 1.10.1, released 2013/08/26
R:\SpecialProjects\tsunami work\>gdalinfo -mm 
NETCDF:pola_cascadiad_sift2.nc:max_speed
Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
Driver: netCDF/Network Common Data Format
Files: pola_cascadiad_sift2.nc
Size is 1810, 1100
Coordinate System is `'
Origin = (-118.30250118700,33.7766194)
Pixel Size = (0.92592592580,-0.92592592587)
Metadata:
  lat#long_name=latituded
  lat#units=degrees_north
  lon#long_name=longitude
  lon#units=degrees_east
  max_speed#long_name=Maximum Current Speed
  max_speed#missing_value=-1e+034
  max_speed#units=centimeters/second
Corner Coordinates:
Upper Left  (-118.3025000,  33.780)
Lower Left  (-118.3025000,  33.6781481)
Upper Right (-118.1349074,  33.780)
Lower Right (-118.1349074,  33.6781481)
Center  (-118.2187037,  33.7290741)
Band 1 Block=1810x1 Type=Float64, ColorInterp=Undefined
Computed Min/Max=0.000,1.#IO
  NoData Value=-9.9995e+033
  Metadata:
long_name=Maximum Current Speed
missing_value=-1e+034
NETCDF_VARNAME=max_speed
units=centimeters/second

-Original Message-
From: neteler.os...@gmail.com [mailto:neteler.os...@gmail.com] On Behalf Of 
Markus Neteler
Sent: Sunday, March 30, 2014 2:39 PM
To: Moskovitz, Bob@DOC
Cc: grass-user@lists.osgeo.org
Subject: Re: [GRASS-user] Set rasters with inf values to nodata

On Wed, Mar 19, 2014 at 8:18 PM, Moskovitz, Bob@DOC 
 wrote:
> Hello all,
>
> I have a set of geotiffs that have inf (infinity) values that I would 
> like to set to nodata.  How can I do that in GRASS or gdal?

what does

gdalinfo -mm yourrastermap.tif

report as Min/Max?

Markus


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Re: [GRASS-user] Set rasters with inf values to nodata

[Moskovitz, Bob@DOC]

This is the first time I've seen such Max values too.  But I have been able to 
solve my problem by converting these files to Arc/info ASCII Grids and replace 
all "1.#INF" to my nodata value.  I'll see if anyone on the GDAL ML can help me.

-Bob

-Original Message-
From: neteler.os...@gmail.com [mailto:neteler.os...@gmail.com] On Behalf Of 
Markus Neteler
Sent: Tuesday, April 01, 2014 12:39 PM
To: Moskovitz, Bob@DOC
Cc: grass-user@lists.osgeo.org
Subject: Re: [GRASS-user] Set rasters with inf values to nodata

On Tue, Apr 1, 2014 at 6:43 PM, Moskovitz, Bob@DOC 
 wrote:
> Thanks for looking into this.
>
> Here it is:
> R:\SpecialProjects\tsunami work>gdalinfo --version GDAL 1.10.1, 
> released 2013/08/26 R:\SpecialProjects\tsunami work\>gdalinfo -mm 
> NETCDF:pola_cascadiad_sift2.nc:max_speed
> Warning 1: No UNIDATA NC_GLOBAL:Conventions attribute
> Driver: netCDF/Network Common Data Format
> Files: pola_cascadiad_sift2.nc
...
> Band 1 Block=1810x1 Type=Float64, ColorInterp=Undefined
> Computed Min/Max=0.000,1.#IO

Honestly, I have never seen such Max value.

Maybe a GDAL or dataset issue?

>   NoData Value=-9.9995e+033
>   Metadata:
> long_name=Maximum Current Speed
> missing_value=-1e+034
> NETCDF_VARNAME=max_speed
> units=centimeters/second


Markus


___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user