[GRASS-user] Landsat SLC-Off gap filling
Hi - Is there anyone out there who has worked on an algorithm to fill gaps in Landsat ETM+ SLC-Off data? Now that the Landsat 7 archive is freely available I am looking for an open source approach to fill the gaps in the SLC-Off data and GRASS seems like a good option. I'm most interested in multi-date approaches but at this stage I'm looking at single-date approaches as well. Any ideas? Ned ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Landsat SLC-Off gap filling
+1 I had just downloaded new Landsat 7 SLC-Off data. It still useable for land use change analysis. A grass tool would be very helpful. cheers, maning On 10/20/08, Ned Horning [EMAIL PROTECTED] wrote: Hi - Is there anyone out there who has worked on an algorithm to fill gaps in Landsat ETM+ SLC-Off data? Now that the Landsat 7 archive is freely available I am looking for an open source approach to fill the gaps in the SLC-Off data and GRASS seems like a good option. I'm most interested in multi-date approaches but at this stage I'm looking at single-date approaches as well. Any ideas? Ned ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user -- |-|--| | __.-._ |Ohhh. Great warrior. Wars not make one great. -Yoda | | '-._7' |Freedom is still the most radical idea of all -N.Branden| | /'.-c |Linux registered user #402901, http://counter.li.org/ | | | /T |http://esambale.wikispaces.com| | _)_/LI |-|--| ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Landsat SLC-Off gap filling
Ned Horning wrote: Hi - Is there anyone out there who has worked on an algorithm to fill gaps in Landsat ETM+ SLC-Off data? Now that the Landsat 7 archive is freely available I am looking for an open source approach to fill the gaps in the SLC-Off data and GRASS seems like a good option. I'm most interested in multi-date approaches but at this stage I'm looking at single-date approaches as well. Any ideas? Some ideas here: http://grass.osgeo.org/wiki/MODIS#Removing_holes both single scene interpolation (r.fillnulls) and rolling composite using a time series (r.patch) methods are explained. Hamish __ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: examples of Web GIS interaction between citizens and governments
Marco: If an owner wants to offer (enter) a parcel, he/she also can do it using the web interface. In the same web (http://www.bantegal.com/sitegal ), the fourth option on the left is ofrecer parcelas (offer parcels). Once selected, this option guides you through several steps to locate the parcel (you could try: just select a province, municipality and parish, click on seguir --go on-- and a webmap will appear on which you can interactively click on the parcel you would like to enter). A few more clicks on seguir will show some cadastre information, a reference rent to be perceived (a minimum of 20 euros a year or 10 euros per hectare-year), and other. For a parcel to be correctly introduced in the system, the owner must register (a registration page would appear near the end of the process). Once the parcel is introduced, the administrative process is started and several days later the owner would be required to sign a contract on paper. By the way, some of the parcels already introduced in the system were managed by the administration, and these were introduced by public officials. 2008/10/18 M. Fioretti [EMAIL PROTECTED] On Fri, Oct 17, 2008 10:07:17 AM +0200, Eduardo corbelle wrote: Marco: Although it is not exactly what you asked for, I think the Galician Land Information System ( http://www.bantegal.com/sitegal ) recently created in Spain deserves a look (unfortunately it doesn't have an English version). Eduardo, yes, thanks, this is an interesting use case. How do people or offices who want to _lease_ unused lands to farmers enter those lands into the system, and who are they? Are they only public officers? Thanks, Marco -- Your own civil rights and the quality of your life heavily depend on how software is used *around* you:http://digifreedom.net/node/84 ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
RE: [GRASS-user] How to install Grass into Ubuntu
RE: [GRASS-user] How to install Grass into Ubuntu sudo apt-get install grass The offical Ubuntu repositories are a bit outdated. A more updated repo is the one from Jachym: http://www.les-ejk.cz/ubuntu/ The Gutsy one has Grass 6.3.0, updated to May 2008. Im working with Ubuntu 7.10 Gutsy and in this repositories not is available the version 6.3.0, the laster version in this repositories is 6.2.2 (http://www.les-ejk.cz/ubuntu/). The people who working in Ubuntu and want to work with other version have to compile it. It's not a so hard task. http://grass.osgeo.org/wiki/Compile_and_Install If you want newer versions... you need to compile it. It's not a so hard task, and we can give advices. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user _ ¿Sigue el calor? Consulta MSN El tiempo http://eltiempo.es.msn.com/___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: grass-user Digest, Vol 30, Issue 46
(r.fillnulls) and rolling composite using a time series (r.patch) methods are explained. Hamish __ Do You Yahoo!? Tired of spam? Yahoo! Mail has the best spam protection around http://mail.yahoo.com -- Message: 8 Date: Mon, 20 Oct 2008 09:59:03 +0200 From: Eduardo corbelle [EMAIL PROTECTED] Subject: [GRASS-user] Re: examples of Web GIS interaction between citizensand governments To: M. Fioretti [EMAIL PROTECTED] Cc: grass-user@lists.osgeo.org Message-ID: [EMAIL PROTECTED] Content-Type: text/plain; charset=iso-8859-1 Marco: If an owner wants to offer (enter) a parcel, he/she also can do it using the web interface. In the same web (http://www.bantegal.com/sitegal ), the fourth option on the left is ofrecer parcelas (offer parcels). Once selected, this option guides you through several steps to locate the parcel (you could try: just select a province, municipality and parish, click on seguir --go on-- and a webmap will appear on which you can interactively click on the parcel you would like to enter). A few more clicks on seguir will show some cadastre information, a reference rent to be perceived (a minimum of 20 euros a year or 10 euros per hectare-year), and other. For a parcel to be correctly introduced in the system, the owner must register (a registration page would appear near the end of the process). Once the parcel is introduced, the administrative process is started and several days later the owner would be required to sign a contract on paper. By the way, some of the parcels already introduced in the system were managed by the administration, and these were introduced by public officials. 2008/10/18 M. Fioretti [EMAIL PROTECTED] On Fri, Oct 17, 2008 10:07:17 AM +0200, Eduardo corbelle wrote: Marco: Although it is not exactly what you asked for, I think the Galician Land Information System ( http://www.bantegal.com/sitegal ) recently created in Spain deserves a look (unfortunately it doesn't have an English version). Eduardo, yes, thanks, this is an interesting use case. How do people or offices who want to _lease_ unused lands to farmers enter those lands into the system, and who are they? Are they only public officers? Thanks, Marco -- Your own civil rights and the quality of your life heavily depend on how software is used *around* you:http://digifreedom.net/node/84 -- next part -- An HTML attachment was scrubbed... URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20081020/c19ef400/attachment-0001.html -- Message: 9 Date: Mon, 20 Oct 2008 10:10:32 +0200 From: Jhon Ortiz [EMAIL PROTECTED] Subject: RE: [GRASS-user] How to install Grass into Ubuntu To: grass-user@lists.osgeo.org, [EMAIL PROTECTED], [EMAIL PROTECTED] Message-ID: [EMAIL PROTECTED] Content-Type: text/plain; charset=iso-8859-1 RE: [GRASS-user] How to install Grass into Ubuntu sudo apt-get install grass The offical Ubuntu repositories are a bit outdated. A more updated repo is the one from Jachym: http://www.les-ejk.cz/ubuntu/ The Gutsy one has Grass 6.3.0, updated to May 2008. Im working with Ubuntu 7.10 Gutsy and in this repositories not is available the version 6.3.0, the laster version in this repositories is 6.2.2 (http://www.les-ejk.cz/ubuntu/). The people who working in Ubuntu and want to work with other version have to compile it. It's not a so hard task. http://grass.osgeo.org/wiki/Compile_and_Install If you want newer versions... you need to compile it. It's not a so hard task, and we can give advices. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user _ �Sigue el calor? Consulta MSN El tiempo http://eltiempo.es.msn.com/ -- next part -- An HTML attachment was scrubbed... URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20081020/1507b616/attachment.html -- ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user End of grass-user Digest, Vol 30, Issue 46 ** This message was sent using IMP, the Internet Messaging Program. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to install Grass into Ubuntu
Hi Matt, it depeds who are these users. If you run grass as you, you can manage mapsets in your home folders without having to run it with sudo... If you want to work on others' homes you need the sudo permissions, as any other operation. 2008/10/20 Matt B [EMAIL PROTECTED]: Hello list, I'm using Grass on ubuntu 8.04. Weird thing is that unless I run it sudo grass I get a bunch of errors (and it's unusable) related to grass not having the correct permissions to write to where I keep the data (ie the users home directory). It's not a show stoppper at the moment and only slightly annoying. I'm not sure why this is but I have some ideas to hunt down when I get time. I used the standard ubuntu install and assume its something to do with the ubuntu setup script. Matt On Mon, Oct 20, 2008 at 7:10 PM, Jhon Ortiz [EMAIL PROTECTED] wrote: RE: [GRASS-user] How to install Grass into Ubuntu sudo apt-get install grass The offical Ubuntu repositories are a bit outdated. A more updated repo is the one from Jachym: http://www.les-ejk.cz/ubuntu/ The Gutsy one has Grass 6.3.0, updated to May 2008. Im working with Ubuntu 7.10 Gutsy and in this repositories not is available the version 6.3.0, the laster version in this repositories is 6.2.2 (http://www.les-ejk.cz/ubuntu/). The people who working in Ubuntu and want to work with other version have to compile it. It's not a so hard task. http://grass.osgeo.org/wiki/Compile_and_Install If you want newer versions... you need to compile it. It's not a so hard task, and we can give advices. ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user Ahora llévate lo mejor de MSN y Windows Live, en tu móvil ___ 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 mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Re: v.digit-grass64 (Martin Landa)
hello Grass users-group To get contibution areas for the sewer system of the city I have converted the base-map (DXF-DWG format) into grass63 and digitized it (75 km2) with v.digit tcktk. every line is green with a blue centroide. Then export them to esri shape-files and import the shape-files in InfoWorks. every thing seems OK. The roofs the streets etc. looks like areas island. and line are closed polygons I suppost v.digit (tcltk grass63) is working well (not buggy). Joop Goedbloed ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] getting the number of occurences of different classes in a moving window neighbourhood
Hi list, I'm looking for a way to extract the occurence (number of pixels of a certain class in a moving window) of a certain class in a moving window neighbourhood and write this to a map. r.neighbors only provides summarized versions of what I need. Is there a way to calculate this basic statistic without resorting to too much coding or exporting my maps into R or matlab? Kind regards, Koen ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Geotiff export error (extra row and column)
Hi list, I've recently noticed that when exporting maps as Geotiff GDAL adds an extra row and column to the map. Within grass this data is not there. For example my map that measures 4509 rows and 4359 columns gets exported as a tiff of size 4510 by 4360. The added row and column are empty and can be removed easily but it still is an inconvenience and it should not be there in the first place. Any idea what causes this error and how it can be solved? Kind regards, Koen ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Geotiff export error (extra row and column)
Hi, it sounds a little bit like the problem which I had with PNG export http://trac.osgeo.org/grass/ticket/138 but I have no idea what is the problem, sorry. Alex Hufkens Koen schrieb: Hi list, I've recently noticed that when exporting maps as Geotiff GDAL adds an extra row and column to the map. Within grass this data is not there. For example my map that measures 4509 rows and 4359 columns gets exported as a tiff of size 4510 by 4360. The added row and column are empty and can be removed easily but it still is an inconvenience and it should not be there in the first place. Any idea what causes this error and how it can be solved? Kind regards, Koen ___ 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] generalize a map of points
2008/10/19 Wolf Bergenheim [EMAIL PROTECTED]: José María, If possible, I think it would be useful feature the following: - Where two labels overlap, comparing an attribute to decide which one stays For example: on a map of cities you can use the number of inhabitants. If the city A has more inhabitants than the city B, maintaining the label of the city A. I have just now committed a change in v.label.sa that will do exactly this. See the new overlap parameter. I would appreciate it if you could test it, and give some feedback on how it is working for you. Thank you, --Wolf Yes, I am doing the tests right now. It's very good for me to working with you in some way. So, thank you. José María -- :3 ) Wolf Bergenheim ( 8: ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
[GRASS-user] Correcting Junctions
Dear all: I have two vector maps, one with road segments, and one with junction points. But I know the Junctions aren't 100% correct. Sometimes there are intersecting road segments. I want to find out which nodes in the junctions map are unnecessary (if any), and which are missing. How can I accomplish that? Any help appreciated. My approaches so far were to use v.net (nodes) on the road segments to create correct segments and nodes in a new layer. My understanding is that the new layer now contains correct nodes and line segments. Then I tried to use v.overlay to find nodes in my nodes map that overlap nodes in the newly created map. But the result was empty. I also tried to use a SQL Query on the new map to select all cat that aren't in the junction map, but that didn't work either. I'm quite new to grass, and unfortunately I'm stuck now. Could you please give some advice? Thanks in advance Dennis ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] getting the number of occurences of different classes in a moving window neighbourhood
Hello, You can use r.mapcalc to perform any calculations based on cell surrounding cell values/presence. See more information in r.mapcalc help page [1] or take a look at r.lake [2] help page to see how cell value is calculated based on neighborhood cell values with mapcalc. WBR, Maris. 1. http://grass.osgeo.org/grass64/manuals/html64_user/r.mapcalc.html 2. http://grass.osgeo.org/grass64/manuals/html64_user/r.lake.html 2008/10/20, Hufkens Koen [EMAIL PROTECTED]: Hi list, I'm looking for a way to extract the occurence (number of pixels of a certain class in a moving window) of a certain class in a moving window neighbourhood and write this to a map. r.neighbors only provides summarized versions of what I need. Is there a way to calculate this basic statistic without resorting to too much coding or exporting my maps into R or matlab? Kind regards, Koen ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] How to install Grass into Ubuntu
Matt B wrote: I'm using Grass on ubuntu 8.04. Weird thing is that unless I run it sudo grass I get a bunch of errors (and it's unusable) related to grass not having the correct permissions to write to where I keep the data (ie the users home directory). It's not a show stoppper at the moment and only slightly annoying. I'm not sure why this is but I have some ideas to hunt down when I get time. I used the standard ubuntu install and assume its something to do with the ubuntu setup script. Note that GRASS won't let you select a mapset as the current mapset (where new files are stored) unless you own it. Write permission isn't sufficient. If you are creating a location which is to be shared by multiple users, you either need to create a mapset directory for each user, owned by the user, or grant all such users write permission on the location directory so that they can create their own mapset directory (which they will own). -- Glynn Clements [EMAIL PROTECTED] ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] Re: v.digit-grass64 (Martin Landa)
Hi, 2008/10/20 Joop Goedbloed [EMAIL PROTECTED]: To get contibution areas for the sewer system of the city I have converted the base-map (DXF-DWG format) into grass63 and digitized it (75 km2) with v.digit tcktk. every line is green with a blue centroide. Then export them to esri shape-files and import the shape-files in InfoWorks. every thing seems OK. The roofs the streets etc. looks like areas island. and line are closed polygons I suppost v.digit (tcltk grass63) is working well (not buggy). as far as I remember from your screenshot: tcl/tk v.digit --- green color for boundaries which shares two areas, orange for boundaries which shares only one area (left/right) wx vdigit -- green for one area / orange for two shared areas in tcl/tk v.digit I would not expect that outline boundaries will be green (i.e. that shares two areas) in wx digitizer I would expect that boundaries between area and isle will be green not orange (i.e. shares only one valid area) Anyway it's just the way how vector features are visually interpreted based on their topology. When you export date e.g. to ESRI Shapefile result will be the same. 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
Re: [GRASS-user] generalize a map of points
2008/10/19 Wolf Bergenheim [EMAIL PROTECTED]: José María, If possible, I think it would be useful feature the following: - Where two labels overlap, comparing an attribute to decide which one stays For example: on a map of cities you can use the number of inhabitants. If the city A has more inhabitants than the city B, maintaining the label of the city A. I have just now committed a change in v.label.sa that will do exactly this. See the new overlap parameter. I would appreciate it if you could test it, and give some feedback on how it is working for you. After some drawbacks to update GRASS, I have achieved some results. - The command d.labels does not show anything. It seems to run normally, without error, but the active monitor does not show anything. - Working with the GIS Manager , the labels are displayed well. However, after exporting the Display Map to raster labels disappear from the image (JPG, PPM, and others). - I got maps with labels using ps.map. See files attached. labels_sa_new.pdf : show labels obtained with the command: $ v.label.sa map=ciudades_region labels=ciudades_v_label_sa_new column=LOCALIDAD font=arial overlap=inhab labels.pdf : show labels obtained with the command: $ v.label map=ciudades_region labels=ciudades_v.label column=LOCALIDAD font=arial The files in $MAPSET/paint/labels have the same number of tags. No label is removed by v.label.sa. Please tell me what else I can do, or what other information I can provide it. Thank you José María Thank you, --Wolf -- :3 ) Wolf Bergenheim ( 8: labels.pdf Description: Adobe PDF document labels_sa_new.pdf Description: Adobe PDF document ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: Re: [GRASS-user] grass
I am getting the exact same error. Windows XP, SP2 and Grass 6.3. 4GB memory. 2147 rows, 1984 columns. All I did was add the raster, draw it and try and run terraflow, both MFD and SFD give the same error. Thanks for any help On Tue, 2008-08-26 at 16:09 +0200, Christian Herold wrote: Assertion failed: nrows * ncols == str-stream_len(), file grass2str.h, line 144 Any suggestions? Yes: tell us which GRASS version and which operating system... Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ -- View this message in context: http://www.nabble.com/Re%3A-grass-tp19217641p19974001.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] Re: Contents of grass-user digest... HOW TO INSTALL GRASS
Hi! (Assuming you work with Ubuntu Hardy Heron 8.04) I would remove everything related with grass and re-install again (using Jachym's repository [1]). Open a terminal and... 1. execute the following command: sudo apt-get purge grass 2. According to [1] edit the sources.list file which is under /etc/apt: sudo gedit /etc/apt/sources.list 3. Add the 2 lines in the end of the file for example (although you can add them anywhere else) as procided in [1]: deb http://les-ejk.cz/ubuntu hardy multiverse deb-src http://les-ejk.cz/ubuntu hardy multiverse 4. Save and close the sources.list file 5. Get Jachym's gpg signature wget -q http://les-ejk.cz/pgp/jachym_cepicky-gpg.pub -O - | sudo apt-key add - 6. Update and Upgrade sudo apt-get update sudo apt-get upgrade 7. install grass with: sudo apt-get install grass 8. Launch grass with: grass -gui Good Luck, Nikos ... [1] http://www.les-ejk.cz/ubuntu/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
r.terraflow error - Re: Re: [GRASS-user] grass
Hi Laura, due you have any clue about below error? thanks Markus On Mon, Oct 20, 2008 at 7:35 PM, Tim Wade [EMAIL PROTECTED] wrote: I am getting the exact same error. Windows XP, SP2 and Grass 6.3. 4GB memory. 2147 rows, 1984 columns. All I did was add the raster, draw it and try and run terraflow, both MFD and SFD give the same error. Thanks for any help On Tue, 2008-08-26 at 16:09 +0200, Christian Herold wrote: Assertion failed: nrows * ncols == str-stream_len(), file grass2str.h, line 144 Any suggestions? Yes: tell us which GRASS version and which operating system... Markus ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user ___ -- View this message in context: http://www.nabble.com/Re%3A-grass-tp19217641p19974001.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 -- Open Source Geospatial Foundation http://www.osgeo.org/ http://www.grassbook.org/ ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: [GRASS-user] generalize a map of points
On 20.10.2008 20:31, José María Michia wrote: After some drawbacks to update GRASS, I have achieved some results. - The command d.labels does not show anything. It seems to run normally, without error, but the active monitor does not show anything. - Working with the GIS Manager , the labels are displayed well. However, after exporting the Display Map to raster labels disappear from the image (JPG, PPM, and others). - I got maps with labels using ps.map. See files attached. labels_sa_new.pdf : show labels obtained with the command: $ v.label.sa map=ciudades_region labels=ciudades_v_label_sa_new column=LOCALIDAD font=arial overlap=inhab Hmm that is really strange... Could you send me that city map (off list), so that I can test a bit. I did try it and it worked for me. Also caopy.paste all commands you used to generate the maps. Thanks for testing, --Wolf -- :3 ) Wolf Bergenheim ( 8: ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user
Re: r.terraflow error - Re: Re: [GRASS-user] grass
Markus Neteler wrote: due you have any clue about below error? I am getting the exact same error. Windows XP, SP2 and Grass 6.3. 4GB memory. 2147 rows, 1984 columns. All I did was add the raster, draw it and try and run terraflow, both MFD and SFD give the same error. Assertion failed: nrows * ncols == str-stream_len(), file grass2str.h, line 144 It might be an EOL issue. In 6.x, r.terraflow/Makefile has its own linking rules, and they don't include $(FMODE_OBJ). The file is opened with open() followed by fdopen(). The fdopen() has the b modifier, but I don't know if that works for fdopen() (the MSDN documentation implies that it should, but that may just be because the fopen() documentation was copied verbatim). In 7.0, r.terraflow uses the normal linking rules from Compile.make (the trade-off is that it doesn't build the short version, only the float version). -- Glynn Clements [EMAIL PROTECTED] ___ grass-user mailing list grass-user@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/grass-user