区域降水实况栅格图API

From CaiyunWiki
Revision as of 10:31, 14 December 2018 by Liuzhonghua (talk | contribs)
Jump to: navigation, search

开发者token需要额外开通权限,请联系商务:api@caiyunapp.com

1.API概述 :

  • 中国区域降水实况栅格图 API 提供三种文件:降水强度图、可用区域mask图、雨雪类型mask图;
  • 文件名为格林威治时间,提供最近十帧(45分钟历史+当前最新),每5分钟更新1次(每小时12次,从00H:00M开始);
  • 等经纬度投影,栅格分辨率 (6000 * 7500) , 纬经度分辨率近似为(0.009, 0.0105),(经度分辨率0.0105非准确值);
  • 数据第1维 (3.9079N ~ 57.9079N)纬度,第2维 (71.9282E ~ 150.6026E)经度;
  • 栅格数据索引和经纬度之间的转换详见下面示例代码;
  • 提供的数据如下:
   1.降水强度图,数据uint8存储,除以255为彩云降水强度(等于0~1.0之间);
   2.降水类型图,数据uint8存储,255代表此像素为雪(固态),100代表此像素为雨(液态),0代表不可用;
   3.可用区域图,数据uint8存储,255代表可用(不缺测),0代表不可用。

2.API调用方式:

2.1 API接口:http://api.caiyunapp.com/v1/radar/cndata?token=$YOUR_TOKEN
2.2 调用API 获取十帧图的url:降水强度图(共10张)、降水类型图(共10张)、可用区域图(共10张)
2.3 降水强度图示例:http://cdn.caiyunapp.com/weather/radar/cnmap_precp_201810250540_xxxxxxxxxx.png4.
2.4 降水类型图示例:http://cdn.caiyunapp.com/weather/radar/cnmap_precptype_201810250540_xxxxxxxxxx.png5.
2.5 可用区域图示例:http://cdn.caiyunapp.com/weather/radar/cnweight_raw_201810250540_xxxxxxxxxx.png

3.API返回示例:

URL: http://api.caiyunapp.com/v1/radar/cndata?token=$YOUR_TOKEN

返回示例如下:

{

status*: "ok",

cndata*: [

         {

          ## 降水强度图(UTC+0 格林威治时间)  

          data_url*: "cdn.caiyunapp.com/weather/radar/ cnmap_precp_201810250540_xxxxxxxxxx.png",

          ## 降水类型图(UTC+0 格林威治时间) 

          data_type_url*: "cdn.caiyunapp.com/weather/radar/ cnmap_precptype_201810250540_xxxxxxxxxx.png",

          ## 可用区域图(UTC+0 格林威治时间)

          weight_url*: "cdn.caiyunapp.com/weather/radar/cnweight_raw_201810250540_xxxxxxxxxx.png",

          datatime*: 1540446000, # 降雨图的时间戳(UTC+0 格林威治时间)

          lat*: [3.9079,57.9079], # 纬度范围

          lng*: [71.9282,150.6026], # 经度范围

          unit*: "grey scale" # 灰度单位,降雨图数据除以255为彩云降水强度;

                             #   可用范围数据中,255代表可用,0代表不可用

          },

          ...

        ]

4.Python参考代码:

##根据用户经纬度,获取彩云降水(雪)强度
##lat_interval, lng_interval = 0.009, 0.0105
from scipy.misc import imread
VALID = 255
NORM = 255
lats = [3.9079, 57.9079]
lngs = [71.9282, 150.6026]
lat_size, lng_size = 6000, 7500
def latlng2xy(lat, lng):
y = int((lng - lngs[0]) / (lngs[1] - lngs[0]) * lng_size)
x = int((lats[1] - lat) / (lats[1] - lats[0]) * lat_size)
return x, y
data = imread('cnmap_precp_201810250540_7bdef9f9ac.png')
weight = imread('cnweight_raw_201810250540_7bdef9f9ac.png')
lat_size, lng_size = data.shape ##(6000, 7500)
lat, lng = 40.35, 116.31 ##查询用户位置降水情况
x, y = latlng2xy(lat, lng)
if weight[x, y] == VALID:
result = float(data[x,y]) / NORM