<html>
<head>
<meta content="text/html; charset=utf-8" http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<font size="+1">Hi Irantzu</font>,<br>
I am just trying again with this dataset:<br>
<a 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>
</body>
</html>