[R-sig-Geo] Need to plot all 5 spatial maps together in a single window

Thu Jun 9 11:52:05 CEST 2022

Thank you Prof. Bivand for correcting my views.
Your insights and books are really helpful for us !

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 à 12:18, Roger Bivand <Roger.Bivand using nhh.no> a écrit :

> 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 using 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 using 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 using 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 using nhh.no
> https://orcid.org/0000-0003-2392-6140
> https://scholar.google.no/citations?user=AWeghB0AAAAJ&hl=en

	[[alternative HTML version deleted]]