我在美国大陆的km网格上有一组数据值,这些列是“纬度”,“经度”和“观测”,例如:
"lat" "lon" "yield"
25.567 -120.347 3.6
25.832 -120.400 2.6
26.097 -120.454 3.4
26.363 -120.508 3.1
26.630 -120.562 4.4
或者,作为R数据帧:
mydata <- structure(list(lat = c(25.567, 25.832, 26.097, 26.363, 26.63),
lon = c(-120.347, -120.4, -120.454, -120.508, -120.562),
yield = c(3.6, 2.6, 3.4, 3.1, 4.4)), .Names = c("lat",
"lon", "yield"), class = "data.frame", row.names = c(NA, -5L))
(完整的数据集可以在此处以csv的形式下载)
数据从30 km x 30 km的作物模型(预期在)上输出(来自Miguez等,2012)。
如何使用GIS相关的元数据(例如地图投影)将它们转换为栅格文件?
理想情况下,该文件将是文本(ASCII?)文件,因为我希望它与平台和软件无关。
作为CSV,这已经是 ASCII中的“文本文件”。另外,由于它根本不使用任何投影,因此可能几乎没有要添加的相关元数据(主要是基准)。您能否更具体地说明您要寻求什么样的输出以及打算如何使用它?
—
ub
我想通过各种地图绘制软件(ArcGIS,Google Maps,Grass,R等)使某人尽可能轻松地使用数据,从而方便重用,例如不需要其他转换步骤。根据GIS文件格式的Wikipedia页面,我推断1)“光栅”文件应具有带有纬度的行名和经度的列名(如图片),以及2)元数据应包括地理信息(拐角位置,覆盖区域)通过数据)。
—
安倍
@Nandini不知道正确的proj4string是什么,我怀疑lambert保形:
—
2014年
proj +proj=lcc +lat_1=50.0 +lat_2=50.0 +units=km +lon_0=-145.5 +lat_0=1.0
。对于另一个csv文件,我不确定您要的是什么-与问题中链接到的csv文件有什么不同,或者由接受的答案产生?