o
    Wid                     @   sF   d dl mZ d dlZd dlZdddZdd Zdd	d
ZdddZdS )    )gdalNmouza_map.pointsc                 C   s   ddddddddddddddddddddg}t |d,}|d	 |D ]}||d
  d|d  d|d  d|d  d q'W d    |S 1 sNw   Y  |S )Nd   g9#JV@g镒6@)pixelXpixelYmapXmapYi  g]CV@gǺ6@wz%mapX,mapY,pixelX,pixelY,enable,order
r   ,r   r   r   z,1,0
)openwrite)
image_pathoutput_pathgcpsfZgcp r   &/var/www/html/qgis-automation/utils.pygenerate_points_file   s   
4
r   c           	      C   s   g }t | d@}|D ]4}|ds|drq
| d}t|dk r%q
tt|d d \}}}}|t	||d|| q
W d    |S 1 sJw   Y  |S )Nr#r   r
      r   )
r   
startswithstripsplitlenmapfloatappendr   ZGCP)	Zpoints_filer   r   linepartsr   r   r   r   r   r   r   	load_gcps   s   
		r    mouza_map_georef.tifc                 C   s0   t |}d}tj|| |d tj||dd |S )Nztemp_with_gcps.tif)ZGCPsT)Ztps)r    r   Z	TranslateZWarp)r   points_pathr   r   Z	temp_pathr   r   r   georeference_image    s
   r#   tiles_outputc              	   C   s   t ddddd| |g |S )Nzgdal2tiles.pyz-zz14-18z-wZgoogle)
subprocessrun)Zgeotiff_path	tiles_dirr   r   r   generate_tiles'   s   r(   )r   )r!   )r$   )Zosgeor   r%   osr   r    r#   r(   r   r   r   r   <module>   s    

