代码妈妈杯大数据数学建模.pdf

Python绘绘图图实实现现台台风风路路径径可可视视化化代代码码实实例例

台风是重大灾害性天气，台风引起的直接灾害通常由方面造成，狂风、暴雨、风暴潮，除此以外台风的这些灾害极易诱发城市内涝、房屋倒塌、山洪、

泥石流等次生灾害。正因如此，台风在科研和业务工作中是研究的重点。希望这次台风路径可视化可以给予大家一点点帮助。

台台风风路路径径的的获获取取

中国气象局(CMA)

中国气象局(CMA)的台风最佳路径数据集(BST)，BST是之后对历史台风路径进行校正后发布的，其经纬度、强度、气压具有更高的可靠性，但是时间分

辨率为6小时，部分3小时，这一点不如观测数据。下载地址：

/

温州台风网

温州台风网的数据是实时发布数据的记录，时间分辨率最高达1小时，对于台风轨迹具有更加精细化的表述。下载地址：

/

示例

导入模块并读取数据，使用BST的201年台风路径数据作为示例，已经将原始的txt文件转换为xls文件。

importos,glob

importpandasaspd

importnumpyasnp

importshapely.geometryassgeom

importmatplotlib.pyplotasplt

frommatplotlib.imageimportimread

frommatplotlib.animationimportFuncAnimation

importmatplotlib.linesasmlines

importcartopy.crsasccrs

importcartopy.featureascfeat

fromcartopy.mpl.tickerimportLongitudeFormatter,LatitudeFormatter

importcartopy.io.shapereaderasshpreader

importcartopy.io.img_tilesascimgt

fromPILimportImage

importwarnings

warnings.filterwarnings(ignore)

df=pd.read_csv(./201typhoon.csv)

定义等级色标

defget_color(level):

globalcolor

iflevel==热带低压orlevel==热带扰动:

color=#FFFF00

eliflevel==热带风暴:

color=#6495ED

eliflevel==强热带风暴:

color=#3CB371

eliflevel==台风:

color=#FFA500

eliflevel==强台风:

color=#FF00FF

eliflevel==超强台风:

color=#DC143C

returncolor

定义底图函数

defcreate_map(title,extent):

fig=plt.figure(figsize=(12,))

ax=fig.add_subplot(1,1,1,projection=ccrs.PlateCarree())

url=/wmts-geo/wmts.cgi

layer=BlueMarble_ShadedRelief

ax.add_wmts(url,layer)

ax.set_extent(extent,crs=ccrs.PlateCarree())

gl=ax.gridlines(draw_labels=False,linewidth=1,color=k,alpha=0.5,linestyle=--)

gl.xlabels_top=gl.ylabels_right=False

ax.set_xticks(np.arange(extent[0],extent[1]+5,5))

ax.set_yticks(np.arange(extent[2],extent[3]+5,5))

ax.xaxis.set_major_formatter(LongitudeFormatter())

ax.xaxis.set_minor_locator(plt.MultipleLocator(1))

ax.yaxis.set_major_formatter(LatitudeFormatter())

ax.yaxis.set_minor_locator(plt.Multi