如何使用Pandas,GeoPandas和ArcGIS实现数据处理和绘图

发布时间：2021-11-30 14:49:44 作者：小新
来源：亿速云阅读：836

这篇文章主要介绍如何使用Pandas,GeoPandas和ArcGIS实现数据处理和绘图，文中介绍的非常详细，具有一定的参考价值，感兴趣的小伙伴们一定要看完！

第一部分是学习使用pandas进行数据数据。pandas中有两大数据结构，分别是series序列和DataFrame数据框。使用pandas的实验对象是之前电脑本地的一些气象数据，存储格式为.xlsx。实验数据仅仅是实验数据而已，用它来练习pandas还是不错的，至少目前。

pandas数据分析的大概步骤也不复杂，导入库—>读取数据—>数据操作(切片，筛选，组合，排序，批处理等)—>写出数据到csv或xlsx或txt纯文本。

pandas的数据读取很容易：

import pandas as pddata = pd.read_excel(r'D:\inuyasha\codeX\Python\气象数据1960-2018\1960.xlsx')
series = data['平均气温'] #序列，就是1,2,3这样排列下去的series_to_dataframe=pd.DataFrame(series) #把series数据结构转为数据框结构，pd两大数据结构series_to_dataframe

如何使用Pandas,GeoPandas和ArcGIS实现数据处理和绘图

重命名，增加索引操作：

datare = data.rename(columns={'平均气温':'平均温'}) #使用rename方法重命名列的名字，字典，1是原名，2是更换后的名字datare['INDEX'] = datare.index #原基础上增加索引，.index方法，名INDEXdatare['均温'] = datare['平均温'] #没有替代原来的那一列，而是新增一列，一模一样的，末尾添加datare.drop(columns='INDEX') #删除原来列datare

如何使用Pandas,GeoPandas和ArcGIS实现数据处理和绘图

删除某一行，也可以删除多行，index是行索引，axis指针对所有列：

df1 = datare.drop(index=0,axis=1)
据索引删除多行，index是行索引，axis指针对所有列df2 = datare.drop(index=[0,1,2,3],axis=1)

定位和提取数据：

datare['区站号'].isin([50136]) #isin函数判断这一列中是否存在这个数，数字这不要加引号，否则就成字符串了datare[datare['区站号'].isin([50136])] #在原始数据中取出定位到的数据datare[-datare['区站号'].isin([50136])] #在原始数据中取出定位到的数据,前面加一个横线表示反选

如果想一次就提取多列呢，那就把需要的列以list形式放进去

data_multicolu2 = datare[['区站号','纬度','经度','年','月','日','平均温']] #列表形式把需要的多列切片提取出来data_multicolu2

如何使用Pandas,GeoPandas和ArcGIS实现数据处理和绘图

排序：

datare.sort_values(by='平均温',ascending=False)

还有更多功能，这里就不一一写出来了。多用几次就熟悉了。之前想的一个问题，如何提取多个xlsx里面的某些数据，这不就发挥作用了嘛哈哈。

第二部分是Geopandas读取shp数据+绘图展示+添加经纬度标识(EWSN)

脚本如下：

import geopandas as gpdimport numpy as npimport matplotlib.pyplot as pltimport cartopy.crs as ccrsfrom cartopy.mpl.ticker import LongitudeFormatter, LatitudeFormatterimport cmaps#import arcpy

gd_boundary = gpd.read_file('D:\ArcGISgeobase\guangdata\gd_boundary.shp')gd_road = gpd.read_file('D:\ArcGISgeobase\guangdata\gd_road.shp')#fig,ax = plt.subplots(figsize=(16,12)) #设置画布大小plt.figure(figsize=(16, 12))ax = plt.subplot(projection=ccrs.PlateCarree())#全局设置plt.rcParams['font.family'] = 'Times New Roman'plt.rcParams['font.size'] = 16
gd_boundary.plot(ax=ax,facecolor='none',edgecolor='k',linewidth=0.99) #边界设置gd_road.plot(ax=ax,edgecolor='b',cmap=cmaps.posneg_2,lw=0.3,legend=False) #路网设置
ax.set_yticks(np.arange(20, 26, 1.0),crs=ccrs.PlateCarree())ax.set_xticks(np.arange(109, 118, 1.0),crs=ccrs.PlateCarree())
ax.xaxis.set_major_formatter(LongitudeFormatter()) #经纬度标识显示，EWSNax.yaxis.set_major_formatter(LatitudeFormatter())

#ax.set_axis_on() #打开矩形边框ax.set_title('Guangdong_Roads',fontsize=16) #设置标题
plt.grid() #设置网格

展示结果如下图，发现了连经纬度标识都没有，虽然可以根据数据变化方向识别，不过增加上经纬度标识还是要好一些。然后古哥浏览器找了一下看看有没有解决方法，结果没找到，看到很多教程都没有添加这个经纬度标识，于是自己动手：

如何使用Pandas,GeoPandas和ArcGIS实现数据处理和绘图

自己动手添加经纬度标识：

plt.figure(figsize=(16, 10))ax = plt.subplot(projection=ccrs.PlateCarree())ax.set_yticks(np.arange(20, 26, 1.0),crs=ccrs.PlateCarree())ax.set_xticks(np.arange(109, 118, 1.0),crs=ccrs.PlateCarree())
ax.xaxis.set_major_formatter(LongitudeFormatter()) #经纬度标识显示，EWSNax.yaxis.set_major_formatter(LatitudeFormatter())

现在来看看效果吧，啊哈，这样的话，经纬度标识不就来了嘛

如何使用Pandas,GeoPandas和ArcGIS实现数据处理和绘图