Tag Info

Hot answers tagged

3

Use the gdalwarp utility. gdalwarp -s_srs epsg:4326 -t_srs epsg:3857 input.tif output.tif


2

I suspect the issue is the GDAL 1.11 bindings you're fetching from gohlke do not match the internal GDAL inside of Arc*. Your bindings need to be compiled against ESRI's GDAL to work reliably. It's possible they can be made to work, but it is going to be a lot of headache. In short, ESRI needs to provide the gdal_i.lib stub file that GDAL generates as part ...


2

I would suggest you to use this alternative way to do the same: To export raster as TIFF, find the oid of the raters’ lob(large object) using the query: SELECT oid, lowrite(lo_open(oid, 131072), png) As num_bytes FROM ( VALUES (lo_create(0), ST_Astiff( (SELECT rast FROM raster_table WHERE rid = 1) ) ) ) As v(oid,png); Use the PostgreSQL \lo_export command ...


1

The easiest solution is to use union, cascaded_unionor unary_union. All the lines are split at the points of intersection: from shapely.geometry import LineString line1 = LineString([(0, 0), (2, 2),(3,1)]) line2 = LineString([(2, 0), (2, 1),(1,2)]) print line1.intersection(line2) POINT (1.5 1.5) for line in line1.union(line2): print line LINESTRING (0 ...


1

There are two different ways that Tiles can be addressed, the TMS standard specifies that tile coordinates start at the bottom left, but in practice most software is using a coordinate system with the Y axis reversed from TMS: it starts at the top left instead of the bottom right. I don't know that there is really an official name for this scheme, some ...


1

You can use gdal_merge which will combine the rasters. Have a look at the description of gdal_merge: This utility will automatically mosaic a set of images. All the images must be in the same coordinate system and have a matching number of bands, but they may be overlapping, and at different resolutions. In areas of overlap, the last image will be copied ...


1

I believe what you are experiencing is more or less a copy of this question. The coordinates in the rainfall data are in longitude/latitude, but with values ranging from 0 to 360, instead of -180 to 180 (as your political boundaries are). See the GPCC spatial note here (emphasis mine): Spatial Coverage: 0.5 degree latitude x 0.5 degree ...


1

After all I found solution, but not sure if that is right and elegant. I used proj4 string that SeaWifs level 3 have: s_srs="+proj=eqc +lat_ts=0 +lat_0=0 +lon_0=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs" and overwrite georeferenced bounds by the values from SeaWifs hdf file a_ullr =c(-20037508.343,10018754.171,20037508.343,-10018754.171) after ...


1

Gdal PYthon binding works with conjuration with gdal. You need to have installed GDAL along with the binding to work. A quick skim over the errors you provided they suggest that while you installed gdal-python you haven't install the gdal by itself: A possible solution chain is: Install GDAL (http://www.gisinternals.com/sdk/) Append the installation ...


1

Couchbase is not Apache CouchDB compatible, hence the OGR driver doesn't work. AFAIK there's currently no OGR driver for Couchbase.


1

you can loop on your files and append them to your list command = [sys.executable,gmerge,'-o','C:\\r.tif','-of','GTiff'] images = glob.glob("D:\\*.tif") for image in images: command.append(image) subprocess.call(command)


1

Your output image will have more pixels than the sum of your input images, but this does not explain the large difference. I suggest that you look at the characteristics of your images based on gdalinfo in order to see what compression is used and check that the extents are correct. (assuming your input images have the same size, it makes 20000 * 12000 ...


1

From the CSV driver documentation: Starting with GDAL 1.8.0, for files structured as CSV, but not ending with .CSV extension, the 'CSV:' prefix can be added before the filename to force loading by the CSV driver. Either rename DGM5_BE.txt to DGM5_BE.csv or change the <SrcDataSource> element to: ...


1

No, as JPG does not support the GDAL create option but only GDAL create copy: you can check it using this code: from osgeo import gdal format = "Jpeg" driver = gdal.GetDriverByName( format ) metadata = driver.GetMetadata() if metadata.has_key(gdal.DCAP_CREATE) and metadata[gdal.DCAP_CREATE] == 'YES': print 'Driver %s supports Create() method.' % format ...


1

Use EPSG:3857. That is the official EPSG code number for Google mercator projection. 900913 was incorproated in GDAL some toime ago, but is now dropped in favour of 3857.


1

If you're scared or unsure about using command line tools like GDAL_Translate then there are other ways to do this. If you have Esri products then you can use Raster to Other Format or Copy Raster, you can even export the image from ArcMap just by right clicking on the layer and select export then fill in the blanks and select GeoTiff. If you want to use ...


1

Manage to solve my problem. THis can be done by instantiating OGR and openup a shapefile database, layers etc Afterwhich, create a OGRPoint object with the neccessary lat lon input OGRPoint* pt = new OGRPoint(lat, lon); set a spatial filter based on this point object on the layer next extract the feature based on the filter resu;lt Finally extract the ...



Only top voted, non community-wiki answers of a minimum length are eligible