地理信息系统

我正在用Python开发，并使用OSGEO的GDAL来处理栅格和shapefile并与之交互。 我想获取一个具有点特征的shapefile并将其插入到表面栅格中。现在，我使用的是'RasterizeLayer'方法，该方法将点要素中的值刻录到栅格中（使用所有nodata值进行设置），但将所有未触及的像素保留为'nodata'值。因此，我剩下的是棋盘式栅格。 使用RasterizeLayer之后，我得到了什么： 我想要的最终产品是： 我相信我要寻找的功能在arcgisscripting导入中称为'Spline_sa（）'。 GDAL是否具有类似的功能，或者有其他方法来获得所需的输出？

20 python  raster  gdal  interpolation 

我们已经下载了本地区的OSM数据，并使用osm2pgsql将其加载到PostGIS中。到目前为止，一切都很好。最上面的是GeoServer，我们想使用WMS公开OSM数据。我们希望像“标准” OpenStreetMaps（例如osm.org）中那样对数据进行样式设置，因为这种制图方法设计得很好。我们试图下载用于OSM数据的som标准SLD，但是使用这些SLD似乎存在各种问题（XML解析错误，无效内容等）。 有什么好主意如何重用OSM SLD通过WMS设置OSM样式？

20 openstreetmap  sld 

给定大量（约100万个）不均匀分布点的样本-是否可以生成包含指定最小n点数量的不规则网格（大小，但形状是否可能不规则？）？ 对于我来说，如果这样的网格的“细胞”精确地包含n个点或至少n个点，那么对我来说就不那么重要了。 我知道诸如ArcGIS中的genvecgrid或QGIS / mmgis中的Create Grid Layer之类的解决方案，但是它们都会创建常规网格，从而导致输出的单元格为空（较小的问题-我可以简单地丢弃它们）或具有点数的单元格小于n（更大的问题，因为我需要一种解决方案来聚合这些单元格，可能使用此处的一些工具？）。 我一直闲逛无济于事，并且对商业（ArcGIS和扩展）或免费（Python，PostGIS，R）解决方案都开放。

20 r  clustering  quadtree 

我想将海拔标签定向到山上。 那可能吗？ 我目前正在使用QGIS的标签工具（不是图层设置中的那个）。我希望标签弯曲并面向上坡（并且我想定义数字显示的频率）。全部用于打印地图。

20 qgis  labeling  contour 

在ArcGIS 10和Python中，我想获取shapefile中每个多边形的范围（xmax，ymax，xmin，ymin）信息。 我可以使用获取整个shapefile的范围 file=r"D:\SCRATCH\ARCGIS\100k_trc_tiles_TVM.shp" desc=arcpy.Describe(file) print desc.extent.Xmax 394551.52085039532 但是我似乎无法弄清楚如何为数据集中的每一行获取相同的信息。 rows = arcpy.SearchCursor("100k_trc_tiles_TVM") for row in rows: print row 打印数据集中的31行，但是 for row in rows: desc=arcpy.Describe(row) print desc.extent.Xmax 给出一个错误。 运行时错误：对象：描述输入值无效的类型 我当时正在考虑使用“计算几何”将范围值添加到表中，但这仅给出了质心。然后，我猜我们可以使用诸如row.GetValue（“ xmax”）之类的东西。 话虽这么说，我知道我们可以使用http://www.ian-ko.com/free/free_arcgis.htm中的函数创建X / Y，max / min，但是最好避免添加字段，尤其是在ArcPy可以获取这些值的情况下。 基本上，我需要获取范围以输入到裁剪工具中，以裁剪出30个数据区域（根据1：100,000地图图）进行地理处理，因为拆分工具由于数据集的大尺寸而失败（请参见为什么相交会给错误999999：错误执行功能无效拓扑[太多lineseg端点]）。我要自动执行此操作，因为它在许多数据集上重复出现。 ===工作脚本=== # Emulates Arc Info SPLIT tool by using Clip but # Requires a …

20 arcgis-desktop  arcpy  arcgis-10.0 

您将如何从建筑脚印生成建筑阴影？ 我有建筑物高度的属性值。 因此，最好使阴影的大小基于此而变化。 使用ArcGIS找到的最佳建议是阴影：http : //blogs.esri.com/Support/blogs/mappingcenter/archive/2008/10/21/3d-effect-for-a-building-footprints.aspx

20 arcgis-desktop  arcmap  polygon-creation  shadow 

我已经在Ubuntu 11.10上安装了PostGIS 1.5.3的PostgreSQL 9.1。 按照此说明，我还安装了pgadmin3 1.14.0-beta1 。 如果我没记错的话，当我在Windows上安装PostGIS时，pgadmin自动安装了Shapefile GUI Loader插件。但是，在Ubuntu上，插件下拉菜单为空。 我发现的说明似乎是Windows特定的。 有没有办法在Ubuntu上安装/启用它？

20 postgis  postgresql  ubuntu  pgadmin-3 

