Overlay shapefile on multiband raster using ggRGB and ggplot2 in R The 2019 Stack Overflow...
What are the motivations for publishing new editions of an existing textbook, beyond new discoveries in a field?
It's possible to achieve negative score?
"Riffle" two strings
What does "rabbited" mean/imply in this sentence?
Manuscript was "unsubmitted" because the manuscript was deposited in Arxiv Preprints
What are the advantages and disadvantages of running one shots compared to campaigns?
Why Did Howard Stark Use All The Vibranium They Had On A Prototype Shield?
What is this 4-propeller plane?
I see my dog run
Should I use my personal or workplace e-mail when registering to external websites for work purpose?
If Wish Duplicates Simulacrum, Are Existing Duplicates Destroyed?
How to create dashed lines/arrows in Illustrator CC
Is it possible to build an equivalent function just looking at the input and output of the original function?
What is the use of option -o in the useradd command?
What is the steepest angle that a canal can be traversable without locks?
Which Sci-Fi work first showed weapon of galactic-scale mass destruction?
What is the motivation for a law requiring 2 parties to consent for recording a conversation
Is "plugging out" electronic devices an American expression?
Can I write a for loop that iterates over both collections and arrays?
What could be the right powersource for 15 seconds lifespan disposable giant chainsaw?
Can't find the latex code for the ⍎ (down tack jot) symbol
What do the Banks children have against barley water?
Dual Citizen. Exited the US on Italian passport recently
How to deal with fear of taking dependencies
Overlay shapefile on multiband raster using ggRGB and ggplot2 in R
The 2019 Stack Overflow Developer Survey Results Are InChoropleth in R: custom breaks and plottingUsing stat_density2d with coord_map in ggplot2Fortifying re-projected shapefile for ggplot2 in RPlotting a certain range in azimuthal projection using ggplot2How can I create a choropleth from the 2010 census ZCTA shapefile using R and ggplot2?Rasterizing getCover not working in R?Map with ggplot2 and bubblesCombining postal codes into sales territories by value using R and ggplot2?Reducing processing time of multiband satellite images using RSToolbox in RGgplot2 map with colored background and coord_map()
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty{ margin-bottom:0;
}
I am trying to use ggplot2 to map a polygon feature class over a raster brick (created from multiple Landsat band images). I'm able to map the polygons and the raster separately, but when I try to have them together (following the last code chunk at http://bleutner.github.io/RStoolbox/rstbx-docu/ggRGB.html RStudio crashes (even after restarting the computer).
ggRGB attempt:
fallowFields <- rgdal::readOGR(dsn = fgdb, layer = 'Fallowing2018')
fallowFieldsWGS84 <- spTransform(fallowFields, crs(L8_20180420_brick))
fallowFieldsDataFrame <- fallowFieldsWGS %>% fortify
landsatAprilCropped <- raster::crop(L8_20180420_brick, raster::extent(fallowFields))
p <- ggplot()
p + ggRGB(img = landsatAprilCropped,
r = 3,
g = 2,
b = 1,
stretch = 'hist',
ggLayer = TRUE) +
geom_polygon(col = 'Black',
fill = 'grey70',
data = fallowFieldsDataFrame,
aes(x = long, y = lat, group = group))
replacing the data frame (fallowFieldsDataFrame) with the shapefile (fallowFieldsWGS84) works, but it warps the Landsat image.
The cropped image and shapefile are in the same projection, but the shapefile seems to be warping the raster. Is there a way to plot the shapefile over the raster without warping it?
It doesn't get warped with plotRGB, but I was hoping to use ggRGB so that I could make a more elegant map with gglot2.
plotRGB(landsatAprilCropped,
r = 3, g = 2, b = 1,
stretch = 'hist', axes = TRUE,
main = 'Fallow Fields Overlaid on Landsat Image: 4/20/2018')
plot(bardFallowFields, col = 'transparent', border = 'red', add = TRUE)
raster coordinate-system r ggplot2
asked Aug 27 '18 at 21:50
IracambiIracambi
13110
bumped to the homepage by Community♦ 19 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 |
I am trying to use ggplot2 to map a polygon feature class over a raster brick (created from multiple Landsat band images). I'm able to map the polygons and the raster separately, but when I try to have them together (following the last code chunk at http://bleutner.github.io/RStoolbox/rstbx-docu/ggRGB.html RStudio crashes (even after restarting the computer).
ggRGB attempt:
fallowFields <- rgdal::readOGR(dsn = fgdb, layer = 'Fallowing2018')
fallowFieldsWGS84 <- spTransform(fallowFields, crs(L8_20180420_brick))
fallowFieldsDataFrame <- fallowFieldsWGS %>% fortify
landsatAprilCropped <- raster::crop(L8_20180420_brick, raster::extent(fallowFields))
p <- ggplot()
p + ggRGB(img = landsatAprilCropped,
r = 3,
g = 2,
b = 1,
stretch = 'hist',
ggLayer = TRUE) +
geom_polygon(col = 'Black',
fill = 'grey70',
data = fallowFieldsDataFrame,
aes(x = long, y = lat, group = group))
replacing the data frame (fallowFieldsDataFrame) with the shapefile (fallowFieldsWGS84) works, but it warps the Landsat image.
The cropped image and shapefile are in the same projection, but the shapefile seems to be warping the raster. Is there a way to plot the shapefile over the raster without warping it?
It doesn't get warped with plotRGB, but I was hoping to use ggRGB so that I could make a more elegant map with gglot2.
plotRGB(landsatAprilCropped,
r = 3, g = 2, b = 1,
stretch = 'hist', axes = TRUE,
main = 'Fallow Fields Overlaid on Landsat Image: 4/20/2018')
plot(bardFallowFields, col = 'transparent', border = 'red', add = TRUE)
raster coordinate-system r ggplot2
asked Aug 27 '18 at 21:50
IracambiIracambi
13110
bumped to the homepage by Community♦ 19 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 |
I am trying to use ggplot2 to map a polygon feature class over a raster brick (created from multiple Landsat band images). I'm able to map the polygons and the raster separately, but when I try to have them together (following the last code chunk at http://bleutner.github.io/RStoolbox/rstbx-docu/ggRGB.html RStudio crashes (even after restarting the computer).
ggRGB attempt:
fallowFields <- rgdal::readOGR(dsn = fgdb, layer = 'Fallowing2018')
fallowFieldsWGS84 <- spTransform(fallowFields, crs(L8_20180420_brick))
fallowFieldsDataFrame <- fallowFieldsWGS %>% fortify
landsatAprilCropped <- raster::crop(L8_20180420_brick, raster::extent(fallowFields))
p <- ggplot()
p + ggRGB(img = landsatAprilCropped,
r = 3,
g = 2,
b = 1,
stretch = 'hist',
ggLayer = TRUE) +
geom_polygon(col = 'Black',
fill = 'grey70',
data = fallowFieldsDataFrame,
aes(x = long, y = lat, group = group))
replacing the data frame (fallowFieldsDataFrame) with the shapefile (fallowFieldsWGS84) works, but it warps the Landsat image.
The cropped image and shapefile are in the same projection, but the shapefile seems to be warping the raster. Is there a way to plot the shapefile over the raster without warping it?
It doesn't get warped with plotRGB, but I was hoping to use ggRGB so that I could make a more elegant map with gglot2.
plotRGB(landsatAprilCropped,
r = 3, g = 2, b = 1,
stretch = 'hist', axes = TRUE,
main = 'Fallow Fields Overlaid on Landsat Image: 4/20/2018')
plot(bardFallowFields, col = 'transparent', border = 'red', add = TRUE)
raster coordinate-system r ggplot2
asked Aug 27 '18 at 21:50
IracambiIracambi
13110
I am trying to use ggplot2 to map a polygon feature class over a raster brick (created from multiple Landsat band images). I'm able to map the polygons and the raster separately, but when I try to have them together (following the last code chunk at http://bleutner.github.io/RStoolbox/rstbx-docu/ggRGB.html RStudio crashes (even after restarting the computer).
ggRGB attempt:
fallowFields <- rgdal::readOGR(dsn = fgdb, layer = 'Fallowing2018')
fallowFieldsWGS84 <- spTransform(fallowFields, crs(L8_20180420_brick))
fallowFieldsDataFrame <- fallowFieldsWGS %>% fortify
landsatAprilCropped <- raster::crop(L8_20180420_brick, raster::extent(fallowFields))
p <- ggplot()
p + ggRGB(img = landsatAprilCropped,
r = 3,
g = 2,
b = 1,
stretch = 'hist',
ggLayer = TRUE) +
geom_polygon(col = 'Black',
fill = 'grey70',
data = fallowFieldsDataFrame,
aes(x = long, y = lat, group = group))
replacing the data frame (fallowFieldsDataFrame) with the shapefile (fallowFieldsWGS84) works, but it warps the Landsat image.
The cropped image and shapefile are in the same projection, but the shapefile seems to be warping the raster. Is there a way to plot the shapefile over the raster without warping it?
It doesn't get warped with plotRGB, but I was hoping to use ggRGB so that I could make a more elegant map with gglot2.
plotRGB(landsatAprilCropped,
r = 3, g = 2, b = 1,
stretch = 'hist', axes = TRUE,
main = 'Fallow Fields Overlaid on Landsat Image: 4/20/2018')
plot(bardFallowFields, col = 'transparent', border = 'red', add = TRUE)
raster coordinate-system r ggplot2
raster coordinate-system r ggplot2
asked Aug 27 '18 at 21:50
IracambiIracambi
13110
asked Aug 27 '18 at 21:50
IracambiIracambi
13110
edited Aug 27 '18 at 22:02
Iracambi
asked Aug 27 '18 at 21:50
IracambiIracambi
13110
asked Aug 27 '18 at 21:50
IracambiIracambi
13110
asked Aug 27 '18 at 21:50
IracambiIracambi
13110
13110
bumped to the homepage by Community♦ 19 mins ago
This question has answers that may be good or bad; the system has marked it active so that they can be reviewed.
bumped to the homepage by Community♦ 19 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 |
add a comment |
1 Answer
1
active
oldest
votes
You miss a very important part. First, an example:
library(raster)
library(ggplot2)
library(RStoolbox)
data(lsat)
lsat
poly <- as(extent(622999,625000,-417000,-415000),'SpatialPolygons')
poly <- fortify(poly)
names(poly)[1:2] <- c('x','y')
e <- extent(lsat)
p <- ggplot() + ggRGB(img = lsat,
r = 3,
g = 2,
b = 1,
stretch = 'hist',
ggLayer = T) +
xlim(e[1:2]) + ylim(e[3:4]) +
geom_polygon(aes(x,y), alpha = 0.9, data = poly)
Using ggRBG() output is:
ggRGB(img = lsat,
r = 3,
g = 2,
b = 1,
stretch = 'hist')
With your procedure I reproduced the same issue:
p
The solution is to add coord_equal()
p + coord_equal()
I recomend you to add ylim() and xlim() with the extent of raster object
answered Aug 30 '18 at 12:00
aldo_tapiaaldo_tapia
7,93451742
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%2f294151%2foverlay-shapefile-on-multiband-raster-using-ggrgb-and-ggplot2-in-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
You miss a very important part. First, an example:
library(raster)
library(ggplot2)
library(RStoolbox)
data(lsat)
lsat
poly <- as(extent(622999,625000,-417000,-415000),'SpatialPolygons')
poly <- fortify(poly)
names(poly)[1:2] <- c('x','y')
e <- extent(lsat)
p <- ggplot() + ggRGB(img = lsat,
r = 3,
g = 2,
b = 1,
stretch = 'hist',
ggLayer = T) +
xlim(e[1:2]) + ylim(e[3:4]) +
geom_polygon(aes(x,y), alpha = 0.9, data = poly)
Using ggRBG() output is:
ggRGB(img = lsat,
r = 3,
g = 2,
b = 1,
stretch = 'hist')
With your procedure I reproduced the same issue:
p
The solution is to add coord_equal()
p + coord_equal()
I recomend you to add ylim() and xlim() with the extent of raster object
answered Aug 30 '18 at 12:00
aldo_tapiaaldo_tapia
7,93451742
add a comment |
You miss a very important part. First, an example:
library(raster)
library(ggplot2)
library(RStoolbox)
data(lsat)
lsat
poly <- as(extent(622999,625000,-417000,-415000),'SpatialPolygons')
poly <- fortify(poly)
names(poly)[1:2] <- c('x','y')
e <- extent(lsat)
p <- ggplot() + ggRGB(img = lsat,
r = 3,
g = 2,
b = 1,
stretch = 'hist',
ggLayer = T) +
xlim(e[1:2]) + ylim(e[3:4]) +
geom_polygon(aes(x,y), alpha = 0.9, data = poly)
Using ggRBG() output is:
ggRGB(img = lsat,
r = 3,
g = 2,
b = 1,
stretch = 'hist')
With your procedure I reproduced the same issue:
p
The solution is to add coord_equal()
p + coord_equal()
I recomend you to add ylim() and xlim() with the extent of raster object
answered Aug 30 '18 at 12:00
aldo_tapiaaldo_tapia
7,93451742
add a comment |
You miss a very important part. First, an example:
library(raster)
library(ggplot2)
library(RStoolbox)
data(lsat)
lsat
poly <- as(extent(622999,625000,-417000,-415000),'SpatialPolygons')
poly <- fortify(poly)
names(poly)[1:2] <- c('x','y')
e <- extent(lsat)
p <- ggplot() + ggRGB(img = lsat,
r = 3,
g = 2,
b = 1,
stretch = 'hist',
ggLayer = T) +
xlim(e[1:2]) + ylim(e[3:4]) +
geom_polygon(aes(x,y), alpha = 0.9, data = poly)
Using ggRBG() output is:
ggRGB(img = lsat,
r = 3,
g = 2,
b = 1,
stretch = 'hist')
With your procedure I reproduced the same issue:
p
The solution is to add coord_equal()
p + coord_equal()
I recomend you to add ylim() and xlim() with the extent of raster object
answered Aug 30 '18 at 12:00
aldo_tapiaaldo_tapia
7,93451742
You miss a very important part. First, an example:
library(raster)
library(ggplot2)
library(RStoolbox)
data(lsat)
lsat
poly <- as(extent(622999,625000,-417000,-415000),'SpatialPolygons')
poly <- fortify(poly)
names(poly)[1:2] <- c('x','y')
e <- extent(lsat)
p <- ggplot() + ggRGB(img = lsat,
r = 3,
g = 2,
b = 1,
stretch = 'hist',
ggLayer = T) +
xlim(e[1:2]) + ylim(e[3:4]) +
geom_polygon(aes(x,y), alpha = 0.9, data = poly)
Using ggRBG() output is:
ggRGB(img = lsat,
r = 3,
g = 2,
b = 1,
stretch = 'hist')
With your procedure I reproduced the same issue:
p
The solution is to add coord_equal()
p + coord_equal()
I recomend you to add ylim() and xlim() with the extent of raster object
answered Aug 30 '18 at 12:00
aldo_tapiaaldo_tapia
7,93451742
answered Aug 30 '18 at 12:00
aldo_tapiaaldo_tapia
7,93451742
answered Aug 30 '18 at 12:00
aldo_tapiaaldo_tapia
7,93451742
answered Aug 30 '18 at 12:00
aldo_tapiaaldo_tapia
7,93451742
7,93451742
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%2f294151%2foverlay-shapefile-on-multiband-raster-using-ggrgb-and-ggplot2-in-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
