<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<font size="+1">Hello all,<br>
after about 10 hours of processing (the computer got to sleep
state over night) still progress 0%. Is it normal on such region
size and core i5 PC?<br>
T. Brunclik<br>
</font><br>
<div class="moz-cite-prefix">Dne 1.6.2015 v 11:44 Tomáš Brunclík
napsal(a):<br>
</div>
<blockquote cite="mid:556C2995.7030401@atlas.cz" type="cite">
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
<font size="+1">Hi Irantzu</font>,<br>
I am just trying again with this dataset:<br>
<a moz-do-not-send="true" class="moz-txt-link-freetext"
href="https://dl.dropboxusercontent.com/u/24062833/ES2003_1m_mb_contour.zip">https://dl.dropboxusercontent.com/u/24062833/ES2003_1m_mb_contour.zip</a><br>
It is cosiderably smaller than yours, with 1m resolution the
region is 5676x4142 pixels, but running r.surf.contour for 15
minutes and still 0% done. Hope it is able to finish before I have
to leave, I have common office PC. In any case I will let you know
what was the result, you can try with the same data for
comparison.<br>
Looking at the data, which are of some inland water body, I see
the most of the region area is covered by land, which is not
covered by contours (there are contours for elevation <= 0
only). Would see how the command copes with that, this could
possibly be source of your problems, should your data also contain
large areas not covered by the contours. It such case a mask may
help.<br>
Regards,<br>
Tomas<br>
<br>
For the record, the commands I used:<br>
<br>
#First created new location using the data as source of its
projection settings and imported the file (from startup GUI),
then:<br>
<br>
g.region vector=ES2003_1m_mb_contour@PERMANENT<br>
<br>
g.region res=1<br>
<br>
g.region -p<br>
<font color="#993399">projection: 1 (UTM)<br>
zone: 10<br>
datum: wgs84<br>
ellipsoid: wgs84<br>
north: 4079792.873656<br>
south: 4074116.9543482<br>
west: 608344.55566446<br>
east: 612486.8266235<br>
nsres: 0.99998578<br>
ewres: 1.00006542<br>
rows: 5676<br>
cols: 4142<br>
cells: 23509992<br>
</font><br>
v.to.rast input=ES2003_1m_mb_contour@PERMANENT
output=contours_rast attribute_column=CONTOUR<br>
<br>
r.surf.contour input=contours_rast output=dmt_rast<br>
<br>
<br>
<br>
<div class="moz-cite-prefix">Dne 31.5.2015 v 20:21 Irantzu Alvarez
napsal(a):<br>
</div>
<blockquote
cite="mid:CAAa2YsX+NSDT2PKV258tH8AbP3nK-RJeLuSdGKcstHjPo-DFcQ@mail.gmail.com"
type="cite">
<div dir="ltr">Hi Markus and Tomas,<br>
Thank you very much for all your support.<br>
<br>
I am using Grass 7.0 in linux and I am working in a high
capacity computer.<br>
<br>
My bathymetry input file contains vector lines, and all of
them have an elevation value, from -2 to -116, with an
interval of 1 m. The file have 300 MB.<br>
<br>
After executing the v.to.rast command, I get a raster, and it
seems that is correct (pixels have the value of the
elevation).<br>
I am trying to change te resolution in order to get the
interpolation. If I choose a resolution of 1 meter, the output
file is exactly the same as the previous one, but if I change
it (for example a resolution of 100 m), it works, but with a
very bis pixexl size. The question is that I need a high
resolution DEM.<br>
This are the commands I have used:<br>
<br>
v.to.rast input=iso1@ras type=line output=ras_ir use=attr
attribute_column=CONTOUR
<br>
r.surf.contour --overwrite input=ras_ir@ras output=mdt_ras
<br>
<br>
Regards,<br>
<br>
Irantzu </div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2015-05-26 10:50 GMT+02:00 Tomáš
Brunclík <span dir="ltr"><<a moz-do-not-send="true"
href="mailto:brunclik@atlas.cz" target="_blank">brunclik@atlas.cz</a>></span>:<br>
<blockquote class="gmail_quote" style="margin:0 0 0
.8ex;border-left:1px #ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000"> <font size="+1">PS:
Just another wild guess ;) <br>
- Does the input vector bathymetry file contain vector
lines? Because when I have tried my test, what I
downloaded to play with contained polygons, not lines
and I had to convert it to lines first (something
like: v.type input=bathymetry output=bathymetry_lines
from_type=boundary to_type=line). That would also
explain your difficulties..</font><br>
Tomas<br>
<br>
<div>Dne 26.5.2015 v 10:38 Tomáš Brunclík napsal(a):<br>
</div>
<div>
<div class="h5">
<blockquote type="cite"> <font size="+1">Hello
Irantzu,<br>
As I said, for anyone to be able to help, please
post</font> the exact commands used (all from
the beginning, not just the last one which fails)
and GRASS version and operating system, just like
Markus did in his post. If you are using GUI
tools, there is the copy button in the command
dialogs, that will allow you to paste the command
parameters to e-mail. Also, address the reply back
to the conference - more people may have ideas.<br>
It occurs to me, that if there is no error, and
the interpolation throws unchanged input raster on
you, that perhaps the input raster contains zero
values instead of no-data (null) values or
something - then there would be nothing to
interpolate... But without seeing the data or the
commands used it is just wild guessing..<br>
Regards,<br>
Tomas<br>
<br>
<div>Dne 25.5.2015 v 13:43 Irantzu Alvarez
napsal(a):<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div>
<div>
<div>
<div>
<div>Hello Tomas,<br>
</div>
As you said, I think that the problem
could be in the volume of data. My
contours are of 300 MB and if I set a
region with resolution of 1 m (the one
I would like to perform) the output
raster has over 3 billion cells.
Although I have a high capacity
computer, grass doen't compute the
surface. With a resolution of 10 m my
region setting are:<br>
g.region
-p
<br>
projection: 1 (UTM)<br>
zone: 30<br>
datum: etrs89<br>
ellipsoid: grs80<br>
north: 4818930<br>
south: 4791480<br>
west: 486510<br>
east: 601280<br>
nsres: 10<br>
ewres: 10<br>
rows: 2745<br>
cols: 11477<br>
cells: 31504365<br>
<br>
</div>
I am trying to generalize contours, in
order to reduce the contour file, but I
don't know which algorithm is better in
that case.<br>
</div>
After that I will try to use diferrent
resolutions to see if that is the problem.<br>
<br>
</div>
All the best,<br>
</div>
irantzu<br>
<br>
<br>
<br>
</div>
<div class="gmail_extra"><br>
<div class="gmail_quote">2015-05-25 10:12
GMT+02:00 Tomáš Brunclík <span dir="ltr"><<a
moz-do-not-send="true"
href="mailto:brunclik@atlas.cz"
target="_blank">brunclik@atlas.cz</a>></span>:<br>
<blockquote class="gmail_quote"
style="margin:0 0 0 .8ex;border-left:1px
#ccc solid;padding-left:1ex">
<div bgcolor="#FFFFFF" text="#000000"> <font
size="+1">Hello,<br>
I tried to process bathymetry lines
with negative elevation and it just
worked in GRASS 7.0.1svn r65198
@Linux. Maybe try if it works with
coarser resolution or region subset to
test if the data volume is the
problem, or if your version of GRASS
really has problem with negative
elevation values. Also try to look if
the input raster contains expected
values, if you did not check that
already.<br>
If the problem persists, post your
findings with log of the commands you
used, your GRASS version and platform.
<br>
If it seems a problem of the GRASS
module, and you are using a recent
version of GRASS, it would be better
to fill a bug, instructions at <a
moz-do-not-send="true"
href="http://grass.osgeo.org/development/bug-tracking/"
target="_blank">http://grass.osgeo.org/development/bug-tracking/</a>
<br>
Regards,<br>
Tomas Brunclik<br>
</font><br>
<div>Dne 18.5.2015 v 16:47 Irantzu
Alvarez napsal(a):<br>
</div>
<blockquote type="cite">
<div>
<div>
<div dir="ltr">
<div> <span><br>
</span> </div>
<div>
<div>
<p>I am trying to create a
DEM from a bathymetric
contour lines. Contours
are in vector format and
the elevation value is
negative (integers). To
create the surface I have
converted the contours in
raster (<a
moz-do-not-send="true"
href="http://v.to"
target="_blank">v.to</a>
rast) and then I have used
r.surf.contour, but It
doesn't work (there is no
error but the resulting
output raster<span
lang="en"><span></span>
<span>is exactly the
same </span>as the
previous one, without
any interpolation)</span>.
<br>
</p>
<p>This are my region
settings:</p>
<p> </p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">g.region
-p -a
raster=contours@iiaa
nsres=1 ewres=1<span>
</span></span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">projection:
1 (UTM)</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">zone:<span>
</span>30</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">datum:<span>
</span>etrs89</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">ellipsoid:<span>
</span>grs80</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">north:<span>
</span>4818927</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">south:<span>
</span>4791484</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">west:<span>
</span>486512</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">east:<span>
</span>601275</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">nsres:<span>
</span>1</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">ewres:<span>
</span>1</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">rows:<span>
</span>27443</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">cols:<span>
</span>114763</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span
style="font-family:"Cambria","serif"" lang="EN-US">cells:<span>
</span><a
moz-do-not-send="true"
href="tel:3149441009"
value="+13149441009"
target="_blank">3149441009</a></span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><br>
<span
style="font-family:"Cambria","serif""
lang="EN-US"></span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><br>
<span>I need a DEM of very
high resolution of a
wide area. I am using a
computer with large
capacities for it.</span></p>
<p class="MsoNormal"
style="margin-bottom:0.0001pt"><span><br>
</span></p>
Could someone help me
please? </div>
</div>
</div>
<br>
<fieldset></fieldset>
<br>
</div>
</div>
<pre>_______________________________________________
grass-user mailing list
<a moz-do-not-send="true" href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a>
<a moz-do-not-send="true" href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a></pre>
</blockquote>
<br>
</div>
<br>
_______________________________________________<br>
grass-user mailing list<br>
<a moz-do-not-send="true"
href="mailto:grass-user@lists.osgeo.org"
target="_blank">grass-user@lists.osgeo.org</a><br>
<a moz-do-not-send="true"
href="http://lists.osgeo.org/mailman/listinfo/grass-user"
target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</blockquote>
</div>
<br>
</div>
</blockquote>
<br>
<br>
<fieldset></fieldset>
<br>
<pre>_______________________________________________
grass-user mailing list
<a moz-do-not-send="true" href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a>
<a moz-do-not-send="true" href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a></pre>
</blockquote>
<br>
</div>
</div>
</div>
<br>
_______________________________________________<br>
grass-user mailing list<br>
<a moz-do-not-send="true"
href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
<a moz-do-not-send="true"
href="http://lists.osgeo.org/mailman/listinfo/grass-user"
target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</blockquote>
</div>
<br>
</div>
</blockquote>
<br>
<br>
<fieldset class="mimeAttachmentHeader"></fieldset>
<br>
<pre wrap="">_______________________________________________
grass-user mailing list
<a class="moz-txt-link-abbreviated" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="http://lists.osgeo.org/mailman/listinfo/grass-user">http://lists.osgeo.org/mailman/listinfo/grass-user</a></pre>
</blockquote>
<br>
</body>
</html>