:Re: [GRASS-user] Some problem about the grass
Thank you very for your reply! In fact I mean the Categories, and what you said is very helpful for my understanding something about the grass though not much. I think what I am now confusing about the grass is how it describe a map either a raster or a vector (actually the raster is my primacy), in other words what can we get from a raster or a vector map? I had read the manual you mentioned Raster data processing in GRASS GIS, it says The geographic boundaries of the raster map are described by the north, south, east, and west fields”, but how?? Is there a example? I often encounter the following two questions when I use the garss: 1) The “r.out.gdal”. When I try to convert a map to a geotiff, however the image I got is often a black blank! 2) Some command which can run normally in the garss can not in the cmd in windows with the tips: C:\Documents and Settings\Administratorr.shaded.relief map=elevation@permanent shadedmap=ttsh altitude=30 azimuth=270 zmult=1 scale=1 --overwrite '' neither a internal or outer command nor a executable program or batch file I have configured the environment variables, and some other commands can run normally such as r.buffe r.info…… PS: I have just register the maillists of the grass and how can I directly ask or reply in the web pages as in a forum? Thanks! Apachemaven 2009-3-18 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] problem about import of shp file
Dear everyone, Hello. I made my first challenge to import Arc-GIS shp file to GRASS. But I have some trouble. This may be an elementary question...,but could anyone help me, please? What I want to do (final destination) is to 1) Import vector map (Arc-GIS shp file) into location1(UTM) 2) Reproject vector data into location2 (lat-lon). 3) Overlay on raster map in location2. I need area information in vector map (Arc-GIS shp file) . But I have a trouble in step 1, and area (polygon) cannot be shown... Outline of what I did is shown below. Details of commands and error message are shown in attached file. 1) Create UTM new location 47N, but the area was extended to eastward because the area includes a part of 48N This seems to succeed. 2) v.in.ogr dsn=lu2543s.shp output=lu2543s failed. 3) v.in.ogr dsn=lu2543s.shp output=lu2543s -o Succeeded, but had an error message. Seems to failed to Import of area and isles...??? 4) d.vect lu2543s. Some map was shown, but had an error to show area. 5) d.what.vect failed. Problem about topology. 6) v.category input=lu2543s option=report failed * I have 6 another files in addition to shp file to be imported (lu2543s.shp), which is lu2543s.cpg, lu2543s.dbf, lu2543s.prj, lu2543s.sbn, lu2543s.sbx, lu2543s.shx. But I am not sure whether I should use them to import shp file correctly or not... How can I import shp file correctly? I hope someone kindly advice me, please... Thank you. N. Yoshifuji 2009-3-18 [ 1) created location] projection: 1 (UTM) zone: 47 north: 330 south: 0 east:1666000 west:0 e-w res: 50 n-s res: 50 total rows: 66000 total cols: 33320 total cells: 2,199,120,000 [ 2) import shp file] GRASS 6.2.3 (utm47n3):~ v.in.ogr dsn=Thai_GIS/lu2543s.shp output=lu2543s A datum name wgs84 (WGS_1984) was specified without transformation parameters. Note that the GRASS default for wgs84 is towgs84=0.000,0.000,0.000. ERROR: Projection of dataset does not appear to match current location. LOCATION PROJ_INFO is: name: UTM datum: wgs84 towgs84: 0.000,0.000,0.000 proj: utm ellps: wgs84 a: 6378137.00 es: 0.0066943800 f: 298.2572235630 zone: 47 Dataset PROJ_INFO is: name: Transverse Mercator proj: tmerc datum: wgs84 a: 6378137.0 es: 0.006694379990141317 lat_0: 0 lon_0: 99 k: 0.999600 x_0: 50 y_0: 0 no_defs: defined You can use the -o flag to v.in.ogr to override this check. Consider to generate a new location with 'location' parameter from input data set. [3) import shp file] GRASS 6.2.3 (utm47n3):~ v.in.ogr dsn=Thai_GIS/lu2543s.shp output=lu2543s -o A datum name wgs84 (WGS_1984) was specified without transformation parameters. Note that the GRASS default for wgs84 is towgs84=0.000,0.000,0.000. Over-riding projection check. Proceeding with import... Layer: lu2543s WARNING: Default driver / database set to: driver: dbf database: $GISDBASE/$LOCATION_NAME/$MAPSET/dbf/ Importing map 89326 features... - Building topology ... 162787 primitives registered Building areas: 100% 15407 areas built 14236 isles built Attaching islands: 100% Attaching centroids: 100% Topology was built. Number of nodes : 88204 Number of primitives: 162787 Number of points: 0 Number of lines : 0 Number of boundaries: 162787 Number of centroids : 0 Number of areas : 15407 Number of isles : 14236 Number of incorrect boundaries : 148487 Number of areas without centroid : 15407 - WARNING: Cleaning polygons, result is not guaranteed! Building topology ... Topology was built. Number of nodes : 88204 Number of primitives: 162787 Number of points: 0 Number of lines : 0 Number of boundaries: 162787 Number of centroids : 0 Number of areas : - Number of isles : - - Break polygons: Registering points ... 4611190assertion n failed: file node.c, line 48 Aborted (core dumped) [4) display imported vector file] GRASS 6.2.3 (utm47n3):~ d.vect lu2543s WARNING: Cannot display areas, topology not available GRASS 6.2.3 (utm47n3):~ v.category input=lu2543s option=report ERROR: Cannot open old vector lu25...@permanent on level 2___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem about import of shp file
Natsuko YOSHIFUJI wrote: Dear everyone, Hello. I made my first challenge to import Arc-GIS shp file to GRASS. But I have some trouble. This may be an elementary question...,but could anyone help me, please? What I want to do (final destination) is to 1) Import vector map (Arc-GIS shp file) into location1(UTM) 2) Reproject vector data into location2 (lat-lon). 3) Overlay on raster map in location2. I need area information in vector map (Arc-GIS shp file) . But I have a trouble in step 1, and area (polygon) cannot be shown... Outline of what I did is shown below. Details of commands and error message are shown in attached file. 1) Create UTM new location 47N, but the area was extended to eastward because the area includes a part of 48N This seems to succeed. 2) v.in.ogr dsn=lu2543s.shp output=lu2543s failed. Projections don't match. The shapefile is in tmerc (Transverse Mercator), not utm (Universal Transverse Mercator). You should create a new location with the projection of the shapefile: Dataset PROJ_INFO is: name: Transverse Mercator proj: tmerc datum: wgs84 a: 6378137.0 es: 0.006694379990141317 lat_0: 0 lon_0: 99 k: 0.999600 x_0: 50 y_0: 0 no_defs: defined 3) v.in.ogr dsn=lu2543s.shp output=lu2543s -o Succeeded, but had an error message. Seems to failed to Import of area and isles...??? Looks like you ran out of memory? This message assertion n failed: file node.c, line 48 means that an attempt to allocate more memory failed. Consequently, cleaning failed, topology was not built and the imported vector is not usable. The solution would be to install more RAM or increase swap space. It is strongly recommended to use grass-6.4 instead of grass-6.2.3 if possible. Markus M ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.net.iso Performance
Hi all, Was unsuccessful in building Grass 6.5 on Windows (new at this!) and thus haven't been able to test v.net.iso. I will try again at some point on an Ubuntu machine, and report back. However, it may be a while! When will 6.5 or newer be available as a package in OSGeo4W? Cheers, Darren Cope http://dmcope.freeshell.org http://bluesignweekly.blogspot.com/ On Fri, Mar 13, 2009 at 3:32 AM, Markus Neteler nete...@osgeo.org wrote: On Thu, Mar 12, 2009 at 2:34 PM, Markus Metz markus.metz.gisw...@googlemail.com wrote: Darren Cope wrote: ... My question--is the fix outlined in the thread above implemented in GRASS 6.4.0RC3-2 via OSGeo4W? Not yet, waiting for more testing. Currently it is only fixed in grass65 and grass7. If so, how long should such a process expect to take? Rather minutes than hours:-) I have made more tests today with LRS http://www.grassbook.org/examples3rd_chapter6.php and the results are identical to the pre-fix times. So also the v.lrs.* which I suppose to use the dblib are ok. And pretty fast now! +1 for backporting of the dglib cache fix. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.net.iso Performance
Hi, 2009/3/18 Martin Landa landa.mar...@gmail.com: 2009/3/18 Darren Cope darrenc...@gmail.com: Ubuntu machine, and report back. However, it may be a while! When will 6.5 or newer be available as a package in OSGeo4W? no, only grass64 is available in OSGeo4W package. sorry I should read more carefully. Well, AFAIK it's not decided. First grass64 should be released, then we can include grass65 / grass7 as devel version. Martin -- Martin Landa landa.martin gmail.com * http://gama.fsv.cvut.cz/~landa ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] v.in.dxf TEXT
Hi, I'm new in grass. I'm trying to import a dxf file to grass using v.in.dxf, but I don't know how to see the text data (entity type TEXT, specificly the text under the group code 1). I'm trying with this http://www.techsoftdesign.com/kt/eng/dxf.htm sample file . For example the t...@15674 text data is Gedeckt 8 When I run: v.in.dxf 'input=/home/mingumg/data/sample.dxf' output=text_sample layers=1 --overwrite I could not see the Gedeckt 8 anywhere. By the way if I run: v.out.dxf input=text_sample output=/home/toopath/data_prova/text_sample.dxf This information has disappeared... Please could anyone help me? Thanks Mingu -- View this message in context: http://n2.nabble.com/v.in.dxf-TEXT-tp2496298p2496298.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] problem about import of shp file
Natsuko YOSHIFUJI wrote: 3) v.in.ogr dsn=lu2543s.shp output=lu2543s -o Succeeded, but had an error message. Seems to failed to Import of area and isles...??? I forgot, it may help to close all other applications to get as much memory as possible. BTW, the procedure to break polygons should use about 20% less memory in grass65 than in grass64 and something between 20% and 40% less memory in grass7 than in grass64. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.in.dxf TEXT
mingumg wrote: Hi, I'm new in grass. I'm trying to import a dxf file to grass using v.in.dxf, but I don't know how to see the text data (entity type TEXT, specificly the text under the group code 1). To the best of my knowledge, v.in.dxf does not import any attribute information. -- Micha I'm trying with this http://www.techsoftdesign.com/kt/eng/dxf.htm sample file . For example the t...@15674 text data is Gedeckt 8 When I run: v.in.dxf 'input=/home/mingumg/data/sample.dxf' output=text_sample layers=1 --overwrite I could not see the Gedeckt 8 anywhere. By the way if I run: v.out.dxf input=text_sample output=/home/toopath/data_prova/text_sample.dxf This information has disappeared... Please could anyone help me? Thanks Mingu ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem about import of shp file
Dear Mr. Markus Metz , Thank you very much for your quick reply! I misunderstood that tmerc is the same as UTM... I tried again following your advice, but I have additional questions... Could you kindly reply me, please? Projections don't match. The shapefile is in tmerc (Transverse Mercator), not utm (Universal Transverse Mercator). You should create a new location with the projection of the shapefile: Dataset PROJ_INFO is: name: Transverse Mercator proj: tmerc datum: wgs84 a: 6378137.0 es: 0.006694379990141317 lat_0: 0 lon_0: 99 k: 0.999600 x_0: 50 y_0: 0 no_defs: defined 3) v.in.ogr dsn=lu2543s.shp output=lu2543s -o Succeeded, but had an error message. Seems to failed to Import of area and isles...??? Looks like you ran out of memory? This message assertion n failed: file node.c, line 48 means that an attempt to allocate more memory failed. Consequently, cleaning failed, topology was not built and the imported vector is not usable. The solution would be to install more RAM or increase swap space. I created new tmerc location. v.in.ogr dsn=lu2543s.shp output=lu2543s worked, but I had same error message ssertion n failed: file node.c, line 48. You told me that this is because of shortage of RAM. So, I tried to use another PC (Grass6.2.2), created tmerc location, and tried to import the same shp file. Then, I received an error message saying that projections don't match! In this time, data set projection was reported as follows; Dataset PROJ_INFO is: name: Universe Transverse Mercator proj: utm datum: wgs84 a: 6378137.0 es: 0.006694379990141316 zone: 47 no_defs: defined The data supplier told me that the projection of this shp file is UTM. So, the message of new PC seems to be correct, and previous message may be wrong. Anyway, it is very confusing... In addition, import of shp file was failed, too, probably because of shortage of RAM... It is strongly recommended to use grass-6.4 instead of grass-6.2.3 if possible. Can such confusing situation be resolved if I install grass-6.4? I forgot, it may help to close all other applications to get as much memory as possible. BTW, the procedure to break polygons should use about 20% less memory in grass65 than in grass64 and something between 20% and 40% less memory in grass7 than in grass64. In grass HP (http://wgrass.media.osaka-cu.ac.jp/grassh/download/index.php), I can find only GRASS-6.4. Where can I get grass-6.5 and grass7? Now I use grass 6.2.2. Does grass-6.4 use smaller memory than grass 6.2.2? Thank you in advance. Natsuko Yoshifuji ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Plot 3D Drillholes in Grass GIS
Casey, I was very interested in your original email, so I tried to track down the GRASS modules mentioned in the paper you cited. As it turns out the authors have written several related papers. While I was able to obtain the other papers, I could not locate the software. So, I emailed the authors (Shinji Masumoto) and just got a response. Unfortunately, the authors are unwilling to release the software into the public domain. I have not specifically tried NVIZ visualization as you describe, but I think the 3rd Ed. of Open Source GIS: A GRASS GIS Approach covers how to create 3D vectors which can be brought into NVIZ to produce what you need. Regards, Tom Casey Vandenberg wrote: I was wondering if anybody is familiar with ways of plotting drill hole data in Grass (i.e. to visualize in NVIZ or to manage attributes, query etc...) based on point features with known dip, strike and azimuth angles. I have come across a paper (Development of Visualization Tools for Geologic Information on Grass GIS) which seems to do just that, however as far as I can tell it has not been released or provided to the Grass Community. If anybody is familiar with anything else, it would be greatly appreciated. Thanks, Casey ___ 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
Re: [GRASS-user] problem about import of shp file
Natsuko YOSHIFUJI wrote: I created new tmerc location. v.in.ogr dsn=lu2543s.shp output=lu2543s worked, but I had same error message ssertion n failed: file node.c, line 48. You told me that this is because of shortage of RAM. I guess. I hope it's not a bug in grass, but to my knowledge this error has not been reported before and I tested v.in.ogr recently (in 6.4, 6.5, and 7) and also did not get such an error. So, I tried to use another PC (Grass6.2.2), created tmerc location, and tried to import the same shp file. Then, I received an error message saying that projections don't match! In this time, data set projection was reported as follows; Dataset PROJ_INFO is: name: Universe Transverse Mercator proj: utm datum: wgs84 a: 6378137.0 es: 0.006694379990141316 zone: 47 no_defs: defined The data supplier told me that the projection of this shp file is UTM. So, the message of new PC seems to be correct, and previous message may be wrong. Anyway, it is very confusing... This is confusing indeed. There should be a file called lu2543s.prj . This is a plain text file with the projection information for lu2543s.shp, you can open it with any text editor, see what's in there and create a location according to this information. If v.in.ogr still complains that projections don't match, you can use v.in.ogr -o, but something may be wrong with gdal and/or proj4. Can such confusing situation be resolved if I install grass-6.4? grass-6.4 has lots of bug fixes and improvements over previous versions, I would always recommend to use grass-6.4 instead of earlier versions. And also update gdal and proj4. Where can I get grass-6.5 and grass7? grass-6.5 and grass-7 are only available as source code via svn. If you are not familiar with compiling from source, this can be a bit tricky not to say it could cost you a day or more. Information on how to download source is here: http://trac.osgeo.org/grass/wiki/DownloadSource There is a wiki page on how to compile and install: http://grass.osgeo.org/wiki/Compile_and_Install It is easier to install grass-6.4 binaries. Does grass-6.4 use smaller memory than grass 6.2.2? Not for breaking polygons. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Plot 3D Drillholes in Grass GIS
I was trying to find related software or a way to work with boreholes and geological cross sections in grass, not just for the drawing part but a way to use it for geological analysis. It's very sad the fact that these modules are not for public use, even with a stricter licensing than GPL. Regards, Vasilis On Wed, Mar 18, 2009 at 1:45 PM, Thomas Adams thomas.ad...@noaa.gov wrote: Casey, I was very interested in your original email, so I tried to track down the GRASS modules mentioned in the paper you cited. As it turns out the authors have written several related papers. While I was able to obtain the other papers, I could not locate the software. So, I emailed the authors (Shinji Masumoto) and just got a response. Unfortunately, the authors are unwilling to release the software into the public domain. I have not specifically tried NVIZ visualization as you describe, but I think the 3rd Ed. of Open Source GIS: A GRASS GIS Approach covers how to create 3D vectors which can be brought into NVIZ to produce what you need. Regards, Tom Casey Vandenberg wrote: I was wondering if anybody is familiar with ways of plotting drill hole data in Grass (i.e. to visualize in NVIZ or to manage attributes, query etc...) based on point features with known dip, strike and azimuth angles. I have come across a paper (Development of Visualization Tools for Geologic Information on Grass GIS) which seems to do just that, however as far as I can tell it has not been released or provided to the Grass Community. If anybody is familiar with anything else, it would be greatly appreciated. Thanks, Casey ___ 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 mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] v.net.iso Performance
On Wed, Mar 18, 2009 at 9:59 AM, Darren Cope darrenc...@gmail.com wrote: Hi all, Was unsuccessful in building Grass 6.5 on Windows (new at this!) and thus haven't been able to test v.net.iso. Colin and me have started to write up a step by step compilation document (special thanks to Colin!): http://trac.osgeo.org/grass/wiki/CompileOnWindows which is still evolving but hopefully already helpful. Comment or direct modifications welcome. I will try again at some point on an Ubuntu machine, and report back. However, it may be a while! When will 6.5 or newer be available as a package in OSGeo4W? That's unlikely unless we find a volunteer. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem about import of shp file
On Wed, Mar 18, 2009 at 1:29 PM, Markus Metz markus.metz.gisw...@googlemail.com wrote: Natsuko YOSHIFUJI wrote: ... Where can I get grass-6.5 and grass7? grass-6.5 and grass-7 are only available as source code via svn. If desired, I can set up cronjobs to generate also the grass-6.5 and grass-7 binaries. For grass-7, there is the problem that I don't manage to compile CAIRO = 1.5.8 on Fedora FC4 which is used on grass.osgeo.org. Help (i.e., 32bit RPM) welcome. Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] truncating filds of type text during export from grass to shp
Hi It is well-known that ogr2ogr truncate character field to 80 chars during export to dbf. But I have a problem with exporting data where atributes are stored in postgree SQL database and one field is type text with more than 255 chars. It seems it cannot be exported with r.out.ogr so I use ogr2ogr, according to its help: ogr2ogr -f ESRI Shapefile -sql select *, cast(branch as text) from cieki_glowna rivers.shp grassdata/Rzeki/glowna/vector/cieki_glowna/head but it also failed So it is a way to do it, without exporting columns form pgSQL as text and attaching them in ESRI externaly?? kindly regards Jarek ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem about import of shp file
Dear Mr. Markus Metz, So, I tried to use another PC (Grass6.2.2), created tmerc location, and tried to import the same shp file. Then, I received an error message saying that projections don't match! In this time, data set projection was reported as follows; Dataset PROJ_INFO is: name: Universe Transverse Mercator proj: utm datum: wgs84 a: 6378137.0 es: 0.006694379990141316 zone: 47 no_defs: defined The data supplier told me that the projection of this shp file is UTM. So, the message of new PC seems to be correct, and previous message may be wrong. Anyway, it is very confusing... This is confusing indeed. There should be a file called lu2543s.prj . This is a plain text file with the projection information for lu2543s.shp, you can open it with any text editor, see what's in there and create a location according to this information. If v.in.ogr still complains that projections don't match, you can use v.in.ogr -o, but something may be wrong with gdal and/or proj4. The contents of lu2543s.prj was as follows. PROJCS[WGS_1984_UTM_Zone_47N,GEOGCS[GCS_WGS_1984,DATUM[D_WGS_1984,SPHEROID[WGS_1984,6378137.0,298.257223563]],PRIMEM[Greenwich,0.0],UNIT[Degree,0.0174532925199433]],PROJECTION[Transverse_Mercator],PARAMETER[False_Easting,50.0],PARAMETER[False_Northing,0.0],PARAMETER[Central_Meridian,99.0],PARAMETER[Scale_Factor,0.9996],PARAMETER[Latitude_Of_Origin,0.0],UNIT[Meter,1.0]] There are two projection information: PROJCS[WGS_1984_UTM_Zone_47Nand PROJECTION[Transverse_Mercator]. I wonder correct projection is UTM47N, isn't it? Where can I get grass-6.5 and grass7? grass-6.5 and grass-7 are only available as source code via svn. If you are not familiar with compiling from source, this can be a bit tricky not to say it could cost you a day or more. Information on how to download source is here: http://trac.osgeo.org/grass/wiki/DownloadSource There is a wiki page on how to compile and install: http://grass.osgeo.org/wiki/Compile_and_Install It is easier to install grass-6.4 binaries. Does grass-6.4 use smaller memory than grass 6.2.2? Not for breaking polygons. Do you mean that installing grass6.4 will not contribute to solve memory-shortage problem? But, I am not familiar with compiling from source... If I can prepare smaller shape file whose area is a part of whole area of lu2543s.shp, is there a possibility that I can avoid memory- shortage problem? Thank you and with regards, Natsuko YOSHIFUJI ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem about import of shp file
On 2009/03/18, at 22:07, Markus Neteler wrote: On Wed, Mar 18, 2009 at 1:29 PM, Markus Metz markus.metz.gisw...@googlemail.com wrote: Natsuko YOSHIFUJI wrote: ... Where can I get grass-6.5 and grass7? grass-6.5 and grass-7 are only available as source code via svn. If desired, I can set up cronjobs to generate also the grass-6.5 and grass-7 binaries. For grass-7, there is the problem that I don't manage to compile CAIRO = 1.5.8 on Fedora FC4 which is used on grass.osgeo.org. Help (i.e., 32bit RPM) welcome. Thank you very very much! I hope binaries of grass-6.5 or grass7. I am not so familiar with PC, so I try to tell you my PC environment anyway. There are 3 type of PC around me. If possible, binaries for Mac is preferable because mac has the biggest memory, but binaries for linux and windows+cygwin is available. PC1 (note pc) Windows XP professional + cygwin 1.5GB RAM Grass - 6.2.3 is already installed * I used this pc to import shp file and failed with an error message n failed node.c, line48 Aborted (core dumped). * Grass on this pc reported that lu2543s.shp is tmerc, but it may not be true... PC2 Fedora Core 5 622MB RAM Grass-6.2.2 is already installed * import of shp file was forced to stop, probably because of shortage of memory...? PC1 Mac OS 10.5.6 2.1GHz Intel Core 2 Duo 3GB 667 MHz DDR2 SDRAM Grass-6.3 is already installed * import of shp file was failed with an error message as follows: Break polygons: Registering points ... 17107648v.in.ogr(1388) malloc: *** mmap(size=624242688) failed (error code=12) *** error: can't allocate region *** set a breakpoint in malloc_error_break to debug ERROR: G_realloc: out of memory With best regards, Natsuko YOSHIFUJI ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] raster merging
Hi you all, I am working with geotiff raster data served on a webgis (mapserver + pmapper) and I would like to serve a merged raster (DEM + its shaded relief) instead of serving two distinct raster/layer with my mapfile (one for the dem and one for the shaded relief). Should I go with r.mapcalc in order to obtain a suitable GTiff for mapserver? Regards raffaele ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] raster merging
On Wed, 2009-03-18 at 16:17 +0100, Cassiel wrote: Hi you all, I am working with geotiff raster data served on a webgis (mapserver + pmapper) and I would like to serve a merged raster (DEM + its shaded relief) instead of serving two distinct raster/layer with my mapfile (one for the dem and one for the shaded relief). Should I go with r.mapcalc in order to obtain a suitable GTiff for mapserver? Regards raffaele Raffaele, perhaps r.blend is what you are looking for. After some trial and error (concerning the Percentage weight of first map for color blending) you'll get a nice looking shaded DEM that you can use (export). Kind regards, Nikos ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] problem about import of shp file
Natsuko YOSHIFUJI wrote: The contents of lu2543s.prj was as follows. PROJCS[WGS_1984_UTM_Zone_47N,GEOGCS[GCS_WGS_1984,DATUM[D_WGS_1984,SPHEROID[WGS_1984,6378137.0,298.257223563]],PRIMEM[Greenwich,0.0],UNIT[Degree,0.0174532925199433]],PROJECTION[Transverse_Mercator],PARAMETER[False_Easting,50.0],PARAMETER[False_Northing,0.0],PARAMETER[Central_Meridian,99.0],PARAMETER[Scale_Factor,0.9996],PARAMETER[Latitude_Of_Origin,0.0],UNIT[Meter,1.0]] There are two projection information: PROJCS[WGS_1984_UTM_Zone_47N and PROJECTION[Transverse_Mercator]. I wonder correct projection is UTM47N, isn't it? You said the data provider said it should be UTM, so I would use UTM. Do you mean that installing grass6.4 will not contribute to solve memory-shortage problem? Probably not, I'm afraid. From your other mail I could see that there is quite a bit of memory on the Mac, but that runs out of memory too. Can't you increase swap space on the Mac or for Windows XP? Just curious, what is the size of lu2543s.shp? And, is it in the public domain so that you can give me a link and I can get it myself from somewhere? If I can prepare smaller shape file whose area is a part of whole area of lu2543s.shp, is there a possibility that I can avoid memory-shortage problem? Yes, but there is no need to prepare a smaller shape file, there are two different options for that in v.in.ogr, either -r or spatial=xmin,ymin,xmax,ymax, see manual. Best, Markus PS: Please don't call me Mr. Markus Metz, Markus is really enough :-) ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: :Re: [GRASS-user] Some problem about the grass
apachemaven: Thank you very for your reply! In fact I mean the Categories, and what you said is very helpful for my understanding something about the grass though not much. I think what I am now confusing about the grass is how it describe a map either a raster or a vector (actually the raster is my primacy), in other words what can we get from a raster or a vector map? Maybe all this is known to you, but I am writing it anyway. These are generic questions. You might want to read about data models in general (basically raster and vector data). A nice and practical view, through GRASS, about all this is the GRASS book [*]. You will get an answer as well about what categories are with respect to raster maps. You can think rasters as matrices and vectors as points, lines or closed polygons. Depending on the nature of what you want to represent in a map, you either choose raster (like for temperature of precipitation maps) or vector (location of weather stations, roads or burned areas). Of course there are data that can be represented in both raster and vector maps. It's always a question of what you have and what you want to achieve. I had read the manual you mentioned Raster data processing in GRASS GIS, it says The geographic boundaries of the raster map are described by the north, south, east, and west fields”, but how?? Is there a example? Absolutely. I (can) set the geographic boundaries for Peloponnese, South Greece for example like that: g.region s=36 e=23.6 n=38.5 w=21 Note that the above values are geographic (=degrees) and NOT projected coordinates. We can do something similar in a projected coordinate system with values given in metres. You know that there are many GRASS-related tutorials [1]. And even more, GFOSS related [2] ? 1) The “r.out.gdal”. When I try to convert a map to a geotiff, however the image I got is often a black blank! That one might NOT be a problem if you just have not set your region correctly. It can be a problem sometimes related with color table rendering. You could try to remove the color table (r.colors -r YourMap) before exporting. There are some tips in the manual about this [3]. And there is a long discussion in GRASS-trac [4] 2) Sorry but I have no working experience with GRASS under Operating Systems (or environments) other than Linux. PS: I have just register the maillists of the grass and how can I directly ask or reply in the web pages as in a forum? You can register your self in GRASS' wiki [5] and add/edit almost anything. The wiki is from and for the community. Cheers, Nikos --- [*] http://www.grassbook.org/ [1] http://grass.osgeo.org/gdp/tutorials.php [2] http://wiki.osgeo.org/wiki/Educational_Content_Inventory [3] http://grass.itc.it/grass64/manuals/html64_user/r.out.gdal.html [4] http://trac.osgeo.org/grass/ticket/73 [5] http://grass.osgeo.org/wiki/Main_Page ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: :Re: [GRASS-user] Some problem about the grass
On Wed, 2009-03-18 at 18:51 +0100, Nikos Alexandris wrote: (like for temperature of precipitation maps) a small but perhaps an important typo: ...like for temperature OR precipitation maps... ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Building GRASS on Windows
Hi all, I'm following along with the instructions here: http://trac.osgeo.org/grass/wiki/CompileOnWindows and am getting stuck at Step 5 under PDCurses. When running: make -f mingwin32.mak DLL=Y all The following error occurs: $ make -f mingwin32.mak DLL=Y all gcc -c -O2 -Wall -I.. -DPDC_DLL_BUILD ../pdcurses/addch.c In file included from ../curspriv.h:25, from ../pdcurses/addch.c:14: ../curses.h:50:57: stdio.h: No such file or directory In file included from ../curspriv.h:25, from ../pdcurses/addch.c:14: ../curses.h:928: error: syntax error before '*' token ../curses.h:1010: error: syntax error before FILE ../curses.h:1027: error: syntax error before FILE make: *** [addch.o] Error 1 Any thoughts? I'd really like to get a new version built so I can play with v.net.iso improvements, but am struggling! Darren Cope http://dmcope.freeshell.org http://bluesignweekly.blogspot.com/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] r.colors- transparent, HSB color
Hi, I have been using r.mapcalc to set various combinations of colors in a 13-color DRG to null to allow the map to be displayed as a series of overlays on top of other maps. It occurs to me that it would be useful to be able to accomplish this using r.colors by defining the color of specific values (or ranges of values) as transparent. Unfortunately, I have not been able to find a way to do this. Is it possible to define a color as transparent without changing the pixels to null pixels? If not, would this be technically feasible in a future version of GRASS? Is there a different approach that would allow me to generate multiple examples of partially transparent maps without the overhead of duplicating the raster map multiple times? Secondly, is it possible to define colors within GRASS using HSB (or HSL) instead of RGB values? There are some kinds of calculations, for example reducing the saturation of all colors by a given amount, that are easier to accomplish in the HSB domain. Thanks, -Dwight ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] r.colors- transparent, HSB color
Dwight Needels wrote: I have been using r.mapcalc to set various combinations of colors in a 13-color DRG to null to allow the map to be displayed as a series of overlays on top of other maps. It occurs to me that it would be useful to be able to accomplish this using r.colors by defining the color of specific values (or ranges of values) as transparent. Unfortunately, I have not been able to find a way to do this. Is it possible to define a color as transparent without changing the pixels to null pixels? No. If not, would this be technically feasible in a future version of GRASS? The main issue is that each file which reads colours would have to be extended to understand transparency. AFAICT, there are 38 modules and 3 libraries which would be affected. Is there a different approach that would allow me to generate multiple examples of partially transparent maps without the overhead of duplicating the raster map multiple times? For integer maps, you can create a reclass map using r.reclass. This has much lower overhead than creating a complete new map with r.mapcalc. Secondly, is it possible to define colors within GRASS using HSB (or HSL) instead of RGB values? There are some kinds of calculations, for example reducing the saturation of all colors by a given amount, that are easier to accomplish in the HSB domain. This isn't supported directly. Are you talking about simply specifying colours in HSB, and keeping the internal data in RGB, or having interpolation performed in HSB space? The former wouldn't be particularly involved; just an extension to G_parse_color_rule(). OTOH, you can obtain this functionality by pre-processing the rules with a script. -- Glynn Clements gl...@gclements.plus.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user