import pandas as pd
import pygrib as pg
import os
import json
import cartopy.crs as ccrs
import cartopy.feature as cfeature
import matplotlib.pyplot as plt
import metpy.calc as mpcalc
from metpy.units import units
import numpy as np
import xarray as xr
import urllib.request
import datetime
from cma.music.DataQueryClient import DataQueryClient
import cartopy.io.shapereader as shpreader
import glob
from cartopy.mpl.gridliner import \
LONGITUDE_FORMATTER, LATITUDE_FORMATTER
# from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from dateutil.relativedelta import relativedelta
import numpy as np
import sys
from hjn.mkNCHJN import mkNCCommonUni,dataClass,clip,clipLat,clipLon,envelope
areaInfoZL={"area":"sc","evn":envelope(33,28,104,108)}
# 获取雨量数据方法,参数:开始时间,结束时间
def CountRainData(startTime, endTime):
client = DataQueryClient()
# 2. 调用方法的参数定义,并赋值
# 2.1 用户名&密码
userId = "xxxxx"
pwd = "xxxx"
# 2.2 接口ID
interfaceId = "getRadaFileByTimeRangeAndStaId"
# 2.3 服务节点ID
serverId = "NMIC_MUSIC_CMADAAS"
# 2.4 接口参数,多个参数间无顺序
# 必选参数 (1)资料:中国地面逐小时资料; (2)检索要素:站号、站名、小时降水、气压、相对湿度、能见度、2分钟平均风速、2分钟风向;
# (3)检索时间;(4)排序:按照站号从小到大;(5)返回最多记录数:10。
params = {'dataCode': "RADA_L2_FMT", \
'timeRange': "(%s,%s]" % (startTime, endTime), \
'staIds': "xxxxx", \
'dataFormat': "json"}
# 3. 调用接口
# result = client.callAPI_to_array2D(userId, pwd, interfaceId, params)
result = client.callAPI_to_serializedStr(userId, pwd, interfaceId, params, serverId)
print(result)
root = json.loads(result)
if root['returnCode'] == str(0):
CountList = []
CountUrlList = []
for row in root['DS']:
countfile = row['FILE_NAME']
countfileurl = row['FILE_URL']
print(countfile)
CountList.append(countfile)
CountUrlList.append(countfileurl)
# for list_ragion in CountList:
# str1 = list_ragion[0] + ' ' + list_ragion[1] + ' ' + list_ragion[2] + ' ' + list_ragion[3] + ' ' + list_ragion[4] + ' ' + list_ragion[5] + '\n'
# with open('result.txt', 'a+') as file_object:
# file_object.write(str1)
return CountList,CountUrlList
list_region,listUrl= CountRainData('20230727131000', '20230727152000') # [20210807040000,20210808160000]
print(listUrl[0])
print(listUrl)
for i in range(1,len(list_region)):
output = "D:\\xiaoshui\\20230727\\%s" %list_region[i]
if os.path.exists(output ):
continue
os.makedirs(os.path.dirname(output), exist_ok=True)
if not os.path.exists(output):
urllib.request.urlretrieve(listUrl[i], output)