I have multiple shapefiles, and I want to import them as image tiles.  How do I convert a shapefile into image format and then convert that single image file into an image pyramid such as that of OpenStreetMap (OSM)?  I have worked with GeoServer for backend and OpenLayers to render the map but now I am asking about converting my own custom shapefiles into OSM tile pyramid.
3177
I want to start working on a web map at work to map some of our statistics by city.    I have an idea, but I'm not exactly sure where to start.   Are there any definitive resources that provide step by step guidance on how to do this?
119192
So I'm working on a project that's been in development for over a year, but I only recently got involved. My task involves converting roughly 2,500 KML files into another format, importing these files into a Geoserver, and then copying over the associated style information for each file. I realize that Geoserver allows you to import KML files directly, however, our files never seem to render correctly, and my boss would like to move away from KML format for philosophical reasons.   Using ogr2ogr, I've tried to convert these files into ESRI Shapefiles, however the conversion process only worked on roughly a third of the files and erased some of the metadata called for in the style information. The files that didn't convert code for polygons and points, which I've read can be an error for other programs like QGIS.  I've also tried to put all of the KML files into a PostGIS database and convert them like that. However, that seems to delete all of the data leaving essentially blank tables for import.  TL;DR I need to find a way to convert 2,500 KML files into a format usable by Geoserver, and will work with XML style information. I'm a complete GIS noob and this is my first time posting, so I'm sorry if this seems a bit ridiculous.
72233
I have an image of size 1GB (.tif), with the width and height 94000x71680. I would like to chunk this image into 20000X20000 tiles so that I can process them.   How can I do this?
25789
I'd like to generate a 3d terrain model moving from a DEM dataset (currently in shp and GeoTIFF format).   As far as I understand it, the first step of the better path would be generating a TIN, but I can't find the way to do.   I found a lot of references to ArcGIS 3D Analyst, but no tools in MacOSX/Unix world.   Are there any?
81329
I've been working with loads of shapefiles projected according to OSGB 1936 (CRS 27700). The project CRS is set to EPSG:27700 and "On the fly reprojection" to EPSG:27700 is set (I keep this as a default as I work with OSGB 1936 95% of my time).  I decided to download OSM layers from the OpenLayers plugin. When I load this, the EPSG in the bottom right of the screen goes from EPSG:27700 to EPSG:3857   Question 1: should this be happening anyway as it would suggest that the whole project has been re-projected to EPSG:3857?   As soon as the OSM layer is loaded it's way out of alignment but a quick pan seems to bring it altogether quiet nicely.   What I have noted is that (Question 2), if I try to re-specify the CRS to 27700, then the shapefiles "move" by a few tens of meters with respect to the OSM. Does this suggest that QGIS is doing a second re-projection (from what to what?)?   Also (Question 3), does it mean that if I import another shapefile that would be originally a dataset digitised with EPSG:27700 it would be re-projected to EPSG:3857 because the CRS in the bottom right of the screen was set to EPSG:3857 after the OSM import. This is all a bit confusing and logic would suggest that the import of a non-ESPG:27700 (project default CRS) should not go and change the EPSG setting in the bottom right of the screen. Anita Graser's book says on page 34 "Your project CRS will automatically be changed to Pseudo Mercator when you load a background map using the OpenLayers plugin option".   Does that mean there is no way of keeping ESPG:27700 active in the project when importing basemaps from OpenLayers? If so, what alternatives are there?  What I am trying to bottom out is to understand the correct way of importing basemaps (usually Google, Bing or OSM through the OpenLayers plugin) into my projects (always EPSG:27700) without losing alignment and getting layers to correctly re-project. I'm trying to support greater use of QGIS in my company but don't want to lose credibility by promoting the use of a tool that is misaligning datasets. Don't get me wrong. QGIS is a great tool and I am very fond of it but importing layers with different CRS's is quite common these days and getting this right from the start is half the job done.  I did quite an extensive Google search before posting and didn't come up with a recent, clear cut answer to my questions above. There is a bit of content out there on problems with older versions of QGIS, PROJ4, etc. but nothing recent for QGIS 2.0 so worth an update.
95660
Is it possible to load a raster file (e.g. geotiff or erdas img) directly into a Leaflet map? If so, can you guide me to the correct documentation or method to do this? I am trying to do this with the imageOverlay() function in Leaflet currently.   var imageUrl = 'http://localhost/test/my_raster.tif',         imageBounds = [[30.430224, -87.193133], [30.435867, -87.185944]];         L.imageOverlay(imageUrl, imageBounds).addTo(map);    I am new to Leaflet.
75830
Alright, I'm toying around with converting a PIL image object back and forth to a numpy array so I can do some faster pixel by pixel transformations than PIL's PixelAccess object would allow.  I've figured out how to place the pixel information in a useful 3D numpy array by way of:  pic = Image.open("foo.jpg") pix = numpy.array(pic.getdata()).reshape(pic.size[0], pic.size[1], 3)   But I can't seem to figure out how to load it back into the PIL object after I've done all my awesome transforms.  I'm aware of the putdata() method, but can't quite seem to get it to behave.
23399
I am trying to export images from the QGIS print composer, the exported image has a basemap that is loaded from an XML file following . The export doesn't look at all like the image I see in the print composer see below and I'm trying to understand why.   What I'm trying to do is simply make a decent looking map without obvious tiling with readable labels. Perhaps I have to use a satellite only image and remake with my own labels. The map is a HERE nokia hybrid map   This is a jpeg at 150 dpi, note the tiny labels in comparison to the screenshot at the bottom .   This is a png export at 300 dpi, which is standard for the print composer, note the labels, and tiles.    Finally, this is a screenshot of the map as it looks in the print composer with the labels readable at this zoom scale.
57359
How would one save a WMS view to a local raster from the python command line?  I've seen the answer to  and that works for me but I'd like to know how to replicate that process exactly from the console.  ie How to progamatically set the layer, the extent to save, and then save that extent for that layer to a tif?  All of the searches I did threw up lots of results of people suggesting that downloading via wfs was the best solution. Or the linked solution via the gui. But none pointed me to the OWSLib solution.    Here is the working code I made on the back of the answer from @nmtoken I got.  from owslib.wms import WebMapService from osgeo import gdal, osr import processing  url='http://&lt;mapserverurl&gt;/MapServer/WMSServer?service=WMS' folder = '~/output' wms = WebMapService(url, version='1.1.1') wms_layers= list(wms.contents)   for layerno in range(0,len(wms_layers)):     raw_tiff=os.path.join(folder,'_raw'+str(layerno)+'.tif')     georeferenced_tiff=os.path.join(folder,'_georeferenced'+str(layerno)+'.tif')     defaultepsg=4326      #Note that available epsgs can be determined with:     print wms[wms_layers[0]].crsOptions     srs_string='EPSG:'+str(defaultepsg)     lonmin=1     lonmax=2     latmin=4     latmax=5     Xpix=400     Ypix=300     #Note output options can be determined with:     print wms.getOperationByName('GetMap').formatOptions     img = wms.getmap( layers=[wms_layers[layerno]], styles=['default'], srs=srs_string, bbox=( lonmin,latmin , lonmax,latmax), size=(Xpix, Ypix), format='image/tiff' )     out = open(raw_tiff, 'wb')     out.write(img.read())     out.close()     out = None      #Georeferencing     src_ds = gdal.Open(raw_tiff)     format = "GTiff"     driver = gdal.GetDriverByName(format)     dst_ds = driver.CreateCopy(georeferenced_tiff, src_ds, 0)     gt = [lonmin, (lonmax-lonmin)/Xpix, 0, latmax, 0, -(latmax-latmin)/Ypix]     dst_ds.SetGeoTransform(gt)     srs = osr.SpatialReference()     srs.ImportFromEPSG(defaultepsg)     dest_wkt = srs.ExportToWkt()     dst_ds.SetProjection(dest_wkt)     dst_ds = None     src_ds = None
60420
I am trying to convert a KML file to Shapefile using GDAL (ogr2ogr). Regular KML files convert fine but KML files containing MultiGeometry features (eg: Geometry containing both Point and Polygon) do not get transformed.   The output is simply a blank shapefile containing no shapes or attributes.  I have tried this command:  ogr2ogr -overwrite -f "ESRI Shapefile" -where OGR_GEOMETRY='MultiGeometry' &lt;ouputfolder&gt; &lt;input kml&gt;   For an example, KML file which has MultiGeometry, you can see the us_states.kml hosted by Google.    Other info: I am using 'GDAL 1.10.1, released 2013/08/26' for windows.
29700
I'm looking for a shapefile containing an 1000m or 2000m UTM grid ("graticule"). All I was about to find on the web is a 10000m UTM grid (here: ) and I am aware that a 1000m UTM grid for the whole world might be quite big.  That's why I'm now wondering if and how it would be possible to generate such a shapefile for a part of the world. I would prefer to use open source tools like gdal. Others output formats than a shapefile might be ok too.   Has anyone an idea?
48371
I need to identify the roofs of buildings using a satellite image, something like what is available in OSM  (something like the image).  Is there anyway to do this in Qgis, R, Arcmap?
32557
I have tried reading up on WMS on wiki and other sources but haven't managed to get my head around a couple of points.   1) Using WMS is it possible to get a local copy of the raster data (say a tiff image tile). The reason, I ask is, I need to use his raster data in a secondary process that needs a local copy of the file.   2) I am primarily interested in allowing users the option of downloading orthoimagery and elevation from the USGS repository. I've tried using the Seamless viewer but I am thinking more along the lines of user specifying a coordinate (or four corner coords) and getting the orthoimage and elevation data corresponding to that. I thought WMS (Geoserver being my first option) is the way to go for this, but not certain. So wanted to get your advice.  Thanks in advance.  Sincerely, ST.
28218
I am putting together a GIS server using GeoServer and need some help with importing my layer styles from ArcGIS 10 to  format.    I have one layer in ArcMap properly symbolized and labeled, how can I get this layers symbolization and labeling into a SLD file?
47330
Is there a tool like  available for Android or Windows systems? I'd like to  turn images into polygones/triangles like shown on the app site.
109294
What I want to do is to fill a polygon (country) with a graphic file in any format (jpg., bmp., png.).  It should look like this:     I'm aware that GIS soft is not proper for such a operations but maybe is there any option to do it only in QGIS.
61183
I've made a noob mistake and set up a whole scene with empty image objects.  Is it possible to convert them so they would show in the render? I know I can delete everything and start importing everything from scratch as planes, but that would take too much time to start all over again.