SAS中文论坛

标题: 持续高温天气，大家纷纷表示很热（从7月30日持续分析到8月5日） [打印本页]

作者: shiyiming 时间: 2010-10-22 13:32
标题: 持续高温天气，大家纷纷表示很热（从7月30日持续分析到8月5日）
From sxlion's blog on Sohu

<div style="FONT-SIZE: 14px; LINE-HEIGHT: 160%">
<div style="FONT-SIZE: 14px; LINE-HEIGHT: 160%">
     <img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1842.img.pp.sohu.com.cn/images/blog/2010/8/4/13/20/12aeed116aeg215.jpg" border="0" />    子弹图（Bullet graph/chart）是由斯蒂芬（Stephen Few）从条形图演化而来。据称灵感来自温度计图和仪表盘中的进度条。子弹图提供了一个用于仪表盘可选图形。相对于表示KPI指标的常用仪表盘的gagues和meters，它只需很小的空间，简洁和提供更多的信息，在这里我用它在表示夏天天气的舒适程度。温湿指数是用温度和湿度为变量用来描述人对天气舒适程度的感觉指标。显然，像处于武汉，长沙，上海等南方城市，仅用温度这一指标来衡量天气是否适宜是不可靠的。如果全面考虑人对热环境的感受而言，影响因素很多，除了生理因素，气象因素包括温度、湿度、风速、风向、气压和辐射等等。但是考虑过多的因素使得指标复杂，难以实施使用。于是本博用易于得到的温度和湿度两个变量，使用温湿指数公式：THI= 温度-0.55*(1-相对湿度)*(温度-14.5)。其中温度采用摄氏温度，单位为摄氏度，相对湿度为干湿球温度计读得,单位为%。    
     上图数据测试地点是上海虹桥机场，时间是2010.08.03早上7：30，温度为33度，相对湿度为75%，THI达到30.5。温湿指数THI反映了群体的人对环境的热感受。当夏季的THI小于21.1时，几乎无人有不舒适感；THI>=21.1时，随着THI的增大，具有不舒适感的人增多；THI达到23.9时，约有50%的人感到不舒适；THI达到26.6时，几乎无人感到舒适；特别是THI达到29.5时，酷热难忍。许多国家政府甚至将这一值作为工厂企业等停业下限指标（文献来自1981年）。              <img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://512.img.pp.sohu.com.cn/images/blog/2010/8/4/13/20/12aeeec84e1g215.jpg" border="0" />
由于用温湿指数信息有限，这里用汽车仪表盘来表示风速大小。显然如果有风的话，感觉会舒服一些。
 
<a href="http://pp.sohu.com/photoview-329082938-36248359.html" target="_blank"><img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://512.img.pp.sohu.com.cn/images/blog/2010/8/4/14/26/12aef291418g214.jpg" border="0" /></a>
用曲线图的形式来展现室外的温度和风速，温度基本保持了30度以上，最高温度达到39度，近40度。Y轴上是表示的风速大小。处于断点的地方表示没有风。怪不得我昨天早上起那么早，原来温度在30度，并且还没有风。 <img alt="郁闷" src="http://js3.pp.sohu.com.cn/ppp/images/emotion/base/gloomy.gif" /> 
 
<a href="http://pp.sohu.com/photoview-329084819-36248359.html" target="_blank"><img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1842.img.pp.sohu.com.cn/images/blog/2010/8/4/15/14/12aef3a0ce7g214.jpg" border="0" /></a>
     从上图可以看出，昨天（8月3日）上海几乎全天都处于非常热的状态，有一半的时间处于停业下限以上。看来今年的热时真的很热啊。不知道参观世博的老人和小孩怎么样？因为他们都是对热耐受能力比较差的群体。
