Vjylku

Can a Hydra make multiple opportunity attacks at once?

Multiplying elements of a list

How does the spell Slow affect freefall?

Can you wish for more wishes from an Efreeti bound to service via an Efreeti Bottle?

Is corrosion inhibitor paste conductive?

If I have Haste cast on me, does it reduce the casting time for my spells that normally take more than a turn to cast?

Is it ethical to apply for a job on someone's behalf?

Have any astronauts or cosmonauts died in space?

Is there a way to pause a running process on Linux systems and resume later?

A cancellation property for permutations

How do I avoid the "chosen hero" feeling?

Coworker is trying to get me to sign his petition to run for office. How to decline politely?

Are encryption algorithms with fixed-point free permutations inherently flawed?

Boss asked me to sign a resignation paper without a date on it along with my new contract

Is there any danger of my neighbor having my wife's signature?

Are all power cords made equal?

Given the mapping a = 1, b = 2, ... z = 26, and an encoded message, count the number of ways it can be decoded

Have the UK Conservatives lost the working majority and if so, what does this mean?

Multiple null checks in Java 8

SQL Server 2017 crashes when backing up because filepath is wrong

What if you do not believe in the project benefits?

Why are `&array` and `array` pointing to the same address?

Do error bars on probabilities have any meaning?

What's the function of the word "ли" in the following contexts?

GDAL: Using gdallocationinfo on a rotated SAR image to get pixel value at Coordinates

Understanding the Map Tile creation propertiesHow can I fix badly specified geographic coordinate system with units in arc-seconds?Rotating rasters using GDAL Geotransform (and ArcGIS Desktop)?Geotransformation for polar stereographic?How do I convert a GeoTiff's neatline metadata from vector format to pixels using GDAL?GeoTiff image not aligned after using gdal_translate and gdalwarpCheck if raster (part of) lies within a pixel - python/gdalExtract a cutline from VRT file in world coordinates?How to calculate the image size knowing its coordinates and pixel size?Transforming Lambert Conformal Conic to EPSG:3857

0

Some information: I must not re-sample the image, the pixels are rotated, and the image opens in Arcmap just fine.

now, there are 15 GCPs on the raster (which make the pixels rotate)

Progress:

1.When I first try to use gdallocationinfo I get:

ERROR 6: No translation for an empty SRS to PROJ.4 format is known.

2.I use gdal_edit -ro -a_srs epsg:4326 "Image", this creates a .aux.xml file with the same name as my image, gdal now reads them both together, image still works on Arcmap.

then I use gdallocationinfo again and get:

ERROR 1: Cannot get geotransform

3.I use gcps2wld "Image" to get the geotransformation parameters, and insert it to the .aux.xml file via python GDAL dataset.SetGeoTransform

now Arcmap can't read the image.
and when I use gdalinfo, I get the following:

