CGIAR SRTM to xyz using GRASS GIS, with sampling The Next CEO of Stack OverflowGDAL: gdallocationinfo of CGIAR mosaic elevation dataError in using/installing v.digit: GRASS GISCorrect ASTER with SRTM or Topographic MapWhy is this SRTM TIFF not visible in GRASS GIS 6.4.3?Importing SRTM Data into GrassCleaning shapefiles with islands in GRASS GISProblem with SRTM 30 m data while creating watershedUnable to open shapefile with GRASS GIS?Problem with sampling points vectors on r.sim.water Grass GISDelineating coastal subbasins using SAGA GIS or GRASS in QGIS?Problem with GRASS GIS addons

I believe this to be a fraud - hired, then asked to cash check and send cash as Bitcoin

Bartok - Syncopation (1): Meaning of notes in between Grand Staff

Won the lottery - how do I keep the money?

Do I need to write [sic] when a number is less than 10 but isn't written out?

Can a Bladesinger Wizard use Bladesong with a Hand Crossbow?

Why isn't the Mueller report being released completely and unredacted?

Writing differences on a blackboard

How to count occurrences of text in a file?

I want to delete every two lines after 3rd lines in file contain very large number of lines :

Where do students learn to solve polynomial equations these days?

How to check if all elements of 1 list are in the *same quantity* and in any order, in the list2?

A Man With a Stainless Steel Endoskeleton (like The Terminator) Fighting Cloaked Aliens Only He Can See

Is it ever safe to open a suspicious HTML file (e.g. email attachment)?

How to invert MapIndexed on a ragged structure? How to construct a tree from rules?

Is micro rebar a better way to reinforce concrete than rebar?

Why is my new battery behaving weirdly?

Proper way to express "He disappeared them"

Why this way of making earth uninhabitable in Interstellar?

Is it okay to majorly distort historical facts while writing a fiction story?

Why doesn't UK go for the same deal Japan has with EU to resolve Brexit?

Why is the US ranked as #45 in Press Freedom ratings, despite its extremely permissive free speech laws?

Math-accent symbol over parentheses enclosing accented symbol (amsmath)

Is French Guiana a (hard) EU border?

Is it my responsibility to learn a new technology in my own time my employer wants to implement?



CGIAR SRTM to xyz using GRASS GIS, with sampling



The Next CEO of Stack OverflowGDAL: gdallocationinfo of CGIAR mosaic elevation dataError in using/installing v.digit: GRASS GISCorrect ASTER with SRTM or Topographic MapWhy is this SRTM TIFF not visible in GRASS GIS 6.4.3?Importing SRTM Data into GrassCleaning shapefiles with islands in GRASS GISProblem with SRTM 30 m data while creating watershedUnable to open shapefile with GRASS GIS?Problem with sampling points vectors on r.sim.water Grass GISDelineating coastal subbasins using SAGA GIS or GRASS in QGIS?Problem with GRASS GIS addons










0















I have used CGIAR interface to select the tiles, like so:
enter image description here



This resulted in a bunch of *.zip-files, and I have extracted them all into
a directory.
Next, I followed the instructions as https://grasswiki.osgeo.org/wiki/HOWTO_import_SRTM_elevation_data to turn everything into mosaic,
like so:



 #create mosaik (optionally reproject on the fly with -t_srs)
gdalwarp *.tif srtm_mosaik.tif
# import
r.in.gdal input=srtm_mosaik output=srtm_mosaik location=srtm_mosaic_cgiar


Finally, I "connect" to the location srtm_mosaic_cgiar using grass-gui and export the tiles as xyz elevation data, using r.out.xyz (I use the File->Export... drop-down menu).
Now my question is, since this operation appears to be hard on my laptop,
and I can imagine the elevation data of Europe (as I have selected)
is going to be huge, how to "restrict" r.out.xyz to sample (somehow) the
exported data, or how to ascertain beforehand the size of the output?



My final goal is rather simple. I would like to have a file with x,y,z, where x,y refer to a "locus" somehwere in Europe, all the x,ys are mutually sufficiently apart (say, 300-500km), and z means the "elevation"
of the respective point x,y. If there is a misconception on my part at this
stage already, please feel free to correct me. I am new to GIS.










share|improve this question






















  • You can set region in GRASS to limit the operations to a smaller rectangular part of your map. See g.region command or Settings/Computational region from the GRASS GUI.

    – Zoltan
    yesterday















0















I have used CGIAR interface to select the tiles, like so:
enter image description here



This resulted in a bunch of *.zip-files, and I have extracted them all into
a directory.
Next, I followed the instructions as https://grasswiki.osgeo.org/wiki/HOWTO_import_SRTM_elevation_data to turn everything into mosaic,
like so:



 #create mosaik (optionally reproject on the fly with -t_srs)
