o
    Ffk                     @  s:  d Z ddlmZ ddlmZ ddlmZmZ ddlZddlm	Z	 ddl
ZddlmZmZmZmZmZ ddlmZmZ dd	lmZ dd
lmZ ddlZddlmZmZmZ ddlmZ ddl m!  m"  m#Z$ ddl%m&Z& e	rxddl'm(Z(m)Z)m*Z* edZ+edZ,d ddZ-d!ddZ.G dd dZ/G dd de&ej0Z1dS )"a  
Read SAS7BDAT files

Based on code written by Jared Hobbs:
  https://bitbucket.org/jaredhobbs/sas7bdat

See also:
  https://github.com/BioStatMatt/sas7bdat

Partial documentation of the file format:
  https://cran.r-project.org/package=sas7bdat/vignettes/sas7bdat.pdf

Reference for binary data compression:
  http://collaboration.cmc.ec.gc.ca/science/rpn/biblio/ddj/Website/articles/CUJ/1992/9210/ross/ross.htm
    )annotations)abc)datetime	timedeltaN)TYPE_CHECKING)read_double_with_byteswapread_float_with_byteswapread_uint16_with_byteswapread_uint32_with_byteswapread_uint64_with_byteswap)Parserget_subheader_index)cast_from_unit_vectorized)EmptyDataError)	DataFrame	Timestampisna)
get_handle)
ReaderBase)CompressionOptionsFilePath
ReadBufferz
1970-01-01z
1960-01-01sas_datetimefloatunitstrc                 C  sR   t | rtjS |dkrtdddt| d S |dkr%tdddt| d S td)Ns     )Zsecondsd)Zdayszunit must be 'd' or 's')r   pdZNaTr   r   
ValueError)r   r    r"   6lib/python3.10/site-packages/pandas/io/sas/sas7bdat.py_parse_datetimeA   s   r$   sas_datetimes	pd.Seriesreturnc                 C  sl   t t d}|dkr#t| jddd}|d| }tj|| jddS t	j
| dd| }tj|d	| jdd
S )a  
    Convert to Timestamp if possible, otherwise to datetime.datetime.
    SAS float64 lacks precision for more than ms resolution so the fit
    to datetime.datetime is ok.

    Parameters
    ----------
    sas_datetimes : {Series, Sequence[float]}
       Dates or datetimes in SAS
    unit : {'d', 's'}
       "d" if the floats represent dates, "s" for datetimes

    Returns
    -------
    Series
       Series of datetime64 dtype or datetime.datetime.
    r   ms)r   Zout_unitzM8[ms]FindexcopyzM8[D]dtypezM8[s]r-   r*   r+   )_sas_origin_unix_originZas_unitr   Z_valuesviewr    Seriesr*   npZarray)r%   r   ZtdZmillisZdt64msZvalsr"   r"   r#   _convert_datetimesO   s   r4   c                   @  sH   e Zd ZU ded< ded< ded< ded< ded< ded	< dddZdS )_Columnintcol_idstr | bytesnamelabelformatbytesctypelengthr'   Nonec                 C  s(   || _ || _|| _|| _|| _|| _d S N)r7   r9   r:   r;   r=   r>   )selfr7   r9   r:   r;   r=   r>   r"   r"   r#   __init__u   s   

