|
首页 >> 地学 >> 查看话题

GIS、python、ENVI、IDL中分别实现点值提取总结

kainan001 举报此信息
如题,最近在用这个功能,总结如下:
1、GIS中实现点值提取: http://www.sunzx.net/archive/1883.html
2、python for gis中实现: http://resources.arcgis.com/zh-cn/help/main/10.1/index.html#//009z0000002t000000
3、ENVI中实现: http://blog.sina.com.cn/s/blog_764b1e9d0101dqte.html
4、IDL+ENVI中实现,参考《遥感二次开发语言IDL》,徐永明著,代码如下:
pro read_index_layer
fn=dialog_pickfile(title='选择采样点文件',get_path=work_dir)
cd,work_dir
data=read_csv(fn,count=nsta,header=header)
Sta_name=data.(0)
Lat=data.(1)
Lon=data.(2)
fn=dialog_pickfile(title='选择图层文件')
envi_open_file,fn,r_fid=fid
map_info=envi_get_map_info(fid=fid)
i_proj=envi_proj_create(/geographic,datum='D_WGS_1984')
o_proj=map_info.proj
envi_convert_projection_coordinates,lon,lat,i_proj,xmap,ymap,o_proj
envi_convert_file_coordinates,fid,xf,yf,xmap,ymap
xf=floor(xf) & yf=floor(yf)
altitude=fltarr(nsta)
for i=0,nsta-1 do begin
  dims=
  altitude =envi_get_data(fid=fid,dims=dims,pos=0)
endfor
o_fn=dialog_pickfile(title='结果保存为')+'.csv'
header=
data=create_struct(data,'Field4',altitude)
write_csv,o_fn,data,header=header
end