>>> from liblas import srs >>> from liblas import point >>> from liblas import header >>> import liblas >>> from liblas import file >>> f = file.File('../test/data/1.2_3.las',mode='r') >>> s = f.header.srs >>> s.proj4 '' >>> del f >>> def new_offset(old_scale, new_scale, old_offset, x): ... return (new_scale*(x - old_offset) - old_scale*x)/(-1.0*old_scale)