如何给数据附加坐标信息

指甲钳

        我的数据是存数在TXT文档上的，就是文本格式的。然后通过lines = asciiread(fpath,-1,"string")
delim = ";"
   lon  = tofloat(str_get_field(lines,2,delim))
  lat  = tofloat(str_get_field(lines,3,delim))
  data =tofloat(str_get_field(lines,4,delim))




这样的方式读取，读取出来后就是三个一维数组了，

lon(),lat(),data()

我现在想坐标信息，就是经纬度信息附加给数据
就是变成 data(lon,lat) 这样的格式？
请问要怎么做啊？

woyaofeixiang

data!0      = "lon"
date&lon  =lon
lon@units = "degrees_east"
data!1     = "lat"
date&lat  =lat
lat@units = "degrees_north"

LZ可以试试上面你的命令

你可以试试这个

longlivehj

直接按SF的做法不行。
首先要对一维的data用reshape整成二维，相应的lon和lat要去重。然后再按SF说的加坐标变量。

指甲钳

longlivehj 发表于 2014-4-3 12:27
直接按SF的做法不行。
首先要对一维的data用reshape整成二维，相应的lon和lat要去重。然后再按SF说的加坐 ...

谢谢！！这位大大，能说的具体些吗？怎么reshap，怎么去重啊？

指甲钳

woyaofeixiang 发表于 2014-4-3 11:25
data!0      = "lon"
date&lon  =lon
lon@units = "degrees_east"

谢谢答复，我尝试过，不可行啊，应为data本来是一个一维数组

longlivehj

指甲钳 发表于 2014-4-3 14:11
谢谢！！这位大大，能说的具体些吗？怎么reshap，怎么去重啊？

把数据发一个上来吧！

指甲钳

longlivehj 发表于 2014-4-3 14:17
把数据发一个上来吧！

G6825;112.8738889;23.05722222;1931.2
G6832;113.0583333;23.13611111;2020.6
G2271;112.9477778;23.06194444;1553.5
G2263;113.0163889;22.92277778;2017.4
G6840;113.0022222;23.21777778;1834.1
G6837;113.1833333;23.12361111;2006.4
G6843;112.9544444;23.20805556;1925.2
G6826;112.9744444;23.03694444;2086.5
G2272;113.1666667;23.07861111;1925.2
G6819;112.9952778;22.87222222;1937.7
G6868;112.9666667;22.93611111;2069.1
G6833;113.0663889;23.10083333;2012.6
G2270;113.0483333;23.06194444;1851.2
G2274;113.1866667;23.01861111;2050.7
G6847;113.133333;23.05;1912
G6842;112.9361111;23.15;1759.5
G6827;112.9244444;22.99694444;1739.9
G2265;113.01;22.83305556;1979.2
G2261;113.1144444;23.11527778;1926.8
G6829;113.095;23.275;2033.3
G6816;113.2494444;23.03694444;2154
G6830;113.125;23.2;2201.6
G2213;113.0086;23.1447;2149.2
G6822;112.8869444;22.93083333;1978.3
G6835;113.1194444;23.075;1827.4
G6841;113.1083333;23.16388889;1857.4
G6818;112.9372222;22.85305556;1610.5
G2268;113.0669444;23.23333333;1918.3
G6817;113.2294444;22.99777778;2067.1
G2276;113.0669444;23.18666667;1916.5
G6828;112.9022222;23.11777778;1931.6
G6821;112.9822222;22.97388889;1933.4
G2206;113.1569444;23.21555556;1752.7
G6844;113.0197222;23.07611111;1911.6
G6839;112.9777778;23.175;1816.7


第二列是经度，三列是纬度，四列是数据

longlivehj

指甲钳 发表于 2014-4-3 14:39
G6825;112.8738889;23.05722222;1931.2
G6832;113.0583333;23.13611111;2020.6
G2271;112.9477778;23.0 ...

这个不完整吧！能发个完整的不？哪怕下载要钱。
因为，我看不出经纬度的规律。如果没有规律的化，还不一定能按你的要求转成二维。

指甲钳

longlivehj 发表于 2014-4-3 14:49
这个不完整吧！能发个完整的不？哪怕下载要钱。
因为，我看不出经纬度的规律。如果没有规律的化，还不一 ...

没有规律的啊~~~那怎么办啊？因为我想对这些数据在NCL里面进行异常值剔除，而且如果转换成二维数组做插值也容易很多~~

指甲钳

longlivehj 发表于 2014-4-3 14:49
这个不完整吧！能发个完整的不？哪怕下载要钱。
因为，我看不出经纬度的规律。如果没有规律的化，还不一 ...

就是都在一个范围112~114,22~24的经纬度范围