z_Column.__init__N)r7   r6   r9   r8   r:   r8   r;   r8   r=   r<   r>   r6   r'   r?   )__name__
__module____qualname____annotations__rB   r"   r"   r"   r#   r5   m   s   
 r5   c                   @  sJ  e Zd ZU dZded< ded< 								dYdZddZd[ddZd[ddZd[dd Zd\d!d"Z	d\d#d$Z
d]d&d'Zd^d*d+Zd_d,d-Zd`d/d0Zdad2d3Zd\d4d5Zdbd6d7Zd\d8d9Zd\d:d;Zdcd<d=Zdcd>d?Zdcd@dAZdcdBdCZdcdDdEZdcdFdGZdcdHdIZdcdJdKZdddedMdNZdOdP Zd]dQdRZdSdT ZdfdWdXZ dS )gSAS7BDATReadera  
    Read SAS files in SAS7BDAT format.

    Parameters
    ----------
    path_or_buf : path name or buffer
        Name of SAS file or file-like object pointing to SAS file
        contents.
    index : column identifier, defaults to None
        Column to use as index.
    convert_dates : bool, defaults to True
        Attempt to convert dates to Pandas datetime values.  Note that
        some rarely used SAS date formats may be unsupported.
    blank_missing : bool, defaults to True
        Convert empty strings to missing values (SAS uses blanks to
        indicate missing character variables).
    chunksize : int, defaults to None
        Return SAS7BDATReader object for iterations, returns chunks
        with given number of lines.
    encoding : str, 'infer', defaults to None
        String encoding acc. to Python standard encodings,
        encoding='infer' tries to detect the encoding from the file header,
        encoding=None will leave the data in binary format.
    convert_text : bool, defaults to True
        If False, text variables are left as raw bytes.
    convert_header_text : bool, defaults to True
        If False, header text, including column names, are left as raw
        bytes.
    r6   _int_lengthzbytes | None_cached_pageNTinferpath_or_bufFilePath | ReadBuffer[bytes]convert_datesboolblank_missing	chunksize
int | Noneencoding
str | Noneconvert_textconvert_header_textcompressionr   r'   r?   c
           
   	   C  s   || _ || _|| _|| _|| _|| _|| _d| _d| _g | _	g | _
g | _g | _g | _d | _g | _g | _g | _d| _d| _d| _t|dd|	d| _| jj| _| j| j| j| j| j| j| j| jd g	| _ z| !  | "  W d S  t#yx   | $   w )Nzlatin-1    r   rbF)Zis_textrV   )%r*   rM   rO   rP   rR   rT   rU   default_encodingrV   column_names_rawcolumn_namescolumn_formatscolumns%_current_page_data_subheader_pointersrI   _column_data_lengths_column_data_offsets_column_types_current_row_in_file_indexZ_current_row_on_page_indexr   handlesZhandle_path_or_buf_process_rowsize_subheader_process_columnsize_subheader_process_subheader_counts_process_columntext_subheader_process_columnname_subheader#_process_columnattributes_subheader_process_format_subheader_process_columnlist_subheader_subheader_processors_get_properties_parse_metadata	Exceptionclose)
rA   rK   r*   rM   rO   rP   rR   rT   rU   rV   r"   r"   r#   rB      sT   
zSAS7BDATReader.__init__
np.ndarrayc                 C     t j| jt jdS )z5Return a numpy int64 array of the column data lengthsr,   )r3   asarrayr_   int64rA   r"   r"   r#   column_data_lengths      z"SAS7BDATReader.column_data_lengthsc                 C  rs   )z0Return a numpy int64 array of the column offsetsr,   )r3   rt   r`   ru   rv   r"   r"   r#   column_data_offsets   rx   z"SAS7BDATReader.column_data_offsetsc                 C  s   t j| jt ddS )zj
        Returns a numpy character array of the column types:
           s (string) or d (double)
        ZS1r,   )r3   rt   ra   r-   rv   r"   r"   r#   column_types   s   zSAS7BDATReader.column_typesc                 C  s   | j   d S r@   )rc   rq   rv   r"   r"   r#   rq      s   zSAS7BDATReader.closec                 C  s  | j d | j d| _| jdttj tjkrtd| tj	tj
}|tjkr:d| _d| _tj| _tj| _nd| _tj| _tj| _d| _| tjtj}|tjkrYtj}nd}| tjtj}|dkrqd	| _tjd
k| _n	d| _tjdk| _| tjtjd }|tj v rtj | | _!| j"dkr| j!| _"nd| d| _!t#ddd}| $tj%| tj&}|t'j(|dd | _)| $tj*| tj+}|t'j(|dd | _,| -tj.| tj/| _0| j | j0d }|  j|7  _t| j| j0krtd| -tj1| tj2| _3d S )Nr   i   z'magic number mismatch (not a SAS file?)T   F      <Zbig>littlerJ   zunknown (code=)r   r   r   )r   z*The SAS7BDAT file appears to be truncated.)4rd   seekreadrI   lenconstmagicr!   _read_bytesZalign_1_offsetZalign_1_lengthZu64_byte_checker_valueU64rH   Zpage_bit_offset_x64_page_bit_offsetZsubheader_pointer_length_x64_subheader_pointer_lengthZpage_bit_offset_x86Zsubheader_pointer_length_x86Zalign_2_offsetZalign_2_lengthZalign_1_checker_valueZalign_2_valueZendianness_offsetZendianness_length
byte_ordersys	byteorderneed_byteswapZencoding_offsetZencoding_lengthZencoding_namesZinferred_encodingrR   r   _read_floatZdate_created_offsetZdate_created_lengthr    Zto_timedeltaZdate_createdZdate_modified_offsetZdate_modified_lengthZdate_modified
_read_uintZheader_size_offsetZheader_size_lengthZheader_lengthZpage_size_offsetZpage_size_length_page_length)rA   bufZalign1Zepochxr"   r"   r#   rn      sd   





