[GRASS-user] v.rast.stats
Hi When I'm running the command below I'm getting some errors, I pasted them below v.rast.stats -ec vector=catchments raster=rainfall colprefix=rain_cols percentile=90 DBMI-DBF driver error: SQL parser error: syntax error, unexpected NULL_VALUE processing 'NULL' in statement: UPDATE catchments SET rain_cols_=-NULL WHERE cat=1 Error in db_execute_immediate ERROR: Error while executing: 'UPDATE catchments SET rain_cols_=-NULL WHERE cat=1' -- Kind Regards TS Gumede Student Researcher Meraka Institute, CSIR office : 012 841 2606 mobile : 072 258 1650 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Projecting a GCS onto a computer screen
Hi, I was just wondering how the GRASS code displays coordinates in a geographic coordinate system on a computer screen. Do they get projected in some way? Or is there another trick to it? Thanks Hanlie ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Help with NVIZ crash
Good morning all, I'm looking for some help with NVIZ which keeps dropping out on my Ubuntu Hardy / GRASS 6.4 set up when loading a raster with the following message - Error in startup script: Togl: X server has no OpenGL GLX extension while executing togl $Nv_(TOP).canvas -rgba true -double true -depth true (procedure Nv_makeGUI line 56) invoked from within Nv_makeGUI $Nv_(AREA) (file /usr/lib/grass64/etc/nviz2.2/scripts/nviz2.2_script line 1052) Any suggestions welcome! Kind regards Rebecca ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] 3D peat modeling
Hello all, I'm in GRASS for a week and I try to model develope phases of peat from a drillhole data. I have seen previous posts about the same problem (geological) but I always get some errors repeating the steps. I have boundary of peat imported from Microstation with v.in.ogr and then I created 2D mask. Then I have drillhole data (5 vector point files) also imported from MicroStation with v.in.dxf Afterwards I made 5 rasters of peat (bottom-spag_rst_mask , gyttja-gytia..., peat1-mech1.., peat2-oles..., top-mech2) withv.surf.rst with the mask I mentioned above. I tried to have r.to.rast3 but I got this error http://monki.us/~kamil/grass/1.txt LINK . I added the layes from the bottom to top. I dont know maybe some 3D regional problem? Another approach was with r.to.rast3elev but I got an error to http://monki.us/~kamil/grass/2.txt LINK the sequence of adding layers from the bottom. My aim is to get something like this http://opengeology.pl/files2/Image/3d/brom_isosur.png AIM I mean each 3D layer(phase) with different colour but I have no idea how to do this. v.vol.rst ??? It would be nice to have shown drillholes (x,y) doesn't matter if I get the results in nviz or paraview, I have both. Just like I said before I'm beginer in grass and I realise I could have basic errors doing my job. I appreciate the respond as quick as you can. -- View this message in context: http://osgeo-org.1803224.n2.nabble.com/3D-peat-modeling-tp5138619p5138619.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] Problem importing netCDF
Hi,
I'm using Win XP and GRASS 6.4RC6 to import a TRMM netCDF file using
the command:
-
r.in.gdal -o input =c:/data/3B42.000202.0.6.nc output=test band=precipitation
-
and am getting the error:
-
G_set_window(): Illegal latitude for North
-
gdalinfo for the file outputs:
-
Driver: netCDF/Network Common Data Format
Files: C:\data\3B42.000202.0.6.nc
C:\data\3B42.000202.0.6.nc.aux.xml
Size is 512, 512
Coordinate System is `'
Metadata:
NC_GLOBAL#Conventions=CF-1.0
Subdatasets:
SUBDATASET_1_NAME=NETCDF:C:\data\3B42.000202.0.6.nc:precipitation
SUBDATASET_1_DESC=[400x1440] precipitation (32-bit floating-point)
SUBDATASET_2_NAME=NETCDF:C:\data\3B42.000202.0.6.nc:relativeError
SUBDATASET_2_DESC=[400x1440] relativeError (32-bit floating-point)
Corner Coordinates:
Upper Left (0.0,0.0)
Lower Left (0.0, 512.0)
Upper Right ( 512.0,0.0)
Lower Right ( 512.0, 512.0)
Center ( 256.0, 256.0)
-
In Panoply (http://www.giss.nasa.gov/tools/panoply/), the information
for the file displays as:
-
netcdf file:/C:/data/3B42.000202.0.6.nc {
dimensions:
lat = 400;
lon = 1440;
variables:
float precipitation(lat=400, lon=1440);
:long_name = precipitation;
:units = mm hr-1;
:_FillValue = -.9f; // float
float relativeError(lat=400, lon=1440);
:long_name = relativeError;
:units = mm hr-1;
:_FillValue = -.9f; // float
float lat(lat=400);
:standard_name = latitude;
:units = degrees_north;
:long_name = latitude;
:_CoordinateAxisType = Lat;
float lon(lon=1440);
:standard_name = longitude;
:units = degrees_east;
:long_name = longitude;
:_CoordinateAxisType = Lon;
:Conventions = CF-1.0;
}
-
In the Panoply array display, the coordinates range from 49.875° to
-49.875° (lattitude) and -179.875° to +179.875° longitude. These are
coordinates for the centres of 0.25° grid cells.
I have set my region according to the documentation for the dataset to:
-
projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 50N
south: 50S
west: 180W
east: 180E
nsres: 0:15
ewres: 0:15
rows: 400
cols: 1440
cells: 576000
-
Can someone perhaps spot my error?
Thanks
Hanlie
___
grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Help with NVIZ crash
Great thanks for the heads up Glenn :) From: Glynn Clements gl...@gclements.plus.com To: Rebecca Bennett rabenn...@ymail.com Cc: grass-user@lists.osgeo.org Sent: Fri, 4 June, 2010 10:47:52 Subject: Re: [GRASS-user] Help with NVIZ crash Rebecca Bennett wrote: I'm looking for some help with NVIZ which keeps dropping out on my Ubuntu Hardy / GRASS 6.4 set up when loading a raster with the following message - Error in startup script: Togl: X server has no OpenGL GLX extension What it says: your X server doesn't support OpenGL. You should ask on an Ubuntu forum for advice on fixing this. -- 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] use of v.surf.rst on big datasets with low memory
Hi, I'm trying to build a regional DEM based on a contour line map (3D vector) using v.surf.rst but the process fails with an memory allocation error when I set my region to a 50m resolution. As I'm using GRASS 6.4svn (OSGEO4W) on a Windows XP x86 with 2Go of RAM, I would like to know if it is gonna be a constant error due to those hardware/software specs. (Thu Jun 03 13:12:31 2010) v.surf.rst input=mnt_i2g...@mnt layer=0 elev=rst_mnt_i2g_r...@mnt tension=40 segmax=40 npmin=300 dmin=25 dmax=125 Pourcentage complet : Lignes de lecture de la carte vectorielle... certains points sont en dehors de la région -- ignorés... Reading nodes from vector map ... il y a des points en-dehors de la région 2D/3D spécifiée -- 2 points ignorés 28763365 points ignorés -- trop dense Processing all selected output files will require 26129572 bytes of disk space for temp files Pourcentage complet : la recherche de points pour l'interpolation prend trop de temps -- veuillez changer la régionpour la zone où se trouvent vos points. Les calculs se poursuivent... The number of points from vector map is 32808058 The number of points outside of 2D/3D region 2 The number of points being used is 4044691 [...] ERROR: G_malloc: unable to allocate 16992 bytes at segmen2d.c:196 (Thu Jun 03 16:12:58 2010) Command finished (10826 sec) Regards, MORREALE Jean Rof ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] gdal
Hi I've got a geotiff file which covers the area of South Africa, from that tif file I want to take a portion of Cape Flats (Western Cape) and run it on GRASS. Which command should I use to do this? -- Kind Regards TS Gumede Student Researcher Meraka Institute, CSIR office : 012 841 2606 mobile : 072 258 1650 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.rast.stats
On Fri, 4 Jun 2010, Sandile Gumede wrote: in statement: UPDATE catchments SET rain_cols_=-NULL WHERE cat=1 Do you really mean '-NULL'? I've never seen NULL negated before. Rich ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Problem with GRASS 6.4RC6 and wxPython 2.8.11 on CentOS 5.5 x86_64
Hello, I'm working on the RPM packaging of GRASS 6.4 for Enterprise Linux (http://wiki.osgeo.org/wiki/Enterprise_Linux_GIS). We rely on the EPEL repository for many libraries and since their upgrade of wxGTK / wxPython to 2.8.11 the wx UI of GRASS is broken. We get the stacktrace below: ... File /usr/lib64/grass-6.4.0RC6/etc/wxpython/gui_modules/wxgui_utils.py, line 73, in __init__ super(LayerTree, self).__init__(parent, id, pos, size, style=style, ctstyle=ctstyle) File /usr/lib64/python2.4/site-packages/wx-2.8-gtk2-unicode/wx/lib/mixins/treemixin.py, line 488, in __init__ super(DragAndDrop, self).__init__(*args, **kwargs) TypeError: __init__() got an unexpected keyword argument 'ctstyle' (complete stack at the end of the mail) This may also be due to the update of other component when updating to CentOS 5.5, but I rather suspect the wx update. Anyhow the wx UI was working on CentOS 5.4 with wxGTK 2.8.9 / wxPython 2.8.7 . Can you please help me analyze this issue further? Thanks in advance! Mathieu Welcome to GRASS 6.4.0RC6 (2010) GRASS homepage: http://grass.osgeo.org/ This version running thru: Bash Shell (/bin/bash) Help is available with the command: g.manual -i See the licence terms with: g.version -c If required, restart the GUI with: g.gui wxpython When ready to quit enter:exit GRASS 6.4.0RC6 (EasternBalkans2010):~ WARNING: Vector digitizer is not available (No module named grass6_wxvdigit). Note that the vector digitizer is currently not working under MS Windows (hopefully this will be fixed soon). Please keep an eye out for updated versions of GRASS. Traceback (most recent call last): File /usr/lib64/grass-6.4.0RC6/etc/wxpython/wxgui.py, line 1819, in ? sys.exit(main()) File /usr/lib64/grass-6.4.0RC6/etc/wxpython/wxgui.py, line 1812, in main app = GMApp(workspaceFile) File /usr/lib64/grass-6.4.0RC6/etc/wxpython/wxgui.py, line 1737, in __init__ wx.App.__init__(self, False) File /usr/lib64/python2.4/site-packages/wx-2.8-gtk2-unicode/wx/_core.py, line 7978, in __init__ self._BootstrapApp() File /usr/lib64/python2.4/site-packages/wx-2.8-gtk2-unicode/wx/_core.py, line 7552, in _BootstrapApp return _core_.PyApp__BootstrapApp(*args, **kwargs) File /usr/lib64/grass-6.4.0RC6/etc/wxpython/wxgui.py, line 1755, in OnInit workspace = self.workspaceFile) File /usr/lib64/grass-6.4.0RC6/etc/wxpython/wxgui.py, line 185, in __init__ self.NewDisplay(show=False) File /usr/lib64/grass-6.4.0RC6/etc/wxpython/wxgui.py, line 1236, in NewDisplay auimgr=self._auimgr, showMapDisplay=show) File /usr/lib64/grass-6.4.0RC6/etc/wxpython/gui_modules/wxgui_utils.py, line 73, in __init__ super(LayerTree, self).__init__(parent, id, pos, size, style=style, ctstyle=ctstyle) File /usr/lib64/python2.4/site-packages/wx-2.8-gtk2-unicode/wx/lib/mixins/treemixin.py, line 488, in __init__ super(DragAndDrop, self).__init__(*args, **kwargs) TypeError: __init__() got an unexpected keyword argument 'ctstyle' ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Importing Dems with r.in.xyz
Hi, I've been following the procedure at http://grass.osgeo.org/wiki/Import_XYZ to import DEMs into GRASS. It has worked for two DEMS, but I have a problem with the third one at the step where one verifies that the number of rows in the ASCII file corresponds to the number of cells in the enlarged region. At this point, g.region reports 1146474 cells in the region, while I have 1146370 lines of coordinates in my file. The output of g.region is: - projection: 99 (Transverse Mercator) zone: 0 datum: ** unknown (default: WGS84) ** ellipsoid: wgs84 north: -49312.49 south: -74587.5 west: -3015862.5 east: -2987512.5 nsres: 25.0989 ewres: 25 rows: 1011 cols: 1134 cells: 1146474 - The first three and last three coordinates in my text file are: - -74575.00,-3015850.00,1548.83 -74575.00,-3015825.00,1548.33 -74575.00,-3015800.00,1547.50 . . . -49324.99,-2987575.01,1510.98 -49324.99,-2987550.01,1511.24 -49324.99,-2987525.01,1511.47 - Can someone help me to figure out what's going on? Thanks Hanlie ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.rast.stats
On 06/04/2010 09:53 AM, Sandile Gumede wrote: Hi When I'm running the command below I'm getting some errors, I pasted them below v.rast.stats -ec vector=catchments raster=rainfall colprefix=rain_cols percentile=90 DBMI-DBF driver error: SQL parser error: syntax error, unexpected NULL_VALUE processing 'NULL' in statement: UPDATE catchments SET rain_cols_=-NULL WHERE cat=1 Error in db_execute_immediate ERROR: Error while executing: 'UPDATE catchments SET rain_cols_=-NULL WHERE cat=1' Hi Sandile Is there raindata for all the catchments? Do they cover the same region? i.e. maybe you can post the output of each of these: g.region -p v.info -g catchments r.info -g rainfall -- Kind Regards TS Gumede Student Researcher Meraka Institute, CSIR office : 012 841 2606 mobile : 072 258 1650 This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver Arava Development Co. +972-52-3665918 http://www.surfaces.co.il ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] gdal
On 06/04/2010 03:28 PM, Sandile Gumede wrote: Hi I've got a geotiff file which covers the area of South Africa, from that tif file I want to take a portion of Cape Flats (Western Cape) and run it on GRASS. Which command should I use to do this? I'd just run gdal_translate -projwin ulx uly lrx lry src.tif capeflats.tif where ulx is upper left x coord, uly is the upper left y coord etc... That will give you a clipped tif to import to GRASS. -- Kind Regards TS Gumede Student Researcher Meraka Institute, CSIR office : 012 841 2606 mobile : 072 258 1650 This mail was received via Mail-SeCure System. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver Arava Development Co. +972-52-3665918 http://www.surfaces.co.il ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Importing Dems with r.in.xyz
On 06/04/2010 04:20 PM, Hanlie Pretorius wrote: Hi, I've been following the procedure at http://grass.osgeo.org/wiki/Import_XYZ to import DEMs into GRASS. It has worked for two DEMS, but I have a problem with the third one at the step where one verifies that the number of rows in the ASCII file corresponds to the number of cells in the enlarged region. At this point, g.region reports 1146474 cells in the region, while I have 1146370 lines of coordinates in my file. Assuming you already did the r.in.xyz -s ... step to get and set your region to match the input data. So it looks like there are about 100 coordinates missing from the ASCII file. Maybe holes in the data? One way to work around this might be to import the point data as a 3D vector, then run the usual v.surf.rst interpolation. v.in.xyz -z in=ascii.txt out=elev_pts z=3 fs=, then g.region vect=elev_pts res=choose appropriate resolution v.surf.rst elev_pts elev=dem layer=0 The output of g.region is: - projection: 99 (Transverse Mercator) zone: 0 datum: ** unknown (default: WGS84) ** ellipsoid: wgs84 north: -49312.49 south: -74587.5 west: -3015862.5 east: -2987512.5 nsres: 25.0989 ewres: 25 rows: 1011 cols: 1134 cells: 1146474 - The first three and last three coordinates in my text file are: - -74575.00,-3015850.00,1548.83 -74575.00,-3015825.00,1548.33 -74575.00,-3015800.00,1547.50 . . . -49324.99,-2987575.01,1510.98 -49324.99,-2987550.01,1511.24 -49324.99,-2987525.01,1511.47 - Can someone help me to figure out what's going on? Thanks Hanlie ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user This mail was received via Mail-SeCure System. -- Micha Silver Arava Development Co. +972-52-3665918 http://www.surfaces.co.il ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
