如何读取R中的Mapinfo文件

时间:2013-11-20 16:14:34

标签: r mapinfo

法国国家研究所(Insee)以MapInfo格式提供地理数据(两个文件.mid和.mif以及一个dbf文件)。如何在R中读取这些文件?

这是一个example

2 个答案:

答案 0 :(得分:8)

MapInfo文件有一个OGR驱动程序(包rgdal):

R> library("rgdal")
R> ogrDrivers()[28, ]
           name write
28 MapInfo File  TRUE

但是文件/几何体存在问题,readOGR会显示错误消息:

R> ogrListLayers("R02_rfl09_UTM20N1000.mid")
[1] "R02_rfl09_UTM20N1000"

R> readOGR("R02_rfl09_UTM20N1000.mid", layer="R02_rfl09_UTM20N1000")
OGR data source with driver: MapInfo File 
Source: "R02_rfl09_UTM20N1000.mid", layer: "R02_rfl09_UTM20N1000"
with 967 features and 4 fields
Feature type: wkbPolygon with 2 dimensions
Error in stopifnot(is.list(srl)) : ring not closed

但是,我能够使用GRASS GIS读取文件,这些文件可以从R(包spgrass6)编写脚本:

v.in.ogr dsn=R02_rfl09_UTM20N1000.mid output=R02_rfl09_UTM20N1000 snap=1e-08

GRASS screenshot

答案 1 :(得分:5)

这有点难以说,因为你的pdf只定义了.mid结构。

这取决于你想对日期做什么,但看看它,.mid文件有每个区域的SW合作,并检查.mif文件,每个区域是1000m2,所以你可以只计算区域(对于此数据样本)而不是加载它们。

所以这是加载它的一种方法,但它取决于你想用数据做什么

首先将.csv文件复制到您的工作目录中,然后

coords<-read.csv(file="R02_rfl09_UTM20N1000.mid", header=FALSE)
colnames(coords)<-c("SW.E","SW.N","ind","indXYNE1")
# add the co-ords for the area
coords$SE.N=coords$SW.N
coords$SE.E=coords$SW.E+1000
coords$NW.N=coords$SW.N+1000
coords$NW.E=coords$SW.E
coords$NE.N=coords$SW.N+1000
coords$NE.E=coords$SW.E+1000

head(coords)

这会给你:

    SW.E    SW.N ind indXYNE1    SE.N   SE.E    NW.N   NW.E    NE.N   NE.E
1 690000 1636000 241        6 1636000 691000 1637000 690000 1637000 691000
2 690000 1637000 414        3 1637000 691000 1638000 690000 1638000 691000
3 690000 1638000 240        6 1638000 691000 1639000 690000 1639000 691000
4 690000 1640000   8        0 1640000 691000 1641000 690000 1641000 691000
5 691000 1634000 142        0 1634000 692000 1635000 691000 1635000 692000
6 691000 1635000 216        5 1635000 692000 1636000 691000 1636000 692000
....

这是每个区域的四个边界点,加上ind和indXYNE1,我猜你正在寻找什么?然后,您可以使用SW点(或新的派生密钥)转换数据作为每个区域的参考。

希望有所帮助!取决于你想要对数据做什么。