To project, or not to project? Extracting raster values with RProjecting shp from GK3 to UTM32 with...
Crop Image to Circle
Use GPLv3 library in a closed system (no software distribution)
Meaning of “Bulldog drooled courses through his jowls”
Cutting a 4.5m long 2x6 in half with a circular saw
What is gerrymandering called if it's not the result of redrawing districts?
Raise Error Concatenation in SQL Server
Why is this a missed win?
Is it allowed to let the engine of an aircraft idle without a pilot in the plane. (For both helicopters and aeroplanes)
How to pronounce correctly [b] and [p]? As well [t]/[d] and [k]/[g]
Is the phrase “You are requested” polite or rude?
How would a race of humanoids with tails design [vehicle] seats?
Why doesn't hot charcoal glow blue?
Minimum number of turns to capture all pieces in Checkers
Why it is a big deal whether or not Adam Schiff talked to the whistleblower?
Is it okay to request a vegetarian only microwave at work ? If, yes, what's the proper way to do it?
Why doesn't English employ an H in front of the name Ares?
I don't want my ls command in my script to print results on screen
Grade changes with auto grader
Rules on "Pets on shoulder"
Can I use Oko's ability targetting a creature with protection from green?
Did I Traumatize My Puppy?
How were Kurds involved (or not) in the invasion of Normandy?
Confronted about an Amazon review
Will I be allowed to enter the US after living there illegally then legally in the past?
To project, or not to project? Extracting raster values with R
Projecting shp from GK3 to UTM32 with BeTA2007How do I make this into a working CRS in ArcGIS?QGIS doesn't create the correct projections when creating new shp filesExtract values from raster using point data with buffer in RProcessing vector to raster faster with RCRS Projection problems prevent spatial merge in RCompute a distance matrix in R with geosphereTransforming map with “custom” CRS to WGS84?
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty{
margin-bottom:0;
}
I have a raster from which I need to extract values at point locations using raster::extract(raster_stack, points_sf).
The raster is at 30 arc-second resolution and came from WorldClim with proj4string: "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"
i.e. it's coordinates are geographic on an ellipsoid. My point location data has the same lon-lat WGS84 CRS.
I've noticed that about 3/4 of the values I extract differ depending on whether or not I've transformed both the raster and the points to a planar projection with proj4string:
"+proj=laea +lat_0=43.15268 +lon_0=-70.30744 +x_0=12019341.4
+y_0=1730136 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs".
My question is, should I expect points to move between raster cells this much as a result of projection? And which method is more accurate? In case spatial scale matters, the bounding box for the points is:
xmin ymin xmax ymax
-73.68631 41.01861 -67.89930 45.23919
Please let me know in the comments if this question needs a reproducible example.
raster coordinate-system r proj
asked 8 hours ago
GregoryGregory
7932 gold badges11 silver badges27 bronze badges
add a comment
|
I have a raster from which I need to extract values at point locations using raster::extract(raster_stack, points_sf).
The raster is at 30 arc-second resolution and came from WorldClim with proj4string: "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"
i.e. it's coordinates are geographic on an ellipsoid. My point location data has the same lon-lat WGS84 CRS.
I've noticed that about 3/4 of the values I extract differ depending on whether or not I've transformed both the raster and the points to a planar projection with proj4string:
"+proj=laea +lat_0=43.15268 +lon_0=-70.30744 +x_0=12019341.4
+y_0=1730136 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs".
My question is, should I expect points to move between raster cells this much as a result of projection? And which method is more accurate? In case spatial scale matters, the bounding box for the points is:
xmin ymin xmax ymax
-73.68631 41.01861 -67.89930 45.23919
Please let me know in the comments if this question needs a reproducible example.
raster coordinate-system r proj
asked 8 hours ago
GregoryGregory
7932 gold badges11 silver badges27 bronze badges
add a comment
|
I have a raster from which I need to extract values at point locations using raster::extract(raster_stack, points_sf).
The raster is at 30 arc-second resolution and came from WorldClim with proj4string: "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"
i.e. it's coordinates are geographic on an ellipsoid. My point location data has the same lon-lat WGS84 CRS.
I've noticed that about 3/4 of the values I extract differ depending on whether or not I've transformed both the raster and the points to a planar projection with proj4string:
"+proj=laea +lat_0=43.15268 +lon_0=-70.30744 +x_0=12019341.4
+y_0=1730136 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs".
My question is, should I expect points to move between raster cells this much as a result of projection? And which method is more accurate? In case spatial scale matters, the bounding box for the points is:
xmin ymin xmax ymax
-73.68631 41.01861 -67.89930 45.23919
Please let me know in the comments if this question needs a reproducible example.
raster coordinate-system r proj
asked 8 hours ago
GregoryGregory
7932 gold badges11 silver badges27 bronze badges
I have a raster from which I need to extract values at point locations using raster::extract(raster_stack, points_sf).
The raster is at 30 arc-second resolution and came from WorldClim with proj4string: "+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0"
i.e. it's coordinates are geographic on an ellipsoid. My point location data has the same lon-lat WGS84 CRS.
I've noticed that about 3/4 of the values I extract differ depending on whether or not I've transformed both the raster and the points to a planar projection with proj4string:
"+proj=laea +lat_0=43.15268 +lon_0=-70.30744 +x_0=12019341.4
+y_0=1730136 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs".
My question is, should I expect points to move between raster cells this much as a result of projection? And which method is more accurate? In case spatial scale matters, the bounding box for the points is:
xmin ymin xmax ymax
-73.68631 41.01861 -67.89930 45.23919
Please let me know in the comments if this question needs a reproducible example.
raster coordinate-system r proj
raster coordinate-system r proj
asked 8 hours ago
GregoryGregory
7932 gold badges11 silver badges27 bronze badges
asked 8 hours ago
GregoryGregory
7932 gold badges11 silver badges27 bronze badges
asked 8 hours ago
GregoryGregory
7932 gold badges11 silver badges27 bronze badges
asked 8 hours ago
GregoryGregory
7932 gold badges11 silver badges27 bronze badges
asked 8 hours ago
GregoryGregory
7932 gold badges11 silver badges27 bronze badges
7932 gold badges11 silver badges27 bronze badges
add a comment
|
add a comment
|
1 Answer
1
active
oldest
votes
Reprojecting rasters is usually a bad thing to do. It involves a non-reversible transformation from one grid system to another grid system that can have a non-linear relationship to the first. Hence the value in a cell of the new system can end up being some average of whichever grid cells in the source raster it overlapped.
If you have a raster and points in the same coordinate system, work with that.
If you have a raster and points in a different coordinate system, reproject the points to the raster system (since that is a precise and reversible transformation with no information loss).
Only if you have two rasters with different coordinate systems and you want to do some kind of operation between them should you think about reprojecting one of them.
answered 6 hours ago
SpacedmanSpacedman
29.4k3 gold badges38 silver badges56 bronze badges
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/4.0/"u003ecc by-sa 4.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%2f338522%2fto-project-or-not-to-project-extracting-raster-values-with-r%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
Reprojecting rasters is usually a bad thing to do. It involves a non-reversible transformation from one grid system to another grid system that can have a non-linear relationship to the first. Hence the value in a cell of the new system can end up being some average of whichever grid cells in the source raster it overlapped.
If you have a raster and points in the same coordinate system, work with that.
If you have a raster and points in a different coordinate system, reproject the points to the raster system (since that is a precise and reversible transformation with no information loss).
Only if you have two rasters with different coordinate systems and you want to do some kind of operation between them should you think about reprojecting one of them.
answered 6 hours ago
SpacedmanSpacedman
29.4k3 gold badges38 silver badges56 bronze badges
add a comment
|
Reprojecting rasters is usually a bad thing to do. It involves a non-reversible transformation from one grid system to another grid system that can have a non-linear relationship to the first. Hence the value in a cell of the new system can end up being some average of whichever grid cells in the source raster it overlapped.
If you have a raster and points in the same coordinate system, work with that.
If you have a raster and points in a different coordinate system, reproject the points to the raster system (since that is a precise and reversible transformation with no information loss).
Only if you have two rasters with different coordinate systems and you want to do some kind of operation between them should you think about reprojecting one of them.
answered 6 hours ago
SpacedmanSpacedman
29.4k3 gold badges38 silver badges56 bronze badges
add a comment
|
Reprojecting rasters is usually a bad thing to do. It involves a non-reversible transformation from one grid system to another grid system that can have a non-linear relationship to the first. Hence the value in a cell of the new system can end up being some average of whichever grid cells in the source raster it overlapped.
If you have a raster and points in the same coordinate system, work with that.
If you have a raster and points in a different coordinate system, reproject the points to the raster system (since that is a precise and reversible transformation with no information loss).
Only if you have two rasters with different coordinate systems and you want to do some kind of operation between them should you think about reprojecting one of them.
answered 6 hours ago
SpacedmanSpacedman
29.4k3 gold badges38 silver badges56 bronze badges
Reprojecting rasters is usually a bad thing to do. It involves a non-reversible transformation from one grid system to another grid system that can have a non-linear relationship to the first. Hence the value in a cell of the new system can end up being some average of whichever grid cells in the source raster it overlapped.
If you have a raster and points in the same coordinate system, work with that.
If you have a raster and points in a different coordinate system, reproject the points to the raster system (since that is a precise and reversible transformation with no information loss).
Only if you have two rasters with different coordinate systems and you want to do some kind of operation between them should you think about reprojecting one of them.
answered 6 hours ago
SpacedmanSpacedman
29.4k3 gold badges38 silver badges56 bronze badges
answered 6 hours ago
SpacedmanSpacedman
29.4k3 gold badges38 silver badges56 bronze badges
answered 6 hours ago
SpacedmanSpacedman
29.4k3 gold badges38 silver badges56 bronze badges
answered 6 hours ago
SpacedmanSpacedman
29.4k3 gold badges38 silver badges56 bronze badges
29.4k3 gold badges38 silver badges56 bronze badges
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%2f338522%2fto-project-or-not-to-project-extracting-raster-values-with-r%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
