地理信息系统

下面以平方度给我了多边形的面积： SELECT ST_Area(ST_GeomFromText('POLYGON((871325.790874952 6105405.3261047,871418.748307692 6105359.72944624,871346.22022442 6105215.141258,871254.85408906 6105261.72007212,871325.790874952 6105405.3261047))',4326)) As sqm; 结果：16586.2319335938 如何获得以平方米为单位的面积？我尝试了以下方法： SELECT ST_Area(ST_Transform(ST_GeomFromText('POLYGON((871325.790874952 6105405.3261047,871418.748307692 6105359.72944624,871346.22022442 6105215.141258,871254.85408906 6105261.72007212,871325.790874952 6105405.3261047))',4326)),31467) As sqm; 但是似乎ST_Transform不起作用。我收到以下错误消息： ERROR: function st_transform(geometry) does not exist LINE 1: SELECT ST_Area(ST_Transform(ST_GeomFromText('POLYGON((871325... ^ HINT: No function matches the given name and argument types. You might need to add explicit type …

10 postgis  area  wgs84 

我正在尝试在WFS查询中使用BBOX（或INTERSECTS）过滤器从一小块区域获得要素。我应该只期望一小部分功能或单个功能，但是查询返回的功能相当丰富。 我有带有Oracle数据存储的GeoServer版本2.2.2。我已禁用“宽松的bbox”选择，因此这不应成为问题的原因。 这是带有BBOX过滤器的查询： <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.1.0" outputFormat="json" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"> <wfs:Query typeName="LiVi:LIIKENNE_ELEMENTTI" srsName="EPSG:3067" xmlns:LiVi="http://172.17.14.211:8080/LiVi"> <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc"> <ogc:BBOX> <ogc:PropertyName>GEOMETRY</ogc:PropertyName> <gml:Envelope xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:3067"> <gml:lowerCorner>316600 6838986</gml:lowerCorner> <gml:upperCorner>327696 6844298</gml:upperCorner> </gml:Envelope> </ogc:BBOX> </ogc:Filter> </wfs:Query> </wfs:GetFeature> 使用INTERSECTS过滤器，问题仍然存在。多边形的大小约为10 x 10米。 <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.1.0" outputFormat="json" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"> <wfs:Query typeName="LiVi:LIIKENNE_ELEMENTTI" srsName="EPSG:3067" xmlns:LiVi="http://172.17.14.211:8080/LiVi"> <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc"> <ogc:Intersects> <ogc:PropertyName>GEOMETRY</ogc:PropertyName> …

10 geoserver  wfs  filter 

我目前正在开发一个用于评估OpenStreetMap数据的小程序，以计算弯道或弯道的状况，以便生成摩托车骑士或其他驾驶爱好者可能感兴趣的道路地图：https : //github.com/adamfranco /曲率/维基 我发现的主要限制是缺乏路面数据，该数据指示道路是沥青，混凝土，碎石路面还是最小铺装/未铺装路面。是否有整个美国或特定州都包含此信息的非专有数据集？据我所知，TIGER数据集似乎并不包含路面信息。

10 data 

我有5个Landsat TM5场景的图像拼接。在场景重叠的地方，只有左侧和右侧有一个暗带，但顶部和底部是无缝的。在Erdas 2011中更正此问题的最佳方法是什么？

10 erdas-imagine  landsat  filling 

如何在弹出窗口中显示图表？我使用的是Leaflet JS，它支持raphael插件http://dynmeth.github.com/RaphaelLayer/。是否可以在弹出窗口中创建div？我在考虑http://softwarebyjosh.com/raphy-charts/ 如果使用Leaflet无法实现，那么我只能使用Raphael地图解决方案。 谢谢

10 javascript  web-mapping  leaflet 

我在qgis-1.7.0中具有一些多边形形状，并希望确保每个多边形都是矩形。参见下图。 我有什么办法可以在QGIS（或任何其他f / oss）中做到这一点？或者如何改善多边形绘图，以确保例如每座建筑物都由矩形组成（所有天使= 90°）？ 任何提示表示赞赏！谢谢！

10 qgis  polygon  editing  digitizing  open-source-gis 

考虑阿拉斯加的NLCD2001土地覆盖数据集（下载链接）。我需要对该数据集进行重新分类，以便仅保留值41、42和43的像素；所有其他像素值应变为NoData（如果需要，则为0）。 这似乎是一项简单的任务，只需要调用一次“重新分类”工具即可。不幸的是，每次调用都会导致模糊而无益的错误消息： Executing: Reclassify "D:\ak_nlcd_2001_land_cover_3-13-08_se5.img" Value "0 40 0;41 41;42 42;43 43;44 255 0;NODATA 0" "D:\alaska_reclassified.tif" DATA Start Time: Thu Jan 03 09:23:13 2013 ERROR 999998: Unexpected Error. Failed to execute (Reclassify). Failed at Thu Jan 03 09:23:13 2013 (Elapsed Time: 0.00 seconds) 如何对该栅格数据集进行重新分类？我正在使用启用了Spatial Analyst扩展的ArcCatalog 10.0，Build 4000。

10 arcpy  raster  arcgis-10.0  reclassify  error-999998 

我刚刚更新到ArcGIS 10.1 for Desktop，当我打开一个属性表时，在任何地方都看不到它。 我想知道它是否在我的屏幕上某处，或者它是否很小，我看不到它。 我正在使用双显示器，但都没有弹出。 有没有办法可以找到我的属性表？

10 arcgis-desktop  arcgis-10.1  attribute-table 

我用Garmin设备创建了许多.gpx数据，但是由于数据量大，我多次从设备中导入了映射数据（也是因为我不想丢失它）... 结果，在将所有gpx文件另存为shapefile然后合并它们之后，我现在具有相同点的两倍甚至六倍（包括轨迹，我也将其显示为点，以便获得高程）。 有没有一种方法可以使用QGIS自动擦除重复项，或者如果没有，可以使用在Ubuntu 12.4上运行的任何其他免费软件自动擦除重复项？ 非常感谢您的帮助！

10 qgis  point 

我在互联网上以及在从ArcGIS for Desktop导出符号系统到Microstation或QGIS上进行了详尽的搜索。但没有找到满意的答案。希望从这里得到专家的解答。 假设我已经设置了一个多边形shapefile的样式（例如：边框宽度= 0.26，边框颜色=绿色，红色的简单填充，标签显示等）。当我将样式保存在QGIS中时，它以XML格式保存，该格式仅特定于QGIS。其他软件以自己的格式存储它的情况也是如此。这意味着我无法在不同软件之间交换样式信息。 尽管OGC SE标准已经到位，但是我还是很困惑为什么没有任何软件来实现它。QGIS 1.8支持SLD，但是将样式存储为QGIS特定的.qml文件的目的是什么？ 另外，由于SLD标准用于Web服务，在这种情况下（在不同软件之间交换数据）可以使用该标准吗？

10 symbology  ogc  style 

我想使用下面示例中所示的地质图，它只是地图的一个很小的部分，因此最终似乎没有任何意义。 我的想法是先将其转换为仅留下基岩和构造断层的矢量，然后再有可能选择使用“清理的”栅格文件。 关于“如何做”的任何提示

10 qgis  raster  vectorization 

我想学习如何使用Python中的国家和居住区的数据集从使用OGR http://www.naturalearthdata.com/downloads/50m-cultural-vectors/。我试图使用过滤器和缓冲区在命名国家/地区的指定缓冲区（从ne_50m_admin_0_countries.shp中的要素类ADMIN过滤）中查找点（ne_50m_populated_places.shp）。问题似乎是我不了解用于buffer（）的单位。在脚本中，我只是使用任意值10来测试脚本是否有效。该脚本会运行，但会从加勒比海地区返回命名国家为“安哥拉”的人口稠密的地方。理想情况下，我希望能够指定一个缓冲距离，例如500km，但是由于我的理解是buffer（）使用的是country.shp单位（将采用wgs84纬度/经度格式），因此无法解决该问题。对此方法的建议将不胜感激。 # import modules import ogr, os, sys ## data source os.chdir('C:/data/naturalearth/50m_cultural') # get the shapefile driver driver = ogr.GetDriverByName('ESRI Shapefile') # open ne_50m_admin_0_countries.shp and get the layer admin = driver.Open('ne_50m_admin_0_countries.shp') if admin is None: print 'Could not open ne_50m_admin_0_countries.shp' sys.exit(1) adminLayer = admin.GetLayer() # open ne_50m_populated_places.shp and get …

10 python  ogr  buffer 

我有2部智能手机（android银河和iphone），我想根据GPS坐标在地图上显示位置。 是否可以下载任何应用程序或可以注册的任何服务？每次我启动该应用程序时，智能手机的坐标都会发送到某些Web服务，几秒钟后，实际位置就会显示在地图上。用户应该通过普通的Internet浏览器从台式机访问/查看地图。 编辑：在任何已知的开源项目中，您都可以设置自己的“跟踪系统”并对其进行扩展/修改？

10 web-mapping  gps  tracking 

我有一个shapefile，其中包含大约6,500个点，我正在尝试使用WMS与Geoserver 2.2.1一起使用。直到我在使用传单的客户端应用程序上实现了过滤功能之前，一切都很好。当我在WMS查询中附加CQL_FILTER（属性过滤器，而不是空间过滤器）时，我注意到缩小时缺少功能。当我放大时，它们有时会重新出现，但并非总是如此。见下图- 在左侧的缩放级别，不显示亚特兰大。当我放大时，它是。但是，有时甚至在左侧的缩放级别上也不会显示坦帕中的点。如果再缩小3个级别，将不会显示任何点。我不确定问题是CQL_FILTER参数，因为有6,500个点很难在全球范围内注意到一些缺失点，但是我在这里作为示例展示的特定过滤器仅过滤了3个特征，何时根据缩放级别，其中有1至3个缺失，这一点特别明显。 我可以使用高级地图选项中的CQL过滤器在geoserver自己的图层预览中重新创建相同的行为，因此，我可以肯定它不是我的客户端应用程序。 至于调试此问题的步骤，我尝试使用默认的Points样式/ sld文件排除自己的图层样式。我已经关闭了所有我知道的缓存。我再次检查了我的投影是否正确-我在WGS_1984_Web_Mercator_Auxiliary_Sphere中使用WGS_1984_Web_Mercator_Auxiliary_Sphere作为投影创建了shapefile，并将图层设置为geoserver中的EPSG：3857，我认为这是等效的。我还从geoserver 2.2升级到2.2.1，并且在这两个版本中都有相同的问题。我还删除了geoserver的空间索引文件（.qix）并重新创建它，因为我在Arc中看到了类似的问题，即空间索引已损坏，但显然也不起作用。 这是Geoserver自己的图层预览的快照，其中启用了CQL过滤器，并放大到了如上所示的相同区域。红色圆圈大约是我应该看到的另一点（亚特兰大）的位置。 我尝试过调整所有我能想到的其他设置，但是没有运气。我还查看了地理服务器日志，并启用了详细日志记录，没有看到任何错误/异常。在chrome的开发人员工具中，我也看不到任何不良要求。 如果我缺少任何重要信息，我将尽力提供，但这是一个内部/非公共应用程序。

10 geoserver  wms  cql-filter 

我已将.gpx文件导入到qgis中，并使用“矢量”>“ GPS”>“ GPS工具”分离了航路点。全面查看时，我可以看到原始的航迹和零部件航路点。但是，当我放大时，轨迹线消失了，但航路点仍然存在。我努力了： 检查依赖于轨迹层比例的渲染是否为OFF 将跟踪图层导出为.shp并添加到地图 将轨迹层移动到层堆栈的顶部 切换地图渲染 但放大时仍然无法查看轨迹。有什么想法吗？ 谢谢， 托尼 [编辑]我已经上传了Pic1，该图片显示了缩放到图层范围（“ GPX跟踪日志导出”）。Pic2是使用鼠标滚轮放大并且shapefile的一部分已消失的情况。请注意，这些是作为shapefile导出的gpx文件跟踪日志。

10 qgis  vector  gpx