【IDL求助】 风矢作图程序出错，求大神指点！

Crazy_Merlin

各位大神，我写了一个程序，但是在画风矢的时候出现图中（见图）的错误，不知道为啥~请大家帮忙看看~谢谢！在线等回复！麻烦啦~
为了检验，我把画风矢所需的信息都输出到屏幕上了，分别是U,V,X,Y。不知道为啥程序会报错呢？

下面给出我写的程序（是整个大程序中的一个小事件）：
Pro VectorProcedure,event
result=dialog_message("Draw Vector!",title='Message')
Close,10,20
OPENR,10,"number.txt"
OPENR,20,"AMVs.txt"
Readf,10,number
Print,number
;
WV1=fltarr(7,number)
U=fltarr(number)
V=fltarr(number)
X=fltarr(number)
Y=fltarr(number)
Lon=fltarr(number)
Lat=fltarr(number)
Readf,20,WV1
;print,size(WV1)
;;Draw Vectors
U=WV1(4,*)
V=WV1(5,*)
X=WV1(0,*)
Y=WV1(1,*)
Lon=WV1(2,*)
Lat=WV1(3,*)


for i=0,number-1 do begin
  Print,U(i),V(i),X(i),Y(i)
endfor
;
Graphic1=VECTOR(U,V,X,Y,$
  /BUFFER,/CURRENT,$
  VECTOR_STYLE=1,$
  COLOR="blue",$
  OVERPLOT=1)
;;
Close,10,20
End

程序所需的数据在附件中给出！
（说明一下：那个AMVs.txt中，从左到右分别是，X,Y,Longitude,Latitude,U,V,Pressure，number.txt是用来存储风矢的数目的）


————————————————————————————————————————————————上面的问题已经解决，具体的是重新建一个数组，对其再赋值。给出代码：
for i=0,number-1 do begin
  u1(i)=U(i)
  v1(i)=V(i)
  x1(i)=X(i)
  y1(i)=Y(i)
endfor
————————————————————————————————————————————————

{:soso_e154:}{:soso_e154:}{:soso_e154:}{:soso_e154:}
【求助】那位大神能帮我把这个问题解决了，或者指点一下~跪谢！！
另外，还有一个问题想要请教，我想把风矢直接叠加到云图中，而不是弹出一个窗口（有X,Y坐标轴）显示，该怎么写程序呢？
能不能有哪位高手帮我解决这个问题~~~~~谢谢啊！
云图是一个1900（宽）*1300（高）的等经纬度云图。上面的X,Y就是代表像素的位置。详见图2

（这是我写的那个大程序，主要是用来读取云图，并且将计算的风矢画到云图中对应的像素。第二个button对应的时间就是我上面给的程序。）


大体意思就是将风矢一个个画到云图中对应的位置~

Crazy_Merlin

大家帮帮忙~求回复！

大吉

庆哥我来了~

Crazy_Merlin

跪求高手啊！！！谢谢了~指点一下！

爱听雨的青蛙

你把画云图的代码贴出来看看啊，是风云awx 5km分辨率的吧。。。

Crazy_Merlin

爱听雨的青蛙 发表于 2013-5-14 09:41

                               
登录/注册后可看大图

你把画云图的代码贴出来看看啊，是风云awx 5km分辨率的吧。。。

height=1300
width=1900
;;1300hang,1900lie
WV0=PTR_NEW(intarr(width,height))
WV1=PTR_NEW(intarr(width,height))
WV2=PTR_NEW(intarr(width,height))
;
AWX=strarr(3)
AWX=DIALOG_PICKFILE(/READ,FILTER='*.AWX',$
  PATH="F:\IDLworkspace\clearskyAMVs\AWX\",$
  TITLE='Please Choose 3 IR1 Image',$
  /MULTIPLE_FILES)
AWXfile,AWX(0),*WV0,height,width,ilongitudeofwest,ilongitudeofeast,ilatitudeofnorth,ilatitudeofsouth,iyear0,imonth0,iday0,ihour0,iminute0
AWXfile,AWX(1),*WV1,height,width,ilongitudeofwest,ilongitudeofeast,ilatitudeofnorth,ilatitudeofsouth,iyear1,imonth1,iday1,ihour1,iminute1
AWXfile,AWX(2),*WV2,height,width,ilongitudeofwest,ilongitudeofeast,ilatitudeofnorth,ilatitudeofsouth,iyear2,imonth2,iday2,ihour2,iminute2
;为了解决调用问题，将文件名存入到一个文件中，方便CalProcedure调用
Close,10
OPENW,10,"Temp.000"
printf,10,AWX(0:2)
Close,10
;
TV,*WV1