IDL基础

时间：2022-10-14 21:48:03
先列后行

arr=indgen(3,4)

SIZE(arr,/dimensions)

 print ,where(arr gt7)

print,arr[where(arr gt 7)]

print,where(arr gt 7,cout,COMPLEMENT=com,ncomplement=ncom)

array[WHERE(array GT 5, /NULL)] = 5

arr=indgen(2,8);

arrReform=reform(arr,1,8);维数一定,改变个数

arrRebin=REBIN(arr,2,8);维数一定，改变每一维的倍数

 arrCongrid=congrid(arr,1,1)

arrRotate5=ROTATE(ARR,5)

arr=indgen(8)

print,arr[sort(arr)];sort()返回排序后的索引

print,uniq(arr,sort(arr));返回arr中的唯一值，即，不重复，如果用uniq(arr)则返回的是将相邻的相同值只取一个，不相邻的却去不了

print,arr1 # arr2;数组相乘

print,arr1 ## arr2;矩阵相乘

print,strmid(str,strlen(str)+1,strlen(str),/reverse_offset);反向截取字符串

file=dialog_pickfile()

print,file

hdr=STRMID(file,pointPosition+1,strlen(file)-pointPosition-1)

help,strtrim(string(a),1);转换为string并trim掉一边的空格

结构体

 str={v1:5,v2:'dfdf'};定义匿名结构体

print,n_tags(str);结构体个成员个数

print,tag_names(str);结构体成员 名称

print,str.v1;某一成员的值

str3={two,inherits one,v4:'1'};two结构体继承自one,并向two添加v4变量

动态扩展

str4=create_struct(str2,'v5',findgen(5));创建结构体，并将str2中的数据也copy过来

指针

data=indgen(5);定义数组

ptr=ptr_new(data);定义指针------ptr_new()

Data=!null;制空

print,*ptr;打印指针,-------*ptr

ptr=ptr_new(data,/no_copy);创建指针后将data制空（undefined）

print,*pt=5;动态赋值还可以动态改变

print,*pt=indgen(5)

链表

obj_destroy,list;销毁list

Obj_valid(list);结果如果为0则说明销毁对象list

逻辑运算符

&&与

||或

~非

位运算符

AND

NOT

OR

XOR异或

if (keyword_set(xxx)) then begin;keyword_set()方法用于判断是否输入xxx

关键字继承

Plotex.por文件

pro plotex,x,y,_extra=_extra

plot,x,y,_extra=_extra			;_strict_extra=_extra严格关键字继承

end

然后在控制台中键入

IDL> x=findgen(100)/10

IDL> plotex,x,sin(x)

IDL> plotex,x,sin(x),thick=2,color=255

效果如下

pro plotex,x,y,_ref_extra=_extra;_ref_extra定义之后plotex只接受下面_extra=[]中指定的变量

plot,x,y,_extra=['thick','color']

End

pro plotex,x,y,_strict_extra=_extra

plot,x,y,_extra=_extra	        ;_strict_extra=_extra严格关键字继承

end

算法优化

a=!null;结果：A               UNDEFINED = !NULL

delvar,a;结果：A               UNDEFINED = <Undefined>

IDL> a=4

IDL> b=5

IDL> c=temporary(a)+temporary(b)

IDL> help,c,a,b

C               INT       =        9

A               UNDEFINED = <Undefined>

B               UNDEFINED = <Undefined>

IDL> 

时间控制

TIC:程序开始

TOC:程序结束

EX:

pro plotex,x,y,_ref_extra=_extra

tic

plot,x,y,_extra=['thick','color']

toc

End

结果：% Time elapsed: 0.0069999695 seconds.
IDL基础

IDL> z=shift(dist(40),20,20)

% Compiled module: DIST.

IDL> z=exp(-(z/10)^2)

IDL> surface,z

IDL>

输出为：
IDL基础

数据输入输出：

print，read，reads,string()

IDL> str='123 45 idl'

IDL> a=0

IDL> b=''

IDL> reads,str,a,b

IDL> help,a,b

A               INT       =      123

B               STRING    = ' 45 idl'

Format---输出格式控制

[n]FC[+][-][width]

Opner	只读操作

Openw	写文件

Openu

File_search()	搜索特定的文件

Dialog_pickfile()

Fstat()

Eof()	判断是否到最后

Close	关闭逻辑设备号

Free_Lun	释放逻辑设备号

;按行读取数据

