weixin_44188622的博客

要想在R里玩转空间数据，是要先搞清楚坐标系怎么操作。接下来一起来学学Coordinate System in R网页的内容吧。查看坐标系方法1：library(sp)st_crs()方法2：library(raster)crs()两种命令分别返回.shp和栅格的坐标系。定义坐标系有EPGS编码和proj4的坐标系字符串两种方法。并非所有坐标系都有EPGS编码，但他们都有proj4字符串。理解Proj4语法每一个描述前都有一个"+"，描述其坐标系、时区、大地水准面、椭球体、单位等信息。

LOL colors

geom_sf()绘制地图坐标转化坐标转换在上一篇设置ylim中有提及，一些网站也附在里面。总体来说，用sf::st_read()读入数据后，首先用st_crs()看看有无空间参考坐标信息。如果没有，若不想用coord_sf()的自定义设置，可以用st_set_crs()，进行设置。如果有，则用st_transform_crs(x,crs)进行转换。crs参数直接用epgs码即可，在上一篇内容附...

用geom_sf()绘制地图使用geom_polygon()绘制地图，每次都会提示该功能已经deprecated，而sf包配合ggplot2的geom_sf()函数则能快速地绘制地图，且无须设置long,lat,group=id这些参数。代码如下：ws <- ggplot() + geom_sf(data = fortify(wsdiv1), colou...

R中转换日期遇到的若干问题：如何将字符串转换为日期？dt <- as.Date("2020-01-03")如何将Day of the year转换为Date？dt <- as.Date(100,origin='2020-01-01')其中origin是必须参数。还可以设置format参数，以进一步format日期。如何更改日期格式？format(dt,"%b%d%Y"...

ggplot2默认绘图边框、网格白色，背景浅灰。一般统计图也还勉强算作美观，但是绘制中国地图时则非常碍事。而有时在ggplot后设置theme(）这一操作会莫名其妙失灵。而众多解决方案则直接移除了所有边框和网格。倒腾半天，终于找到了设置四周边框黑色，背景白色的最简便方法。绘制出了美观的中国1级流域地图。ggplot() + geom_polygon(data = fortify(wsdiv...

写在前面：本文上接上一篇博文，是Emmanuel Paradis《R for beginners》的笔记、部分翻译和实践代码的整理，部分例子有所修改。虽然本书已经足够精简，但是便于查阅和复习，便有了该系列。由于笔者之前对R有过简单了解和使用，部分内容有所精简。Data with R1. Objects所有的对象有两个独有的属性：mode和length。类型和长度。有numeric, ch...

《R for beginners》|R语言基本概念PreambleR 语言和Python, MATLAB等高级语言一样，具有强大的数据分析和可视化能力。而对于科研狗来说，更是不可不掌握的必备技能。作为一个科研入门阶段的本科小白，搜罗多方信息，对比之后，选择了这本短小精悍的**《R for beginners》**进行一个系统的学习入门。这本《R for beginners》作者是Emmanu...