24小时热门版块排行榜    

CyRhmU.jpeg
查看: 3405  |  回复: 5

kainan001

铁杆木虫 (正式写手)

[交流] GIS、python、ENVI、IDL中分别实现点值提取总结已有3人参与

如题,最近在用这个功能,总结如下:
1、GIS中实现点值提取:http://www.sunzx.net/archive/1883.html
2、python for gis中实现:http://resources.arcgis.com/zh-c ... /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=[-1,xf,xf,yf,yf]
  altitude=envi_get_data(fid=fid,dims=dims,pos=0)
endfor

o_fn=dialog_pickfile(title='结果保存为')+'.csv'
header=[header,'Altitude']
data=create_struct(data,'Field4',altitude)
write_csv,o_fn,data,header=header
end
回复此楼

» 收录本帖的淘帖专辑推荐

学习

» 本帖已获得的红花(最新10朵)

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dannyjoyride

木虫 (正式写手)

LZ真善于总结
Nothing is impossible!
2楼2016-12-21 13:34:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wzwvk

金虫 (初入文坛)

送红花一朵
非常棒!
心中有猛虎,细嗅蔷薇!
3楼2016-12-21 14:14:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kainan001

铁杆木虫 (正式写手)

引用回帖:
3楼: Originally posted by wzwvk at 2016-12-21 14:14:17
非常棒!

哈哈,谢谢小红花,共同学习!

发自小木虫Android客户端
4楼2016-12-21 14:22:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

kainan001

铁杆木虫 (正式写手)

引用回帖:
2楼: Originally posted by dannyjoyride at 2016-12-21 13:34:56
LZ真善于总结

共同学习,一起进步

发自小木虫Android客户端
5楼2016-12-21 14:22:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

阿毛虫

银虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
楼主,能不能共享一份arcgis pro 软件
学习吧
6楼2020-01-06 21:24:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 kainan001 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见