docs: def _add_prj_file(original_gis_file, new_gis_file):
          """
          Adds projection file
          """
          out_prj_file = "{0}.prj".format(os.path.splitext(new_gis_file)[0])
          if original_gis_file.endswith(".shp"):
              dataset = ogr.Open(original_gis_file)
              layer = dataset.GetLayer()
              spatial_ref = layer.GetSpatialRef()
              spatial_ref.MorphToESRI()
              spatial_ref_str = spatial_ref.ExportToWkt()
          else:
              dataset = gdal.Open(original_gis_file)
              spatial_ref_str = dataset.GetProjection()
  
          with open(out_prj_file, 'w') as prj_file:
              prj_file.write(spatial_ref_str)