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

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

            QGIS export composer to PDF scale the map [closed] Planned maintenance scheduled April 23, 2019 at 23:30 UTC (7:30pm US/Eastern) Announcing the arrival of Valued Associate #679: Cesar Manara Unicorn Meta Zoo #1: Why another podcast?Print Composer QGIS 2.6, how to export image?QGIS 2.8.1 print composer won't export all OpenCycleMap base layer tilesSave Print/Map QGIS composer view as PNG/PDF using Python (without changing anything in visible layout)?Export QGIS Print Composer PDF with searchable text labelsQGIS Print Composer does not change from landscape to portrait orientation?How can I avoid map size and scale changes in print composer?Fuzzy PDF export in QGIS running on macSierra OSExport the legend into its 100% size using Print ComposerScale-dependent rendering in QGIS PDF output

            PDF-ში გადმოწერა სანავიგაციო მენიუproject page