pro readAndWriteFile

file=dialog_pickfile();打开文件

openr,lun,file,/get_lun;读取文件

temp='';定义中间变量用来存储读取的数据

openw,lunOpen,'F:\envi\Temp\openw1.txt',/get_lun

while(~eof(lun))do begin

  readf,lun,temp

  printf,lunOpen,temp

Endwhile

free_lun,lun

free_lun,lunOpen

end

;分块读取数据

file=dialog_pickfile()

openr,lun,file,/get_lun

data1=strarr(20)

data2=findgen(6,30)

readf,lun,data1

readf,lun,data2

print,data1

print,data2

;向导式读取数据

pro ascii_templatemethond

file=dialog_pickfile()

if ~file_test(file) then return

template=ascii_template(file);存为二进制文件

if size(template,/type) eq 2 then return;判断读进来的数据是否为整形

data=read_ascii(file,template=template);读取二进制文件

p=plot(data.(1));用读到的数据的第二列制图

help,data

print,data.(0)

End

;规则数据读取

pro array

file=dialog_pickfile()

openr,lun,file,/get_lun

readf,lun,firstLine

colum=0

row=0

type=0

reads,firstLine,colum,row,type

arr=make_array(colum,row,type=type)

readf,lun,temp1

readf,lun,temp2

;忽略两行

readf,lun,temp3

print,temp3

end

二进制文件读取

Readu

writeu

风云数据读取

file=dialog_pickfile()

openr,lun,file,/get_lun

headline=intarr(3)

 point_lun,lun,20;跳过前面20个数据

readu,lun,headline

fyarr=bytarr(1201,1201)

point_lun,lun,2*1201;跳过两行

tvscl,fyarr

读图片

file=dialog_pickfile()

read_jpeg,file,data

tvscl,data,/true;true为1，从help data-->DATA            BYTE      = Array[3, 600, 400],3在第一个位置，

 write_jpeg,"桌面\1.jpg",tvrd(true=1),quality=75,true=1

随机推荐

  1. 项目中libevent几个问题

    几个问题: .libevent到底用的是select还是iocp,然后是如何突破64限制的 typedef struct fd_set { u_int fd_count; /* how many ar ...

  2. DirectX截图黑屏的解决办法

    好久没有更新博客了,今天开始继续耕耘. 生活要继续 工作要继续 梦想也一定要继续! 之前写过一篇关于DirectX截屏的文章,其中有网友留言提到了截图黑屏的问题,于是这些日子研究了一下,与大家一同分享 ...

  3. hdu2044java递推

    一只小蜜蜂... Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 65536/32768 K (Java/Others)Total Su ...

  4. redis在spring中的配置及java代码实现

    1.建一个redis.properties属性文件 # Redis Setting redis.addr = 127.0.0.1 redis.port = 6379 redis.auth = mast ...

  5. pycharm&plus;QT4的helloworld

    # -*- coding: utf-8 -*- from PyQt4 import QtCore, QtGui try: _fromUtf8 = QtCore.QString.fromUtf8 exc ...

  6. 350 - Pseudo-Random Numbers

     Pseudo-Random Numbers  Computers normally cannot generate really random numbers, but frequently are ...

  7. loadrunner中如何将MD5加密的值转换为大写

    上篇博客中写过如何将MD5加密,但是我们在实际的测试过程中可能需要将加密的结果进行大小写转换.我在这次的测试过程中就遇见了这样的问题, 我在测试时发现开发人员代码传的sign值是大写,而我加密出来的s ...

  8. sublime使用技巧之集成VI

    熟悉开发工具,减少多余的操作流程有助于提高开发效率,而Sublime Text 2是sublime产品的经典版本,因此本文基于Sublime Text 2讲解sublime的使用技巧. VI的主要作用 ...

  9. Chapter4 复杂度分析&lpar;下&rpar;&colon;浅析最好&comma;最坏&comma;平均&comma;均摊时间复杂度

    四个复杂度分析: 1:最好情况时间复杂度(best case time complexity) 2:最坏情况时间复杂度(worst case time complexity) 3:平均情况时间复杂度( ...

  10. excel导出的时候从程序后台写到excel里的是文本，所以无法在excel中计算怎么办？

    文章引用自:http://www.cnblogs.com/rayray/p/3414452.html excel导出的时候从程序后台写到excel里的是文本,所以无法在excel中计算怎么办? 需要导 ...

相关文章