A Geospatial Data Abstraction Library (GDAL) image mosaicing, reprojection, and warping utility

learn more… | top users | synonyms

1
vote
0answers
12 views

Convert MODIS HDF file in sinusoidal projection into GeoTiff with

I'm looking for a way to convert MODIS LST data from HDF to Tiff format, and reproject from sinusoidal to WGS84 - all using python GDAL. An example file is here (I only care about band 0) : ...
0
votes
0answers
8 views

Resampling using gdalwarp gives me no data

I have an image with the following information. Driver: GTiff/GeoTIFF Files:test_DEM.tiff Size is 91, 73 Coordinate System is: GEOGCS["WGS 84", DATUM["WGS_1984", ...
0
votes
0answers
22 views

Resampling geotiff images? [closed]

I am trying to resample an image (geotiff) to a 5 m x 5m resolution using python. I am new to gdal and not sure how to go about it. Using gdal.ReprojectImage is what I should be doing, It would be ...
3
votes
0answers
21 views

How to blend edges of aerial imagery with gdalwarp

I need to fade/blend the edges of a series of aerial images with gdal. The closest thing that I've been able to accomplish is using the -cutline and -cblend switches of gdalwarp to get the edges to ...
0
votes
0answers
30 views

How to convert raster in UTM to lat/lon using GDAL (or other tool)

I have a .tiff file given with UTM coordinates that I'd like to convert to a .nc (NetCDF) file with latitude/longitude coordinates. I've tried with gdal_translate and gdal_warp but these did not work. ...
2
votes
1answer
28 views

Corners of GeoTiff's not lining up after doing gdalwarp EPSG:3006 -> EPSG:4326

I have some GeoTiff's that are in EPSG:3006 (SWEREF99 TM) atm, and I'm trying to convert them to EPSG:4326 (WGS84) because some other software needs it in this projection. But after doing the ...
1
vote
0answers
56 views

Slow performance clipping raster with GDAL

I have a few large LiDAR DEMs that I want to merge. Each DEM has some glitches around the borders which need to be removed in order to get a seamless final product. I manually traced around each DEM ...
4
votes
1answer
59 views

Georeferencing Himawari-8 in GDAL (or other)

Does anyone know the appropriate project information to georeference (and hence reproject) Himawari-8 AHI data in the IR? The images are 5500x5500 pixels. Currently I try these: ulx=-5570248.832537 ...
1
vote
0answers
65 views

How to make NoData transparent?

I have been using python to turn a netcdf into a geotiff. Once I have the geotiff, I am using gdalwarp to transfer it into the right projection. Then use gdal_translate to get the tif into a vrt ...
1
vote
0answers
31 views

Gdalwarp UTM reprojection causes image distortions

I'm trying to re-project Sentinel-2 images to a common UTM zone in order to create a mosaic. For this I'm using gdalwarp with the following syntax: C:/OSGeo4W64/OSGeo4W.bat; gdalwarp -t_srs ...
1
vote
0answers
35 views

gdalwarp won't run with “cubic” resampling method

I am having an issue with running the gdalwarp command from command prompt. Basically I have a raster WGS84 file, which I would like to reproject to spatial reference which is using the Azimuthal ...
1
vote
1answer
36 views

gdalwarp output shifted 15 KM

This is similar to this question regarding shifts in raster output, but the difference I'm experiencing is fairly large-- the output raster is shifted 15 KM to the south. Here's an image: When ...
1
vote
1answer
14 views

Error using gdalwarp in EPSG:4326

I tried to change the projection of AMSR-2 in Python using: 'gdalwarp -geoloc -t_srs EPSG:4326 C:\GW1AM2_201301311114_050A_L1SGBTBR_1110110.vrt C:\test100.tif' It said there are something wrong ...
1
vote
1answer
33 views

gdalwarp causing shift in pixels

Can somebody help me with gdalwarp issue? I am having an input raster and would like to clip it with polygons. After running the below code, while overlaying output raster on input raster, i see there ...
0
votes
0answers
31 views

How to improve poor image quality of “gdalwarp” operation result?

I'm trying to improve the result of gdalwarp operation. The image width and height are 6732x4721 GCP's are: -gcp 0 0 38.005301407753 55.048725423557 -gcp 6732 0 38.184943751144 ...
0
votes
0answers
35 views

gdalwarp leaves black boundary of cropped geotiff with cutline

This is a new question based on a response on the thread white page border when converting geoPDF to geoTiff with gdal_translate I am using gdalwarp as described to do a cutline of geotiff images ...
0
votes
0answers
25 views

intersect geotiff with bbox

In a project using gdal I need to crop a geotiff file with bbox in pixel or coordinates epsg:4326. This step is ok. The next step is to manage cases where the image size is "lower" than the bbox size ...
2
votes
1answer
33 views

Using gdalwarp to remove one row of pixels or add one row of pixels

I am trying to use GDAL_calc but my two rasters are off by one pixel in the y dimension. One is 256005 the other is 256004 so GDAL_calc fails with a multiplication. Is it possible in gdalwarp to ...
1
vote
1answer
71 views

GRIB to NetCDF with GDAL, specifying dimensions

I have a GRIB dataset representing temperature over North America, in Lambert Conformal Conic (with specified central meridian, latitude of origin, etc.). It has what should be understood as a ...
4
votes
3answers
33 views

gdalwarp to geographic coordinate system

I'm getting the folllowing Error when trying to re-project a raster from NAD_1983_Albers to GCS_North_American_1983. I used the proj4 defenition from going to Properties in QGIS of the layer whose ...
1
vote
1answer
51 views

GDAL: gdalwarp wrong transformation to “South_Pole_Lambert_Azimuthal_Equal_Area” EPSG:102020

I am trying to reproject a raster that is defined in EPSG:4326 to EPSG:102020. Command used gdalwarp -ot Float32 -s_srs EPSG:4326 -t_srs EPSG:102020 -r near -of GTiff -co COMPRESS=DEFLATE -co ...
0
votes
1answer
52 views

Reprojecting and retiling a large set of tiles

I am trying to understand the process to re-project and re-tile a set of raster (GeoTiff) but can't wrap my head around how to use the different GDAL tools to achieve my goal. More concretely, I ...
7
votes
1answer
203 views

gdalwarp producing empty rasters when called from python script, but not via command line

I'm trying to aggregate some raster data (all .tif) up to a coarser resolution (from .05 degree to .25 degree) using gdalwarp in python, but the command is not working. Instead of getting an output ...
2
votes
1answer
72 views

Raster interpolation with gdalwarp, and replicating with python API

My problem is similar to this: Replicating result of gdalwarp using gdal Python bindings. I want to resample a raster using bilinear interpolation, which I've done successfully with gdalwarp by adding ...
1
vote
0answers
33 views

Retaining resolution of MODIS data while reprojecting using gdalwarp

Can anyone help me to use gdal(gdalwarp specifically) for reprojection? I am using it to reproject MODIS data : I originally have an hdf file - MCD43A4.A2005001.h08v05.005.2007351192333.hdf which has ...
10
votes
0answers
257 views

PyQGIS save raster as rendered image, then use GDAL tools on it

I have tried to use the code from QGis Save Raster as Rendered Image. i is a raster layer: pipelayer = i pipeextent = pipelayer.extent() pipewidth, pipeheight = (pipelayer.width(), ...
3
votes
1answer
74 views

Crop raster from vector problem

So I have a shapefile of a vector for Philadelphia's city lines and a raster image of Philadelphia county that I got from the lovely people from USDA:NRCS. I converted the raster from a sid file to ...
1
vote
1answer
158 views

Transform an image with projection EPSG:4326 to EPSG:3857

I have an image with known boundary and projection EPSG:4326 and want to overlay it to openstreetmap using Leaflet.js. I have tried using function imageoverlay, but the image is not aligned with the ...
1
vote
0answers
34 views

getting started designing satellite projection raster

How do I get started designing a map with a satellite projection like this:http://bl.ocks.org/joyrexus/6619251. So far using d3 I've created a map of North Africa w/ topojson and tiles from Mapbox but ...
3
votes
1answer
114 views

How to ‘invert’ or ‘reverse’ a vector cutline when gdalwarp'ing a raster?

I have a raster of North America and I'd like to remove the larger lakes from it, replacing them with the same no-data value as the original raster. So I downloaded Natural Earth's 10 meter lakes ...
1
vote
1answer
13 views

Create single map image from two images with different color tables

Is it possible to merge the color tables of two images when using gdal to create a single image? Or is there another way to merge said images not using a color table?
3
votes
0answers
62 views

Downsampling for beginners using QGIS

I apologise from the outset. I need a full run through of working gdalwarp to downsample a raster to match another i.e. How to resample GeoTIFF images to the same resolution?. Using QGIS 2.12. Using ...
2
votes
1answer
56 views

gdalwarp leaves big hole at 180° longitude when warping from 4326 to custom AEqD

I'm trying to reproject the Natural Earth I 1:50m raster into an Azimuthal Equidistant projection centered on my hometown. As you may know, the source raster covers from -180° to 180° latitude, and ...
1
vote
0answers
83 views

gdalwarp UTM Zone 1 to WGS84 - fails to transform

I have a global dataset consisting of tiles in UTM projection that I want to mosaic and therefore need to reproject to a common projection. Can I solve my reprojection error with options to gdalwarp ...
1
vote
2answers
126 views

gdalwarp - merge two images, but only the no-data values

I am preparing tiff elevation files. For one particular country I have a 10 meter elevation model, but the data is present only for the specific country - so when the image is sliced into smaller ...
0
votes
0answers
21 views

how to change whole raster folder?

I used to gdalwarp to reproject raster tiff file, but I did it one by one How can I manage to do all raster in my folder (about 70 file) ? Is it possible? I am using ubuntu-64bit
0
votes
0answers
68 views

How to transform any GTiff to Lambert Azimuthal Equal Area with GDAL?

I am trying to transform a GTiff EPSG:4326 +proj=longlat +datum=WGS84 +no_defs and a GTiff EPSG:31982 +proj=utm +zone=22 +south +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs both to ...
2
votes
2answers
128 views

projection conversion from local to UTM

I have some images at a local coordinate system with the following specifications: TIFF header file: -gdalinfo im.tif Driver: GTiff/GeoTIFF Coordinate System is: PROJCS["Eastern Ridge Project Grid ...
1
vote
1answer
32 views

resample etopo.grd with gdal_warp gives many null values

I'm having a problem with resampling ETOPO NetCDF GRD to a bigger pixelsize for using it with Generic Mapping Tools. I'm using gdal_warp for the resampling (just for testing with pixelsize of 1°): ...
3
votes
1answer
72 views

Mass reprojection of raster files in a loop [Linux]

I have a folder with numerous .asc raster files that each accompanied by a corresponding .prj file. The file basename is the same for both files. The projection however can vary between all of them ...
1
vote
0answers
68 views

GDAL - replacing rgb value (0,0,0) with (nodata,nodata,nodata)

I have a set of rasters in which I want to replace any rgb value of (0,0,0) with (nodata,nodata,nodata) I've tried using gdal_translate -a_nodata "0,0,0" in.tif out.tif But this seems to convert ...
1
vote
0answers
23 views

Transform stereographic projection using GDAL

How can I read/wrap a JPG image without any embedded projection/coordinate system to an autoreferenced GEOTIFF with Mercator or similar coordinate system using GDAL? The image is 500x500 px and I ...
1
vote
2answers
38 views

gdalwarp change spatial resolution

I just realized that I have Geo-TIFF's in the samme coordinatesystem but different resolution. DOP20 (20cm=1px) and DOP10 (10cm=1px) mix. Now wanted to use gdalwarp to transform the DOP20 TIFF into ...
2
votes
2answers
209 views

How to re-project a netcdf file?

I want to re-project and resample a netcf file but I am getting a usage problem: gdalwarp -t_srs EPSG:3410 -ts 1383 586 NETCFD:"C:\data\product.nc":gsfc output.tiff
1
vote
0answers
30 views

Changes in MOD09A1 quality mask as processed with GDAL

I'm analyzing MOD09A1 spectral data using GDAL/Python. The quality mask is a uint32 dataset and want to pull out the value of the first bit to mask the spectral data (1 == poor data; 0 == good data). ...
1
vote
1answer
57 views

gdal_warp some geotiffs change SRS

Hello I#m quite new to GDAL. I have some Geotiffs that I want reproject. I tried that with one file and it works great, but because the TIFF has to be slightly rotated for the conversion I have some ...
1
vote
0answers
91 views

how to reproject an image programatically in c# using gdal

I have to reproject an image programmatically in c# using gdal library. For the moment i'm a little bit lost in gdal and i don't really now which methods i have to use. I think i will have to ...
2
votes
1answer
81 views

Reproject “unknown” coordinate system to LCC 2SP across lon 180

Starting with a grib file with limited projection information, I am trying to import it as a GeoTiff in Geoserver 2.7.x. The original gdalinfo is: Size is 218, 81 Coordinate System is: ...
0
votes
1answer
114 views

Use lamber 93 imageOverlay with leaflet on top of OSM tiles layer

I'm trying to use a lambert 93 image overlay on top of an Open Street Map title layer. var map = L.map('map'); L.tileLayer('http://{s}.tile.osm.org/{z}/{x}/{y}.png', { }).addTo(map); var coords = { ...
0
votes
0answers
79 views

Define band specific metadata of a GeoTIFF file using GDAL (gdalwarp/gdal_translate)

Does anybody know how to change or set the "Description" option/tag (band specific metadata) of a GeoTIFF file using GDAL? To specify what I mean, this is an example of gdalinfo return from a GeoTIFF ...