[Proj] lla file format

Richard Greenwood richard.greenwood at gmail.com
Tue Feb 26 23:23:49 EST 2008


That's a bit more than I was looking for. I'm just trying to decypher
the .lla file. I assume that an .lla file is a combination of .las and
.los files (that the lat and long shifts are in a single file) and
that the shifts are in angular units. But that's about as far as I
have gotten.

Thanks for the reply,
Rich


On Tue, Feb 26, 2008 at 12:19 PM, Clifford J Mugnier <cjmce at lsu.edu> wrote:
>
>
>
> Is this what you need?
>
> Cliff Mugnier
> LOUISIANA STATE UNIVERSITY
>
> c       program convert.ftn  -- cpfs csm golden,co
> c
> c       this routine converts ascii files of x,y,z1,z2,...zn
> c       data to the center format for binary file storage
> c       the interactive nature allows for creation of the
> c       appropriate headers and file formats
> c
> c----------------------------------------------------------------
>       program convdam
>       parameter(nfile=50,nx=6000,ny=10)
>       character rident*56,pgm*8,filin(nfile)*30,
>      1          filout*30,formin*56
>       character gname(nfile)*8
>       dimension data(ny,nx),zero(ny)
>       common /inout/ iit,iot
>       common /input/ nv,nr,formin,rident,pgm,filout
> c      data iit,iot/5,6/,in,iout/30,31/,nc/-1/
> c      data pgm/'convertd'/
>       data x0,dx,y0,dy,angle/5*0.0/,zero/ny*0.0/
>  iit = 5
>  iot = 6
>  in = 30
>  iout = 31
>  nc = -1
>  pgm = 'convertd'
> c
> c......................................................................
> 1000  format(///,5x,'this program converts a formatted data file to',/,
>      &           5x,'an unformatted binary file conforming to the  ',/,
>      &           5x,'center format',//,' ** warning - the input files',
>      &              ' are assumed to have a uniform format;',
>      &            /,' all the groups are to have the same number',
>      &              ' and type of variables',//,
>      &           5x,'each group will come from one input file.',//,
>      &           5x,'program is dimensioned to read',i3,' files of',/,
>      &           8x,i6,' variables per sample',/,
>      &           8x,i6,' samples per group',//)
> 1005  format(' enter number of input-files/output-groups',
>      &       ' (eg. flightlines):')
> 1010  format(' enter filename for input file ',i3,' :')
> 1020  format(/,' enter filename for output file:')
> 1030  format(a30)
> 1040  format(///,5x,'enter header information,up to 56 characters',/,
>      &           5x,'this should include data identifier',/,
>      &           5x,'and coordinate system used !!')
> 1050  format(a56)
> 1070  format(/,5x,'enter the number of variables at each location',/,
>      &         5x,'this includes the x and y locations')
> 1080  format(/,5x,'enter the format of the input file records',/,
>      &         5x,'e.g. (3f20.4)    ')
> 1100  format(' enter an groupname (8 char) for input file ',a30,)
> 1110  format(a8)
> 1120  format('there will be',i3,' sample sets on each binary record')
> 1130  format(' there are ',i6,' samples in group: ',a8)
> 1300  format(' cannot open file ',a30,/,' reenter correct name:')
> 1400  format(' error in read; current format: ',a56,//,
>      &       ' fix file and/or reenter correct format:')
> c......................................................................
> c
> c...write intro to screen
>       write(iot,1000)nfile,ny,nx
> c...get in and out filenames
>
>       write(iot,1005)
>       read(iit,*)nr
>       do 50 i = 1,nr
>          write(iot,1010)i
>          read(iit,1030)filin(i)
> 50    continue
>       write(iot,1020)
>       read(iit,1030)filout
> c
> c...get input information from user and set pgm for header
> c
>       write(iot,1040)
>       read(iit,1050)rident
>       write(iot,1070)
>       read(iit,*)nv
>       write(iot,1080)
>       read(iit,1050)formin
>       do 60 i = 1,nr
>          write(iot,1100)filin(i)
>          read(iit,1110)gname(i)
> 60    continue
> c
> c...check input information and allow for changes
> c
>       call go(nfile,filin,gname)
> c
> c...now figure number of samples that will fit in 256 bytes
> c
>       nrec=64/nv
>       write(iot,1120)nrec
> c
> c...open output file
> c
>       open(iout,file=filout,access='direct',
>      *     form='unformatted',status='unknown',recl=256)
> c      rewind(iout)
> c
> c...write header to output file
> c
>       nrd = 1
>       write(iout,rec=nrd)rident,pgm,nc,nr,nv,x0,dx,y0,dy,angle
>       nrd = nrd + 1
> c
> c...proceed to open input data files
> c
>       do 70 if=1,nr
> 1     open(in,file=filin(if),status='old',err=2000)
>       rewind(in)
> c...read in file/group data and recover the number of samples (nsamp)
>
>       j = 1
> 80    read(in,fmt=formin,end=100,err=3000) (data(k,j),k=1,nv)
>       j = j + 1
>       goto 80
> 100   nsamp = j - 1
> c...write group header to binary file
>       write(iout,rec=nrd)gname(if),nsamp
>       nrd = nrd + 1
>       write(iot,1130)nsamp,gname(if)
> c...check to see if last binary record is filled
> c   (use integer division)
>
>       numrec = nsamp / nrec
>       nsleft = nsamp - (nrec*numrec)
>       if(nsleft.eq.0)then
>          numwri = numrec
>       else
>          numwri = numrec + 1
>       endif
> c...see how many zeros are in remainder part of record
>       nremain = 64 - nv * nrec
>
> c...write this data to records (extra space in records are zeros)
>
>       do 110 k=1,numwri
>          write(iout,rec=nrd)((data(j,i),j=1,nv),i=(k-1)*nrec+1,k*nrec),
>      1                                         (zero(ir),ir=1,nremain)
>          nrd = nrd + 1
> 110   continue
> c...zero data array for next file read
> c
>       do 120 kk=1,nx
>          do 130 kkk=1,ny
>             data(kkk,kk) = 0.0
> 130      continue
> 120   continue
>
> c...loop back up and write next file
>       close(in)
>       go to 70
> 2000  write(iot,1300)filin(if)
>       read(iit,1030)filin(if)
>       goto 1
> 3000  write(iot,1400)formin
>       read(iit,1050)formin
>       goto 80
> 70    continue
> c
>       close(iout)
>       stop
>       end
> c-------------------------------------------------------------------
>       subroutine go(nfile,filin,gname)
> c
> c  purpose: allows alteration of input parameters
> c
>       character rident*56,pgm*8,filin(nfile)*30,ac*1,
>      1     filout*30,formin*56,gname(nfile)*8,
>      1     buffer*200
>       common /inout/ iit,iot
>       common /input/ nv,nr,formin,rident,pgm,filout
> c....................................................................
> 100   format(/,' there are ',i3,' input data files:',/)
> 102   format(' file ',i3,' is :',a30)
> 103   format(/,'  enter filename number to change: ',/,
>      1         '  (enter 0 if no changes desired)')
> 104   format(/,' file ',i3,' is presently :',a30,/,
>      1         ' please enter correct filename :')
> 105   format(a30)
> 110   format(/,' there are ',i3,' group names:',/)
> 111   format(' group ',i3,' is :',a8)
> 112   format(/,'  enter groupname number to change: ',/,
>      1         '  (enter 0 if no changes desired)')
> 113   format(/,' group ',i3,' is presently :',a8,/,
>      1         ' please enter correct groupname (8 char):')
> 114   format(a8)
> 300   format(/,' okay, enter changes ... (c/r twice when finished)',/)
> 301   format(a1)
> c....................................................................
> 200     write(iot,400)nr,nr,formin,filout,rident,pgm,nr,nv
> 400     format(/////,
>      1  ' the following are the control parameters:',//,
>      1  ' 1. enter  1  to change any of the ',i3,' input files',/,
>      1  ' 2. enter  2  to change any of the ',i3,' group names',/,
>      1  ' 3. input files format (assumed the same): ',/,4x,a56,/,
>      1  ' 4. output filename: ',/,4x,a30,/,
>      1  ' 5. output file header comment:',/,4x,a56,/,
>      1  ' 6. program title (8 CHARACTERS): ...................',A8,/,
>      1  ' 7. NUMBER OF DATA Groups in output file: ...........',i3,/,
>      1  ' 8. number of variables per sample: .................',i3,//)
>         write(iot,*)' do you wish to change anything ? (y/n) :'
>         read(iit,301)ac
>             if(ac.eq.'n'.or.ac.eq.'n')return
>             write(iot,300)
> c    an 80 character buffer is used to input the changed parameters
> c    from the crt screen.  the number identifying the parameter is
> c    read first then the parameter is read.  changes are accordingly
> c    made.
> c
> 201           read(iit,fmt='(a)',end=200)buffer
>                 read(buffer(1:2),403)chng
> 403             format(f2.0)
> 410             format(f2.0,1x,a8)
> 420             format(f2.0,1x,a30)
> 430             format(f2.0,1x,a56)
> 440             format(f2.0,1x,i3)
>               if(chng.eq.1.0)goto 500
>               if(chng.eq.2.0)goto 600
>               if(chng.eq.3.0)read(buffer,430)chng,formin
>               if(chng.eq.4.0)read(buffer,420)chng,filout
>               if(chng.eq.5.0)read(buffer,430)chng,rident
>               if(chng.eq.6.0)read(buffer,410)chng,pgm
>               if(chng.eq.7.0)read(buffer,440)chng,nr
>               if(chng.eq.8.0)read(buffer,440)chng,nv
>               if(chng.eq.0.0)goto 200
>               goto 201
> c
> 500   write(iot,100)nr
>       do 501 i=1,nr
>          write(iot,102)filin(i)
> 501   continue
>       write(iot,103)
>       read(iit,*)ichng
>       if(ichng.eq.0)goto 200
>       write(iot,104)ichng,filin(ichng)
>       read(iit,105)filin(ichng)
>       goto 500
> 600   write(iot,110)nr
>       do 601 i=1,nr
>          write(iot,111)gname(i)
> 601   continue
>       write(iot,112)
>       read(iit,*)ichng
>       if(ichng.eq.0)goto 200
>       write(iot,113)ichng,gname(ichng)
>       read(iit,114)gname(ichng)
>       goto 600
> c
>       end
>
>  ________________________________
>  From: proj-bounces at lists.maptools.org on behalf of Eric Miller
> Sent: Tue 26-Feb-08 12:00
> To: PROJ.4 and general Projections Discussions
> Subject: Re: [Proj] lla file format
>
>
>
>
>
> As far as I can tell, it is only documented in the Fortran source code for
> nadgrd.  Get the sources and read "readme.grd" for some help...
>
> """
> The program NADCON reads binary grids that contain NAD 27 to NAD 83
> latitude and longitude shifts.  The first record in a grid file consists of
> header information.  All the other records consist of FORTRAN REAL*4
> numbers.  The grid files are unformatted and direct access.
>
> ...
> """
>
> >>> On 2/16/2008 at 5:30 PM, "Richard Greenwood"
> <richard.greenwood at gmail.com>
> wrote:
> > Can anyone point me to documentation of the *.lla file format?
>
>
> _______________________________________________
> Proj mailing list
> Proj at lists.maptools.org
> http://lists.maptools.org/mailman/listinfo/proj
>
> _______________________________________________
>  Proj mailing list
>  Proj at lists.maptools.org
>  http://lists.maptools.org/mailman/listinfo/proj
>



-- 
Richard Greenwood
richard.greenwood at gmail.com
www.greenwoodmap.com


More information about the Proj mailing list