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

Thu May 12 13:04:42 CEST 2022

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!"*

	[[alternative HTML version deleted]]