zSAS7BDATReader._get_propertiesr   c                 C  s(   | j | jpdd}|jr|   t|S )Nr   )nrows)r   rP   emptyrq   StopIteration)rA   dar"   r"   r#   __next__A  s
   zSAS7BDATReader.__next__offsetwidthc                 C  sN   | j d usJ |dkrt| j || jS |dkrt| j || jS |   td)Nr|   r{   zinvalid float width)rI   r   r   r   rq   r!   rA   r   r   r"   r"   r#   r   I  s   

zSAS7BDATReader._read_floatc                 C  s~   | j d usJ |dkr| |dd S |dkrt| j || jS |dkr+t| j || jS |dkr7t| j || jS |   td)Nr   r      r|   r{   zinvalid int width)rI   r   r	   r   r
   r   rq   r!   r   r"   r"   r#   r   X  s"   


zSAS7BDATReader._read_uintr>   c                 C  sB   | j d usJ || t| j kr|   td| j |||  S )NzThe cached page is too small.)rI   r   rq   r!   rA   r   r>   r"   r"   r#   r   l  s
   zSAS7BDATReader._read_bytesr8   c                 C  s   |  | ||dS )N     )_convert_header_textr   rstripr   r"   r"   r#   _read_and_convert_header_texts  s   z,SAS7BDATReader._read_and_convert_header_textc                 C  sV   d}|s)| j | j| _t| jdkrd S t| j| jkr!td|  }|rd S d S )NFr   z2Failed to read a meta data page from the SAS file.)rd   r   r   rI   r   r!   _process_page_meta)rA   Zdoner"   r"   r#   ro   x  s   zSAS7BDATReader._parse_metadatac                 C  sZ   |    tjtjtjg }| j|v r|   | jtjk}| jtjk}t|p+|p+| j	g kS r@   )
_read_page_headerr   page_meta_typesZpage_amd_typepage_mix_type_current_page_type_process_page_metadatapage_data_typerN   r^   )rA   ptZis_data_pageZis_mix_pager"   r"   r#   r     s   
z!SAS7BDATReader._process_page_metac                 C  s^   | j }tj| }| |tjtj@ | _tj| }| |tj| _	tj
| }| |tj| _d S r@   )r   r   Zpage_type_offsetr   Zpage_type_lengthZpage_type_mask2r   Zblock_count_offsetZblock_count_lengthZ_current_page_block_countZsubheader_count_offsetZsubheader_count_length_current_page_subheaders_count)rA   
bit_offsetZtxr"   r"   r#   r     s   



z SAS7BDATReader._read_page_headerc                 C  s  | j }t| jD ]}tj| }|| j|  }| || j}|| j7 }| || j}|| j7 }| |d}|d7 }| |d}|dksG|tjkrHq| 	|| j}	t
|	}
| j|
 }|d u r|tjdfv }|tjk}| jrx|rx|rx| j||f q|   td|	 ||| qd S )Nr   r   zUnknown subheader signature )r   ranger   r   Zsubheader_pointers_offsetr   r   rH   Ztruncated_subheader_idr   r   rm   Zcompressed_subheader_idZcompressed_subheader_typerV   r^   appendrq   r!   )rA   r   ir   Ztotal_offsetZsubheader_offsetZsubheader_lengthZsubheader_compressionZsubheader_typeZsubheader_signatureZsubheader_indexZsubheader_processorf1f2r"   r"   r#   r     s<   





z%SAS7BDATReader._process_page_metadatac                 C  s   | j }|}|}| jr|d7 }|d7 }n|d7 }|d7 }| |tj|  || _| |tj|  || _| |tj|  || _	| |tj
|  || _tj| }| || || _| |d| _| |d| _d S )Ni  i  ib  iz  r   )rH   r   r   r   Zrow_length_offset_multiplierZ
row_lengthZrow_count_offset_multiplier	row_countZcol_count_p1_multipliercol_count_p1Zcol_count_p2_multipliercol_count_p2Z'row_count_on_mix_page_offset_multiplierZ_mix_page_row_count_lcs_lcp)rA   r   r>   int_lenZ
lcs_offsetZ
lcp_offsetZmxr"   r"   r#   re     s4   

