; ; input the wedges ; data will be in dar[9] {} ; baseD='/share/megs2_u1/phil/data/leica_p50/200311/' scanlist=file_search(baseD + "Scan*.las") a=stregex(scanlist,".*Scan([0-9]*).las",/extr,/sub) scanNum=reform(long(a[1,*])) ii=where((scanNum ge 10) and (scanNum le 18),nscans) ; a={scanNum:0l,$ npnts :0l,$ pd :ptr_new()$ ; points to d[npnts] } dar=replicate(a,nscans) for i=0,nscans-1 do ptr_free,dar[i].pd for i=0,nscans-1 do begin &$ j=ii[i] &$ print,"start "+ scanList[j] &$ npnts=blminpfile(scanList[j],hdr,d) &$ if i eq 0 then begin &$ a={scanNum:0l,$ npnts :0l,$ hdr :hdr,$ pd :ptr_new() } &$ dar=replicate(a,nscans) &$ endif &$ dar[i].scanNum=scanNum[j] &$ dar[i].npnts=npnts &$ dar[i].hdr =hdr &$ dar[i].pd =ptr_new(d,/no_copy) &$ endfor end