On Thu, 9 Jun 2022, Ranjeet Kumar Jha wrote:
Hello Roger, Thanks for your kind response! I have provided the reproducible data and code now in this shared link. https://drive.google.com/drive/folders/1dNmGTI8_c9PK1QqmfIjnpbyzuiCXgxFC?msclkid=4816425dd04111ec912cf8b4175aa93c Please now look at it! I appreciate your help.
No, now it is up to you to take the reproducible example, and modify the alternative you choose to suit your needs. You can keep your existing code but convert each sf plot to a grob in turn, use tmap, or use mapsf, but if I do that, you will not learn from the process. That is why reproducible examples are crucial; they provide a shared platform for learning. Providing access to your data may be sensible in an off-list setting when no reproducible example permits a bug to be found, but preferably not otherwise. Roger
Regards and Thanks, Ranjeet On Thu, Jun 9, 2022 at 2:48 PM Roger Bivand <Roger.Bivand at nhh.no> wrote:
On Thu, 9 Jun 2022, J. Marius RAKOTONDRAMANGA wrote:
Hi Ranjeet, Have a look to this function cowplot::plot_grid(), from
Marius,
Thanks for trying to help, but cowplot is designed to combine plots
created using grid, not base, graphics. To use it (or gridExtra), sf
plots, which are base graphics, must be converted to grid objects (grobs),
for example using the gridGraphics package.
Since Ranjeet did not provide a reproducible example, this may not
correspond well, but given the lack of information, it's hard to be more
precise:
library(sf)
nc <- st_read(system.file("gpkg/nc.gpkg", package="sf"))
plot(nc["BIR74"])
creates output on the default device, but no object that can be
manipulated. Using:
par(mfrow=c(2,2))
plot(nc["BIR74"])
gives exactly the same output, because sf's plot() itself takes control of
the graphics device, so ignoring the request for multiple panels.
par(mfrow=c(1,1))
plot(nc["BIR74"])
gridGraphics::grid.echo()
library(grid)
g <- grid.grab()
grid.newpage()
gridExtra::grid.arrange(g, b, ncol=2)
gives two copies of our initial output after conversion to grid graphics.
However, it is easier to choose alternatives:
library(tmap)
tm_shape(nc) + tm_fill(c("BIR74", "BIR79"))
and probably library(mapsf) too, but it, like sf, uses base graphics; it
respects par(mfrow=):
library(mapsf)
par(mfrow=c(1,2))
mf_map(x=nc, var=c("BIR74"), type="choro")
mf_map(x=nc, var=c("BIR79"), type="choro")
Hope this clarifies,
Roger
https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwilkelab.org%2Fcowplot%2Farticles%2Fplot_grid.html&data=05%7C01%7CRoger.Bivand%40nhh.no%7C01ce113052f9492ab51d08da49ec9930%7C33a15b2f849941998d56f20b5aa91af2%7C0%7C0%7C637903577780613163%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=vrH2iUZeaQw6TgzuYAYLXS6Yv%2Bt0hVXgu6PPDYcriEg%3D&reserved=0
It has mainy features and permit you to save also the combined plots : here
https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fsearch.r-project.org%2FCRAN%2Frefmans%2Fcowplot%2Fhtml%2Fsave_plot.html&data=05%7C01%7CRoger.Bivand%40nhh.no%7C01ce113052f9492ab51d08da49ec9930%7C33a15b2f849941998d56f20b5aa91af2%7C0%7C0%7C637903577780613163%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=gRj%2BPj%2B29W0mrOtT0Vp5HNEigvBQUsOL4VAdF%2B64fe0%3D&reserved=0
or
https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fwilkelab.org%2Fcowplot%2Freference%2Fggsave2.html&data=05%7C01%7CRoger.Bivand%40nhh.no%7C01ce113052f9492ab51d08da49ec9930%7C33a15b2f849941998d56f20b5aa91af2%7C0%7C0%7C637903577780613163%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=Idsy6TN5bxrqZtqgi6CGpqzRKjcofl71cR9EBTc9jvY%3D&reserved=0
Regards, Marius, --------- J. Marius RAKOTONDRAMANGA, MSc PhD Student in Biostatistics/Biomathematics ED 393 Pierre Louis de Sant? Publique : Epid?miologie et Sciences de l'Information Biom?dicale Sorbonne Universit? (ex. UPMC) Le jeu. 9 juin 2022 ? 08:29, Ranjeet Kumar Jha <
ranjeetjhaiitkgp at gmail.com>
a ?crit :
Hi, I like to plot all my 5 spatial maps of drought indices in a single
window
using the sf plot function in R. I am able to plot all maps separately
but
unable to plot altogether. All indexes are showing drought indices.
I really appreciate your help!
Here is the code:
#merging using merge function
merge_drought_indices <- merge(shp_dist, drought_indices)
view(merge_drought_indices)
#plot in a single window
library(sf)
par(mfrow=c(2,2))
plot(merge_drought_indices["r_count_lt_p05_harvest_end"],
pal = colorRampPalette(c("red", "white", "blue","green")),
main = "End harvest index variation in India during drought year 2000",
key.pos = 4, axes = TRUE, key.width = lcm(1.3), key.length = 1.0)
plot(merge_drought_indices["r_count_lt_p05_harvest"],
pal = colorRampPalette(c("red", "white", "blue","green")),
main = "Harvest variation index in India during drought year 2000",
key.pos = 4, axes = TRUE, key.width = lcm(1.3), key.length = 1.0)
plot(merge_drought_indices["r_count_lt_p05_main"],
pal = colorRampPalette(c("red", "white", "blue","green")),
main = "Main growth stage index variation in India during drought year
2000",
key.pos = 4, axes = TRUE, key.width = lcm(1.3), key.length = 1.0)
plot(merge_drought_indices["r_count_lt_p05_plant.start"],
pal = colorRampPalette(c("red", "white", "blue","green")),
main = "Plant start index variation in India during drought year 2000",
key.pos = 4, axes = TRUE, key.width = lcm(1.3), key.length = 1.0)
plot(merge_drought_indices["r_count_lt_p05_plant"],
pal = colorRampPalette(c("red", "white", "blue","green")),
main = "main plant stage index variation in India during drought year
2000",
key.pos = 4, axes = TRUE, key.width = lcm(1.3), key.length = 1.0)
Regards and Thanks,
Ranjeet
--
*"Simple Heart, Humble Attitude and Surrender to Supreme Being make our
lives beautiful!"*
[[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org
https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-geo&data=05%7C01%7CRoger.Bivand%40nhh.no%7C01ce113052f9492ab51d08da49ec9930%7C33a15b2f849941998d56f20b5aa91af2%7C0%7C0%7C637903577780613163%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=awwczk%2Bpse%2B5tTKYB1ecwtuNE3C6AjgtgvzbE7Aga6w%3D&reserved=0
[[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo at r-project.org
https://eur02.safelinks.protection.outlook.com/?url=https%3A%2F%2Fstat.ethz.ch%2Fmailman%2Flistinfo%2Fr-sig-geo&data=05%7C01%7CRoger.Bivand%40nhh.no%7C01ce113052f9492ab51d08da49ec9930%7C33a15b2f849941998d56f20b5aa91af2%7C0%7C0%7C637903577780613163%7CUnknown%7CTWFpbGZsb3d8eyJWIjoiMC4wLjAwMDAiLCJQIjoiV2luMzIiLCJBTiI6Ik1haWwiLCJXVCI6Mn0%3D%7C3000%7C%7C%7C&sdata=awwczk%2Bpse%2B5tTKYB1ecwtuNE3C6AjgtgvzbE7Aga6w%3D&reserved=0
-- Roger Bivand Emeritus Professor Department of Economics, Norwegian School of Economics, Postboks 3490 Ytre Sandviken, 5045 Bergen, Norway. e-mail: Roger.Bivand at nhh.no https://orcid.org/0000-0003-2392-6140 https://scholar.google.no/citations?user=AWeghB0AAAAJ&hl=en
Roger Bivand Emeritus Professor Department of Economics, Norwegian School of Economics, Postboks 3490 Ytre Sandviken, 5045 Bergen, Norway. e-mail: Roger.Bivand at nhh.no https://orcid.org/0000-0003-2392-6140 https://scholar.google.no/citations?user=AWeghB0AAAAJ&hl=en