IDL编程学习之筛选数据对的心法

15195775117

在工作和研究中经常遇到这种情况：
现在有3列数据：经度-纬度-测值，要求在某经纬度范围内统计测值的规律。
这时，就需要过滤下数据。我们使用where函数获得符合要求的经度列的索引，和纬度列的索引，这两个索引的共同元素就是符合经纬度要求的行。
案例：附近txt文件有2列数据，要求筛选出第一列大于0，第二列小于0的行。
代码：
a=read_ascii('C:\Users\Administrator\Desktop\x.txt',data_start=0)
x=a.field1
p1=where(x[0,*] gt 0,count1)
p2=where(x[1,*] lt 0,count2)
foreach element,p1 do begin
  p=where(p2 eq element,count)
  if(count eq 1)then begin
    print,x[*,element]
    endif
endforeach

ybmy001

你写复杂了，害我下载你这个txt文件，浪费我5个金币啊。


  file = 'D:\x.txt'
  
  data = fltarr(2, file_lines(file))
  
  openr, lun, file, /get_lun
    readf, lun, data
  free_lun, lun
  
  lines = where(data[0, *] GT 0 and data[1, *] LT 0)
  
  PRINT, lines

ybmy001

还有，IDL这语言好奇葩，用了5年了。我现在Python和IDL都在用。在逐步转向Python

15195775117

ybmy001 发表于 2018-2-13 18:02
还有，IDL这语言好奇葩，用了5年了。我现在Python和IDL都在用。在逐步转向Python

毕竟IDL资源少