gdalwarp *.tif srtm_mosaik.tif
# import
r.in.gdal input=srtm_mosaik output=srtm_mosaik location=srtm_mosaic_cgiar


Finally, I "connect" to the location srtm_mosaic_cgiar using grass-gui and export the tiles as xyz elevation data, using r.out.xyz (I use the File->Export... drop-down menu).
Now my question is, since this operation appears to be hard on my laptop,
and I can imagine the elevation data of Europe (as I have selected)
is going to be huge, how to "restrict" r.out.xyz to sample (somehow) the
exported data, or how to ascertain beforehand the size of the output?



My final goal is rather simple. I would like to have a file with x,y,z, where x,y refer to a "locus" somehwere in Europe, all the x,ys are mutually sufficiently apart (say, 300-500km), and z means the "elevation"
of the respective point x,y. If there is a misconception on my part at this
stage already, please feel free to correct me. I am new to GIS.










share|improve this question






















  • You can set region in GRASS to limit the operations to a smaller rectangular part of your map. See g.region command or Settings/Computational region from the GRASS GUI.

    – Zoltan
    yesterday













0












0








0








I have used CGIAR interface to select the tiles, like so:
enter image description here



This resulted in a bunch of *.zip-files, and I have extracted them all into
a directory.
Next, I followed the instructions as https://grasswiki.osgeo.org/wiki/HOWTO_import_SRTM_elevation_data to turn everything into mosaic,
like so:



 #create mosaik (optionally reproject on the fly with -t_srs)
gdalwarp *.tif srtm_mosaik.tif
# import
r.in.gdal input=srtm_mosaik output=srtm_mosaik location=srtm_mosaic_cgiar


Finally, I "connect" to the location srtm_mosaic_cgiar using grass-gui and export the tiles as xyz elevation data, using r.out.xyz (I use the File->Export... drop-down menu).
Now my question is, since this operation appears to be hard on my laptop,
and I can imagine the elevation data of Europe (as I have selected)
is going to be huge, how to "restrict" r.out.xyz to sample (somehow) the
exported data, or how to ascertain beforehand the size of the output?



My final goal is rather simple. I would like to have a file with x,y,z, where x,y refer to a "locus" somehwere in Europe, all the x,ys are mutually sufficiently apart (say, 300-500km), and z means the "elevation"
of the respective point x,y. If there is a misconception on my part at this
stage already, please feel free to correct me. I am new to GIS.










share|improve this question














I have used CGIAR interface to select the tiles, like so:
enter image description here



This resulted in a bunch of *.zip-files, and I have extracted them all into
a directory.
Next, I followed the instructions as https://grasswiki.osgeo.org/wiki/HOWTO_import_SRTM_elevation_data to turn everything into mosaic,
like so:



 #create mosaik (optionally reproject on the fly with -t_srs)
gdalwarp *.tif srtm_mosaik.tif
# import
r.in.gdal input=srtm_mosaik output=srtm_mosaik location=srtm_mosaic_cgiar


Finally, I "connect" to the location srtm_mosaic_cgiar using grass-gui and export the tiles as xyz elevation data, using r.out.xyz (I use the File->Export... drop-down menu).
Now my question is, since this operation appears to be hard on my laptop,
and I can imagine the elevation data of Europe (as I have selected)
is going to be huge, how to "restrict" r.out.xyz to sample (somehow) the
exported data, or how to ascertain beforehand the size of the output?



My final goal is rather simple. I would like to have a file with x,y,z, where x,y refer to a "locus" somehwere in Europe, all the x,ys are mutually sufficiently apart (say, 300-500km), and z means the "elevation"
of the respective point x,y. If there is a misconception on my part at this
stage already, please feel free to correct me. I am new to GIS.







grass srtm






share|improve this question













share|improve this question











share|improve this question




share|improve this question










asked yesterday









IlonpilaajaIlonpilaaja

1105




1105












  • You can set region in GRASS to limit the operations to a smaller rectangular part of your map. See g.region command or Settings/Computational region from the GRASS GUI.

    – Zoltan
    yesterday

















  • You can set region in GRASS to limit the operations to a smaller rectangular part of your map. See g.region command or Settings/Computational region from the GRASS GUI.

    – Zoltan
    yesterday
















You can set region in GRASS to limit the operations to a smaller rectangular part of your map. See g.region command or Settings/Computational region from the GRASS GUI.

– Zoltan
yesterday





You can set region in GRASS to limit the operations to a smaller rectangular part of your map. See g.region command or Settings/Computational region from the GRASS GUI.

