欢迎来到奇趣后备箱! 手机访问:pm25分布图pm25分布图

休闲八卦

pm25分布图pm25分布图

吕云来自:陕西省 安康市 白河县 时间:2019-01-10 20:14 影响: 281162人

我们找到第1篇与pm25分布图pm25分布图有关的信息,分别包括:

以下是的一些我们精选的pm25分布图pm25分布图

@Conda 博文已更新,新文章见http://blog.sina.cn/dpool/blog/s/blog_403aa80a0102vrti.html?vt=4
更多R语言文章:http://r.lovedatasci.com
@Conda其实这段代码是很久以前的了,大概是10月底的代码,只是因为那段时间忙着考试、写论文以及足球赛的事情,就一直耽搁下来了,前两天这些事情全部告一段落,想着把这些代码跟大家分享一下。
废话不多说,先上图:
用R语言抓取pm2.5数据绘制全国pm2.5分布图
这张图是前天的全国pm2.5的情况,也就是我们决赛夺冠那天的情况,可惜长沙那天是一块偏红的大圈圈,真是苦了我们兄弟们。这张图里面颜色越偏红,圈越大表示pm2.5越高。
绘制这张图需要以下几个条件:
1、计算机能联网;
2、要安装RgoogleMaps、RCurl、XML这三个R包。
数据来源是http://www.cnpm25.cn,代码主要分两部分,第一部分是负责上http://www.cnpm25.cn抓取数据,第二部分是负责利用google的placeAPI获取IP地址并绘图。
下面是代码,具体代码的意思就请各位看官自行解决,中间有些冗余的代码,不过我想不影响大家的使用,我也没有花太多时间去修饰,毕竟这是个兴趣爱好,而且个人事情也太多。欢迎大家前来讨论交流。
用R语言抓取pm2.5数据绘制全国pm2.5分布图
library(RgoogleMaps)
library(RCurl)
library(XML)
Start = Sys.time()
day = 13
skip = 0
website = 'http://www.cnpm25.cn'
pm = getURL(website,.encoding='UTF-8')
pmweb = htmlParse(pm,encoding='UTF-8')
pmtotal = getNodeSet(pmweb,'//div[@class='warp']//a')
pmurl = sapply(X=pmtotal,FUN=xmlGetAttr,name='href')
pmcityA = getNodeSet(pmtotal[[1]],'//strong')
pmcity = sapply(X=pmcityA,FUN=xmlValue)
pmcity = gsub(pattern='[A-Z]',replacement='',x=pmcity)
remove.160 = function(x)
{
x = strsplit(x=x,split='')[[1]]
index = sapply(X=x,FUN=utf8ToInt)
x = x[index!=160&index!=32]
return(paste(x,collapse=''))
}
pmcity = sapply(X=pmcity,FUN=remove.160)
pmcity = as.character(pmcity[pmcity!=''])
pmurl = pmurl[!grepl(pattern='article',x=pmurl)]
index =match(c('黄山风景区','菏泽','九华山风景区','恩施','黄冈','随州'),table=pmcity)
pmcity = pmcity[-index]
pmurl = pmurl[-index]
result = NULL
for( i in pmcity)
{
tmp = China[grepl(pattern=i,x=China$County),]
if(nrow(tmp)>1)
tmp = tmp[1,]
result = rbind(result,tmp)
rm(tmp)
}
data = NULL
date = NULL
for( i in pmurl)
{
#i = pmurl[45]
city = pmcity[match(i,table=pmurl)]
cat(paste('Get data for',city,'... ...'))
cityweb =getURL(url=paste(website,i,sep='/'),.encoding='UTF-8')
cityweb = strsplit(x=cityweb,split='flashvalue')[[1]]
cityweb = cityweb[grepl(pattern='set name',x=cityweb)]
cityweb = cityweb[grepl(pattern='月',x=cityweb)]
pm = NULL
if(length(cityweb) == 0)
{
cat('Date is incorrect... skip...')
next
}
for(j in cityweb[day])
{
j = cityweb[day]
tmp = strsplit(x=j,split='\s')[[1]]
name = tmp[grepl(pattern='name',x=tmp)]
name = strsplit(x=name,split='\'')[[1]][2]
value = tmp[grepl(pattern='value',x=tmp)]
value = as.numeric(strsplit(x=value,split='\'')[[1]][2])
rm(tmp)
pm = cbind(pm,value)
}
if(!is.null(date))
{
if(date != name)
{
cat('Date is incorrect... skip...')
skip = skip + 1
next
}
}else{
date = name
}
data = rbind(data,cbind(city=city,pm))
}
lon = NULL
lat = NULL
for( i in data[,1])
{
tmp = China[grepl(pattern=i,x=China$County),]
tmp = tmp[1,]
lon = c(lon,tmp[1,4])
lat = c(lat,tmp[1,5])
}
#col.base = c('green','orange','red')
pm = as.numeric(data[,2])
#col =ifelse(test=pm<100,yes=1,ifelse(pm>200,yes=3,no=2))
#col = col.base[col]
tp.pm = GetMap.bbox(lonR=lon,latR=lat)
pm.max = max(pm)
red = seq(from=0,to=1,length.out=pm.max)
green = seq(from=1,to=0,length.out=pm.max)
blue = seq(from=0,to=0,length.out=pm.max)
col = rgb(red=red,green=green,blue=blue,alpha=0.5)
col = col[pm]
cex = pm/100 + 5
suppressWarnings(expr=PlotOnStaticMap(MyMap=tp.pm,lat=lat,lon=lon,pch=20,col=col,cex=cex,mar=c(0,0,2,0)))
title(main=paste(date,'全国',length(pm),'个主要城市PM2.5分布情况',sep=''))
text.col = rgb(red=1,green=1,blue=1,alpha=0.9)
text(x=210,y=-250,labels=paste('数据来源:',website,sep=''),cex=0.8,col=text.col)
text(x=210,y=-230,labels='R语言 RgoogleMaps@Conda',cex=0.8,col=text.col)
Stop = Sys.time()
Cost = Stop - Start
cat(paste('It cost',Cost,'minutes.'))
cat(paste('Gather',length(pm),'cities pm 2.5 data, andskip',skip,'city for incorrect date.'))
更多R语言文章:http://r.lovedatasci.com
上一篇:诺基亚7373主题
下一篇:没有了
  • 魔禁同人游戏:十大影响力极高的Galgame排行,要优雅不要污

    ::::
  • www.10xo.com_www10xocom

    我们找到第1篇与www.10xo.com_www10xocom有关的信息,分别包括:

    以下是的一些我们精选的www.10xo.com_www10xocom

    “什么事儿?”
    陈丰没料到拉拉会这态度,楞了一下说:“没什么特别的事。”
    拉拉冷冰冰地说:“那你想说什么?”
    陈丰越发觉着不对劲,赔着小心说:“怎么了?要不我先听WWW)10XO_CoM你说吧。”
    拉拉没好气地嗔怪道:“是你打给我的,你让我说啥?”
    陈丰解释说:“我就找你随便聊聊天。”
    拉拉话中带刺道:“那你想了解什么呢?”
    陈WWW)10XO_CoM丰辩白说:“我真没有什么想了解的啊。你今天怎么了?”
    拉拉根本不信,她不耐烦起来,不觉地声音就像刀片划过玻璃那样刺人耳朵:“你到底什么WWW)10XO_CoM事儿吧?”WWW)10XO_CoM
    陈丰也急了:“没事儿就不能找你聊天吗?”
    拉拉怀疑地说:“那WWW)10XO_CoM好吧,你想聊什么?”
    陈丰见不WWW)10XO_CoM是个事,就说:“拉拉,我怎么觉得你对我有意见?是我做错了什么?”
    ......
  • dnf魔道加点模拟器

    ::::
  • 99crw.com_99crwcom

    我们找到第1篇与99crw.com_99crwcom有关的信息,分别包括:

    以下是的一些我们精选的99crw.com_99crwcom

    WWW_99crw_CoM的们置位分得实好最,的们餐用湾岸在“
    。仆女小分得实她和莎尔艾分得实城丁马及以,的们莎罗妹妹聊上网分得实后身林罗在跟。林罗,的们哥三分得实他,的们人熟分得实杰罗聊上网然竟人分得实首为而。来人个几进走又外厅大,的们候时分得实躇踌而字文丽华分得实单菜满着对正菲苏在就
    。量数分得实钱掏方对WWW_99crw_CoM心关只们他,的们里这在少至,的们个哪聊上网论无。祖世二分得实粹纯聊上网就么要,的们贯万缠腰分得实真聊上网么要,的们人惊么那都费消分得WWW_99crw_CoM实童门给赏打人个两这,的们辣神眼然当者侍。花鲜分得实香淡种各聊上网满围周,的们水流山假有侧一,的们一之置位分得实雅优最厅大个整聊上网里这,的们桌小人双处一WWW_99crw_CoM到来杰罗着带者侍
    ?步一WWW_99crw_CoMWWW_99crwWWW_99crw_WWW_99crw_CoMCoM_CoM进更人两明说以可聊上网不聊上网这,的们何如论无,的们来起心开然突却她,的们后之秒几在就但
    ......
  • 8月份适合去哪里旅游_8月份适合哪里旅游

    ::::
  • 快乐美美往前冲_快乐美美往前

    我们找到第1篇与快乐美美往前冲_快乐美美往前有关的信息,分别包括:

    以下是的一些我们精选的快乐美美往前冲_快乐美美往前

      最近,关于湖南卫视芒果台的《快乐美美往前冲》的一位超大胸MM在闯关过程中不幸掉落水中,而由于冲击力太大(不排除救人男恶意行为)其上身紧身衣“缩水”,整个胸部几乎全部外露!此情形被现场的观众拍照,并迅速发布到网上,网友们戏称为“落水门”事件。
      《快乐美美向前冲》是湖南经视的一档大型户外竞技秀节目。湖南经视联袂湖北经视特别制作的大型户外竞技秀节目—快乐美美向前冲。在快乐美美向前冲0627期间,发生了一件意想不到的事情:湖北选手廖慧敏在闯关的过程中掉入水中,由于冲击力太大,其上身紧身衣“缩水”,整个胸部几乎全部外露!在网络上落水视频一下子疯传。
    湖南卫视《快乐美美往前冲》廖慧敏落水露胸【落水门】

    钟佳宜遭遇歌迷行凶,衣服被扒当场走光【扒衣门】
    刘嘉亮MV玩双飞,删节版泄露遭疯狂下载。。。
    ......
  • 西班牙足球明星

    ::::
  • 害怕恨会像爱一样浓_害怕恨会像爱一样

    我们找到第1篇与害怕恨会像爱一样浓_害怕恨会像爱一样有关的信息,分别包括:

    以下是的一些我们精选的害怕恨会像爱一样浓_害怕恨会像爱一样

    昨天刚刚从深圳参加完“中国音乐流行榜颁奖晚会”回到北京,我带着我的新歌《害怕恨像爱一样浓》获得了“08-09年度十大新人新歌奖”,在此我要向公司、我的家人、朋友还有支持我的歌迷朋友们表达感谢之情~是你们的帮助、支持和鼓励,使我获得了今天的成绩。新闻网址:http://music.yule.sohu.com/20090929/n267074208.shtml这首《害怕恨像爱一样浓》是我个人非常喜欢的一首歌,歌词带着依恋、无奈还有点小忧伤,希望可以跟年轻的歌迷朋友们一起分享。大家可以在中国音乐流行榜的官网上抢先听哦~
    我的获奖照
    《害怕恨像爱一样浓》获“最佳新人新歌奖”
    ......
  • 积分入户体检:2018年办理积分入户需要体检哪些项目,想入深户的

    :::
  • e71 code_e71code

    我们找到第1篇与e71 code_e71code有关的信息,分别包括:

    以下是的一些我们精选的e71 code_e71code

    教你查询/修改诺基亚手机CODE码 并附CODE码意思说明
    说明:本文内容适用于几乎所有的S60系统手机(含S60v5,如5230、5530、5800、N97、X6、C6)与Maemo手机(含N900)。
    诺基亚CODE码是什么意思:
      诺基亚手机的CODE代表这台手机的身份,由于就算是同一种型号的手机,在同一地区有时也会通过不同的渠道或者电信运营商来进行销售,诺基亚为了更进一步的区分不同渠道(甚至颜色)的手机,专门设定了一个CODE属性。相同CODE值的手机基本上所有的配置都是一模一样的(甚至包括机身颜色)。
      相同类型标志符(如何查询Nokia手机类型标志符)的手机由于发售国家、渠道(比如是否电信定制机)、机身颜色的不同,也会有不同的CODE值。这些CODE值可通过软件的方式来任意更改。
    更改诺基亚CODE吗的好处:
      由于诺基亚的手机版本升级并不是全球同步更新,部分区域会提前升级。
    ......
  • b7722机_b7722

    我们找到第1篇与b7722机_b7722有关的信息,分别包括:

    以下是的一些我们精选的b7722机_b7722

    一、机前准备
    1、 拆掉手机电池,USB线也不要连接电脑。
    二、运行机软件
    1、 打开机文件夹里的Flash loader 7.2.4 SEC 2.4 Lite And Sle0.4.exe软件,如下图:
    三星GT-B7722机教程2(机步骤)
    2、 出现如下界面:
    三星GT-B7722机教程2(机步骤)
    3、 点SET USB ID按钮,打开文件“usbInstanceID.usbid”
    三星GT-B7722机教程2(机步骤)
    4、 点SET MODEL按钮,打开机软件文件夹中的B7722_Sorrento_Setting_v01.mdl
    ......
  • 陈滢微博:又一TVB女高材生,曾与王儿子相恋,今住千万豪宅、开

    ::::
  • WWW.SMICAP.COM_WWWSMICAPCOM

    我们找到第1篇与WWW.SMICAP.COM_WWWSMICAPCOM有关的信息,分别包括:

    以下是的一些我们精选的WWW.SMICAP.COM_WWWSMICAPCOM


    最初三天。 缓和
    伤害,我以华语抓住一本自己的研究书。 被最近几星期的学校,它
    明显没有我会以拉丁文语制造比一个弱的 D 好的方法,但是
    充足通过了。 我藏了在那一本暴虐外国看书后面WWW*smicap*COM的我的耻辱
    与厚地,在全掩护的黑色中文字体。 我埋葬拉丁文的所有想法进入
    变酸葡萄而且在那里坐了而且学习华人改为!
    中国水手不说拉丁文语 WWW*smicap*COM
    为华语放弃拉丁文是我自己形式的少年的蔑视。 然而,我有
    自从几乎每天在某些方面用了华语。 我对偶然的好奇心承认关于
    那些一位从莱斯利拉丁文分类的正在最近做与他们的
    拉WWW*smicap*COM丁文。
    在夏天期间假期我们去迈阿密海滩叁观我的祖父或祖母。 在一之上
    跌倒,因为比尔叔父从在迈阿密中的火车站到迈阿密海滩驾驶我们,我们通过一
    大群体的游行水手。
    ......
  • 铁甲小宝女主角

    ::::
  • www.000kkk.com_www000kkkcom

    我们找到第1篇与www.000kkk.com_www000kkkcom有关的信息,分别包括:

    以下是的一些我们精选的www.000kkk.com_www000kkkcom

    wWw-000kkk-com
    惟独灰衣老面色极为难看李慕婉是wWw-000kwWw-000kkk-comkk-com他的准儿媳妇此时居然被那人一wWw-000kwWw-000kkk-comkk-com字之下就顺从的走来这让他的老脸立刻有些难堪。
    那老者目露可惜之色。摇头说道:“可wWw-000kkk-com惜。若是无人指点怕是此生无法结丹。”
    “此阵是古wWw-000kkk-com传wWw-000kkk-com送阵当年我师父现此处wWw-000kkk-com洞府时就现了这个阵法此阵虽然因为在石室内的关系没有风吹雨打所以保持完整但开启古传送阵却是需要极品灵石。极品灵石在整个修魔海内邱某尚没有听到wWw-000kkk-com有谁拥有wWw-000kkk-com所以这阵法也就从来没wWw-000kkk-com有开启过。wWw-000kkk-com”邱四平转过头看了那阵法一眼平淡的说道。%DwWw-000kkk-com%A只有wWw-000kkk-com那种在壮年期自愿以生命为代价才可以产生传承。
    随着那些神识一一的自爆四周地飓风wWw-000kkk-com。
    ......
  • pm25分布图pm25分布图相关文章