Questions tagged «coordinate-system»

我们使用形状文件进行数据交换。我看了.prj文件：不同的程序会产生不同的prj文件。 是否有一个标准，必须如何在prj文件中指定参数？ 可以自由指定投影的名称（例如PROJCS [“ ETRS89 / UTM区域32N（zE-N）”，…]）。但是我不相信基准名称（例如DATUM [“ D_ETRS_1989”，…]）也可以自由分配。 我没有找到有关此的信息。即使在ESRI Shapefile技术说明中，也没有有关.prj文件的信息。在ESRI知识库文章14056 Howto：如何创建shape文件的投影元数据（.PRJ）文件说：“你必须使用地图投影和参数对象的预定义的名字。” 但是本文没有指定在哪里可以找到这些预定义名称。每个GIS公司是否都指定自己的名称？还是EPSG指定这些名称？ 以下是QGIS，AutoCAD Map和spatialreference.org之间的.prj文件（EPSG 4647）部分中的一些区别。

15 shapefile  coordinate-system  standards 

墨卡托投影与纬度一起增加了比例。格陵兰岛似乎比澳大利亚大，但实际上要小得多。但是，等面积投影的形状远离赤道。我想不存在同时保留大小和形状的投影。但是，出于比较国家的目的，实际上并不需要完全保留任何财产。变形在地图的任何一点都应该大致相同。最接近完成效果的地图投影是什么？

15 coordinate-system  maps 

如果我希望将某些shapefile的投影从..更改NAD83 (4269)为WGS84 (4326)...，并使用某些第三方应用程序....这是否意味着该.prj文件仅被更新..而不是.shp文件或其他内容的内容关联文件？

15 shapefile  coordinate-system 

我想知道在QGIS中集成NTv2转换的正确方法是什么。 我已经成功地将NTv2转换集成到QGIS中，但是我不确定它是否按照预期的方式完成。 为了在QGIS中集成NTV2网格，我发现了一些信息，即网格名称应存储在srs.db的表tbl_datum_transform的p1列中： https://github.com/qgis/QGIS/pull/1502#issuecomment-50580199 https://dragons8mycat.com/2014/03/03/gis-tips-getting-ostn02-to-work-with-qgis-2-2/ https://web.archive.org/web/20140224182157/http://sourcepole.ch/ntv2-transformations-with-qgis http://isticktoit.net/?p=371 （http://www.gkg-kassel.de/pdf/Praezise_Transformation_im_QGIS_mit_Batchmodus.pdf（德语）的屏幕截图） 使我怀疑是否正确的事实是，该列似乎具有数字字段类型： 七参数转换的第一个参数不是存储在QGIS的tbl_datum_transform表的p1列中吗？ 根据sqlite常见问题解答，可以将字符串存储到sqlite中的数字列中（https://sqlite.org/faq.html#q3），但是在我看来，将字符串存储到数字列中是一种hack或解决方法。这真的是要完成的方式吗？ tl; dr：在QGIS中集成NTv2转换的正确方法是什么？

15 qgis  coordinate-system  proj  epsg 

我无法在QGIS 3.0中捕捉2个多边形。我尝试从同一图层捕捉2个多边形，也从不同图层捕捉2个多边形。 我确保在“设置” >>“选项” >>“数字化”选项卡中启用了捕捉，并增加了捕捉距离。 我在图层级别找不到任何“快照设置”。投影层。

15 qgis  coordinate-system  polygon  qgis-3  snapping 

在R的HistData数据包（https://r-forge.r-project.org/R/?group_id=574）中，我具有与约翰·斯诺（John Snow）1854年伦敦霍乱疫情地图有关的数据集。我相信它们在Walter Tobler的监督下经过精心数字化处理后，具有权威性。John Mackenzie在http://www1.udel.edu/johnmack/frec480/cholera/cholera2.html上描述了有关这些数据集的一些详细信息。 不幸的是，死亡，水泵和街道的坐标使用任意坐标系，而不是适用于R中其他GIS应用程序或地图软件（空间包，ggmap等）的地图坐标。 在http://freakonometrics.hypotheses.org/19213中， Arthur Charpentier将ggmap与来自http://www.rtwilson.com/downloads/SnowGIS_v2.zip的John Snow数据版本一起使用 。Cholera_Deaths.shp但是，该文件仅列出489人死亡，而不是我在中记录的578人HistData::Snow.deaths。 一种想法是找到均值与（x，y）坐标的标准偏差之间的关系并进行线性缩放，但是也许有更好的方法吗？ 到目前为止，这是我尝试过的 > data(Snow.deaths, package="HistData") > D <- Snow.deaths[,2:3] > colMeans(D) x y 13.03312 11.69721 > var(D) x y x 3.8150987 0.3802654 y 0.3802654 2.7213828 读取Cholera_deaths文件 > folder <- "C:/Dropbox/R/data/Snow/SnowGIS_v2/SnowGIS" > library(maptools) > deaths <- readShapePoints(file.path(folder, "Cholera_Deaths")) > head(deaths@coords) coords.x1 …