z)SAS7BDATReader._process_rowsize_subheaderc                 C  sX   | j }||7 }| ||| _| j| j | jkr*td| j d| j d| j d d S d S )Nz Warning: column count mismatch (z + z != z)
)rH   r   column_countr   r   print)rA   r   r>   r   r"   r"   r#   rf     s   
z,SAS7BDATReader._process_columnsize_subheaderc                 C     d S r@   r"   r   r"   r"   r#   rg     s   z(SAS7BDATReader._process_subheader_countsc           	      C  s  || j 7 }| |tj}| ||}|d| d}| j| t| jdkrd}tj	D ]}||v r5|}q-|| _
|| j 8 }|d }| jrI|d7 }| || j}|d}|dkrwd| _|d }| jrg|d7 }| || j}|d| j | _nB|tjkr|d	 }| jr|d7 }| || j}|d| j | _n"| jdkrd| _|d }| jr|d7 }| || j}|d| j | _t| d
r| | j| _d S d S d S )Nr   r   r   rW      r|           (   creator_proc)rH   r   r   Ztext_block_size_lengthr   r   rZ   r   r   Zcompression_literalsrV   r   r   r   r   Zrle_compressionhasattrr   )	rA   r   r>   Ztext_block_sizer   Z	cname_rawZcompression_literalZclZoffset1r"   r"   r#   rh     sT   






"z,SAS7BDATReader._process_columntext_subheaderc                 C  s   | j }||7 }|d|  d d }t|D ]Q}|tj|d   tj }|tj|d   tj }|tj|d   tj }| |tj}	| |tj	}
| |tj
}| j|	 }||
|
|  }| j| | qd S )Nr      r{   r   )rH   r   r   Zcolumn_name_pointer_lengthZ!column_name_text_subheader_offsetZcolumn_name_offset_offsetZcolumn_name_length_offsetr   Z!column_name_text_subheader_lengthZcolumn_name_offset_lengthZcolumn_name_length_lengthrZ   r[   r   r   )rA   r   r>   r   Zcolumn_name_pointers_countr   Ztext_subheaderZcol_name_offsetZcol_name_lengthidx
col_offsetZcol_lenZname_rawZcnamer"   r"   r#   ri      sB   
z,SAS7BDATReader._process_columnname_subheaderc           
      C  s   | j }|d|  d |d  }t|D ]Y}|| tj ||d   }|d|  tj ||d   }|d|  tj ||d   }| ||}	| j|	 | |tj	}	| j