– Zoltan
yesterday










1 Answer
1






active

oldest

votes


















0














Since you have started with gdal, why no do the whole process using gdal utilities. I don't see the need for GRASS here (unless you need something more than what you have asked).



I would not create a new raster for all of Europe, rather use the Virtual Raster capability buildvft, then the gdallocationinfo locationinfo should get you the elevation at each X-Y locations. i.e.:
(after downloading and unzipping the tiles as you described)



gdalbuildvrt srtm.vft srtm*.tif
gdallocationinfo srtm.vft <your x coord> <your y coord>


That should do it. You could, of course, write a loop to extract Z values for a list of coord locations.






share|improve this answer























    Your Answer








    StackExchange.ready(function()
    var channelOptions =
    tags: "".split(" "),
    id: "79"
    ;
    initTagRenderer("".split(" "), "".split(" "), channelOptions);

    StackExchange.using("externalEditor", function()
    // Have to fire editor after snippets, if snippets enabled
    if (StackExchange.settings.snippets.snippetsEnabled)
    StackExchange.using("snippets", function()
    createEditor();
    );

    else
    createEditor();

    );

    function createEditor()
    StackExchange.prepareEditor(
    heartbeatType: 'answer',
    autoActivateHeartbeat: false,
    convertImagesToLinks: false,
    noModals: true,
    showLowRepImageUploadWarning: true,
    reputationToPostImages: null,
    bindNavPrevention: true,
    postfix: "",
    imageUploader:
    brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
    contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
    allowUrls: true
    ,
    onDemand: true,
    discardSelector: ".discard-answer"
    ,immediatelyShowMarkdownHelp:true
    );



    );













    draft saved

    draft discarded


















    StackExchange.ready(
    function ()
    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f317101%2fcgiar-srtm-to-xyz-using-grass-gis-with-sampling%23new-answer', 'question_page');

    );

    Post as a guest















    Required, but never shown

























    1 Answer
    1






    active

    oldest

    votes








    1 Answer
    1






    active

    oldest

    votes









    active

    oldest

    votes






    active

    oldest

    votes









    0














    Since you have started with gdal, why no do the whole process using gdal utilities. I don't see the need for GRASS here (unless you need something more than what you have asked).



    I would not create a new raster for all of Europe, rather use the Virtual Raster capability buildvft, then the gdallocationinfo locationinfo should get you the elevation at each X-Y locations. i.e.:
    (after downloading and unzipping the tiles as you described)



    gdalbuildvrt srtm.vft srtm*.tif
    gdallocationinfo srtm.vft <your x coord> <your y coord>


    That should do it. You could, of course, write a loop to extract Z values for a list of coord locations.






    share|improve this answer



























      0














      Since you have started with gdal, why no do the whole process using gdal utilities. I don't see the need for GRASS here (unless you need something more than what you have asked).



      I would not create a new raster for all of Europe, rather use the Virtual Raster capability buildvft, then the gdallocationinfo locationinfo should get you the elevation at each X-Y locations. i.e.:
      (after downloading and unzipping the tiles as you described)



      gdalbuildvrt srtm.vft srtm*.tif
      gdallocationinfo srtm.vft <your x coord> <your y coord>


      That should do it. You could, of course, write a loop to extract Z values for a list of coord locations.






      share|improve this answer

























        0












        0








        0







        Since you have started with gdal, why no do the whole process using gdal utilities. I don't see the need for GRASS here (unless you need something more than what you have asked).



        I would not create a new raster for all of Europe, rather use the Virtual Raster capability buildvft, then the gdallocationinfo locationinfo should get you the elevation at each X-Y locations. i.e.:
        (after downloading and unzipping the tiles as you described)



        gdalbuildvrt srtm.vft srtm*.tif
        gdallocationinfo srtm.vft <your x coord> <your y coord>


        That should do it. You could, of course, write a loop to extract Z values for a list of coord locations.






        share|improve this answer













        Since you have started with gdal, why no do the whole process using gdal utilities. I don't see the need for GRASS here (unless you need something more than what you have asked).



        I would not create a new raster for all of Europe, rather use the Virtual Raster capability buildvft, then the gdallocationinfo locationinfo should get you the elevation at each X-Y locations. i.e.:
        (after downloading and unzipping the tiles as you described)



        gdalbuildvrt srtm.vft srtm*.tif
        gdallocationinfo srtm.vft <your x coord> <your y coord>


        That should do it. You could, of course, write a loop to extract Z values for a list of coord locations.







        share|improve this answer












        share|improve this answer



        share|improve this answer










        answered yesterday









        MichaMicha

        11k1421




        11k1421



























            draft saved

            draft discarded
















































            Thanks for contributing an answer to Geographic Information Systems Stack Exchange!


            • Please be sure to answer the question. Provide details and share your research!

            But avoid


            • Asking for help, clarification, or responding to other answers.

            • Making statements based on opinion; back them up with references or personal experience.

            To learn more, see our tips on writing great answers.




            draft saved


            draft discarded














            StackExchange.ready(
            function ()
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f317101%2fcgiar-srtm-to-xyz-using-grass-gis-with-sampling%23new-answer', 'question_page');

            );

            Post as a guest















            Required, but never shown





















































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown

































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown







            Popular posts from this blog

            រឿង រ៉ូមេអូ និង ហ្ស៊ុយលីយេ សង្ខេបរឿង តួអង្គ បញ្ជីណែនាំ

            Crop image to path created in TikZ? Announcing the arrival of Valued Associate #679: Cesar Manara Planned maintenance scheduled April 17/18, 2019 at 00:00UTC (8:00pm US/Eastern)Crop an inserted image?TikZ pictures does not appear in posterImage behind and beyond crop marks?Tikz picture as large as possible on A4 PageTransparency vs image compression dilemmaHow to crop background from image automatically?Image does not cropTikzexternal capturing crop marks when externalizing pgfplots?How to include image path that contains a dollar signCrop image with left size given

            Romeo and Juliet ContentsCharactersSynopsisSourcesDate and textThemes and motifsCriticism and interpretationLegacyScene by sceneSee alsoNotes and referencesSourcesExternal linksNavigation menu"Consumer Price Index (estimate) 1800–"10.2307/28710160037-3222287101610.1093/res/II.5.31910.2307/45967845967810.2307/2869925286992510.1525/jams.1982.35.3.03a00050"Dada Masilo: South African dancer who breaks the rules"10.1093/res/os-XV.57.1610.2307/28680942868094"Sweet Sorrow: Mann-Korman's Romeo and Juliet Closes Sept. 5 at MN's Ordway"the original10.2307/45957745957710.1017/CCOL0521570476.009"Ram Leela box office collections hit massive Rs 100 crore, pulverises prediction"Archived"Broadway Revival of Romeo and Juliet, Starring Orlando Bloom and Condola Rashad, Will Close Dec. 8"Archived10.1075/jhp.7.1.04hon"Wherefore art thou, Romeo? To make us laugh at Navy Pier"the original10.1093/gmo/9781561592630.article.O006772"Ram-leela Review Roundup: Critics Hail Film as Best Adaptation of Romeo and Juliet"Archived10.2307/31946310047-77293194631"Romeo and Juliet get Twitter treatment""Juliet's Nurse by Lois Leveen""Romeo and Juliet: Orlando Bloom's Broadway Debut Released in Theaters for Valentine's Day"Archived"Romeo and Juliet Has No Balcony"10.1093/gmo/9781561592630.article.O00778110.2307/2867423286742310.1076/enst.82.2.115.959510.1080/00138380601042675"A plague o' both your houses: error in GCSE exam paper forces apology""Juliet of the Five O'Clock Shadow, and Other Wonders"10.2307/33912430027-4321339124310.2307/28487440038-7134284874410.2307/29123140149-661129123144728341M"Weekender Guide: Shakespeare on The Drive""balcony"UK public library membership"romeo"UK public library membership10.1017/CCOL9780521844291"Post-Zionist Critique on Israel and the Palestinians Part III: Popular Culture"10.2307/25379071533-86140377-919X2537907"Capulets and Montagues: UK exam board admit mixing names up in Romeo and Juliet paper"Istoria Novellamente Ritrovata di Due Nobili Amanti2027/mdp.390150822329610820-750X"GCSE exam error: Board accidentally rewrites Shakespeare"10.2307/29176390149-66112917639"Exam board apologises after error in English GCSE paper which confused characters in Shakespeare's Romeo and Juliet""From Mariotto and Ganozza to Romeo and Guilietta: Metamorphoses of a Renaissance Tale"10.2307/37323537323510.2307/2867455286745510.2307/28678912867891"10 Questions for Taylor Swift"10.2307/28680922868092"Haymarket Theatre""The Zeffirelli Way: Revealing Talk by Florentine Director""Michael Smuin: 1938-2007 / Prolific dance director had showy career"The Life and Art of Edwin BoothRomeo and JulietRomeo and JulietRomeo and JulietRomeo and JulietEasy Read Romeo and JulietRomeo and Julieteeecb12003684p(data)4099369-3n8211610759dbe00d-a9e2-41a3-b2c1-977dd692899302814385X313670221313670221