14 r  coordinate-system  scale 

我有一个多边形的shapefile和另一个CSV文件，其中包含（Lat，Lng）对点列表。 我想检查CSV文件中每个（纬度，经度）对属于哪个多边形。 shapefile被投影，proj文件的内容如下： PROJCS["Transverse_Mercator",GEOGCS["GCS_OSGB 1936", DATUM["D_OSGB_1936",SPHEROID["Airy_1830",6377563.396,299.3249646]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse_Mercator"],PARAMETER["latitude_of_origin",49],PARAMETER["central_meridian",-2],PARAMETER["scale_factor",0.9996012717],PARAMETER["false_easting",400000],PARAMETER["false_northing",-100000],UNIT["Meter",1]] 我的计划如下： 使用readShapePolyR MapTools包中的函数读取shapefile 。 将点坐标从CSV文件读取到数据框中，然后将其转换为SpatialPointsDataFrame 使用over函数确定其落入哪个多边形。 为此，我需要proj4string在第1步中加载shapefile 的同时指定时间，并在使用第3步中的spTransform功能之前，先使用函数将CSV文件中的坐标转换为相同的投影系统over，因为它要求点和多边形必须在同一投影系统下。 是否知道上面显示的proj文件内容的正确值是什么？

14 coordinate-system  shapefile  r  proj 

也许是一个奇怪的问题，但是在我提出实际问题之前，让我先简单介绍一下背景知识： 简而言之（据我所知），原子力显微镜（AFM）使研究人员可以在微米和纳米级扫描区域。它通过使用各种探针“扫描”区域来工作。由于我对此没有真正的了解，因此我难以解释的更多内容。我所知道的，引起我好奇的是，结果实际上是“高度”值的“网格”（一个512x512值的矩阵，描述了该点的高度）。 然后我想：好吧，除了比例尺，这实际上是一个数字高程模型！而且，这意味着，如果我可以设法创建GIS工具所理解的DEM文件，则可以对其应用GIS分析！ 实际上，我的其他重要工作是在拥有AFM机器的实验室中进行的，并且正在她的一个项目中使用它。我从她那里得到了一些扫描文件，并使用Python（struct和numpy）进行了管理，以解析这些二进制文件，而我现在得到的是一个numpy数组，大小为512x512，填充有int16值。 我下一步要计划的是什么，而我需要帮助的是“映射到适当的DEM”部分。我对DEMS有一些了解，但是当涉及到DEMS的实际生成时，我还是很新的。 我在想的是我必须以某种方式对数据进行地理参考，为此，我需要一个自定义（平面）坐标系。我设想我的坐标系将以微米或纳米为单位。然后，只需查找用AFM扫描的区域的大小即可（我相信这是在二进制文件中的某个位置，假设已知道这一点）。 更新：我也以不同的分辨率进行了几次扫描，但扫描范围相同。例如我有关于两次扫描的信息： 放大图片： Scan Size: 51443.5 nm X Offset: 0 nm Y Offset: 0 nm 较小（详细）图像： Scan Size: 5907.44 nm X Offset: 8776.47 nm Y Offset: 1486.78 nm 我的想法是，我的自定义坐标系应以0,0为原点，对于较大的图像，我可以为像素0,0分配（0,0）的坐标值，为像素512,512分配坐标值（51443.5、51443.5 ）（想必您已获得所需其他点的图片）。 然后，较大的图像会将像素（0,0）映射到（8776.47，1486.78），将（512,512）映射到（8776.47 + 5907.44，1486.78 + 5907.44） 那么第一个问题是：如何为这样的坐标系创建proj4 def？即：如何将这些“现实世界的坐标”分配给我的自定义坐标系统（或者，如果我遵循野蛮人的建议并使用局部坐标系统并围绕单位说谎（即，将纳米视为千米） 然后，我必须将我的numpy二维数组转换为地理参考DEM文件格式。我当时在考虑使用GDAL（或更确切地说是Python绑定）。 然后，第二个问题是：如何从诸如我的“任意”数据创建地理参考DEM？最好在Python中使用开放源代码库。 剩下的就应该相当容易了，只是使用正确的分析工具即可。问题是，这项任务是由我自己的好奇心驱动的，因此我不确定我应该对纳米级DEM实际做什么。这乞求 第三个问题：如何处理纳米DEM？可以进行什么样的分析，什么是进行DEM分析的合适工具，最后：根据该数据制作具有山体阴影和轮廓线的地图是否可行？ 我欢迎所有建议和指示，但是请记住，我正在寻找免费替代方法，因为这是一个严格的基于兴趣的项目，没有预算或资金（而且我无法使用任何许可的GIS应用程序）。此外，我确实知道出售这些AFM机器的公司Bruker会提供一些软件，但是使用它不会很有趣。

14 python  coordinate-system  dem 

我在地理坐标系中有点，我想将它们转换为瑞士网格（CH1903 +）。 样本数据： id lon lat 2 7.173500 45.86880 3 7.172540 45.86887 4 7.171636 45.86924 5 7.180180 45.87158 6 7.178070 45.87014 7 7.177229 45.86923 8 7.175240 45.86808 9 7.181409 45.87177 10 7.179299 45.87020 预期结果： id E N 2 2579408.2431 1079721.1499 3 2579333.7158 1079729.1852 4 2579263.6502 1079770.1125 5 2579928.0358 1080028.4605 6 …

14 convert  coordinate-system  r 

我安装了proj 4.8.0的QGIS。我想添加Winkel Tripel投影，该投影似乎定义为 +proj=wintri 但是该定义字符串不起作用。这是我的proj版本中的错误，还是正确的定义字符串？

14 qgis  coordinate-system  proj 

我正在使用传单地图来创建环球挑战的表示。我想添加一条折线，该折线从东京向东行驶，然后在地图上出现在南美洲的西部-但我得到一条沿相反方向与地图交叉的线（请参见黄线）。 我认为这可能与日期线和/或坐标系有关，但在细节上有些粗略。谁能解释我需要做些什么才能使它起作用？我正在使用美国国家航空航天局的bluemarble投影： var bluemarble = new L.TileLayer.WMS("http://demo.opengeo.org/geoserver/wms", { layers: 'bluemarble', attribution: "Data © NASA Blue Marble, image service by OpenGeo", minZoom: 2, maxZoom: 5 });

14 coordinate-system  polyline-creation  leaflet 

我对GIS非常陌生。 我正在gdal阅读土地利用/土地覆盖图，我需要选择某些土地覆盖类型的经/纬度以索引到仅以经/纬度表示的不同数据集。不幸的是，我不了解从geotransform给我的x和y坐标的形式，特别是originXand originY以下： geotransform = dataset.GetGeoTransform() originX = geotransform[0] originY = geotransform[3] 打印这些值会给我类似的坐标(447466.693808, 4952570.40529)。这些与原始纬度和经度有何关系？ 编辑： 这是一个简单的python示例，为我提供了我想要的东西： srs = osr.SpatialReference() srs.ImportFromWkt(dataset.GetProjection()) srsLatLong = srs.CloneGeogCS() ct = osr.CoordinateTransformation(srs,srsLatLong) print ct.TransformPoint(originX,originY) 从以下地方窃取：tolatlong.py

14 gdal  coordinate-system 

我正在使用OpenLayers 3与某些地图进行交互。我首先声明我的地图： map = new ol.Map({ target: 'map', layers: [ new ol.layer.Tile({ source: new ol.source.OSM() }) ], view: new ol.View({ center: [0, 0], zoom: 2 }) }); 我有一个事件触发操作，以更改地图的视图中心。这样，（我的坐标是EPSG:4326格式）： function CenterMap(lat, long) { console.log("Lat: " + lat + " Long: " + long); map.setView(new ol.View({ center: ol.proj.transform([lat, long], 'EPSG:3857', 'EPSG:4326'), zoom: …

14 openlayers-2  coordinate-system  coordinates  epsg  openlayers 

我正在尝试使用qgis和“来自太空的世界” 投影http://spatialreference.org/ref/sr-org/6980/（本质上是正交投影）来创建球形视图。ArcGIS正确包装了形状，但是QGIS（2.01）产生了讨厌的工件。 我必须定期制作具有不同角度的地球仪，所以外面的人是否有解决此问题的想法？

14 qgis  arcgis-desktop  shapefile  coordinate-system 

我想知道光栅图像是如何工作的， 我有一组没有aux，tfw或xml的geotiffs 当我检查arccatalog中的空间信息时，它们确实具有投影。并且它们显示正确。 如果投影嵌入到Geotiff中，为什么需要.tfw？ 我需要知道它的作用，如果Geotiff在真实坐标wgs84中并且我想要另一个用于快速显示的投影，它将用于快速转换吗？还是我的地理学者缺乏信息？ 当我定义另一个投影时，arccatalog创建辅助，我是否需要保持该辅助？ 所以。缩短问题，如果信息已经存储在Geotiff中，则aux，tfw有什么用？

14 coordinate-system  geotiff-tiff