transbigdata笔记:数据栅格化

发布时间:2024年01月11日

1?area_to_grid

在边界或形状中生成矩形栅格

1.1 主要使用方法

transbigdata.area_to_grid(
    location,?
    accuracy=500,?
    method='rect',?
    params='auto')

1.2 主要参数

location

(bounds(List) or?shape(GeoDataFrame)

生成栅格的位置。

  • 如果边界为 [lon1, lat1, lon2, lat2](WGS84),其中 lon1 , lat1 是左下角坐标,lon2 , lat2 是右上角坐标
  • 如果是形状,则应为 GeoDataFrame
accuracy栅格尺寸(米)
method?直角、三角或六角(?rect, tri or hexa)

1.3 返回内容

grid?(GeoDataFrameLONCOL 和 LATCOL 是栅格的索引,geometry是各栅格的中心
params栅格参数

1.3.1 举例(rect)

import transbigdata as tbd
bounds = [113.6, 22.4, 114.8, 22.9]
res=tbd.area_to_grid(bounds, accuracy=500)
res[0]

res[1]
'''
{'slon': 113.6,
 'slat': 22.4,
 'deltalon': 0.004872390756896538,
 'deltalat': 0.004496605206422906,
 'theta': 0,
 'method': 'rect',
 'gridsize': 500}
'''

1.3.2 举例(tri)?

import transbigdata as tbd
bounds = [113.6, 22.4, 114.8, 22.9]
res=tbd.area_to_grid(bounds, accuracy=5000,method='tri')
res[0]

?

res[1]
'''
{'slon': 113.6,
 'slat': 22.4,
 'deltalon': 0.048723907568965386,
 'deltalat': 0.044966052064229066,
 'theta': 0,
 'method': 'tri',
 'gridsize': 5000}
'''

?

res[0].plot(figsize=(15,8))

?

2 area_to_params

transbigdata.area_to_params(
    location, 
    accuracy=500, 
    method='rect')

和area_to_grids参数一样,用法上就是返回area_to_grids返回的第二个元素

bounds = [113.6, 22.4, 114.8, 22.9]
res=tbd.area_to_params(bounds, accuracy=500)
res
'''
{'slon': 113.6,
 'slat': 22.4,
 'deltalon': 0.004872390756896538,
 'deltalat': 0.004496605206422906,
 'theta': 0,
 'method': 'rect',
 'gridsize': 500}
'''

3?GPS_to_grid

将 GPS 数据与栅格匹配

transbigdata.GPS_to_grid(lon,?lat,?params)

3.1 主要参数

lon?经度栏
lat?纬度栏
params2 area_to_params返回的内容

3.2 返回内容

  • 矩形栅格

    • [LONCOL,LATCOL]?(list) – 两列 LONCOL 和 LATCOL 一起可以指定栅格。

  • 三角形和六边形栅格

    • [loncol_1,loncol_2,loncol_3]?(list) – 栅格纬度的索引。两列 LONCOL 和 LATCOL 一起可以指定一个栅格。

    • 【这里是不是loncol_3就是LATCOL啊?】

3.3 举例

比如data是这样的内容:

transbigdata 笔记:官方文档案例1(出租车GPS数据处理)-CSDN博客?中清楚异常出租车记录后

data=tbd.clean_taxi_status(data)
data

3.3.1 栅格method为rect

params1=tbd.area_to_params([113.6, 22.4, 114.8, 22.9],accuracy=500)
params1
'''
{'slon': 113.6,
 'slat': 22.4,
 'deltalon': 0.004872390756896538,
 'deltalat': 0.004496605206422906,
 'theta': 0,
 'method': 'rect',
 'gridsize': 500}
'''

tbd.GPS_to_grid(data['Lng'],data['Lat'],params)
'''
[array([51, 51, 51, ..., 77, 77, 77]), array([65, 66, 66, ..., 32, 32, 32])]
'''

4?grid_to_polygon

文章来源:https://blog.csdn.net/qq_40206371/article/details/135507342
本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。