Hot answers tagged

16

For QGIS 3.8 or greater, generating tiles is built in - you can find it in the Processing tool box as Generate XYZ tiles (Directory) and Generate XYZ tiles (MBTiles).


13

SRTM tiles (3601 px * 3601 px in this case) have 1-pixel overlaps in between. When you apply transparency (or reduced opacity) to your hillshade layer, such overlapping pixels stand out. You may have observed this also on the original images, if you apply transparency (see below). ..... original SRTM, Pseudo-color + 60% opacity Anyway, you can avoid ...


7

You're asking two questions here: administrative borders related to territorial disputes, and map label localization. Since you mention OpenStreetMap, I feel obliged to answer: The OpenStreetMap Foundation changed the policy regarding territorial disputes in 2013. You can find the full document here, but let me quote the relevant paragraph: Borders and ...


4

As described by @Kazuhito, the tiles should be combined into a virtual raster. Their method is performed within QGIS. If you have a large number of tiles you may prefer to do this from the commandline. To do so open a shell in the directory containing the SRTM tiles and run gdalbuildvrt combined.vrt *.hgt The resulting virtual raster (combined.vrt) can ...


4

The layer info at the link you provide indicates the MapServer should provide tiles up to level 19, not 20 as you'd like. But that shouldn't create problems at level 17! However, I've played around with zooming in on the layer the QGis map window and it shows up at scales up to 1:6019, but not more detailed than that. So there's something funny with the ...


4