Corner Coordinates:
Upper Left ( 0.0000360, -0.0000277) ( 0d 0' 0.13"E, 0d 0' 0.10"S)

Lower Left ( -0.145, 887589.718) ( 0d 8'43.59"W,Invalid angle)

Upper Right ( -0.068,-1152505.040) ( 0d 4' 3.91"W,Invalid angle)

Lower Right ( -0.213, -264915.321) ( 0d12'47.63"W,Invalid angle)

Center ( -0.107, -132457.661) ( 0d 6'23.75"W,Invalid angle)

and if I Ignore that and use gdallocationinfo, coordinates that are supposed to be in the image give me errors for being outside the image borders.

I tried using gdal_edit -a_ullr to set the corner coordinates, but that just norths up the image (no rotation).

how do I set up the image right (without resampling) so I can get the right pixels in the right coordinates?

raster arcmap gdal ground-control sar

share|improve this question

asked Jan 21 at 9:26

YardoodYardood

62

bumped to the homepage by Community♦ 4 mins ago

This question has answers that may be good or bad; the system has marked it active so that they can be reviewed.

add a comment | 

0

Some information: I must not re-sample the image, the pixels are rotated, and the image opens in Arcmap just fine.

now, there are 15 GCPs on the raster (which make the pixels rotate)

Progress:

1.When I first try to use gdallocationinfo I get:

ERROR 6: No translation for an empty SRS to PROJ.4 format is known.

2.I use gdal_edit -ro -a_srs epsg:4326 "Image", this creates a .aux.xml file with the same name as my image, gdal now reads them both together, image still works on Arcmap.

then I use gdallocationinfo again and get:

ERROR 1: Cannot get geotransform

3.I use gcps2wld "Image" to get the geotransformation parameters, and insert it to the .aux.xml file via python GDAL dataset.SetGeoTransform

now Arcmap can't read the image.
and when I use gdalinfo, I get the following:

Corner Coordinates:
Upper Left ( 0.0000360, -0.0000277) ( 0d 0' 0.13"E, 0d 0' 0.10"S)

Lower Left ( -0.145, 887589.718) ( 0d 8'43.59"W,Invalid angle)

Upper Right ( -0.068,-1152505.040) ( 0d 4' 3.91"W,Invalid angle)

Lower Right ( -0.213, -264915.321) ( 0d12'47.63"W,Invalid angle)

Center ( -0.107, -132457.661) ( 0d 6'23.75"W,Invalid angle)

and if I Ignore that and use gdallocationinfo, coordinates that are supposed to be in the image give me errors for being outside the image borders.

I tried using gdal_edit -a_ullr to set the corner coordinates, but that just norths up the image (no rotation).

how do I set up the image right (without resampling) so I can get the right pixels in the right coordinates?

raster arcmap gdal ground-control sar

share|improve this question

asked Jan 21 at 9:26

YardoodYardood

62

bumped to the homepage by Community♦ 4 mins ago

This question has answers that may be good or bad; the system has marked it active so that they can be reviewed.

add a comment | 

Some information: I must not re-sample the image, the pixels are rotated, and the image opens in Arcmap just fine.

now, there are 15 GCPs on the raster (which make the pixels rotate)

Progress:

1.When I first try to use gdallocationinfo I get:

ERROR 6: No translation for an empty SRS to PROJ.4 format is known.

2.I use gdal_edit -ro -a_srs epsg:4326 "Image", this creates a .aux.xml file with the same name as my image, gdal now reads them both together, image still works on Arcmap.

then I use gdallocationinfo again and get:

ERROR 1: Cannot get geotransform

3.I use gcps2wld "Image" to get the geotransformation parameters, and insert it to the .aux.xml file via python GDAL dataset.SetGeoTransform

now Arcmap can't read the image.
and when I use gdalinfo, I get the following:

Corner Coordinates:
Upper Left ( 0.0000360, -0.0000277) ( 0d 0' 0.13"E, 0d 0' 0.10"S)

Lower Left ( -0.145, 887589.718) ( 0d 8'43.59"W,Invalid angle)

Upper Right ( -0.068,-1152505.040) ( 0d 4' 3.91"W,Invalid angle)

Lower Right ( -0.213, -264915.321) ( 0d12'47.63"W,Invalid angle)

Center ( -0.107, -132457.661) ( 0d 6'23.75"W,Invalid angle)

and if I Ignore that and use gdallocationinfo, coordinates that are supposed to be in the image give me errors for being outside the image borders.

I tried using gdal_edit -a_ullr to set the corner coordinates, but that just norths up the image (no rotation).

how do I set up the image right (without resampling) so I can get the right pixels in the right coordinates?

raster arcmap gdal ground-control sar

share|improve this question

asked Jan 21 at 9:26

YardoodYardood

62

share|improve this question

asked Jan 21 at 9:26

YardoodYardood

62

share|improve this question

asked Jan 21 at 9:26

YardoodYardood

62

asked Jan 21 at 9:26

YardoodYardood

62

asked Jan 21 at 9:26

YardoodYardood

62

1 Answer
1

active

oldest

votes

0

Okay, I found the cause of the problem. I got the order of the geotransformation parameters wrong.

the georeferencing transform takes the following form :

adfGeoTransform[0] /* top left x */

adfGeoTransform[1] /* w-e pixel resolution */

adfGeoTransform[2] /* 0 */

adfGeoTransform[3] /* top left y */

adfGeoTransform[4] /* 0 */

adfGeoTransform[5] /* n-s pixel resolution (negative value) */

I hope this helps anyone!

share|improve this answer

answered Jan 22 at 23:23

YardoodYardood

62

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
});


}
});

draft saved

draft discarded

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

Name

Email

Required, but never shown

StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f309331%2fgdal-using-gdallocationinfo-on-a-rotated-sar-image-to-get-pixel-value-at-coordi%23new-answer', 'question_page');
}
);

Post as a guest

Name

Email

Required, but never shown

0

Okay, I found the cause of the problem. I got the order of the geotransformation parameters wrong.

the georeferencing transform takes the following form :

adfGeoTransform[0] /* top left x */

adfGeoTransform[1] /* w-e pixel resolution */

adfGeoTransform[2] /* 0 */

adfGeoTransform[3] /* top left y */

adfGeoTransform[4] /* 0 */

adfGeoTransform[5] /* n-s pixel resolution (negative value) */

I hope this helps anyone!

share|improve this answer

answered Jan 22 at 23:23

YardoodYardood

62

add a comment | 

0

Okay, I found the cause of the problem. I got the order of the geotransformation parameters wrong.

the georeferencing transform takes the following form :

adfGeoTransform[0] /* top left x */

adfGeoTransform[1] /* w-e pixel resolution */

adfGeoTransform[2] /* 0 */

adfGeoTransform[3] /* top left y */

adfGeoTransform[4] /* 0 */

adfGeoTransform[5] /* n-s pixel resolution (negative value) */

I hope this helps anyone!

share|improve this answer

answered Jan 22 at 23:23

YardoodYardood

62

add a comment | 

Okay, I found the cause of the problem. I got the order of the geotransformation parameters wrong.

the georeferencing transform takes the following form :

adfGeoTransform[0] /* top left x */

adfGeoTransform[1] /* w-e pixel resolution */

adfGeoTransform[2] /* 0 */

adfGeoTransform[3] /* top left y */

adfGeoTransform[4] /* 0 */

adfGeoTransform[5] /* n-s pixel resolution (negative value) */

I hope this helps anyone!

share|improve this answer

answered Jan 22 at 23:23

YardoodYardood

62

adfGeoTransform[0] /* top left x */

adfGeoTransform[1] /* w-e pixel resolution */

adfGeoTransform[2] /* 0 */

adfGeoTransform[3] /* top left y */

adfGeoTransform[4] /* 0 */

adfGeoTransform[5] /* n-s pixel resolution (negative value) */

share|improve this answer

answered Jan 22 at 23:23

YardoodYardood

62

answered Jan 22 at 23:23

YardoodYardood

62

answered Jan 22 at 23:23

YardoodYardood

62