地理信息系统

制图师,地理学家和GIS专业人员的问答



2
QGIS 3继续创建“调色板”无用文件夹
全新安装的QGIS 3.0.x在Windows的c:/中系统地创建一个“ palettes”文件夹。即使您删除它,它也会在每次重新启动/重新安装QGIS后重新创建。 一劳永逸地摆脱此文件夹的技巧吗?或至少要更改其创建路径?
13 qgis  qgis-3 

3
如何在QGIS 3.0中编辑GDAL / OGR控制台调用?
在QGIS 2.XX中,下拉菜单栏工具的底部有一个可编辑的呼叫框,可以在其中手动输入gui中不可用的命令(单击铅笔按钮)。在3.0中,菜单栏工具与处理工具箱中的GDAL工具相同,这更有意义,但是现在我看不到向工具操作添加简单命令的方法。我附带的示例是该warp (reproject)工具,我希望添加命令“ -wo CUTLINE_ALL_TOUCHED=TRUE”。我怎么做?
13 gdalwarp  qgis-3 

4
使用Shapely将单个线段焊接到一个LineString中
我在python中使用Shapely,并给了MultiLineString它一堆Linestring对象。我可以保证所有LineString对象都是只有2个顶点的简单线,并且它们都是一条直线的一部分(无分支)。 我想“连接点”并创建一个LineString。我是否需要为此编写递归焊接方法,还是有更快的方法?
13 python  shapely 

2
如何将geoJSON数据导入MapBox
在尝试将可编辑数据导入mapbox时,我无法以系统接受的格式导入数据。 它拒绝该数据并显示一条错误消息,指出: city.geojson:输入失败。第1行不建议使用老式的crs成员。 在阅读时,我认为文件转换一定有问题,因此我将其与提供给MapBox 的示例文件进行了比较,而没有发现任何差异。 接下来是最令人困惑的部分:尝试导入示例文件时,我从系统中收到非常相似的错误消息: stations.geojson:输入失败。不建议使用老式的crs成员,>此对象等效于默认对象,应在第1行上将其删除。 现在我的问题是:如何将geoJSON文件导入MapBox?
13 geojson  mapbox 


3
用于逻辑运算符和其他功能的gdal_calc栅格计算器语法
在gdal_calc的文档中,使用numpy语法声明了命令行栅格计算器。稍后有几个示例,其中一个示例: gdal_calc.py -A input.tif --outfile = result.tif --calc =“ A *(A> 0)” --NoDataValue = 0- 表示将零以下的值设置为null 不幸的是,没有逻辑运算符的示例,例如: --calc =“ A *(A> 0和A> B)”- 表示如果A大于零且B较大,则保留A,并将其余值设置为null 基于Numpy / Scipy逻辑函数,我希望将逻辑运算符写为: --calc =“ A * logical_and(A> 0,A> B)” 我尝试了一下,它似乎可以工作,但是我想向您保证这是正确的。 以类似的方式,如果您想要最少的A和B: --calc =“ A *(A <= B)+ B *(A> B)” 您可以这样写: --calc =“ minimum(A,B)” 我的问题是我找不到任何食谱来确保我做对了。是否有一些不错的食谱,其中包含有关gdal_calc什么是和不可能的高级示例?
13 python  gdal  numpy  gdal-calc 

1
Mapbox错误:样式未完成加载
我正在尝试基于suncalc结果并使用geojson在mapbox地图上绘制。首先,我尝试创建2个函数,对于要绘制的每一行,一个函数。但是当我这样做时,它只会显示最后一个调用的函数。所以我得出结论,因为我不熟悉mapbox和Leaflet语法,所以我不知道如何管理图层。 使用传单示例,我得到了以下代码: JAVASCRIPT function drawOnMap(sAz){ //sun calculation stuff, doesn't matter for the problem var lt = parseFloat(document.getElementById('lat').value); var ln = parseFloat(document.getElementById('long').value); var R = 6371000; //Earth's radius var d = 10000* Math.sin(1); // Distance var brngSR = deg2rad(sAz); var ltr = deg2rad(lt); var latSR = rad2deg(Math.asin( Math.sin(ltr)*Math.cos(d/R) + Math.cos(ltr)*Math.sin(d/R)*Math.cos(brngSR) )); …

2
有没有一种方法可以为QGIS中的按列分配表的行号?
字段计算器在“记录”下有一个很棒的功能,称为“ $ rownum”,它返回(分配)当前行的编号。 有没有办法使其与排序一起使用?可以说我想按“长度”列对表格进行排序,并按从最短到最长的长度获得行的顺序(顺序)。 我知道我可以轻松地使用一些Office电子表格软件来对此进行按列排序,并用顺序(顺序)填充新列。但是我更愿意在QGIS中完成整个工作。 更新:作为解决方法,还有如何更改shapefile中特征的顺序?-使用mmqgis插件和功能进行修改/排序,然后在新保存的文件上使用“ $ rownum”。

1
在R中创建给定栅格的范围和边界的多边形
我想创建两个多边形。 一是我的栅格的矩形范围。我知道我可以使用来获得栅格的范围 r <- raster("band5.tif") e <- extent(r) plot(e): 如何创建栅格边界的一个多边形,如下所示? 这是我的栅格的样子:
13 raster  r  rgdal 


3
geopandas空间连接极慢
我正在使用下面的代码来查找数百万个GPS点的国家(有时是州)。该代码当前每点大约需要一秒钟,这非常慢。shapefile为6 MB。 我读到geopandas使用rtree进行空间连接,这使它们效率极高,但这在这里似乎不起作用。我究竟做错了什么?我希望每秒能获得一千个积分。 可以在此处(5MB)下载shapefile和csv:https ://www.dropbox.com/s/gdkxtpqupj0sidm/SpatialJoin.zip ? dl =0 import pandas as pd import geopandas as gpd from geopandas import GeoDataFrame, read_file from geopandas.tools import sjoin from shapely.geometry import Point, mapping,shape import time #parameters shapefile="K:/.../Shapefiles/Used/World.shp" df=pd.read_csv("K:/.../output2.csv",index_col=None,nrows=20)# Limit to 20 rows for testing if __name__=="__main__": start=time.time() df['geometry'] = df.apply(lambda z: Point(z.Longitude, z.Latitude), …

3
QGIS Print Composer不会从横向更改为纵向吗?
使用2.10.1,我打开一个新的Print Composer,然后立即将页面设置方向从横向更改为纵向。但是屏幕的方向不会改变。它继续显示风景。 同样,更改“纸张尺寸”不会影响“打印编辑器”窗口的尺寸/方向。 如果我将地图添加到Print Composer并导出为PDF(无论页面设置设为什么),结果始终是横向的。 看来Print Composer仅能使用A4横向。 我想念什么?

2
Leaflet的免费磁贴提供商
我将Leaflet用于我的网站地图。到目前为止,我使用的是Mapbox地图(很酷),但显然,他们每月为拥有100,000张以上地图的网站收取费用。 我正在尝试寻找瓷砖供应商,这些瓷砖的供应商不收取任何费用。最初,我使用OpenCycleMap进行此操作,并替换为 https:// {s} .tiles.mapbox.com / v3 / {id} / {z} / {x} / {y} .png 与 http:// {s} .tile.opencyclemap.org / cycle / {z} / {x} / {y} .png 两个问题: 1)OpenCycleMap是否提供每月没有地图限制的免费图块? 2)有人认识其他免费提供者吗?

By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.