[R-sig-Geo] Need to merge .RDS dataframe and India Shape file, and then plot yield map with tile legend

Thu May 12 17:19:45 CEST 2022

Hi Ranjeet,

I didn’t find out how to do it with sf, but you can get it using the mapsf
package:

library(mapsf)

brks <- mf_get_breaks(merge_yield$mod.yield.all, breaks = "pretty")
nclass <- length(brks) - 1

# set a theme
mf_theme("ink")

# plot yield
mf_map(
    x = merge_yield,
    var = "mod.yield.all",
    type = "choro",
    breaks = brks,
    pal = colorRampPalette(c("red", "white", "blue"))(nclass),
    border = "grey",
    lwd = 0.5,
    leg_val_rnd = 0,
    leg_pos = "topright",
    leg_title = "Yield (tons/ha)"
)

# title etc.
mf_layout(title = "Rice yield by district in India",
          credits = "", scale = FALSE, arrow = FALSE)

Kind regards,
Zivan

On Thu, May 12, 2022 at 3:42 PM Ranjeet Kumar Jha <
ranjeetjhaiitkgp using gmail.com> wrote:

> Hi Zivan, Just a quick question - Can we add legend title - "Yield
> (tons/ha) with this plot function in sf package. I could n't find so.
>
> Thanks,
> Ranjeet
>
> On Thu, May 12, 2022 at 4:34 PM Ranjeet Kumar Jha <
> ranjeetjhaiitkgp using gmail.com> wrote:
>
>> Thank you Zivan for your kind help and support ! Code really worked now
>> and I understood the concept now!
>>
>> I again appreciate your help!
>>
>> Regards and Thanks,
>> Ranjeet
>>
>> On Thu, May 12, 2022 at 4:06 PM Zivan Karaman <zivan.karaman using gmail.com>
>> wrote:
>>
>>> Hi Ranjeet,
>>>
>>> Perhaps you could try the following script:
>>>
>>> #####Yield data
>>> #average rice yield
>>> dist_yield <- readRDS(file = "dist_yield.RDS")
>>>
>>> #India shape file for districts
>>> library(sf)
>>> shp_dist <- read_sf("District_sim.shp")
>>>
>>> #merging using merge function
>>> merge_yield <- merge(shp_dist, dist_yield)
>>>
>>> #plot
>>> nbrks <- 10
>>> plot(merge_yield["mod.yield.all"], nbreaks = nbrks, pal =
>>> colorRampPalette(c("red", "white", "blue"))(nbrks), main = "yield by
>>> district")
>>>
>>> The plotting is done by the plot from the sf package, you can consult
>>> the help page for more options.
>>> You can, of course, change the color palette, the number of breaks,
>>> title, and other decorations as you wish.
>>>
>>> Kind regards,
>>> Zivan
>>>
>>>
>>> On Thu, May 12, 2022 at 9:03 AM Ranjeet Kumar Jha <
>>> ranjeetjhaiitkgp using gmail.com> wrote:
>>>
>>>> Hello Roger,
>>>>
>>>> I really appreciate your support and help! Now I have uploaded all the
>>>> required .shx, .dbf, and .prj files in the link folder.
>>>>
>>>>
>>>> https://drive.google.com/drive/folders/1dNmGTI8_c9PK1QqmfIjnpbyzuiCXgxFC?msclkid=4816425dd04111ec912cf8b4175aa93c
>>>>
>>>> On Wed, May 11, 2022 at 5:42 PM Roger Bivand <Roger.Bivand using nhh.no>
>>>> wrote:
>>>>
>>>> > On Tue, 10 May 2022, Ranjeet Kumar Jha wrote:
>>>> >
>>>> > > Hello Roger,
>>>> > >
>>>> > > Thanks for your kind reply! I am going through this. However,
>>>> sharing the
>>>> > > link of code, shape file and .rds data frame for better
>>>> understanding and
>>>> > > helping through the real files.
>>>> >
>>>> > District_sim.shp must of course be accompanied by District_sim.shx,
>>>> > District_sim.dbf and most likely District_sim.prj, which are absent.
>>>> >
>>>> > Roger
>>>> >
>>>> > >
>>>> > >
>>>> >
>>>> https://drive.google.com/drive/folders/1dNmGTI8_c9PK1QqmfIjnpbyzuiCXgxFC?msclkid=4816425dd04111ec912cf8b4175aa93c
>>>> > >
>>>> > > Regards and Thanks,
>>>> > > Ranjeet
>>>> > >
>>>> > > On Tue, May 10, 2022 at 1:34 PM Roger Bivand <Roger.Bivand using nhh.no>
>>>> > wrote:
>>>> > >
>>>> > >> On Tue, 10 May 2022, Ranjeet Kumar Jha wrote:
>>>> > >>
>>>> > >>> Hello Everyone,
>>>> > >>>
>>>> > >>> I am trying to plot an Indian district shapefile filled with
>>>> modeled
>>>> > >> yield,
>>>> > >>> corresponding to each GDMID. I needed to first merge the Yield
>>>> data
>>>> > with
>>>> > >>> respect to each GDMID from the .RDS file with the district
>>>> shapefile of
>>>> > >>> India. Each district is represented by a GDMID. After merging, the
>>>> > Indian
>>>> > >>> district map filled with modeled yield values needs to be plotted
>>>> > similar
>>>> > >>> to an example attached figure. Shape file has five features, and
>>>> one of
>>>> > >>> them is "GDMID".
>>>> > >>
>>>> > >> The code is not a reproducible example; it needs to run with no
>>>> > reference
>>>> > >> to local directories. The use of tibbles or other misunderstandings
>>>> > >> represents state.ID and GDMID as numeric, they should be
>>>> character, you
>>>> > >> cannot do arithmetic on them. The level 2 boundaries from
>>>> > >> https://geodata.ucdavis.edu/gadm/gadm4.0/gpkg/gadm40_IND.gpkg do
>>>> not
>>>> > use
>>>> > >> GDMID, so unless you provide a download link to your boundaries,
>>>> no-one
>>>> > >> can help. Also, your data.frame has 694 observations, GADM has 666.
>>>> > >>
>>>> > >> My guess is that you have found a blog somewhere, probably very
>>>> > >> out-of-date (you use rgdal::readOGR() not sf::st_read()), and
>>>> which is
>>>> > >> causing lots of problems, because it complicates things.
>>>> > >>
>>>> > >> Use the sf package, avoid tidyverse which can convert
>>>> representations
>>>> > >> without asking, make sure that the IDs are character not numeric
>>>> (to
>>>> > >> preserve any leading zeros), and use tmap or mapsf not ggplot (much
>>>> > >> simpler and more intuitive interface for mapping).
>>>> > >>
>>>> > >> See https://geocompr.robinlovelace.net/ and the sf vignettes.
>>>> Once you
>>>> > >> have control of the IDs (for example try match() to see if they do
>>>> match
>>>> > >> in type and value). sf has a merge method, which gives a good deal
>>>> of
>>>> > >> control.
>>>> > >>
>>>> > >> Hope this helps,
>>>> > >>
>>>> > >> Roger
>>>> > >>
>>>> > >>>
>>>> > >>> I feel there is something wrong with my code for merging and
>>>> plotting
>>>> > by
>>>> > >>> ggplot.
>>>> > >>>
>>>> > >>> I have attached herewith my code, and RDS data-frame.
>>>> > >>>
>>>> > >>> It would be great if you could help me with merging issue and
>>>> plotting
>>>> > >> the
>>>> > >>> map, showing yield value with legend in tile format.
>>>> > >>>
>>>> > >>> Regards and Thanks,
>>>> > >>> Ranjeet
>>>> > >>>
>>>> > >>>
>>>> > >>
>>>> > >> --
>>>> > >> 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
>>>> > >>
>>>> > >
>>>> > >
>>>> > >
>>>> >
>>>> > --
>>>> > 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
>>>> >
>>>>
>>>>
>>>> --
>>>> Ranjeet  Kumar Jha, M.Tech. (IIT Kharagpur), Ph.D. (USA)
>>>> https://www.linkedin.com/in/ranjeet-kumar-jha-ph-d-usa-73a5aa56
>>>> -----------------------------------------------------------
>>>> Email: *ranjeetjhaiitkgp using gmail.com <ranjeetjhaiitkgp using gmail.com>*
>>>>
>>>>
>>>> *"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://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>>>
>>>
>>
>> --
>> Ranjeet  Kumar Jha, M.Tech. (IIT Kharagpur), Ph.D. (USA)
>> https://www.linkedin.com/in/ranjeet-kumar-jha-ph-d-usa-73a5aa56
>> -----------------------------------------------------------
>> Email: *ranjeetjhaiitkgp using gmail.com <ranjeetjhaiitkgp using gmail.com>*
>>
>>
>> *"Simple Heart, Humble Attitude and Surrender to Supreme Being make our
>> lives beautiful!"*
>>
>>
>
> --
> Ranjeet  Kumar Jha, M.Tech. (IIT Kharagpur), Ph.D. (USA)
> https://www.linkedin.com/in/ranjeet-kumar-jha-ph-d-usa-73a5aa56
> -----------------------------------------------------------
> Email: *ranjeetjhaiitkgp using gmail.com <ranjeetjhaiitkgp using gmail.com>*
>
>
> *"Simple Heart, Humble Attitude and Surrender to Supreme Being make our
> lives beautiful!"*
>
>

	[[alternative HTML version deleted]]