| Re: Read an ASCII file in a table based on delimeter line string [message #93457 is a reply to message #93455] |
Thu, 21 July 2016 11:17   |
Craig Markwardt
Messages: 1869
Registered: November 1996
|
Senior Member |
|
|
On Thursday, July 21, 2016 at 11:36:56 AM UTC-4, burkina wrote:
> Hi everyone,
>
> suppose I have a file like this one (there will be no regularity in the real ones!):
...
> I want to read in IDL in separated arrays:
>
> x[0]=[a,b,c]
> x[1]=[d,e,f,g]
> x[2]=[h]
>
> y[0]=[1,2,3]
> y[1]=[4,5,6,7]
> y[2]=[8]
>
> or a table t[2,8], or something equivalent.
> In practice, I would like to use something like readcol but with IDL understanding the reading should stop every time it encounters a delimiter (#### in this case), and then start again separating the arrays for each set.
It's not going to be really easy, or efficient. It's not hard for someone who has programming experience, though.
Step one is deciding on what kind of data structure you want. The IDL arrays you selected do not support varying dimensions so you can't use that type of data structure. IDL lists might be a better approach, if you want to keep similar semantics.
Step two is deciding how to parse the data from an open file unit. Think of how you would do that as a human. The following script may help. I didn't test it in any way, so who knows if it will work.
openr, unit, 'myfile.dat', /get_lun
output = [] ;; Initialize output
while eof(unit) EQ 0 do begin
;; Initialize X and Y to empty arrays
x = []
y = []
s = '' ;; Initialize dummy string
readf, unit, s ;; Read one line from input file
;; s is now one line from the input file
;; Begin parsing
if s EQ '###' then begin ;; It is a delimiter
output = [output, {x: x, y: y}
endif else begin ;; Or else parse it into data
x1 = '' & y1 = 0L ;; Initialize dummy variables
;; NOTE: this will change based on the actual format of your data file
reads, s, x1, y1, format='(A1,I)'
;; Append to existing list
x = [x, x1]
y = [y, y1]
endelse
endwhile
free_lun, unit
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
