Update girds

update cells set grid_id = cast(floor((y(geom)-2575200)/5)*25160+floor((x(geom)-12572400)/5)+1 as integer);

update grid5 set grid_id = cast(floor((y(center)-2575200)/5)*25160+floor((x(center)-12572400)/5)+1 as integer);

gleft = 112.940000000 # 原点x

gbottom = 22.53000000 # 原点y

gright = 114.07000000 # 右边界

gtop = 23.95000000000 # 上边界

step = 12.50000000000 # 网格中心点步距(单位：米)

lng2km = 102.70000000 # 从经度换算到公里

lat2km = 111.70000000 # 从纬度换算到公里

lngstep = step/1000/lng2km # 步距换算成经度 //0.000121713729308666

latstep = step/1000/lat2km # 步距换算成纬度 //0.000111906893464637


ngalaxylngs=((gright-gleft)/lngstep).floor # 9284

ngalaxylats=((gtop-gbottom)/latstep).floor # 12689