当您用谷歌搜索“纬度/经度”时，您获得的结果比“经度/纬度”高10倍。我觉得这很令人困惑，因为“ x / y”的普遍性似乎是“ y / x”的20倍。但是在地图上，纬度在Y轴上，经度在X轴上。 也许我只是在发泄，但也许我错过了使所有这些疯狂变得合乎逻辑的线索。 有什么解释吗？

20 latitude-longitude 

我们是.NET开发人员的小团队。我们拥有丰富的GIS经验，而且我们都不是软件/数据库开发或系统管理的新手。我们拥有技术学位和多年行业经验。我们参加了Esri开发者峰会。 Esri的技术-主要是ArcGIS Server，ArcSDE和ArcObjects-在我们开发的所有软件中起着很小但必不可少的作用。尽管ESRI在我们的技术堆栈中仅占少数地位，但我们仍花费大量时间来解决难以捉摸的错误，制定变通办法，破译其模糊的错误消息，跟踪性能问题以及回收流程。 通常，我们的问题源于真正的错误，不良的异常处理，有限的设计/架构决策，缺乏文档，不稳定或它们的某种组合。（我在这里说的是ESRI堆栈。） 从项目经理的角度来看，我非常担心团队的生产力。这花费了我们很多时间。我们没有时间学习ESRI堆栈的每个特质，但是我们仍然需要完成工作。（不能忍受它，不能没有它。） 对于混合使用ESRI的开发人员，您有什么实用的建议？ 我不是在寻找有关替代技术堆栈的建议。

20 arcgis-platform 

是否有任何教程，讲座或演示文稿讨论使用EC2或其他云计算服务托管GIS服务器？有没有人在使用Amazon EC2上有过任何经验？设置有多难？我知道ESRI的服务器产品有一些版本专门针对EC2，有人使用过它们吗？ 非常感谢你。

20 arcgis-server  postgis  cloud-gis  ec2  remote-hosting 

在PostGIS文档中，它说使用SQL创建空间表有两个步骤： 创建一个普通的非空间表。 使用OpenGIS“ AddGeometryColumn”功能将空间列添加到表中。 如果遵循这些示例，则将创建一个terrain_points像这样的表： CREATE TABLE terrain_points ( ogc_fid serial NOT NULL, elevation double precision, ); SELECT AddGeometryColumn('terrain_points', 'wkb_geometry', 3725, 'POINT', 3 ); 另外，如果我查看pgAdmin III中的现有表，似乎我可以像这样创建相同的表： CREATE TABLE terrain_points ( ogc_fid serial NOT NULL, wkb_geometry geometry, elevation double precision, CONSTRAINT terrain_points_pk PRIMARY KEY (ogc_fid), CONSTRAINT enforce_dims_wkb_geometry CHECK (st_ndims(wkb_geometry) = 3), …

20 postgis  sql  postgresql 

我正在尝试使用Python自动执行ArcGIS Desktop（通常使用ArcMap）中的各种任务，并且一直需要一种将shapefile添加到当前地图的方法。（然后做一些事情，但这是另一个故事）。 到目前为止，我能做的最好的就是使用以下命令将图层文件添加到当前地图（“ addLayer”是图层文件对象）： def AddLayerFromLayerFile(addLayer): import arcpy mxd = arcpy.mapping.MapDocument("CURRENT") df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0] arcpy.mapping.AddLayer(df, addLayer, "AUTO_ARRANGE") arcpy.RefreshActiveView() arcpy.RefreshTOC() del mxd, df, addLayer 但是，我的原始数据始终是shapefile，因此我需要能够打开它们。（等效地：无需打开即可将shapefile转换为图层文件，但我不希望这样做）。

20 arcpy  shapefile  layers 

我在（100，100）周围有一层矢量，我想将它们移动到新坐标（1000，1000）。 我应该如何使用QGIS？

20 qgis  vector  coordinates  coordinate-system  affine-transformation 

我努力说服我组织中的其他人，GIS工具（特别是数字化工具）可以为大众使用，我们应该促进其使用，而不是使用纸质地图。这样可以节省金钱（邮政费）和时间（张贴地图所需的时间），我想这样对环境更好（不使用纸张）。 相反的论点是，人们用笔/铅笔在纸质地图上感觉更自在，并在纸质地图上勾画出区域/多边形等，它们也往往更准确。我的业务要求申请人在绘制背景图边界时必须相当准确。 您认为我该如何说服那些不信的人？:)我们如何为非GIS用户改进草图绘制工具，以便他们在地图上绘制准确的草图。您知道哪些技巧/技术，您自己曾经遇到过这个问题吗？我们是否注定要再使用纸质地图！欢迎您的想法。 PS-GIS应用程序可能是基于Web的，并向公众开放。可能的技术是ESRI Javascript API-具有缓存，优化的地图服务，合理的比例级别和清晰的符号等。

20 gui  digitizing 

是否有任何来源可用于美国主要城市的免费和可分配建筑足迹数据集？ 理想情况下，数据还应包含每个建筑物的楼层数。那也可以吗？ 我将给您一个我想做的事情的想法：我想根据建立平方英尺数和白天/晚上的人口密度来估计地图区域内的人口。

20 data  united-states