df=pd.read_csv(r’C:\Users\admin\Desktop\6-8month_processed_data\shanchu.csv’)
df.columns=[‘n’,‘time’,‘lat’,‘lon’]
density=pd.DataFrame(columns=(‘num’,‘latitude’,‘longitude’))
lamax=df.lat.max()
lomax=df.lon.max()
lamin=df.lat.min()
lomin=df.lon.min()
m=int((lamax-lamin)/0.01)
n=int((lomax-lomin)/0.01)
print(m,n)
lati=df.lat
longi=df.lon
y=0
fre=np.zeros((m,n),dtype=int)
for i in range(0,m):
for j in range(0,n):
for k in range(0,len(lati)):
if (lati[k]>=(lamin+0.01i)) & (lati[k]<(lamin+0.01(i+1))) & (longi[k]>=(lomin+0.01j)) & (longi[k]<(lomin+0.01(j+1))):
fre[i,j]+=1
density.loc[y]=[fre[i,j],lamin+0.01i,lomin+0.01j]
y+=1