R中是否有类似Rasters库的Python库

编程问答 2022-08-12

问题描述

我正在尝试在python中读取和解析.grc / .grd文件，并将某些数据转换为数据帧中的列

有关这些文件类型的一些常规信息： 由MapInfo，Vertical Mapper和各种电信无线规划工具生成的.GRD（网格）和.GRC（分类网格）栅格文件。 这些文件通常用于表示给定蜂窝塔位置的信号强度预测分析的输出。 这些栅格通常会转换为矢量格式，然后在其上执行各种类型的分析

我正在尝试使用gdal库并解析文件，但未找到任何将数据转换为数据框的API python中是否还有其他库可以解析.grc / .grd文件

我的一位同事在R中编写了实现此目的的代码

这是R中的代码段

parse_grc <- function(file_name){
 
tab_grc <- raster(paste0(file_name,".grc"))
grc_data <- as.data.frame(rasterToPoints(tab_grc))
coordinates(grc_data) <- ~x+y
crs(grc_data) <- crs(tab_grc)
grc_data <- spTransform(grc_data,CRS("+proj=longlat +datum=WGS84"))
site_info <- as.data.frame(tab_grc@data@attributes)
 
grc_data %>%
  as.tibble() %>%
  rename(id = contains(".")) %>%
  inner_join(site_info,c("id" = "ID")) %>%
  select(-id) %>%
  rename(sector = category) %>%
  mutate(flname = file_name) %>%
  separate(flname,into = c("type","mkt","band"),sep = "-",remove = TRUE) -> grc_df
 
return(grc_df)
}

解决方法

暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！

如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。

小编邮箱:dio#foxmail.com (将#修改为@）

arcgispython-3.xrrr-raster