原始数据如下:

将站点值插值到1km和10km的网格上,插值的数据格式要求如下:
插值网格图用gis的grid格式
使用的投影:

Alias:
Abbreviation:
Remarks:
Projection: Albers
Parameters:
  False_Easting: 0.000000
  False_Northing: 0.000000
  Central_Meridian: 105.000000
  Standard_Parallel_1: 25.000000
  Standard_Parallel_2: 47.000000
  Latitude_Of_Origin: 0.000000
Linear Unit: Meter (1.000000)
Geographic Coordinate System:
Name: GCS_Krasovsky_1940

处理成ArcGIS可以处理的格式

-999.99 是无效值点

cellsize如何计算呢?

原始数据有两个点经纬度坐标

(21.1437,97.5284)

(29.2437,106.1864)

可以这样认为分别对应左下角点和右上角点,处理数据时使用左下角点.

两点的经度之差除以行数,或者纬度之差除以列数得到的值近似为0.008992

---------

利用ArcToolbox工具

得到没有投影信息的栅格文件,用北京54 Proj一下

按照初始数据的要求自定义投影

北京54就是GCS_Krasovsky_1940

投影后的栅格空间参考信息

利用ArcToolbox 的Resample工具重采样

将图中蓝框里的out cell size修改为1000,实际上数据就是1K+-的由于前面的cellsize设置不准确导致此处有较大误差,修改为10K将重采样为10K的栅格数据


posted on 2010-01-02 22:58  秋波渡  阅读(918)  评论(0编辑  收藏  举报