<html><head><style type="text/css"><!-- DIV {margin:0px;} --></style></head><body><div style="font-family:times new roman, new york, times, serif;font-size:12pt"><DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">Unsubscribe this mail please<BR><BR>
<DIV style="FONT-SIZE: 12pt; FONT-FAMILY: times new roman, new york, times, serif">----- Original Message ----<BR>From: "grassuser-request@grass.itc.it" <grassuser-request@grass.itc.it><BR>To: grassuser@grass.itc.it<BR>Sent: Wednesday, March 21, 2007 3:54:59 PM<BR>Subject: grassuser Digest, Vol 11, Issue 40<BR><BR>
<DIV>Send grassuser mailing list submissions to<BR> grassuser@grass.itc.it<BR><BR>To subscribe or unsubscribe via the World Wide Web, visit<BR> <A href="http://grass.itc.it/mailman/listinfo/grassuser" target=_blank>http://grass.itc.it/mailman/listinfo/grassuser</A><BR>or, via email, send a message with subject or body 'help' to<BR> grassuser-request@grass.itc.it<BR><BR>You can reach the person managing the list at<BR> grassuser-owner@grass.itc.it<BR><BR>When replying, please edit your Subject line so it is more specific<BR>than "Re: Contents of grassuser digest..."<BR><BR><BR>Today's Topics:<BR><BR> 1. Re: v.to.rast "unknown use type" (Markus Neteler)<BR> 2. Re: r.surf.contour (Markus Neteler)<BR> 3. rasterizing points (xrubio@arsystel.com)<BR> 4. Re: project units vs. imported image units ? (Ricardo Oliveira)<BR> 5. Re:
rasterizing points (Markus Neteler)<BR> 6. Inputting vector points in non Lat/Lon cords (Bryan Main)<BR> 7. create a slope vector map (Luca Penasa)<BR> 8. Help the newbie importing data into GRASS (Jorge Echeverri)<BR> 9. Re: Encountering an error importing shapefiles with v.in.ogr<BR> (Thomas Adams)<BR><BR><BR>----------------------------------------------------------------------<BR><BR>Message: 1<BR>Date: Wed, 21 Mar 2007 09:23:33 +0100<BR>From: Markus Neteler <neteler@itc.it><BR>Subject: Re: [GRASS-user] v.to.rast "unknown use type"<BR>To: Seb <spluque@gmail.com><BR>Cc: grassuser@grass.itc.it<BR>Message-ID: <4600EB85.5070603@itc.it><BR>Content-Type: text/plain; charset=ISO-8859-1<BR><BR>Seb wrote on 03/21/2007 04:08 AM:<BR>> Hi,<BR>><BR>> I'm having problems with the following call to rasterize a vector:<BR>><BR>><BR>> v.to.rast
input=test output=test use=z type=point layer=1 value=1 rows=4096<BR>><BR>><BR>> The vector contains 1 minute gridded GEBCO bathymetry data, which I<BR>> imported with v.in.ascii, assigning depth data to the 'z' coordinate.<BR>> However, the call above gives the error message:<BR>><BR>> ---<---------------cut here---------------start-------------->---<BR>> Converted areas: [0] of [0]<BR>><BR>> Converted points/lines: [329427] of [329427]<BR>><BR>> Unknown use type: 4<BR>> ---<---------------cut here---------------end---------------->---<BR>><BR>><BR>> I'd greatly appreciate some clues on what may be going on.<BR>> <BR><BR>I have fixed this yesterday, was running into the same problem.<BR><BR>Markus<BR><BR>------------------<BR>ITC -> dall'1 marzo 2007 Fondazione Bruno Kessler<BR>ITC -> since 1 March 2007 Fondazione Bruno
Kessler<BR>------------------<BR><BR><BR><BR>------------------------------<BR><BR>Message: 2<BR>Date: Wed, 21 Mar 2007 09:25:21 +0100<BR>From: Markus Neteler <neteler@itc.it><BR>Subject: Re: [GRASS-user] r.surf.contour<BR>To: GRASS user list <grassuser@grass.itc.it><BR>Message-ID: <4600EBF1.70007@itc.it><BR>Content-Type: text/plain; charset=ISO-8859-1<BR><BR>Amaginary Studios wrote on 03/21/2007 12:10 AM:<BR>> I am running r.surf.contour and it is ... taking quite a long time...<BR>><BR>> I am just wondering how long the average processing time is with the<BR>> following g.region<BR>><BR>> /projection: 1 (UTM)/<BR>> /zone: 9/<BR>> /datum: nad83/<BR>> /ellipsoid: grs80/<BR>> /north: 6052143.70894/<BR>> /south: 6018782.69306/<BR>>
/west: 644196.986023/<BR>> /east: 657085.826063/<BR>> /nsres: 4.00012181/<BR>> /ewres: 4.00026072/<BR>> /rows: 8340/<BR>> /cols: 3222/<BR>> /cells: 26 871 480/<BR>><BR>Sidenote:<BR>Not sure if a resolution of /4.00012181 is desired.<BR>If you want to get it to 4.0, run<BR><BR>g.region -a -p res=4<BR>and the coordinates box will be slightly extended to match the<BR>resolution precisely.<BR><BR>Markus<BR><BR>/<BR><BR>------------------<BR>ITC -> dall'1 marzo 2007 Fondazione Bruno Kessler<BR>ITC -> since 1 March 2007 Fondazione Bruno Kessler<BR>------------------<BR><BR><BR><BR>------------------------------<BR><BR>Message: 3<BR>Date: Wed, 21 Mar 2007 11:38:35 +0100 (CET)<BR>From:
xrubio@arsystel.com<BR>Subject: [GRASS-user] rasterizing points<BR>To: grassuser@grass.itc.it<BR>Message-ID:<BR> <56237.161.116.25.209.1174473515.squirrel@mail.arsystel.com><BR>Content-Type: text/plain;charset=iso-8859-1<BR><BR><BR>Hi all, I'm trying to raster a vector layer with 3 or 4 points but I don't<BR>know why it is not working.<BR>The command I'm using is:<BR><BR>v.to.rast input=landmarks_v output=landmarks_r value=1 use=val<BR><BR>but it simply generates an empty raster map. I think that grass is trying<BR>to rasterize lines:<BR><BR>Loading vector information ... 0 mins 00 secs<BR>Sorting areas by size ... 0 areas 0 mins 00 secs<BR>Processing lines ... 4
lines 0 mins 00 secs<BR>Writing raster map ... 0 mins 00 secs<BR>Creating support files for raster map ... 0 mins 00 secs<BR><BR>but as the information are points it is not generating any line. I want to<BR>create a raster map with null values except for 3 or 4 cells (where there<BR>are points) with '1'.<BR>Any idea about how it can be done?<BR><BR>Thanks!<BR><BR>Xavi<BR><BR><BR><BR>------------------------------<BR><BR>Message: 4<BR>Date: Wed, 21 Mar 2007 10:51:31 -0000<BR>From: "Ricardo Oliveira" <autopilot@oniduo.pt><BR>Subject: Re: [GRASS-user] project units vs. imported image units ?<BR>To: <grassuser@grass.itc.it><BR>Message-ID:
<CD40D420FBE630488EF0CC85AED081423B271F@RVS2.at.isp><BR>Content-Type: text/plain; charset="iso-8859-1"<BR><BR>Thanks for the quick replies.<BR>I have been testing the solutions:<BR><BR>I can use g.region to create different regions (with specific grid size) so they can match the resolution of the tifs i am importing (this will be useful because the tifs have different resolutions). After that r.region -c will recompute the image to the extents of the current region. This seems to be what i needed :D<BR><BR>Again, thank you<BR>Ricardo<BR><BR>On Tue, 20 Mar 2007 13:19:23 -0700, "Amaginary Studios"<BR><amagine@telus.net> said:<BR>> Hmmm,<BR>> <BR>> I believe you will need to specify the correct "resolution" in your <BR>> "g.region" settings which will ... match the resolution of your .tiff <BR>> files<BR>> <BR>> .tfw files are simply Text files so I suppose you could create your <BR>> own within a text
editor.<BR>> <BR>> for example a .twf file with the following<BR>> <BR>> ****<BR>> <BR>> 1 nresolution<BR>> 0 x-pixel position reference point<BR>> 0 y-pixel position reference point<BR>> -1 sresolution<BR>> 0 x-postion reference<BR>> 0 y-position reference<BR>> <BR>> ****<BR>> <BR>> says, that each pixel in the tiff image is equal to 1x1 (in your <BR>> case probably meters). and the Top left corner is referenced at 0m x <BR>> 0m.<BR>> <BR>> So, it is a matter of matching your image resolution size to the <BR>> correct resolution region settings in GRASS.<BR>> <BR>> I hope
this helps a bit...<BR>> <BR>> Mars<BR>-------------- next part --------------<BR>An HTML attachment was scrubbed...<BR>URL: <A href="http://grass.itc.it/pipermail/grassuser/attachments/20070321/c8510632/attachment-0001.html" target=_blank>http://grass.itc.it/pipermail/grassuser/attachments/20070321/c8510632/attachment-0001.html</A><BR><BR>------------------------------<BR><BR>Message: 5<BR>Date: Wed, 21 Mar 2007 13:55:14 +0100<BR>From: Markus Neteler <neteler@itc.it><BR>Subject: Re: [GRASS-user] rasterizing points<BR>To: grassuser@grass.itc.it<BR>Message-ID: <46012B32.8070806@itc.it><BR>Content-Type: text/plain; charset=ISO-8859-1<BR><BR>xrubio@arsystel.com wrote on 03/21/2007 11:38 AM:<BR>> Hi all, I'm trying to raster a vector layer with 3 or 4 points but I don't<BR>> know why it is not working.<BR>> The command I'm using is:<BR>><BR>> v.to.rast input=landmarks_v output=landmarks_r value=1 use=val<BR>><BR>> but it simply generates
an empty raster map. I think that grass is trying<BR>> to rasterize lines:<BR>><BR>> Loading vector information ... 0 mins 00 secs<BR>> Sorting areas by size ... 0 areas 0 mins 00 secs<BR>> Processing lines ... 4 lines 0 mins 00 secs<BR>> Writing raster map ... 0 mins 00 secs<BR>> Creating support files for raster map ... 0 mins 00 secs<BR>><BR>> but as the information are points it is not generating any line. I want to<BR>> create a raster map with
null values except for 3 or 4 cells (where there<BR>> are points) with '1'.<BR>> Any idea about how it can be done?<BR>> <BR><BR>what does<BR>r.univar landmarks_r<BR>say?<BR><BR>To see the map, try this:<BR><BR>g.region rast=landmarks_r -p<BR>d.erase<BR>d.rast landmarks_r<BR><BR><BR>Markus<BR><BR>------------------<BR>ITC -> dall'1 marzo 2007 Fondazione Bruno Kessler<BR>ITC -> since 1 March 2007 Fondazione Bruno Kessler<BR>------------------<BR><BR><BR><BR>------------------------------<BR><BR>Message: 6<BR>Date: Wed, 21 Mar 2007 10:03:25 -0400<BR>From: "Bryan Main" <bryan.main@gmail.com><BR>Subject: [GRASS-user] Inputting vector points in non Lat/Lon cords<BR>To: grassuser@grass.itc.it<BR>Message-ID:<BR> <33aba5420703210703q178f13deg2c9ee56ba6837988@mail.gmail.com><BR>Content-Type: text/plain; charset=ISO-8859-1; format=flowed<BR><BR>I'm having trouble finding out if I can input vector points in a<BR>different
coordinate system other then Lat/Long. I'm trying to use<BR>the Irish Grid system since all of my information is already in that<BR>format. Is this possible, or will I have to convert all of my points<BR>from the Irish grid to lat/long?<BR><BR>Bryan<BR><BR><BR><BR>------------------------------<BR><BR>Message: 7<BR>Date: Wed, 21 Mar 2007 15:32:17 +0100<BR>From: Luca Penasa <luca.penasa@email.it><BR>Subject: [GRASS-user] create a slope vector map<BR>To: grassuser@grass.itc.it<BR>Message-ID: <460141F1.8020805@email.it><BR>Content-Type: text/plain; charset=ISO-8859-1; format=flowed<BR><BR>i would like to create a vector map with vectors rapresenting the slope...<BR>i need to create these features:<BR>- every vector should point to the max slope line, with an arrow <BR>representig the flow direction<BR>- the lenght should rapresent the slope value in degrees, and the vector <BR>line should start from the center of the square<BR>- i want to
have a vector for every 10x10m square, with an average slope <BR>value rapresented by the vector lenght<BR><BR>I know how to create the slope map, with 10x10m resolution... but.. how <BR>can i create the vector map with all the features i said???<BR><BR>How can i do this??? any suggestions????<BR><BR>Sorry for my english<BR>Thanks, Luca<BR><BR><BR>--<BR>Email.it, the professional e-mail, gratis per te: <A href="http://www.email.it/f" target=_blank>http://www.email.it/f</A><BR><BR>Sponsor:<BR>Analisi di Bilancio: come capire correttamente un bilancio aziendale<BR>* Esempi, modelli pratici, software utili<BR>Clicca qui: <A href="http://adv.email.it/cgi-bin/foclick.cgi?mid=6199&d=21-3" target=_blank>http://adv.email.it/cgi-bin/foclick.cgi?mid=6199&d=21-3</A><BR><BR><BR><BR>------------------------------<BR><BR>Message: 8<BR>Date: Wed, 21 Mar 2007 16:52:49 +0100<BR>From: Jorge Echeverri <jorge.echeverri@Stein.de><BR>Subject: [GRASS-user] Help the newbie importing
data into GRASS<BR>To: grassuser@grass.itc.it<BR>Message-ID: <1174492369.3837.27.camel@linux-a46z.site><BR>Content-Type: text/plain; charset="us-ascii"<BR><BR>Hi,<BR><BR>I need to do the following:<BR><BR>Display a sewer network and generate charts with classification in a<BR>colour scale based on attributes like Age, Material, Diameter, etc.<BR><BR>Data are in a MySQL database:<BR><BR>The network is geometrically defined by the relationship between these<BR>two tables:<BR><BR>tbl_line:<BR>"LINE_ID","FROM_NODE_ID","TO_NODE_ID"<BR>1,1,2<BR>tbl_node<BR>"NODE_ID","X_COORDINATE","Y_COORDINATE","Z_COORDINATE","NET_TYPE_ID"<BR>1,3400630.04,5718486.603,,1,<BR><BR>Assuming that attributes where also stored in the same table: tbl_line,<BR>which is not the case, how would be the better way to perform this task?<BR>Attributes are actually stored in another table, tbl_basedata, which is<BR>in the same database and has a 1:c relationship to tbl_line.<BR><BR>So far I have managed
to connect to the MySQL-database but cannot figure<BR>out how to further proceed.<BR><BR>Previously I had tried following workaround:<BR>1.I performed a MySQL-query on tbl_node and tbl_line to obtain the<BR>couples of points defining a line and stored this in a .csv-file. <BR>2. From this I created a new file with an awk-programm like this:<BR>#'X_COORDINATE, Y_COORDINATE'<BR>L 2 0<BR>3400630.04 5718486.603<BR>3400639.238 5718515.776<BR>L 2 0<BR>3400639.238 5718515.776<BR>3400648.641 5718546.057<BR>L 2 0<BR>3400648.641 5718546.057<BR>3400665.041 5718562.99<BR>L 2 0<BR>3400665.041 5718562.99<BR>3400625.996 5718597.349<BR>L 2 0 <BR>3.I successfully imported this with v.in.ascii and can depict the<BR>geometry.<BR><BR>But I don't think that is the most logical way to do it as GRASS can<BR>connect directly to the DB. Any suggestions???<BR><BR>Thanks in advance<BR>-- <BR>Jorge Echeverri<BR>SP<BR>-------------- next part --------------<BR>An HTML attachment was
scrubbed...<BR>URL: <A href="http://grass.itc.it/pipermail/grassuser/attachments/20070321/59681e75/attachment-0001.html" target=_blank>http://grass.itc.it/pipermail/grassuser/attachments/20070321/59681e75/attachment-0001.html</A><BR><BR>------------------------------<BR><BR>Message: 9<BR>Date: Wed, 21 Mar 2007 11:53:59 -0400<BR>From: Thomas Adams <Thomas.Adams@noaa.gov><BR>Subject: Re: [GRASS-user] Encountering an error importing shapefiles<BR> with v.in.ogr<BR>To: Raffaele Morelli <raffaele.morelli@gmail.com><BR>Cc: GRASS user list <grassuser@grass.itc.it><BR>Message-ID: <46015517.4090205@noaa.gov><BR>Content-Type: text/plain; charset=windows-1252; format=flowed<BR><BR>Raffaele,<BR><BR>After further testing, I do see the same problem you have had, <BR>namely,"(some) areas which are not available after using v.in.ogr." This <BR>is very apparent if I randomly color the polygons — some do not get
<BR>filled. Also, if I do a v.to.rast conversion, I get holes in the newly <BR>created raster.<BR><BR>Does anyone have thoughts on this?<BR><BR>Regards,<BR>Tom<BR><BR><BR>Raffaele Morelli wrote:<BR>><BR>> Cannot display areas, topology not available.<BR>><BR>> I am virtually certain that I was able to successfully import this<BR>> and<BR>> other shapefiles in previous versions of GRASS (GRASS 4.x & GRASS<BR>> 5.x).<BR>> The shapefiles files include the *.dbf, *.shp, and *.shx files,<BR>> which I<BR>> think, should be enough. If I import shapefiles from the US National<BR>> Atlas (created by the USGS), I have no problems. It seems that<BR>> with all<BR>> shapefiles
generated by my office using ArcGIS 9.0, I have this<BR>> problem.<BR>><BR>><BR>> Similar problem with shapefiles generated using ArcGIS, I recently <BR>> posted a topic about (some) areas which are not available after v.in.ogr.<BR>><BR>> Does anyone have suggestions as to what could be going on? <BR>><BR>><BR>> Hamish suggested using 'v.in.ogr -c ' (does not clean polygons) which <BR>> didn't worked for me.<BR>> I partially solved using v.in.ogr without min_area and snap parameters <BR>> but when patching adjacent shapes I do not some polygons lost area.<BR>> Are you experiencing something similar with ArcGIS shapes?<BR>><BR>> regards<BR>> raffaele<BR>> ------------------------------------------------------------------------<BR>><BR>> _______________________________________________<BR>> grassuser mailing list<BR>> grassuser@grass.itc.it<BR>> <A
href="http://grass.itc.it/mailman/listinfo/grassuser" target=_blank>http://grass.itc.it/mailman/listinfo/grassuser</A><BR>> <BR><BR><BR>-- <BR>Thomas E Adams<BR>National Weather Service<BR>Ohio River Forecast Center<BR>1901 South State Route 134<BR>Wilmington, OH 45177<BR><BR>EMAIL: thomas.adams@noaa.gov<BR><BR>VOICE: 937-383-0528<BR>FAX: 937-383-0033<BR><BR><BR><BR>------------------------------<BR><BR>_______________________________________________<BR>grassuser mailing list<BR>grassuser@grass.itc.it<BR><A href="http://grass.itc.it/mailman/listinfo/grassuser" target=_blank>http://grass.itc.it/mailman/listinfo/grassuser</A><BR><BR><BR>End of grassuser Digest, Vol 11, Issue 40<BR>*****************************************</DIV></DIV><BR></DIV></div><br>
<hr size=1>Expecting? Get great news right away with <a href="http://us.rd.yahoo.com/evt=49982/*http://advision.webevents.yahoo.com/mailbeta/newmail_tools.html">email Auto-Check.</a><br>Try the <a href="http://us.rd.yahoo.com/evt=49982/*http://advision.webevents.yahoo.com/mailbeta/newmail_tools.html">Yahoo! Mail Beta.</a></body></html>