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
I have used CGIAR interface to select the tiles, like so:
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,y
s 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
add a comment |
I have used CGIAR interface to select the tiles, like so:
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,y
s 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
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
add a comment |
I have used CGIAR interface to select the tiles, like so:
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,y
s 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
I have used CGIAR interface to select the tiles, like so:
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,y
s 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
grass srtm
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
add a comment |
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
add a comment |
1 Answer
1
active
oldest
votes
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.
add a comment |
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
);
);
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
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.
add a comment |
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.
add a comment |
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.
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.
answered yesterday
MichaMicha
11k1421
11k1421
add a comment |
add a comment |
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.
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
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
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
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
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