据我个人感受上海从上个星期到现在几乎10天一直持续这样的温度（其间少数几天有过阵雨天气），今年真是太热了。
代码:
filename wdata WEBDAV "http://www.wunderground.com/history/airport/ZSSS/2010/08/03/DailyHistory.html?format=1" debug  ; %let day='08.03.2010'; data ex; infile wdata dlm=',' firstobs=3 ; input  @1 Time $ Temperature  dew  RH airp a1 $ a2 $ speed a3 $ a4 $ a5 ~$20. ; dtime=input(cat(&day.,'', time),mdyampm20.); a5=scan(a5,1,'<'); format dtime datetime13.; THI=  temperature-0.55*(1-rh/100)*(temperature-14.5) ; run;
proc print;run;
title h=1 font=bold  "Shanghai 2010.08.03 Outside temperature and wind speed"; footnote  j=l h=0.7 "DATA from <a class="postlink" href="http://www.wunderground.com">www.wunderground.com</a>"&nbsp;&nbsp;&nbsp; ; footnote2 j=l h=0.7 "Code avaliable at sxlion.blog.sohu.com"; ods listing  ; ods graphics / reset width=350px height=250px imagefmt=jpg; proc sgplot CYCLEATTRS NOAUTOLEGEND ;    series x= dtime y=temperature/ LINEATTRS=( color=red) BREAK ;    series x= dtime y=speed/ Y2AXIS LINEATTRS=( color=cx9AC0F1) BREAK ; xaxis  display=(nolabel) ; yaxis  GRID   values =( 0 to 40 by 10 ) label="Temperature (℃)"                MAX=45 min=0 VALUESHINT ; y2axis  values =( 0 to 12 by 3 ) label="Wind speed (km/h)"                MAX=15 min=0 VALUESHINT ; run; title h=1 font=bold  "THI of Shanghai at 2010.08.03 "; title2 h=0.2 "  "; proc sgplot CYCLEATTRS NOAUTOLEGEND ;    series x= dtime y=thi/ LINEATTRS=(color=orange thickness =3)  ; xaxis  display=(nolabel) ; yaxis     values =( 15 to 30 by 5 ) label="THI"                MAX=32 min=15 VALUESHINT  ; REFLINE 21.1 23.8 26.6 29.5 / AXIS= y  LABEL=("Cool" "Hot" "Very hot" "Warning" ) LABELLOC= outSIDE        LINEATTRS =(color=red THICKNESS=0.5); ; run;
 附带列出这几天的THI指数趋势吧。
  <img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1832.img.pp.sohu.com.cn/images/blog/2010/8/5/15/20/12af4813331g214.jpg" border="0" /><img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1832.img.pp.sohu.com.cn/images/blog/2010/8/5/15/20/12af480c8d8g213.jpg" border="0" /><img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1801.img.pp.sohu.com.cn/images/blog/2010/8/5/15/18/12af47f7789g215.jpg" border="0" /><img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1862.img.pp.sohu.com.cn/images/blog/2010/8/5/15/17/12af47e2971g214.jpg" border="0" />
<img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://512.img.pp.sohu.com.cn/images/blog/2010/8/5/15/17/12af47e60a1g215.jpg" border="0" />
<img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://512.img.pp.sohu.com.cn/images/blog/2010/8/5/15/17/12af47e5e8fg213.jpg" border="0" />
<img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1832.img.pp.sohu.com.cn/images/blog/2010/8/5/15/18/12af47edd6bg214.jpg" border="0" />
   连续7天的温湿指数曲线如上面所示，根据我本人的体会，特别是早上5：00-8：00这一段，因为在睡觉，并且没有空调。从上个星期开始，我明显的感觉到，早上越来难越睡觉了，特别是2日，3日，因为太热我很早就爬起来冲凉。3号我5点爬起来，就再也没有睡了，直到昨天晚上我才感觉风吹得特别凉爽，睡眠很好，并且一直睡到8点。从本人的亲身体会可以看出THI还是很能反映情况的，今天8点之前THI都保持在29.5以下，前些天都是6点多就开始越线了。
   另贴出这7天地温度和风速图，气温也也表现同样的规律，今天8点之前都是保持在30度。
   <img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1801.img.pp.sohu.com.cn/images/blog/2010/8/5/15/28/12af48806c4g214.jpg" border="0" /><img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1872.img.pp.sohu.com.cn/images/blog/2010/8/5/15/28/12af4882762g215.jpg" border="0" /><img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1811.img.pp.sohu.com.cn/images/blog/2010/8/5/15/28/12af4881882g213.jpg" border="0" /><img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1872.img.pp.sohu.com.cn/images/blog/2010/8/5/15/28/12af48879b7g214.jpg" border="0" /><img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1872.img.pp.sohu.com.cn/images/blog/2010/8/5/15/28/12af4888dd1g215.jpg" border="0" /><img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://1842.img.pp.sohu.com.cn/images/blog/2010/8/5/15/28/12af488774bg213.jpg" border="0" /><img style="DISPLAY: block; MARGIN: 0px auto 10px; TEXT-ALIGN: center" alt="" src="http://512.img.pp.sohu.com.cn/images/blog/2010/8/5/15/28/12af488df72g214.jpg" border="0" />
参考：</div>
<a href="http://en.wikipedia.org/wiki/Bullet_graph"><a class="postlink" href="http://en.wikipedia.org/wiki/Bullet_graph">http://en.wikipedia.org/wiki/Bullet_graph</a></a>
<a href="http://www.wunderground.com/" target="_blank"><a class="postlink" href="http://www.wunderground.com">http://www.wunderground.com</a></a>
王远飞，沈愈 1998 上海市夏季温湿效应与人体舒适度华东师范大学学报(自然科学版) No3, pp 61-66.</div>

欢迎光临 SAS中文论坛 (https://mysas.net/forum/)