Skip to content

Instantly share code, notes, and snippets.

#!/usr/bin/env python
#buffer.py
import sys
import os
from osgeo import ogr
def buffer(infile,outfile,buffdist):
try:
ds_in=ogr.Open( infile )
lyr_in=ds_in.GetLayer( 0 )
@percursoaleatorio
percursoaleatorio / compressTIFS.bat
Last active December 15, 2015 16:59
GDAL batch file to compress TIFF orthophotos, and place CRS information into the file header.
@echo off
rem “compressTIFS.bat”
rem >> compressTIFS diretoria_input directoria_output
rem nery v2013.03
rem OBJETIVOS: Processar ortofotomapas
rem 1) integrar a informação de georreferenciação do ficheiro TFW no header do ficheiro GeoTiff;
rem 2) forçar a atribuição e escrita da informação sobre o sistema EPSG 3763 (PT-TM06) no header do ficheiro GeoTiff
rem 3) marcar a inexistência de um canal alpha na imagem (requer GDAL 20DEV)
rem 4) comprimir o ficheiro resultante com compressão DEFLATE e utilizando o algoritmo de predição adequado para imagens de satélite ou fotografias aéreas
@percursoaleatorio
percursoaleatorio / CRS_PT
Last active March 5, 2016 13:00
PROJ.4 config file for PT Coordinate Reference Systems. Configuration file for PROJ.4 coordinate transformations. Contains all commonly used CRS and coordinate transformations for Portugal. Includes obsolete CRSs. Includes CRSs for the Madeira and Azores Islands.
# EXAMPLES:
# ogr2ogr -f "ESRI Shapefile" -s_srs "+init=CRS_PT:4274_9615 +wktext" -t_srs "+init=CRS_PT:3763_NIL +wktext" output.shp input.shp
# gdalwarp -s_srs "+init=CRS_PT:3763_NIL +wktext" -t_srs "+init=CRS_PT:3857_NIL +wktext" -r bilinear input.tif output.tif
# SYNTAX: EACH TRANSFORMATION IS IDENTIFIED AS "epsgCoordinateReferenceSystemCode_epsgCoordinateOperationMethod"
# IF, WITHIN A CRS AREA OF USE, DIFFERENT TRANSFORMATIONS ARE AVAILABLE FOR SPECIFIC AREAS THEN EACH TRANSFORMATION IS IDENTIFIED AS "epsgCoordinateReferenceSystemCode_epsgCoordinateOperationMethod_epsgAreaOfUse"
#
### ENUMERATION OF epsg Coordinate Reference System Code
# 4326 "WGS84/Geographic 2D"
# 4258 "ETRF89/Geographic 2D"