Terrain correction / Orthorectification with gdalMODIS Imagery Orthorectification (Open Source)How to decrease GeoTIFF file size by eliminating or negating the blackfill?Re-projecting from Polar projection to EPSG:4326?Cropping EUDEM with GDAL creates zig-zag terrainSubdividing large files due to gdalwarp - too big image sizeprojecting INSAT3D L1B data from HDF5 to GeotiffOrthorectification in QGIS (ORFEO OTB)gdal: apply orthorectification to ogr datasetError in orthorectification of Radarsat-2 scenes with GDAL and RPC + DEM informationBlank Image After Terrain CorrectionOrthorectification with R
Why don't electron-positron collisions release infinite energy?
A Journey Through Space and Time
How is it possible for user to changed after storage was encrypted? (on OS X, Android)
What would the Romans have called "sorcery"?
Why CLRS example on residual networks does not follows its formula?
TGV timetables / schedules?
Do airline pilots ever risk not hearing communication directed to them specifically, from traffic controllers?
What is the command to reset a PC without deleting any files
Why are 150k or 200k jobs considered good when there are 300k+ births a month?
Prevent a directory in /tmp from being deleted
Compute hash value according to multiplication method
What exactly is the parasitic white layer that forms after iron parts are treated with ammonia?
Japan - Plan around max visa duration
I’m planning on buying a laser printer but concerned about the life cycle of toner in the machine
Why don't electromagnetic waves interact with each other?
Book about a traveler who helps planets in need
Modification to Chariots for Heavy Cavalry Analogue for 4-armed race
What is the offset in a seaplane's hull?
Infinite past with a beginning?
What would happen to a modern skyscraper if it rains micro blackholes?
What do you call a Matrix-like slowdown and camera movement effect?
How is this relation reflexive?
What Brexit solution does the DUP want?
Why is this code 6.5x slower with optimizations enabled?
Terrain correction / Orthorectification with gdal
MODIS Imagery Orthorectification (Open Source)How to decrease GeoTIFF file size by eliminating or negating the blackfill?Re-projecting from Polar projection to EPSG:4326?Cropping EUDEM with GDAL creates zig-zag terrainSubdividing large files due to gdalwarp - too big image sizeprojecting INSAT3D L1B data from HDF5 to GeotiffOrthorectification in QGIS (ORFEO OTB)gdal: apply orthorectification to ogr datasetError in orthorectification of Radarsat-2 scenes with GDAL and RPC + DEM informationBlank Image After Terrain CorrectionOrthorectification with R
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;
Since MODIS Imagery Orthorectification (Open Source) had no answer, I will try to get one here with a similar problem:
I got some raster datasets with GCPs in lat/lon, but without height information. When I use gdal to project/re-project the data, I get an offset of approx. 2km, which I think comes from terrain distortion.
Using ESA's SNAP Toolbox, I get the results I want, but they are using the actual SRTM DEM in the background.
Is there a way to provide gdal with the raster, the GCPs and a DEM to get everything done "at once" (or in several steps)?
So far, my approach (which gives the offset) looks like this:
gdal_translate -gcp pix1 line1 x1 y1 -gcp pix2 line2 x2 y2 (...) input.tif output.tif
The data are SAR raster from ERS-1, which might have something to do with the offset due to the viewing angle of the sensor
gdal terrain orthorectification radar
edited Oct 25 '17 at 20:29
PolyGeo♦
53.9k1781245
asked Oct 25 '17 at 9:45
s6heberns6hebern
619210
add a comment |
Since MODIS Imagery Orthorectification (Open Source) had no answer, I will try to get one here with a similar problem:
I got some raster datasets with GCPs in lat/lon, but without height information. When I use gdal to project/re-project the data, I get an offset of approx. 2km, which I think comes from terrain distortion.
Using ESA's SNAP Toolbox, I get the results I want, but they are using the actual SRTM DEM in the background.
Is there a way to provide gdal with the raster, the GCPs and a DEM to get everything done "at once" (or in several steps)?
So far, my approach (which gives the offset) looks like this:
gdal_translate -gcp pix1 line1 x1 y1 -gcp pix2 line2 x2 y2 (...) input.tif output.tif
The data are SAR raster from ERS-1, which might have something to do with the offset due to the viewing angle of the sensor
gdal terrain orthorectification radar
edited Oct 25 '17 at 20:29
PolyGeo♦
53.9k1781245
asked Oct 25 '17 at 9:45
s6heberns6hebern
619210
For orthorectification I would have a try with OSSIM but it probably requires sensor models too. Offset of 2km feels rather big if the view angle is not very low and/or terrain doesn't have great variation in heights.
– user30184
Oct 25 '17 at 14:59
Since the datasets are SAR images, view angles are indeed relaitvely low and elevation differs as well, because my region of interest is at the northern edge of the Alps. I have had no experience with OSSIM whatsoever, do you have any hints for me?
– s6hebern
Oct 26 '17 at 7:10
Start from trac.osgeo.org/ossim, see from trac.osgeo.org/ossim/wiki/sensors that ERS-SAR is a supported sensor. I fear that OSSIM developers are not active in gis.stackexchange but the mailing list is responsive.
– user30184
Oct 26 '17 at 7:16
Cool, thank you. I'll post here if I get any feasible results.
– s6hebern
Oct 26 '17 at 8:33
add a comment |
Since MODIS Imagery Orthorectification (Open Source) had no answer, I will try to get one here with a similar problem:
I got some raster datasets with GCPs in lat/lon, but without height information. When I use gdal to project/re-project the data, I get an offset of approx. 2km, which I think comes from terrain distortion.
Using ESA's SNAP Toolbox, I get the results I want, but they are using the actual SRTM DEM in the background.
Is there a way to provide gdal with the raster, the GCPs and a DEM to get everything done "at once" (or in several steps)?
So far, my approach (which gives the offset) looks like this:
gdal_translate -gcp pix1 line1 x1 y1 -gcp pix2 line2 x2 y2 (...) input.tif output.tif
The data are SAR raster from ERS-1, which might have something to do with the offset due to the viewing angle of the sensor
gdal terrain orthorectification radar
edited Oct 25 '17 at 20:29
PolyGeo♦
53.9k1781245
asked Oct 25 '17 at 9:45
s6heberns6hebern
619210
Since MODIS Imagery Orthorectification (Open Source) had no answer, I will try to get one here with a similar problem:
I got some raster datasets with GCPs in lat/lon, but without height information. When I use gdal to project/re-project the data, I get an offset of approx. 2km, which I think comes from terrain distortion.
Using ESA's SNAP Toolbox, I get the results I want, but they are using the actual SRTM DEM in the background.
Is there a way to provide gdal with the raster, the GCPs and a DEM to get everything done "at once" (or in several steps)?
So far, my approach (which gives the offset) looks like this:
gdal_translate -gcp pix1 line1 x1 y1 -gcp pix2 line2 x2 y2 (...) input.tif output.tif
The data are SAR raster from ERS-1, which might have something to do with the offset due to the viewing angle of the sensor
gdal terrain orthorectification radar
gdal terrain orthorectification radar
edited Oct 25 '17 at 20:29
PolyGeo♦
53.9k1781245
asked Oct 25 '17 at 9:45
s6heberns6hebern
619210
edited Oct 25 '17 at 20:29
PolyGeo♦
53.9k1781245
asked Oct 25 '17 at 9:45
s6heberns6hebern
619210
edited Oct 25 '17 at 20:29
PolyGeo♦
53.9k1781245
edited Oct 25 '17 at 20:29
PolyGeo♦
53.9k1781245
edited Oct 25 '17 at 20:29
PolyGeo♦
53.9k1781245
53.9k1781245
asked Oct 25 '17 at 9:45
s6heberns6hebern
619210
asked Oct 25 '17 at 9:45
s6heberns6hebern
619210
asked Oct 25 '17 at 9:45
s6heberns6hebern
619210
619210
For orthorectification I would have a try with OSSIM but it probably requires sensor models too. Offset of 2km feels rather big if the view angle is not very low and/or terrain doesn't have great variation in heights.
– user30184
Oct 25 '17 at 14:59
Since the datasets are SAR images, view angles are indeed relaitvely low and elevation differs as well, because my region of interest is at the northern edge of the Alps. I have had no experience with OSSIM whatsoever, do you have any hints for me?
– s6hebern
Oct 26 '17 at 7:10
Start from trac.osgeo.org/ossim, see from trac.osgeo.org/ossim/wiki/sensors that ERS-SAR is a supported sensor. I fear that OSSIM developers are not active in gis.stackexchange but the mailing list is responsive.
– user30184
Oct 26 '17 at 7:16
Cool, thank you. I'll post here if I get any feasible results.
– s6hebern
Oct 26 '17 at 8:33
add a comment |
For orthorectification I would have a try with OSSIM but it probably requires sensor models too. Offset of 2km feels rather big if the view angle is not very low and/or terrain doesn't have great variation in heights.
– user30184
Oct 25 '17 at 14:59
Since the datasets are SAR images, view angles are indeed relaitvely low and elevation differs as well, because my region of interest is at the northern edge of the Alps. I have had no experience with OSSIM whatsoever, do you have any hints for me?
– s6hebern
Oct 26 '17 at 7:10
Start from trac.osgeo.org/ossim, see from trac.osgeo.org/ossim/wiki/sensors that ERS-SAR is a supported sensor. I fear that OSSIM developers are not active in gis.stackexchange but the mailing list is responsive.
– user30184
Oct 26 '17 at 7:16
Cool, thank you. I'll post here if I get any feasible results.
– s6hebern
Oct 26 '17 at 8:33
For orthorectification I would have a try with OSSIM but it probably requires sensor models too. Offset of 2km feels rather big if the view angle is not very low and/or terrain doesn't have great variation in heights.
– user30184
Oct 25 '17 at 14:59
For orthorectification I would have a try with OSSIM but it probably requires sensor models too. Offset of 2km feels rather big if the view angle is not very low and/or terrain doesn't have great variation in heights.
– user30184
Oct 25 '17 at 14:59
Since the datasets are SAR images, view angles are indeed relaitvely low and elevation differs as well, because my region of interest is at the northern edge of the Alps. I have had no experience with OSSIM whatsoever, do you have any hints for me?
– s6hebern
Oct 26 '17 at 7:10
Since the datasets are SAR images, view angles are indeed relaitvely low and elevation differs as well, because my region of interest is at the northern edge of the Alps. I have had no experience with OSSIM whatsoever, do you have any hints for me?
– s6hebern
Oct 26 '17 at 7:10
Start from trac.osgeo.org/ossim, see from trac.osgeo.org/ossim/wiki/sensors that ERS-SAR is a supported sensor. I fear that OSSIM developers are not active in gis.stackexchange but the mailing list is responsive.
– user30184
Oct 26 '17 at 7:16
Start from trac.osgeo.org/ossim, see from trac.osgeo.org/ossim/wiki/sensors that ERS-SAR is a supported sensor. I fear that OSSIM developers are not active in gis.stackexchange but the mailing list is responsive.
– user30184
Oct 26 '17 at 7:16
Cool, thank you. I'll post here if I get any feasible results.
– s6hebern
Oct 26 '17 at 8:33
Cool, thank you. I'll post here if I get any feasible results.
– s6hebern
Oct 26 '17 at 8:33
add a comment |
1 Answer
1
active
oldest
votes
Here is the command I use to do it all at once, plus some other helpful settings to help it run a bit quicker. If needed, you can also specify coordinates to chip the image using the -te flag. This command will chip the image, change the output format, and orthorectify using a DEM in a single step. You can download SRTM here: http://srtm.csi.cgiar.org/srtmdata/
gdalwarp -multi -wo NUM_THREADS=ALL_CPUS -co NUM_THREADS=ALL_CPUS --config GDAL_CACHEMAX 1536 -wm 1536 -rpc -to RPC_DEM=$DEM_FILE -of GTiff -te $aoi_coords $inFile $outFile
answered Apr 3 at 15:32
MikeF22MikeF22
11
New contributor
MikeF22 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.
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%2f259604%2fterrain-correction-orthorectification-with-gdal%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
Here is the command I use to do it all at once, plus some other helpful settings to help it run a bit quicker. If needed, you can also specify coordinates to chip the image using the -te flag. This command will chip the image, change the output format, and orthorectify using a DEM in a single step. You can download SRTM here: http://srtm.csi.cgiar.org/srtmdata/
gdalwarp -multi -wo NUM_THREADS=ALL_CPUS -co NUM_THREADS=ALL_CPUS --config GDAL_CACHEMAX 1536 -wm 1536 -rpc -to RPC_DEM=$DEM_FILE -of GTiff -te $aoi_coords $inFile $outFile
answered Apr 3 at 15:32
MikeF22MikeF22
11
New contributor
MikeF22 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.
add a comment |
Here is the command I use to do it all at once, plus some other helpful settings to help it run a bit quicker. If needed, you can also specify coordinates to chip the image using the -te flag. This command will chip the image, change the output format, and orthorectify using a DEM in a single step. You can download SRTM here: http://srtm.csi.cgiar.org/srtmdata/
gdalwarp -multi -wo NUM_THREADS=ALL_CPUS -co NUM_THREADS=ALL_CPUS --config GDAL_CACHEMAX 1536 -wm 1536 -rpc -to RPC_DEM=$DEM_FILE -of GTiff -te $aoi_coords $inFile $outFile
answered Apr 3 at 15:32
MikeF22MikeF22
11
New contributor
MikeF22 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.
add a comment |
Here is the command I use to do it all at once, plus some other helpful settings to help it run a bit quicker. If needed, you can also specify coordinates to chip the image using the -te flag. This command will chip the image, change the output format, and orthorectify using a DEM in a single step. You can download SRTM here: http://srtm.csi.cgiar.org/srtmdata/
gdalwarp -multi -wo NUM_THREADS=ALL_CPUS -co NUM_THREADS=ALL_CPUS --config GDAL_CACHEMAX 1536 -wm 1536 -rpc -to RPC_DEM=$DEM_FILE -of GTiff -te $aoi_coords $inFile $outFile
answered Apr 3 at 15:32
MikeF22MikeF22
11
New contributor
MikeF22 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.
Here is the command I use to do it all at once, plus some other helpful settings to help it run a bit quicker. If needed, you can also specify coordinates to chip the image using the -te flag. This command will chip the image, change the output format, and orthorectify using a DEM in a single step. You can download SRTM here: http://srtm.csi.cgiar.org/srtmdata/
gdalwarp -multi -wo NUM_THREADS=ALL_CPUS -co NUM_THREADS=ALL_CPUS --config GDAL_CACHEMAX 1536 -wm 1536 -rpc -to RPC_DEM=$DEM_FILE -of GTiff -te $aoi_coords $inFile $outFile
answered Apr 3 at 15:32
MikeF22MikeF22
11
New contributor
MikeF22 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.
edited Apr 3 at 15:38
answered Apr 3 at 15:32
MikeF22MikeF22
11
New contributor
MikeF22 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.
answered Apr 3 at 15:32
MikeF22MikeF22
11
answered Apr 3 at 15:32
MikeF22MikeF22
11
11
New contributor
MikeF22 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.
New contributor
MikeF22 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.
MikeF22 is a new contributor to this site. Take care in asking for clarification, commenting, and answering.
Check out our Code of Conduct.
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%2f259604%2fterrain-correction-orthorectification-with-gdal%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
For orthorectification I would have a try with OSSIM but it probably requires sensor models too. Offset of 2km feels rather big if the view angle is not very low and/or terrain doesn't have great variation in heights.
– user30184
Oct 25 '17 at 14:59
Since the datasets are SAR images, view angles are indeed relaitvely low and elevation differs as well, because my region of interest is at the northern edge of the Alps. I have had no experience with OSSIM whatsoever, do you have any hints for me?
– s6hebern
Oct 26 '17 at 7:10
Start from trac.osgeo.org/ossim, see from trac.osgeo.org/ossim/wiki/sensors that ERS-SAR is a supported sensor. I fear that OSSIM developers are not active in gis.stackexchange but the mailing list is responsive.
– user30184
Oct 26 '17 at 7:16
Cool, thank you. I'll post here if I get any feasible results.
– s6hebern
Oct 26 '17 at 8:33