<div dir="ltr"><div><div><div><div>THANKS!<br><br></div>I was downloading the SRTM data from earthexplorer and I choose the dataset that appeared under Digita Elevation/SRTM/SRTM 1 Arc-Second Global.<br></div>The srtmgl1 dataset is "hidden" under NASA LPDAAC Collections / Nasa SRTM (SRTM 3) Collections / NASA SRTM3 SRTMGL1.<br><br></div>Back to downloading the data...<br><br></div>Cheers and thanks once again<br>Daniel<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Jan 15, 2015 at 5:59 AM, Markus Metz <span dir="ltr"><<a href="mailto:markus.metz.giswork@gmail.com" target="_blank">markus.metz.giswork@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><span class="">On Wed, Jan 14, 2015 at 11:48 AM, Daniel Victoria<br>
<<a href="mailto:daniel.victoria@gmail.com">daniel.victoria@gmail.com</a>> wrote:<br>
> Hi Markus,<br>
><br>
> I was using r.terraflow because I had converted my voids to negative values<br>
> while the areas outside the country where set to null.<br>
> But the data is proving to be too big and I'll just skip the gap-filling for<br>
> now. Might do it later on, filling the gaps in the import process.<br>
><br>
> The SRTM v3 you pointed out is very good but it's the 3-arcsec resolution.<br>
> I'm working with the 1-arcsec res.<br>
<br>
</span>Then you could try the gap-filled 1 arc sec version of SRTM v3:<br>
<br>
<a href="https://lpdaac.usgs.gov/products/measures_products_table/srtmgl1" target="_blank">https://lpdaac.usgs.gov/products/measures_products_table/srtmgl1</a><br>
<br>
Markus M<br>
<div class="HOEnZb"><div class="h5"><br>
><br>
> Cheers and thanks<br>
> Daniel<br>
><br>
> On Tue, Jan 13, 2015 at 5:40 PM, Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com">markus.metz.giswork@gmail.com</a>><br>
> wrote:<br>
>><br>
>> On Tue, Jan 13, 2015 at 7:52 PM, Daniel Victoria<br>
>> <<a href="mailto:daniel.victoria@gmail.com">daniel.victoria@gmail.com</a>> wrote:<br>
>> ><br>
>> > Now, the thing is, my main interest is to fill the voids in the SRTM<br>
>> > data.<br>
>> > So I though about using r.terraflow to get the filled elevation. Should<br>
>> > I<br>
>> > continue on that path or would it be better to use r.fill.null on a<br>
>> > large<br>
>> > dataset?<br>
>><br>
>> r.terraflow does not fill voids, it fills sinks. Use SRTM v3 [0], this<br>
>> is gap-filled. Using r.fill.null on SRTM v2 will not give you the same<br>
>> quality. r.terraflow has a limit on the number of rows and columns<br>
>> (max 32,767 each). r.watershed also has a limit on the number of rows<br>
>> and columns, but that is max 2,147,483,647 each. In any case, you<br>
>> should use GRASS 7 for such large datasets.<br>
>><br>
>> Markus M<br>
>><br>
>> [0] <a href="https://lpdaac.usgs.gov/products/measures_products_table/srtmgl3" target="_blank">https://lpdaac.usgs.gov/products/measures_products_table/srtmgl3</a><br>
>><br>
>> ><br>
>> > Thanks<br>
>> > Daniel<br>
>> ><br>
>> > On Tue, Jan 13, 2015 at 1:46 PM, Daniel Victoria<br>
>> > <<a href="mailto:daniel.victoria@gmail.com">daniel.victoria@gmail.com</a>><br>
>> > wrote:<br>
>> >><br>
>> >> Charlie,<br>
>> >><br>
>> >> I just downloaded some SRTM 1arc sec. from EarthExplorer. The data is<br>
>> >> supplied in 3 different file types, GeoTIFF, DTED or BIL and they are<br>
>> >> all in<br>
>> >> Integer values (Int16). No floating point elevation values.<br>
>> >><br>
>> >> Cheers<br>
>> >> Daniel<br>
>> >><br>
>> >> On Tue, Jan 13, 2015 at 12:12 PM, Charlie Shobe <<a href="mailto:chsh5846@colorado.edu">chsh5846@colorado.edu</a>><br>
>> >> wrote:<br>
>> >>><br>
>> >>> Hi Daniel,<br>
>> >>><br>
>> >>> I believe that SRTM data does in fact provide floating point elevation<br>
>> >>> values, so you may want to try working with the DEM as type FCELL from<br>
>> >>> the<br>
>> >>> very beginning. I don't know if this will help solve your problem, but<br>
>> >>> the<br>
>> >>> last time I brought in SRTM (1 arc second) data it was type FCELL and<br>
>> >>> contained several decimal places of precision.<br>
>> >>><br>
>> >>> Good luck,<br>
>> >>><br>
>> >>> Charlie<br>
>> >>><br>
>> >>><br>
>> >>> On Jan 13, 2015, at 4:00 AM, Daniel Victoria<br>
>> >>> <<a href="mailto:daniel.victoria@gmail.com">daniel.victoria@gmail.com</a>><br>
>> >>> wrote:<br>
>> >>><br>
>> >>> Stephan, I'll give r.watershed a try and let it run for a couple of<br>
>> >>> days.<br>
>> >>> Thanks<br>
>> >>><br>
>> >>> Thayer, I used r.recode because the person that sent me the data<br>
>> >>> messed<br>
>> >>> up the null values. So in order to fix that I did:<br>
>> >>> 1) use r.external to bring the data to Grass<br>
>> >>> 2) fix null values with r.recode, which was faster than r.null<br>
>> >>><br>
>> >>> But I gave up on that path and since imported the data (r.in.gdal) and<br>
>> >>> fixed the null values with r.null.<br>
>> >>> I'm also using integer values (CELL type) because from what I heard,<br>
>> >>> SRTM<br>
>> >>> does not provide floating point data.<br>
>> >>><br>
>> >>> Now I tried to run r.terraflow but I got a type error which asked me<br>
>> >>> to<br>
>> >>> use r.terraflow.short. Since I don't have that command in my grass<br>
>> >>> instalatin, I converted the SRTM data to float. But that was to no<br>
>> >>> avail<br>
>> >>> since it's now giving me a dimensions type overflow error and asking<br>
>> >>> me to<br>
>> >>> change the dimension_type and recompile.<br>
>> >>><br>
>> >>> I'll give r.watershed a try. If that does not work, I'll see if I can<br>
>> >>> recompile.<br>
>> >>><br>
>> >>> Thanks<br>
>> >>> Daniel<br>
>> >>><br>
>> >>> On Mon, Jan 12, 2015 at 9:22 PM, Thayer Young <<a href="mailto:thayeray@yahoo.com">thayeray@yahoo.com</a>><br>
>> >>> wrote:<br>
>> >>>><br>
>> >>>> With regards to the size of DEM, if you can get the command to run it<br>
>> >>>> may take several days to finish, this is more than twice the size of<br>
>> >>>> the<br>
>> >>>> Washington State DEM that took the authors 33 hours in 2003 (Arge,<br>
>> >>>> Chase,<br>
>> >>>> Halpin, Toma, Urban, Vitter, Wickremesinghe  Efficient Flow<br>
>> >>>> Computation on<br>
>> >>>> Massive Grid Terrain Datasets Geoinformatica Volume 7 Issue 4,<br>
>> >>>> December 2003<br>
>> >>>> Pages 283 - 313 ).  My mid-2014 computer runs at about twice the<br>
>> >>>> speed<br>
>> >>>> quoted in the paper.<br>
>> >>>><br>
>> >>>> Did you r.recode it to make it smaller?  In flat areas you need the<br>
>> >>>> data<br>
>> >>>> contained after the decimal point to allow the flow a chance to make<br>
>> >>>> it down<br>
>> >>>> hill, otherwise you can get streams going no where near where they do<br>
>> >>>> in<br>
>> >>>> real life.  I have never tried running r.terraflow on an integer<br>
>> >>>> raster, I<br>
>> >>>> know that it works on decimal rasters.<br>
>> >>>><br>
>> >>>> -Thayer<br>
>> >>>><br>
>> >>>><br>
>> >>>><br>
>> >>>> ________________________________<br>
>> >>>><br>
>> >>>> Date: Mon, 12 Jan 2015 09:59:40 -0200<br>
>> >>>> From: Daniel Victoria <<a href="mailto:daniel.victoria@gmail.com">daniel.victoria@gmail.com</a>><br>
>> >>>> To: grass <<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>><br>
>> >>>> Subject: [GRASS-user] r.terraflow on massive DEM<br>
>> >>>> Message-ID:<br>
>> >>>><br>
>> >>>> <CA+irsJjf7xhWgU968D+4+xc2YGJ=_<a href="mailto:n4xJtP0_WeU1ZpoRTt9ng@mail.gmail.com">n4xJtP0_WeU1ZpoRTt9ng@mail.gmail.com</a>><br>
>> >>>> Content-Type: text/plain; charset="utf-8"<br>
>> >>>><br>
>> >>>><br>
>> >>>> Hi list,<br>
>> >>>><br>
>> >>>> I'm trying to run r.terraflow on a very large DEM but I wander if<br>
>> >>>> it's<br>
>> >>>> __too large__.<br>
>> >>>><br>
>> >>>> The region dimensions are:<br>
>> >>>> ncol=141114<br>
>> >>>> nrow=140487<br>
>> >>>> Data type = CELL<br>
>> >>>><br>
>> >>>> The map is actually a reclass map of a raster that I imported using<br>
>> >>>> r.external.<br>
>> >>>><br>
>> >>>> When I run r.terraflow I get:<br>
>> >>>><br>
>> >>>> r.terraflow elevation=srtm_brasil@PERMANENT filled=srtm_fill<br>
>> >>>> direction=flowdir swatershed=sink accumulation=flowacc tci=srtm_tci<br>
>> >>>> directory=E:\terraflow_temp<br>
>> >>>><br>
>> >>>> WARNING: raster srtm_brasil is of type CELL_TYPE --you should use<br>
>> >>>> r.terraflow.short<br>
>> >>>> ERROR: [nrows=140487, ncols=141114] dimension_type overflow -- change<br>
>> >>>> dimension_type and recompile<br>
>> >>>> (Mon Jan 12 09:57:18 2015) Command finished (0<br>
>> >>>> sec)<br>
>> >>>><br>
>> >>>> However, I don't have an r.terraflow.short command<br>
>> >>>> I'm running Grass 7.0.0beta4 from OsGeo4W<br>
>> >>>><br>
>> >>>> Cheers<br>
>> >>>> Daniel<br>
>> >>>> -------------- next part --------------<br>
>> >>>> An HTML attachment was scrubbed...<br>
>> >>>> URL:<br>
>> >>>><br>
>> >>>> <<a href="http://lists.osgeo.org/pipermail/grass-user/attachments/20150112/a27c139f/attachment-0001.html" target="_blank">http://lists.osgeo.org/pipermail/grass-user/attachments/20150112/a27c139f/attachment-0001.html</a>><br>
>> >>>><br>
>> >>>> ------------------------------<br>
>> >>>><br>
>> >>>> _______________________________________________<br>
>> >>>> grass-user mailing list<br>
>> >>>> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
>> >>>> <a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
>> >>>><br>
>> >>>> End of grass-user Digest, Vol 105, Issue 16<br>
>> >>>> *******************************************<br>
>> >>>><br>
>> >>>><br>
>> >>>><br>
>> >>>> _______________________________________________<br>
>> >>>> grass-user mailing list<br>
>> >>>> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
>> >>>> <a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
>> >>><br>
>> >>><br>
>> >>> _______________________________________________<br>
>> >>> grass-user mailing list<br>
>> >>> <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
>> >>> <a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
>> >><br>
>> >><br>
>> ><br>
>> ><br>
>> > _______________________________________________<br>
>> > grass-user mailing list<br>
>> > <a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a><br>
>> > <a href="http://lists.osgeo.org/mailman/listinfo/grass-user" target="_blank">http://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
><br>
><br>
</div></div></blockquote></div><br></div>