2013-10-02 85 views
2

我必须读取表示地图的csv文件。 我的CSV是这样(数据块):/将csv数据绘制为地图 - Python

"","x","y","sim1","sim2","sim3","sim4","sim5","sim6","sim7","sim8","sim9","sim10","sim11","sim12" 
"1",181180,333740,5.56588745117188,6.29487752914429,7.4835410118103,5.75873327255249,6.62183284759521,5.81478500366211,4.85671949386597,5.90418815612793,6.32611751556396,6.99649047851562,6.52076387405396,5.68944215774536 
"2",181140,333700,6.36264753341675,6.5217604637146,6.16843748092651,5.55328798294067,7.00429201126099,6.43625402450562,6.17744159698486,6.72836923599243,6.38574266433716,6.81451606750488,6.68060827255249,6.14339065551758 
"3",181180,333700,6.16541910171509,6.44704437255859,7.51744651794434,5.46270132064819,6.8890323638916,6.46842670440674,6.07698059082031,6.2140531539917,6.43774271011353,6.21923875808716,6.43355655670166,5.90692138671875 

其中第一列显示有多少对x,y坐标我有,X和Y一直经度/纬度硅锰一个属性的值用于模拟#N。

使用D3,我可以快速读取这个CSV和显示它就像一张地图,它代表:

d3.csv("https://dl.dropboxusercontent.com/u/11442023/Geostatistic/sim100.csv", function(error, data) { 
if (error) { //If error is not null, something went wrong. 
    console.log("Error loading sim100 data."); 
    console.log(error); //Log the error. 
} else {  //If no error, the file loaded correctly. Yay! 
    console.log("Data sim100 loaded successfully!"); 
    //Include other code to execute after successful file load here 
    simulationDataset = data; 
    scaleX = d3.scale.linear().domain([d3.min(simulationDataset, function(d) {return d.x;}), d3.max(simulationDataset, function(d) { return d.x; })]); 
    scaleY = d3.scale.linear().domain([d3.min(simulationDataset, function(d) {return d.y;}), d3.max(simulationDataset, function(d) { return d.y; })]); 
    } 
}); 

enter image description here

我已经保存在一个numpy的数组:

np.loadtxt("Data/sim.csv", delimiter=',', usecols=range(1,4), skiprows=1) 

在Python中,我不知道我该怎么做到这一点。此外,除了绘制模拟地图之外,我还要保存这个矩阵,以备将来需要做的工作。但是,因为x和y是经度/纬度,所以我没有像(0,0)(0,1)(0,2)这样的点,就像一个简单的矩阵。但我不知道是否有像上面提到的那样来扩展域。

我想过寻找最大/最小经度和纬度来定义我的矩阵大小,并且他们做了一些微积分来适合每一对经度/纬度到这个结构中。但是这对我来说似乎不是一个好方法,并且已经看到了我看到Python可以做的所有事情,我想这存在一个更好的解决方案。

我最想知道的是,是否有一个“魔术”函数将这些数据保存在一个有意义的矩阵中(对所有经度/纬度索引进行转换并将值保存在此矩阵的索引中),或者如果我必须做我认为的事情,找到最小/最大值,并与翻译工作。

在此先感谢。

+2

提示:通常,当我想弄清楚如何使用'matplotlib'绘制某些东西时,我会转到[gallery](http://matplotlib.org/gallery.html)并寻找类似的东西。 – DSM

+0

你能显示你的R代码吗? –

回答

-2

您应该能够相当容易地编写一些解析该CSV并将数据输入到数组中的数据。 Numpy是科学数据计算最常用的模块之一。从那里,你可以使用matplotlib或mayavi来提供数据的图形表示。

+0

我已经把它放在一个numpy数组中。但除了绘制这个(我目前无法做到这一点),我必须将其保存为一个矩阵。 – pceccon