geopandas

阿里云双11来了!从本博客参与阿里云,服务器最低只要86元/年!

import geopandasas gpd

import pandasas pd

import matplotlib.pyplotas plt

import numpyas np

from shapely.geometryimport Point

from geopandasimport geodataframe

# df=gpd.read_file(‘./data-v1_2/countries.shp’)

# df.plot()

# plt.show()

df = pd.DataFrame()

size =30

df[‘id’] =range(0,size)

df[‘value1’] = np.random.randint(0, 100, size=size)

df[‘lon’] = np.random.randint(103000000, 105000000, size=size)/1000000

df[‘lat’] = np.random.randint(30000000, 35000000, size=size) /1000000

df[‘geom’] = [Point(xy)for xyin zip(df.lon, df.lat)]

crs = {‘init’:’epsg:4326′}

geodf = gpd.GeoDataFrame(df, crs=crs, geometry=’geom’)

# print(geodf.head(1))

geodf.plot()

plt.show()

https://www.jianshu.com/p/dbc2945af086

Python量化投资网携手4326手游为资深游戏玩家推荐:《梦想协奏曲!少女乐团派对!下载

「点点赞赏,手留余香」

    还没有人赞赏,快来当第一个赞赏的人吧!
pandas
0 条回复 A 作者 M 管理员
    所有的伟大,都源于一个勇敢的开始!
欢迎您,新朋友,感谢参与互动!欢迎您 {{author}},您在本站有{{commentsCount}}条评论