You can use UNION to combine the queries. Your server can then concat the binaries before sending them out. I found out that doing this in SQL took to long and increased the compute time to much. You can stack the UNION queries for more objects. Here is an example of the query I used: ( SELECT 'road' as name, ST_AsMVT(q, 'road', 4096, 'geom') as tile FROM ...


4

It works for me if I add a new XYZ Tiles service in the Browser panel. The URL is https://nls.tileserver.com/nls/{z}/{x}/{y}.jpg The tile layer uses EPSG:3857 but if your project CRS is 27700 then QGIS will reproject it on the fly. But be careful of trying to draw the whole world with a project CRS like 27700 because its not valid over the whole globe and ...


2

GDAL can create vector tiles with the MVT driver https://gdal.org/drivers/vector/mvt.html For converting GeoJSON data into vector tiles which are saved into MBTiles database file use a command like ogr2ogr -f MVT -dsco FORMAT=MBTILES -dsco MAXZOOM=10 target.mbtiles source.geojson Another option is to use Tippecanoe the https://github.com/mapbox/tippecanoe....


2

Unfortunately, TileLayerPlugin is compatible with QGIS versions 2.8 - 2.99 So, the other option to add custom tiles stored locally on disk - fill in URL (I use Linux, for Windows I suppose you have to change backslashes): file:///path/to/your/custom/tiles/{z}/{x}/{y}.png


2

As @Houska has noted the map service that you are using only provides tiles down to level 17. However, the data supplied in the map service is the USGS National Hydrography Dataset which is freely available to download in vector form, so if you want to generate tiles at any zoom level, this can be done by using a styled version of the vector NHD data.


2

So I just did a bit of Googling and have managed to get the XYZ URL to render in QGIS. However, I think there is a problem with the coordinate system. So long as you are making your own shapefiles, this shouldn't be a problem. If you're using imported shapefiles, it might be possible to manually move them into place within the Web Mercator projection because ...


2

Raster tiles tend to use a regular grid and each tile on certain zoom level covers as large area. If data is dense then the tile will be bigger as compressed. As uncompressed you would get tiles with uniform size but they would be all big. I fear you are out of luck with current slippy map tiling systems. You are after something like r-tree https://en....


2

Probably the simplest solution for this is to use TMS tiling schema in SAS Planet, since it's natively supported in Leaflet. In SAS Planet you go to Settings > Options > Cache. Then in the Default cache type field select Tile Map service (TMS). The you can declare your TMS layer simply with: var tmsLayer = L.tileLayer('path to cache files/{z}/{x}/{y}.jpg',...


2

gprojector Software, its great for Raster https://www.giss.nasa.gov/tools/gprojector/


2

Advanced projection handling is supposed to take care of the issue (indeed it gets worse when you disable it), but it's really just a set of heuristics. It's likely not accounting for some specific issue in your dataset. You should open a bug report attaching the sample data set, style, and requests to make the issue happen (ideally a preview request of ...


2

You have to set refresh_before or cleanups : remove_before: setting in your seed.yaml. Without them new tiles will not be created.


2

After weeks of research and headache, I finally figured out the correct architecture of the local OSM server. My problem was that the employee who created the local OSM server hasn't actually finished his work (all he has done is cloning the open source OSM project and deploying it on the local server). Probably, the following information won't be helpful ...


2

You are looking for an image pyramid layer, there is a full tutorial in the manual as well.


2

Sentinel 1 products do not use a a tile grid. However, to identify products over your polygon you can use one of several Sentinel 1 search/download resouces. This is not an exhaustive list by any means. ESA offer the SciHub which allows searching by drawing a polygon area of interest (AOI). You have to create an account (free) before you can search and ...


1

Please note that the current recommended method for defining the tiles argument is as follows: # Define a method for displaying Earth Engine image tiles to folium map. def add_ee_layer(self, eeImageObject, visParams, name): map_id_dict = ee.Image(eeImageObject).getMapId(visParams) folium.raster_layers.TileLayer( tiles = map_id_dict['tile_fetcher']....


1

Your matrix is almost certainly not called just 6, if you look in the GetCapabilities response from GeoServer you will see a list of available matrix sets like this: <TileMatrixSetLink> <TileMatrixSet>EPSG:4326</TileMatrixSet> <TileMatrixSetLimits> <TileMatrixLimits> <TileMatrix>EPSG:4326:0</TileMatrix> <...


1

To use this function you need to register (for free) at USGS. Go to https://earthexplorer.usgs.gov/ and log in. Now draw a bounding box of your desired area or upload a shapfile or KML. Note that this is limited to 30 vertex: Add additional criteria such as a time range or the number of visible results (you can download unlimited, see below) if you want ...


1

You can write two .VRT files and use them chained: Create a .VRT that mosaics all your source files into one virtual layer with gdalbuiltvrt https://www.gdal.org/gdalbuildvrt.html Make a new vrt that is referencing this Mosaic as gdal_translate -of vrt mosaic.vrt filtered.vrt Edit filtered.vrt and turn it into KernelFilteredSource as in the answer How to ...


1

One possibility would be to use Leaflet's L.CRS.Simple flat (Cartesian) CRS. Here one pixel corresponds to one map unit. Let's say we want map units that correspond to original rasters map pixels (before it was cut to tiles): one map unit is one pixel on original raster. With tile layer at zoom level 0 whole map is within 256 x 256 pixels, so factor for ...


1

OK, I found the solution. Firstly change tile's src by calling map.getLayers().getArray()[0].getSource().getTile(2,1,-3).image_.src="base64 or new url" Then call changed() function on the tile map.getLayers().getArray()[0].getSource().getTile(2,1,-3).changed() After all to refresh this tile call changed() function on the source object map.getLayers()....


1

The second argument to addLayer is the ID of the later to insert before. That way you can control what appears on top see https://docs.mapbox.com/mapbox-gl-js/example/wms/ Though I think what you'll find better is just creating your own style object when you initialise the map with your raster tiles. Then any other layers you add after will go on top. See ...


1

I have two suggestions Client side My experience has been with leaflet. In the leaflet API you have a choice of the client adding as Vector layers or Raster Layers. It sounds like you are currently using the equivalent of Vector layers. which creates an HTML element for each item you are mapping. Raster layers do the equivalent of a canvas and builds a ...


1

Keep in mind that the plan for Landsat-8 is only for a period of 3 days at most while for Sentinel-2 is a bit better with 10-15 days, so it's difficult to make an accurate prediction beyond that. A bit old question but if someone still interested in the acquisition plan for Sentinel-2 and Landsat-8 you can check planned acquisitions for an AOI through ...


1

Dezoomify You can use dezoomify to download the full image. Instructions Take the original tile URL. In your case: http://fotogrammi.regione.abruzzo.it/ecwp/ImageX.dll?image?cache=true&transparent=true&type=jpg&l=0&tx=9&ty=9&ts=256&fill=FFFFFF&quality=60&layers=/VIGM54/1/48/1686.ecw&srs=EPSG:RAW Replace the values ...


1

It is a generic error but here is good Esri Technical Support article on “ERROR 001369: Failed to create the service” that gives a clue into possible problems/solutions. I would recommend you to go in the servers logs and try to identify the specific error report there.


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