2
TanDEM GeoTIFFs are of type "Pixel is point"
gdalinfo TDM1_DEM__30_N49E010_DEM.tif
Driver: GTiff/GeoTIFF
Files: TDM1_DEM__30_N49E010_DEM.tif
Size is 1201, 1201
Coordinate System is:
GEOGCRS["WGS 84",
DATUM["World Geodetic System 1984",
ELLIPSOID["WGS 84",6378137,298.257223563,
LENGTHUNIT["metre",1]]],
PRIMEM["Greenwich",0,
...
2
I've found an answer on the project-page of exiftool itself.
Hereafter, I'll summarize what I found out with the kind help of the page-author Phil Harvey.
To copy all GeoTiff tags from one file to another, do this:
exiftool -tagsfromfile SRCFILE -GeoTiffDirectory -GeoTiffDoubleParams -GeoTiffAsciiParams DSTFILE
That does the trick.
I assume that those ...
1
You can download the CSVs first and look at it in QGIS
https://dataforgood.fb.com/docs/high-resolution-population-density-maps-demographic-estimates-documentation/#data-format
Then you can move on to GeoTiff which can also be viewed in QGIS. When you have that working you can export the data into a vector format like GeoJSON and use it in CartoDB or another ...
1
Pillow isn't a geospatial library, I recommend looking at rasterio resampling
For example:
import rasterio
from rasterio.enums import Resampling
# Register GDAL format drivers and configuration options with a
# context manager.
with rasterio.Env():
with rasterio.open('/path/to/input.tif') as dataset:
data = dataset.read(1, out_shape=(1200, ...
1
I was able to recover the spatial referencing of the original georefernced tif using gdal translate.
import gdal, osr
original = gdal.Open('path to tif') # the original georeferenced tif
in_tif = 'path to watermarked tif' # the one that lost spatial reference
out_tif = 'path to tif to be created' # the georeferenced output tif
rf = ...
1
gdal_retile -s_srs option sets the source spatial reference system, from your question it is already set as something different. -s_srs is for situations when you know the projection of the data but the data (or associated metadata files) doesn't know or is incorrect. So never use this option if there is a projection set.
If you would like to reproject the ...
1
It looks like the only piece missing from your write call is a transform:
with rasterio.open(imgPath) as src:
win = Window(1000, 1000, 500, 300)
tile = src.read(1, window=win)
win_transform = src.window_transform(win)
with rasterio.open(
resPath,
'w',
driver='GTiff',
width=500,
height=300,
count=1,
dtype=tile.dtype,
...
1
I got the information from here that gdal 3.0.2 does not yet include the "fast overviews for big TIFF" feature It will be included in the 3.1 version
Only top voted, non community-wiki answers of a minimum length are eligible
