Re: [GRASS-user] get values for each point within lines
Milton Cezar Ribeiro wrote: Dear Grass gurus I have a vector map with lines, and I need to get raster values each 5 meters. So I used v.split to spring my lines in segments of 5 meters, and after that I used v.to.points to generate points on each node of my line segments. Now I am trying to overlap my points against a elevation map, but when I run v.what.rast vect=lines_split_05m_points rast=dem col=elevation the numbers of records on my v.db.select is the same of my lines instead of my points. The question is how to get the elevation values and keep the other attributes from my source line. You need to add unique categories for the points in a new layer with v.category option=add. Then transfer the original line attributes from the original layer to the new layer with v.to.db, finally run v.what.rast. HTH, Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Parallel Grass job
I'm new to GRASS. I want to parallelize a program (Grass script), but keep getting the following error. copmpi.c:7:17: error: mpi.h: No such file or directory copmpi.c: In function ‘main’: copmpi.c:177: error: redeclaration of ‘row’ with no linkage copmpi.c:45: error: previous declaration of ‘row’ was here copmpi.c:187: error: ‘MPI_Status’ undeclared (first use in this function) copmpi.c:187: error: (Each undeclared identifier is reported only once copmpi.c:187: error: for each function it appears in.) copmpi.c:187: error: expected ‘;’ before ‘status’ copmpi.c:189: error: ‘MPI_COMM_WORLD’ undeclared (first use in this function) copmpi.c:196: error: ‘MPI_INT’ undeclared (first use in this function) copmpi.c:211: warning: assignment makes integer from pointer without a cast copmpi.c:212: warning: assignment makes integer from pointer without a cast copmpi.c:218: error: ‘MPI_DOUBLE’ undeclared (first use in this function) copmpi.c:224: error: ‘status’ undeclared (first use in this function) copmpi.c:242: error: expected expression before ‘inrast_redchan’ copmpi.c:265: error: subscripted value is neither array nor pointer make: *** [OBJ.x86_64-unknown-linux-gnu/copmpi.o] Error 1 tried other scripts but again myMPIexample.c:1:69: error: mpi.h: No such file or directory But MPI alone (mpi.h) works. and gis.h works alone as well too. That is no problem with libraries. how can I solve this problem??? and what to do? help me please. Thank you !!! ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] copy elements from user-mapset to PERMANENT
Hi, I tried to copy a vector from my user mapset into the PERMANENT mapset of the same location with g.copy. g.copy vect=Treene_GK@Treene,Treene_simple@PERMANENT region=region_treene@Treene,region_Treene@PERMANENT Illegal filename. Character @ not allowed. Treene_simple@PERMANENT is an illegal file name Illegal filename. Character @ not allowed. region_Treene@PERMANENT is an illegal file name (Wed Feb 15 10:57:31 2012) Command finished (0 sec) Maybe that is not the correct command..is there a simple way to perform that copy? Or do I need to export/import the files? So far as I understand I can read the PERMAMENT maps from all mapsets in that location...Thats why I'd like to have certain basic maps in there Best regards... Johannes -- NEU: FreePhone 3-fach-Flat mit kostenlosem Smartphone! Jetzt informieren: http://www.gmx.net/de/go/freephone/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] g.mlist with logical operators?
Thank you all, I'll have a look into regular expressions... sound very promissing to me. /johannes Original-Nachricht Datum: Tue, 14 Feb 2012 20:01:44 + Von: Glynn Clements gl...@gclements.plus.com An: Johannes Radinger jradin...@gmx.at CC: GraSS user list grass-user@lists.osgeo.org Betreff: Re: [GRASS-user] g.mlist with logical operators? Johannes Radinger wrote: I'd like to list (using g.list) all raster maps that start with rast_* but not those which start with rast_A_* OR rast_B_* OR rast_C_*. I tried to use the OR operator and | but without success so far. Are such logical operators generally applicable with g.mlist? No. You may be able to achieve what you want with regexps. Otherwise, you can process the output from g.mlist with additional programs. Writing a regexp which matches the union of multiple patterns (A OR B OR C) is straightforward enough; just use (A|B|C) with extended regexps or '\(A\|B\|C\)' with basic regexps (care needs to be taken to ensure that backslashes are passed verbatim rather than being interpreted by the shell). There isn't an equivalent operator for intersection or difference, but you can achieve the same result by filtering the output with grep or grep -v respectively. -- Glynn Clements gl...@gclements.plus.com -- NEU: FreePhone 3-fach-Flat mit kostenlosem Smartphone! Jetzt informieren: http://www.gmx.net/de/go/freephone/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] copy elements from user-mapset to PERMANENT
Hi, 2012/2/15 Johannes Radinger jradin...@gmx.at: I tried to copy a vector from my user mapset into the PERMANENT mapset of the same location with g.copy. you have write access only for the current mapset. In other words, you need to switch to PERMANENT and call `g.copy` from the target mapset. g.copy vect=Treene_GK@Treene,Treene_simple@PERMANENT region=region_treene@Treene,region_Treene@PERMANENT Illegal filename. Character @ not allowed. Treene_simple@PERMANENT is an illegal file name Illegal filename. Character @ not allowed. region_Treene@PERMANENT is an illegal file name (Wed Feb 15 10:57:31 2012) Command finished (0 sec) g.mapset PERMANENT g.copy vect=Treene_GK@Treene,Treene_simple region=region_treene@Treene,region_Treene Maybe that is not the correct command..is there a simple way to perform that copy? Or do I need to export/import the files? Well, the error message is not clear enough... Martin -- Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] copy elements from user-mapset to PERMANENT
Hi, Hi, 2012/2/15 Johannes Radinger jradin...@gmx.at: I tried to copy a vector from my user mapset into the PERMANENT mapset of the same location with g.copy. you have write access only for the current mapset. In other words, you need to switch to PERMANENT and call `g.copy` from the target mapset. g.copy vect=Treene_GK@Treene,Treene_simple@PERMANENT region=region_treene@Treene,region_Treene@PERMANENT Illegal filename. Character @ not allowed. Treene_simple@PERMANENT is an illegal file name Illegal filename. Character @ not allowed. region_Treene@PERMANENT is an illegal file name (Wed Feb 15 10:57:31 2012) Command finished (0 sec) g.mapset PERMANENT g.copy vect=Treene_GK@Treene,Treene_simple region=region_treene@Treene,region_Treene Thank you... of course I have to change the mapset...don't know why I missed that step before...:( cheers, johannes Maybe that is not the correct command..is there a simple way to perform that copy? Or do I need to export/import the files? Well, the error message is not clear enough... Martin -- Martin Landa landa.martin gmail.com * http://geo.fsv.cvut.cz/~landa -- NEU: FreePhone 3-fach-Flat mit kostenlosem Smartphone! Jetzt informieren: http://www.gmx.net/de/go/freephone/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: grass and palsar processing
Thanks markus for this. I added an initial workflow for the ENVI tutorial I attended my goal is to port this to GRASS. Most of the process is already available in the pdf Alexander provided. However, I'm having difficulty importing images from the KC mosaic product [0]. In the training I attended the mosaic product was converted to GTiff using ASF MapReady [1]. According to ASF MapReady the format is CEOS L1. I tried using r.in.gdal and r.in.bin (based on the pdf) but I cannot import the KC mosaic [0] http://www.eorc.jaxa.jp/ALOS/en/kc_mosaic/kc_mosaic.htm [1] http://www.asf.alaska.edu/downloads/software_tools/mapready/version_history On Sat, Feb 11, 2012 at 7:54 AM, Markus Neteler nete...@osgeo.org wrote: Hi, I am almost there:http://grass.osgeo.org/wiki/Radar_data_processing only imagemagick currently fails to make the thumbnails after debian upgrade.will check asap and upload the missing files. On Fri, Feb 10, 2012 at 7:43 PM, Alexander Muriy amu...@gmail.com wrote: Also there is PDF http://gis-lab.info/forum/download/file.php?id=2370, copied from http://www.safe.iis.u-tokyo.ac.jp/, about SAR GRASS. thanks, added on the Wiki page. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- cheers, maning -- Freedom is still the most radical idea of all -N.Branden wiki: http://esambale.wikispaces.com/ blog: http://epsg4253.wordpress.com/ -- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] How to remove small vectors in a shapefile in GRASS 6.4.1
Hello, I want to remove small vector in a shp. I tried with v.generalize but I only managed to remove segment of a line. The image shows in black the vectors I want to clean. And in green is the result I obtained with v.generalize. Has someone an idea? Thanks Lucien http://osgeo-org.1560.n6.nabble.com/file/n4472852/courbes-niveau.png -- View this message in context: http://osgeo-org.1560.n6.nabble.com/How-to-remove-small-vectors-in-a-shapefile-in-GRASS-6-4-1-tp4472852p4472852.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
Re: [GRASS-user] How to remove small vectors in a shapefile in GRASS 6.4.1
On Wed, Feb 15, 2012 at 5:46 PM, lucien.blandenier lucien.blanden...@unine.ch wrote: Hello, I want to remove small vector in a shp. I tried with v.generalize but I only managed to remove segment of a line. The image shows in black the vectors I want to clean. And in green is the result I obtained with v.generalize. Has someone an idea? Try v.clean tool=rmarea threshold=XXX That will remove all areas smaller than XXX. HTH, Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: grass and palsar processing
On Sat, Feb 11, 2012 at 12:54 AM, Markus Neteler nete...@osgeo.org wrote: I am almost there:http://grass.osgeo.org/wiki/Radar_data_processing only imagemagick currently fails to make the thumbnails after debian upgrade.will check asap and upload the missing files. Fixed by Martin right now, the thumbnails are generated again in the Wiki. Cheers Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Not being able to produce decent-look color-composites
what I normally do is first run r.colors -e color=grey on each band so I get a equalized color composite. Then I run r.composite. So far I've got good results... On Mon, Feb 13, 2012 at 8:16 PM, Glynn Clements gl...@gclements.plus.com wrote: Luisa Peña wrote: I'have 3 layers of Landsat (Green, Red and Near-Infrared) and I'm trying to produce a false-color composite (r.composite( with those but... I'm getting very frustrated because I'm producing very bad look composites. My layers are reflectances from 0-1.8 or something. Shall I rescale them to 0-255? r.composite uses the colour table associated with each map to convert the map's value to an intensity. If 1.8 represents full intensity, then use e.g.: r.colors mapname rules=- EOF 0.0 black 1.8 white EOF If you want a non-linear relationship, you can either contruct a colour table using a piecewise-linear approximation, or apply the transformation to the data with r.mapcalc. -- Glynn Clements gl...@gclements.plus.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] bug in i.topo.corr
Dear list, I am using a current build of grass70 on Ubuntu 11.10. I installed the module 'i.topo.corr' using Trying to use it, I get the following error-message Anyone understands this? I saw that this error message appeared for other modules as well (e.g. http://osgeo-org.1560.n6.nabble.com/python-scripting-changing-region-settings-reiases-error-td3889577.html) and it was reported as a bug. Stefan -- View this message in context: http://osgeo-org.1560.n6.nabble.com/bug-in-i-topo-corr-tp4474772p4474772.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
Re: [GRASS-user] Not being able to produce decent-look color-composites
On Thu, Feb 16, 2012 at 12:54 AM, Daniel Victoria daniel.victo...@gmail.com wrote: what I normally do is first run r.colors -e color=grey on each band so I get a equalized color composite. Then I run r.composite. So far I've got good results... Added to http://grass.osgeo.org/wiki/LANDSAT#Natural_color_composites Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user