[GRASS-user] GRASS 7 on Ubuntu - was Re: New GRASS tutorials for ecologists
(please use an appropriate subject line, i.e. open a new thread for a new topic) On Fri, Nov 16, 2012 at 10:48 PM, Gianluca Massei agr.gianluca.mas...@gmail.com wrote: Hi, I'm trying to install some extension in grass 7 (compiled from svn grass_trunk code in ubuntu 12.10) but I get this message for all my tests: ERROR: MAPSET PERMANENT - permission denied Are you the owner of the mapset? make: *** [v.surf.icw.tmp.html] Errore 1 rm v.surf.icw.tmp.html ERROR: Compilation failed, sorry. Please check above error messages. I can't understand why the permission is denied in MAPSET PERMANENT: I've set all files and folders with chmod 777 * and my $GRASS_ADDON_BASE=/home/gianluca/.grass7/addons. Do you use the wxGUI integrated extension manager or g.extension? Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Problems in adding tables to a map
GRASS user, I'm having some problem in connecting tables to a map. This is the procedure (in order to simplify I used the same table but with different names) : g.copy vect=roads,roads1 g.copy vect=roads,roads2 The two tables are in the table list. db.tables -p v.db.connect map=roads1 table=roads2 layer=2 The connections seem ok: v.db.connect -p map=roads1 Vector map roads1@esercizi is connected by: layer 1 table roads1 in database /home/clerici/GRASSdata/spearfish60/esercizi/dbf/ through driver dbf with key cat layer 2 table roads2 in database /home/clerici/GRASSdata/spearfish60/esercizi/dbf/ through driver dbf with key cat But: d.mon x0 This displays the map with lines and labels d.vect map=roads1 disp=shape,attr attrcol=label layer=1 This doesn't produce any output d.erase d.vect map=roads1 disp=shape,attr attrcol=label layer=2 I tried a check by producing the ASCII file with: v.out.ascii in=roads1 form=standard but after each sequence of coordinates only layer 1 is reported ... L 6 1 590744.10526028 4926835.64296535 590742.5374415 4926835.30309294 590624.60065406 4926810.04019906 590503.0158192 4926784.15943798 590383.85077876 4926762.55284068 590253.15645239 4926731.77335036 1 4 L 2 1 591497.07137115 4926772.322352 591667.99469283 4926778.18725746 1 4 ... Another check with: v.extract in=roads1 out=roads_extr layer=2 list=1-3 produces an empty map. Did I miss something in the connection command? Thanks Aldo Clerici ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] adding a datum to the GRASS database
Hi Ali, On Fri, Nov 16, 2012 at 9:05 PM, Ali Swanson swanson@gmail.com wrote: Hi all - I'm brand new to GRASS (though have basic familiarity with ArcInfo). welcome :) I work in a datum that apparently isn't contained in the basic GRASS library -- from what I've uncovered searching this list the general web, I need to add my datum by going to the $GISBASE/etc/datum.table and appending this text: The Arc 1960 datum; source: EPSG database arc60 Arc 1960 clark80 dx=-160 dy=-6 dz=-302 Which EPSG code is that? EPSG 4210, this one? epsg_tr -pretty_wkt 4210 GEOGCS[Arc 1960, DATUM[Arc_1960, SPHEROID[Clarke 1880 (RGS),6378249.145,293.465, AUTHORITY[EPSG,7012]], TOWGS84[-160,-6,-302,0,0,0,0], AUTHORITY[EPSG,6210]], PRIMEM[Greenwich,0, AUTHORITY[EPSG,8901]], UNIT[degree,0.0174532925199433, AUTHORITY[EPSG,9122]], AUTHORITY[EPSG,4210]] It should be there. You can use the Location Wizard to select this EPSG code. BUT: I just tried and the datum is not recognised: g.proj -w GEOGCS[clark80, DATUM[unknown, SPHEROID[Clarke_1880,6378249.145,293.465]], PRIMEM[Greenwich,0], UNIT[degree,0.0174532925199433]] Do these tickets strike back? http://trac.osgeo.org/grass/ticket/1452 http://trac.osgeo.org/grass/ticket/1748 Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] error in g.extension
I'm sorry I've sent the message too many times in the list! %-( Are you the owner of the mapset? Yes, I'm owner of the mapset. With ls -l -a I've this: drwxrwxrwx 6 gianluca gianluca 4096 nov 6 21:58 UmbriaRoma40 I've used g.extension from command line. I've the same problem in grass 6.4 Thanks for the help 2012/11/17 Markus Neteler nete...@osgeo.org On Fri, Nov 16, 2012 at 10:53 PM, Gianluca Massei agr.gianluca.mas...@gmail.com wrote: Hi, I'm trying to install some extension in grass 7 (compiled from svn grass_trunk code in ubuntu 12.10) but I get this message for all my tests: ERROR: MAPSET PERMANENT - permission denied Are you the owner of the mapset? make: *** [v.surf.icw.tmp.html] Errore 1 rm v.surf.icw.tmp.html ERROR: Compilation failed, sorry. Please check above error messages. I can't understand why the permission is denied in MAPSET PERMANENT: I've set all files and folders with chmod 777 * and my $GRASS_ADDON_BASE=/home/gianluca/.grass7/addons. Do you use the wxGUI integrated extension manager or g.extension (as the subject states)? Markus -- Dott. Gianluca Massei, PhD ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] error in g.extension
If you only want to use this addon, then you just need to download the addon and run the command sudo make MODULE_TOPDIR=/usr/local/grass-7.0.svn/ in the addon folder. dehui -原始邮件- 发件人: Markus Neteler nete...@osgeo.org 发送时间: 2012年11月17日 星期六 收件人: Gianluca Massei agr.gianluca.mas...@gmail.com 抄送: GRASS user list grass-user@lists.osgeo.org 主题: Re: [GRASS-user] error in g.extension On Fri, Nov 16, 2012 at 10:53 PM, Gianluca Massei agr.gianluca.mas...@gmail.com wrote: Hi, I'm trying to install some extension in grass 7 (compiled from svn grass_trunk code in ubuntu 12.10) but I get this message for all my tests: ERROR: MAPSET PERMANENT - permission denied Are you the owner of the mapset? make: *** [v.surf.icw.tmp.html] Errore 1 rm v.surf.icw.tmp.html ERROR: Compilation failed, sorry. Please check above error messages. I can't understand why the permission is denied in MAPSET PERMANENT: I've set all files and folders with chmod 777 * and my $GRASS_ADDON_BASE=/home/gianluca/.grass7/addons. Do you use the wxGUI integrated extension manager or g.extension (as the subject states)? 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
Re: [GRASS-user] error in g.extension
If you only want to use this addon, then you just need to download the addon and run the command sudo make MODULE_TOPDIR=/usr/local/grass-7.0.svn/ in the addon folder. I haven't problem with addons installation from source, but I'd like try g.extension. Thanks ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Does v.kernel have to take 16+ hours?
I have a dataset of just over 700,000 incidents that happened in square-ish Texas county that's about 30 miles on each side. Here's the parameters reported by v.kernel as it's executing: *STDDEV: 1000.00* *RES: 111.419043 ROWS: 458 COLS: 447* * * *Writing output raster map using smooth parameter=1000.00.* * * *Normalising factor=6482635.018778.* I am running this on a Windows 7 x64 machine with 8 GB RAM and an Intel Core i7 Q720 1.6 GHz with 4 physical cores. I notice that it's not multithreaded, only using 1 core. It takes about 16 hours to complete. Is this correct? I'd like to use this on a dataset with closer to 5 million records, and I'm really concerned how long it may take. I tried the v.kernel again with about 10X more cells on each dimension, just to see what would happen. It ran for several hours, and based on the progress meter, I guess it would take less than a day to complete. Unfortunately, it crashed mid-way. See http://trac.osgeo.org/grass/ticket/1800. I posted my question about the 16+ hours at http://gis.stackexchange.com/questions/41058/how-do-i-compute-v-kernel-maps-in-less-than-16-hours/. Bill Huber, who si apparently knowledgeable about kernel density calculations in general, posted a response, and he felt like a kernel density map shouldn't take much time at all. But digging more deeply, turns out he had come up with a kernel density calculation method over a decade ago using Fourier transforms. See http://www.directionsmag.com/features/convolution/129753 and the next two articles linked to it (they are short articles). Apparently this transforms it from an O(n^2) problem to an O(n ln n) complexity problem. I inspected v.kernel's main.c ( http://trac.osgeo.org/grass/browser/grass/trunk/vector/v.kernel/main.c), and looks like v.kernel uses an output-centric method (using Bill's wording) of calculating the output, which seems like O(n^2) complexity. So I guess what I'm getting at is it appears to me that the algorithm behind GRASS GIS's v.kernel is straightforward but is a greedy algorithm ( http://en.wikipedia.org/wiki/Greedy_algorithm), which is fine, but it make take a while to execute. Is this true? Is there not spatial indexing I could add to the dataset? I've done various Google searches on that and can't come up with anything clear. Aren ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user