|	 | |tj}	| j|	dkridnd qd S )Nr   r   r{   r      d   s)rH   r   r   Zcolumn_data_offset_offsetZcolumn_data_length_offsetZcolumn_type_offsetr   r`   r   Zcolumn_data_length_lengthr_   Zcolumn_type_lengthra   )
rA   r   r>   r   Zcolumn_attributes_vectors_countr   Zcol_data_offsetZcol_data_lenZ	col_typesr   r"   r"   r#   rj   A  s,   
z2SAS7BDATReader._process_columnattributes_subheaderc                 C  r   r@   r"   r   r"   r"   r#   rl   [  s   z,SAS7BDATReader._process_columnlist_subheaderc                 C  sx  | j }|tj d|  }|tj d|  }|tj d|  }|tj d|  }|tj d|  }|tj d|  }	| |tj	}
t
|
t| jd }| |tj}| |tj}| |tj}t
|t| jd }| |tj}| |	tj}| j| }| ||||  }| j| }| ||||  }t| j}t|| j| ||| j| | j| }| j| | j| d S )N   r   )rH   r   Z)column_format_text_subheader_index_offsetZcolumn_format_offset_offsetZcolumn_format_length_offsetZ(column_label_text_subheader_index_offsetZcolumn_label_offset_offsetZcolumn_label_length_offsetr   Z)column_format_text_subheader_index_lengthminr   rZ   Zcolumn_format_offset_lengthZcolumn_format_length_lengthZ(column_label_text_subheader_index_lengthZcolumn_label_offset_lengthZcolumn_label_length_lengthr   r]   r5   r[   ra   r_   r\   r   )rA   r   r>   r   Ztext_subheader_formatZcol_format_offsetZcol_format_lenZtext_subheader_labelZcol_label_offsetZcol_label_lenr   Z
format_idxZformat_startZ
format_lenZ	label_idxZlabel_startZ	label_lenZlabel_namesZcolumn_labelZformat_namesZcolumn_formatZcurrent_column_numbercolr"   r"   r#   rk   _  sX   


	z(SAS7BDATReader._process_format_subheaderr   c                 C  s   |d u r| j d ur| j }n|d u r| j}t| jdkr#|   td|dkr0| j| jkr0t S t|| j| j }| j	d}| j	d}t
j||ftd| _t
j|d| ft
jd| _d| _t| }|| |  }| jd urw|| j}|S )Nr   zNo columns to parse from filer   r   r,   r{   )rP   r   r   ra   rq   r   rb   r   r   countr3   r   object_string_chunkZzerosZuint8_byte_chunk_current_row_in_chunk_indexr   r   _chunk_to_dataframer*   Z	set_index)rA   r   Zndnsprsltr"   r"   r#   r     s*   

zSAS7BDATReader.readc                 C  s   g | _ | j| j| _t| jdkrdS t| j| jkr3|   dt| jdd| jdd}t||   | j	t
jv rA|   | j	t
jt
jt
jg vrQ|  S dS )Nr   Tz-failed to read complete page from file (read r   z of z bytes)F)r^   rd   r   r   rI   r   rq   r!   r   r   r   r   r   r   r   _read_next_page)rA   msgr"   r"   r#   r     s,   
zSAS7BDATReader._read_next_pagec                 C  sn  | j }| j}t|| |}i }d\}}t| jD ]}| j| }| j| dkrm| j|d d f j| jd d}	t	j
|	tj|dd||< | jrh| j| tjv rWt|| d||< n| j| tjv rht|| d||< |d7 }q| j| d	krt	j
| j|d d f |dd
||< | jr| jd ur| || j||< |d7 }q|   tdt| j|  t|| j|dd}
|
S )N)r   r   r   r   r,   Fr.   r   r   r   r)   zunknown column type )r]   r*   r+   )r   rb   r   r   r[   ra   r   r1   r   r    r2   r3   Zfloat64rM   r\   r   Zsas_date_formatsr4   Zsas_datetime_formatsr   rT   rR   _decode_stringr   rq   r!   reprr   )rA   nmZixr   ZjsZjbjr9   Zcol_arrZdfr"   r"   r#   r     s2   
 
"
z"SAS7BDATReader._chunk_to_dataframec                 C  s   | | jp| jS r@   )decoderR   rY   rA   br"   r"   r#   r     s   zSAS7BDATReader._decode_stringr   r<   c                 C  s   | j r| |S |S r@   )rU   r   r   r"   r"   r#   r     s   
z#SAS7BDATReader._convert_header_text)NTTNNTTrJ   )rK   rL   rM   rN   rO   rN   rP   rQ   rR   rS   rT   rN   rU   rN   rV   r   r'   r?   )r'   rr   )r'   r?   )r'   r   )r   r6   r   r6   )r   r6   r   r6   r'   r6   )r   r6   r>   r6   )r   r6   r>   r6   r'   r8   )r'   rN   )r   r6   r>   r6   r'   r?   r@   )r   rQ   r'   r   )r   r<   r'   r8   )!rC   rD   rE   __doc__rF   rB   rw   ry   rz   rq   rn   r   r   r   r   r   ro   r   r   r   re   rf   rg   rh   ri   rj   rl   rk   r   r   r   r   r   r"   r"   r"   r#   rG      sP   
 
?




F









+



-
!

7
rG   )r   r   r   r   )r%   r&   r   r   r'   r&   )2r   Z
__future__r   collectionsr   r   r   r   typingr   Znumpyr3   Zpandas._libs.byteswapr   r   r	   r
   r   Zpandas._libs.sasr   r   Zpandas._libs.tslibs.conversionr   Zpandas.errorsr   Zpandasr    r   r   r   Zpandas.io.commonr   Zpandas.io.sas.sas_constantsioZsasZsas_constantsr   Zpandas.io.sas.sasreaderr   Zpandas._typingr   r   r   r0   r/   r$   r4   r5   IteratorrG   r"   r"   r"   r#   <module>   s0    

