[求助帖]全球0.5*0.5栅格数据如何处理成国家-年份

menglhl

地理气象小白求助各位专业大神。

变量描述

变量描述

因为要研究国家GDP与其年平均温度的关系，需要用此0.5*0.5的栅格数据（其实并不太懂是什么意思）估算出全球每个国家的年平均温度。研究了各位大神的经验贴，现在已经由月平均温度得到每个格点的年平均温度，但是不知道下一步该如何用网格数据得出国家数据。
搜了一下各个国家的大致经纬度信息，有些都不在数据描述里的实际变量（lat,lon）范围内,实在不懂地理和气象的专业知识，跪求各位大神指点迷津，最好能浅显易懂，先感谢各位大神了

kongfeng0824

可以输出成GeoTiff后用ArcGIS按国家处理

menglhl

kongfeng0824 发表于 2016-8-22 12:29
可以输出成GeoTiff后用ArcGIS按国家处理

感谢大神回复！
可不可以具体说一下该用ArcGIS的什么功能呢？ 有没有合适的资料推荐学习下，完全没接触过这个软件

MeteoInfo

可以参考这个帖子：MeteoInfoLab脚本示例：计算不同区域平均值
http://bbs.06climate.com/forum.p ... 661&fromuid=106
(出处: 气象家园)

menglhl

MeteoInfo 发表于 2016-8-22 14:07
可以参考这个帖子：MeteoInfoLab脚本示例：计算不同区域平均值
http://bbs.06climate.com/forum.php?mod=v ...

感谢老师回复，根据您帖子里的代码做了改动

1. f = addfile_surfer('E:/Ymeantemp.nc')
   
2. tdata = f['var'][:,:]
   
3. un = shaperead('E:/TEMPERATURES/map/global.shp')
   
4. i = 0
   
5. for rpoly in un.getShapes():
   
6.     name = un.getCellValue('NAME', i)
   
7.     mdata = tdata.maskout(rpoly)
   
8.     tave = mdata.ave()
   
9.     tmin = mdata.min()
   
10.     tmax = mdata.max()
    
11.     print name + ', Ave: %.2f, Min: %.2f, Max: %.2f' %(tave, tmin, tmax)
    
12.     i += 1
    

复制代码

运行后提示错误，老师能不能帮忙看一下哪里的问题。非常感谢！！
File "", line 6
        name = un.getCellValue('NAME', i)

    ^
SyntaxError: expected an indented block

这里是我使用的地图信息
COLUMN                TYPE                        DESCRIPTION

Shape                Polygon                        Country/area border as polygon(s)
FIPS                String(2)                FIPS 10-4 Country Code
ISO2                String(2)                ISO 3166-1 Alpha-2 Country Code
ISO3                String(3)                ISO 3166-1 Alpha-3 Country Code
UN                Short Integer(3)        ISO 3166-1 Numeric-3 Country Code
NAME                String(50)                Name of country/area
AREA                Long Integer(7)                Land area, FAO Statistics (2002)
POP2005                Double(10,0)                 Population, World Polulation Prospects (2005)
REGION                Short Integer(3)         Macro geographical (continental region), UN Statistics
SUBREGION        Short Integer(3)        Geogrpahical sub-region, UN Statistics
LON                FLOAT (7,3)                Longitude
LAT                FLOAT (6,3)                Latitude

MeteoInfo

menglhl 发表于 2016-8-23 13:21
感谢老师回复，根据您帖子里的代码做了改动

运行后提示错误，老师能不能帮忙看一下哪里的问题。非常感 ...

脚本是Python语法，最重要的就是正确缩进（判断、循环中），错误提示显示第6行缩进不正确，注意tab键缩进和4个空格缩进是不一样的。如果是在MeteoInfoLab中写脚本，tab键会自动转为4个空格。建议统一用4个空格进行缩进。

menglhl

MeteoInfo 发表于 2016-8-22 14:07
可以参考这个帖子：MeteoInfoLab脚本示例：计算不同区域平均值
http://bbs.06climate.com/forum.php?mod=v ...

补充一下，用java版本运行脚本仍提示第六行代码有误

rrun script...
  File "<iostream>", line 6
    &#194;&#160; &#194;&#160; code = un.getCellValue('NAME', i)
    ^
SyntaxError: no viable alternative at character '&#160;'

试了很多不同的字符替换’NAME' 都是这样的错误提示，不知道这里应该用的是什么，求解答

menglhl

MeteoInfo 发表于 2016-8-23 14:56
脚本是Python语法，最重要的就是正确缩进（判断、循环中），错误提示显示第6行缩进不正确，注意tab键缩进 ...

谢谢您，我再确认一下。

menglhl

MeteoInfo 发表于 2016-8-23 14:56
脚本是Python语法，最重要的就是正确缩进（判断、循环中），错误提示显示第6行缩进不正确，注意tab键缩进 ...

王老师您好，缩进问题改正以后仍然没有得出结果，已经多次检验使用的格点数据（nc文件）没有错误，运行结果如下：
是因为格点数据文件的格式（.nc）的原因么？数据为三维 lon*lat*t。
感谢您的回复。

MeteoInfo

menglhl 发表于 2016-9-5 19:18
王老师您好，缩进问题改正以后仍然没有得出结果，已经多次检验使用的格点数据（nc文件）没有错误，运行结 ...

建议你仔细看看这里：http://www.meteothinker.com/docs ... /data_tutorial.html