| number formats inputing ascii data [message #43826] |
Fri, 29 April 2005 07:26 |
dave
Messages: 31
Registered: February 1994
|
Member |
|
|
A quick question that I guess follows on to Ed's precision question.
I'm having similar problems when I import an ascii datafile of X,Y,Z
values.
the x and y values are utm meter and the z value is mm's
When I import the data it definitely looks different than it should?
I'm using the transread program from Craig M. to import the data. I've
attached a snipet of the data and the code below...any ideas?
Dave
----data file
LAT LON ELEV
37626408 5352625 668694
37626406 5352080 671058
37626406 5352015 674405
37626405 5351812 682801
37626409 5352791 668146
37626405 5351757 686651
37626405 5351655 694568
37626407 5352339 669945
37626405 5351701 690698
37626405 5351914 677368
37626411 5353354 666662
37626409 5352874 667493
37626410 5353031 666536
37626404 5351470 702865
37626408 5352625 668677
37626407 5352209 670414
37626411 5353456 668980
37626406 5352034 673583
37626412 5353788 666713
37626409 5352920 667493
37626404 5351554 698988
37626409 5352837 667973
37626404 5351406 709303
37626404 5351424 707796
37626404 5351517 700894
37626410 5353003 666463
37626409 5352883 667347
37626409 5352874 667493
37626410 5353003 666474
37626409 5352948 667574
37626408 5352717 668299
37626405 5351868 679227
37626407 5352394 670108
37626404 5351646 695965
37626408 5352569 669591
37626406 5351988 675233
37626407 5352339 670110
37626406 5352172 670095
37626408 5352597 668901
37626407 5352422 670716
37626408 5352560 668695
37626405 5351674 690904
---end data file
pro plane_fit, x, y, z, file=filnam
IF NOT Keyword_set(filnam) THEN BEGIN
filnam=dialog_PICKFILE( PATH=path,TITLE='Select xyz Data for',
/MULTIPLE_FILES, /MUST_EXIST)
ENDIF
print, ''
print,'Filename: ', filnam
print, '
A = 0L & B = 0L & C =0L
FOR i=0, N_ELEMENTS(FILNAM)- 1 DO BEGIN
Aa = 0L & Bb = 0L & Cc =0L
OPENR, lun, filnam[i], /GET_LUN
TRANSREAD, lun, Aa, Bb, Cc,FORMAT='D0.0, D0.0, D0.0', $
FILENAME=filnam[i]
FREE_LUN, lun
xA=[A,Aa]
yB=[B,Bb]
zC=[C,Cc]
ENDFOR
print, 'number of data lines: ', n_elements(xA)
;---Extracting the data from the variables
lon=xA[*]
lat=yB[*]
elev=zC[*]
;-- Creating a 3 column dbl array to house the lon, lat and elev data
combo=fltarr(3,N_ELEMENTS(lon))
;-- Passing the data into the array
combo[0,*] = lon
combo[1,*] = lat
combo[2,*] = elev
data=combo
print, 'first lines of data from original file: ', data(*, 0:2)
;--Display the original data
return
end;
Thanks in advance,
Dave
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
