Questions tagged «coordinate-system»

我是QGIS的新手，我不太清楚网格和缩放的工作方式。我"Project CRS"和我"Grid CRS"应该一样吗？ 我玩了一些设置，经过一段时间后，设法得到了我认为正确的网格（基于对现有地图的一些比较）。但是我仍然不确定如何在我打算制作的其他地图上进行设置，而Google并没有帮助。 我的问题是-我的秤现在没了。例如，我知道萨格勒布到Sisak的距离约为50公里，根据我的地图，该距离约为70/80。其他一些例子甚至更极端。问题是什么？为什么我的秤现在可以工作？ 我完美地对我的图像进行了地理定位-完美地与下面的OpenStreetMap图层匹配-但我的比例尺仍然关闭。投影是世界墨卡托。

10 qgis  coordinate-system  qgis-3 

关于坐标 31.96212, -103.004715 UTM转换器给出的UTM坐标为13/R/FR。 示例转换器在这里：http : //www.rcn.montana.edu/resources/converter.aspx 但是它们很多，对于这些坐标，它们给出相似的答案。 同时，在Sentinel-2数据集中，请访问http://sentinel-s2-l1c.s3-website.eu-central-1.amazonaws.com/#tiles/13/R/ 我找不到FR子目录。 在google中，此位置在这里： 在我看到的Sentinel图片浏览器中找到相同的位置，图块是不同的 其代表13/S/FR即相同的UTM和方形的，但不同的频段。 这怎么可能？ 更新 具有Sentinel-2磁贴的KML还会报告S给定位置的磁贴 更新2 根据这张图片 采取从这里，该FR广场位于一半SUTM区和半R区。显然，大多数自动转换器将此正方形分配给R区域，而Sentinel-2则将其分配给S区域。 这里有什么真相吗？ 更新3 简单的Python代码，从此处获取https://gis.stackexchange.com/a/224994/32207 bandVals = "CDEFGHJKLMNPQRSTUVWXX" lon = 31.96212 lat = -103.004715 zone = int(lat + 186.0) / 6 if (lon >= 84.0): band = 'Y' if (lat < 0.0) else …

10 coordinate-system  utm  sentinel-2  mgrs 

我正在将许多地理数据集导入PostGIS，而它们具有不同SRID的。（有些有EPSG:3857，有些EPSG:4326，有些其他）。 我想另外创建一个geometry column，例如。the_geom_mercator使用SRID EPSG:3857，并将原始geom列保留SRID在其中。 如何使用PostGIS功能做到这一点？

10 postgis  coordinate-system  srid  mercator 

当我打开QGIS，添加图层并通过字段计算器计算shapefile的面积时，我得到的面积与打开QGIS并选中“启用CRS动态转换”并计算面积的面积不同。尽管要确保项目和图层具有相同的坐标系（相同的EPSG编号），但这仍然可行。我究竟做错了什么？ 我有一个shapefile，其中包含使用ArcGIS进行面积计算的信息（不是我自己，数据是交给我的，我不知道使用ArcGIS计算面积的CRS）。形状文件层CRS是EPSG：21781（瑞士）。在QGIS中，如果我不更改OTF设置并将项目CRS保留为EPSG：4326（WGS84），则得到的值与ArcGIS区域值相同。但是，如果在将图层添加到EPSG之前更改了OTF：21781，则会得到不同的面积值。据我了解，这表明ArcGIS Area是使用CRS EPSG：4326计算的。 第一个工作流程： 打开QGIS 项目CRS：EPSG 4326 添加层 项目CRS自动适应，现在为EPSG 21781 使用字段计算器计算$ area 第二个工作流程： 打开QGIS 项目CRS：EPSG 4326 打开OTF，将项目CRS设置为EPSG 21781 添加层 使用字段计算器计算$ area 第一和第二个工作流程的步骤5不要产生相同的区域。

10 qgis  coordinate-system  area 

我想使用python和对栅格进行地理配准GDAL。我当前的方法是调用gdal_translate和gdalwarp使用os.system地面控制点列表。我真的很想在本地进行此操作python。 这是我正在使用的当前过程： import os os.system('gdal_translate -of GTiff -gcp 1251.92 414.538 -7.9164e+06 5.21094e+06 -gcp 865.827 107.699 -7.91651e+06 5.21104e+06 "inraster.tif" "outraster1.tif"') os.system('gdalwarp -r bilinear -tps -co COMPRESS=NONE "outraster2.tif" "outraster3.tif"') 有一个以前的问题和答案在2012年里面写明了gdal_translate可以导入后进行访问gdal。我不确定是否已过时，或者是否错误，但是当我运行时，from osgeo import gdal我不认为这gdal.gdal_translate是一个选择。 我不知道它是否存在，但是如果我能以pythonic的方式翻译和重新投影栅格，我将不胜感激。例如： # translate gcp_points = [(1251.92, 414.538), (-7.9164e+06, 5.21094e+06)] gdal.gdal_translate(in_raster, gcp_points, out_raster1) # warp gdal.gdalwarp(out_raster1, out_raster2, 'bilinear', args*) 这样的方法可行吗？

10 python  raster  gdal  coordinate-system  gdalwarp 

每当我在全球或极地查看Google Maps时，投影都会很麻烦。是否存在具有可缩放世界地图的Web服务，该服务使用等面积投影-最好是伪矩形（例如Mollweide投影）？甚至《国家地理》在线地图也使用墨卡托。 另一种选择是Google Earth，它并没有特别扭曲地球，但是也没有给出全局视图。 另请参阅：为什么Web Mercator（辅助领域）成为Web地图标准？ 编辑2018-08-24： Google Maps不再使用Web Mercator，但已移至看起来是正交投影的位置。

10 coordinate-system  web-mapping 

我在EPSG 4326中有一个世界大陆的shapefile。 我必须将文件重新投影到哪个参考系统上，才能在下面像这样的地球仪进行投影？

10 coordinate-system 

我有一个NetCDF文件，其中包含英国下雨天数的数据。我想将其转换为GeoTIFF，但无法正确获取GDAL地理参考。请参阅下面来自gdalinfo的转储。 我想要一个显示“ raindays_tmean_abs”变量的TIFF乐队。“ lat”和“ lon”变量似乎具有合理的最小/最大值，可以将数据放到整个英国，但是当我使用gdal_translate进行转换时，它最终出现在非洲西海岸。 关于第一个尺寸不是“时间或垂直”尺寸的警告。我读到gdal期望数据按一定顺序排列，所以也许这就是问题所在，但是我不知道如何修改它。 如何转换/重新投影此数据？ gdalinfo NETCDF：RainDays1_1961-1990_LTA_25km.nc Warning 1: dimension #0 (meaning_period) is not a Time or Vertical dimension. Driver: netCDF/Network Common Data Format Files: RainDays1_1961-1990_LTA_25km.nc RainDays1_1961-1990_LTA_25km.nc.aux.xml Size is 39, 52 Coordinate System is `' Origin = (-17.650008352179277,11.329999312758446) Pixel Size = (0.219999614514803,-0.219999998807907) Metadata: axis_0#units= meaning_period#comment=Meaning periods are all …

10 gdal  coordinate-system  netcdf 

我有一个带有纬度经度（NAD27）列的表。我计算另外两个列X和Y，它们分别代表Web Mercator（WGS84）的位置。 目前，我正在使用Arcmap来完成此任务，方法是将研究区域的推荐地理转换（即3参数（地心）地理转换）应用到NAD27到WGS84。 我想完全在Sql Server 2012中执行此操作。据我所知，Sql Server不支持开箱即用的数据转换。有人知道支持这种地理转换的Sql库吗？我想简单地在Arcmap中使用与当前在Arcmap中使用的系数相同的系数。 我还需要将WGS84经纬度投影到网络墨卡托。我看到这个公式是在javascript中实现的，但是如果有人拥有执行此操作的Sql存储过程，那就太好了。

10 sql  sql-server  coordinate-system  datum 

我在ArcGIS中有两个点数据集，这两个数据集都是在WGS84纬度/经度坐标中给出的，并且点分布在整个世界上。我想找到数据集A中与数据集B中每个点最近的点，并以公里为单位获取它们之间的距离。 这似乎是对Near工具的完美使用，但这使我得到了输入点坐标系中的结果：即十进制度。我知道我可以重新投影数据，但是我（从这个问题）收集到（很难（如果不是不可能的话））很难找到一个可以给出全世界准确距离的投影。 该问题的答案建议使用Haversine公式直接使用纬度-经度坐标来计算距离。有没有办法做到这一点，并使用ArcGIS获得以km为单位的结果？如果没有，什么是解决此问题的最佳方法？

10 arcgis-desktop  coordinate-system  r  distance  units 

我有一些地磁数据。也就是说，它们是经度和纬度，但是是指加拿大的磁性北极而不是地理北极。您可能会猜到，我想将它们与其他坐标系中的其他数据进行匹配。 （数据起源于极光椭圆的计算；即，它们起源于地磁，即使在原则上也没有可以参考的地理资源。） 数据将存储在PostGIS中。 这样做的简单方法是找到地磁坐标系的EPSG代码，用该代码标记数据，然后将所有内容神奇地重新投影到我当时需要的任何位置。但是，我在spatialreference.org上找不到包含地磁或磁场的任何EPSG代码。各种Google搜索也都没有结果。 对于我需要的精度水平，任何经修改以具有不同北极的纬度/经度坐标系（例如WGS84）就足够了。因此，如果有一种简单的方法来定义这样的坐标系并将其加载到PostGIS中，那也可以。 澄清度 例如，假设我只是想从UTM 15N转换为WGS84，而不是从地磁坐标转换为地理。我可以将UTM对象未经修改地保存到PostGIS中，并用EPSG代码26915对其进行标记。然后，当我稍后进行查询时，我可以对PostGIS说“请在EPSG代码4326中返回结果”，然后PostGIS将自动转换为WGS84。 同样，我想用一些坐标系标记地磁数据，并让PostGIS在幕后进行重新投影，而无需reproject()手动调用函数。

10 postgis  coordinate-system  epsg 

我想扩展我的知识坐标系，投影和重新投影数据，以期编写一个重新投影引擎。 您能为我提供一些有关坐标系和重新投影数据背后的数学知识的资源（书籍，网站，卷轴）吗？

10 coordinate-system  references  reprojection-mathematics 

我得到了一些没有投影的shapefile。因此，我需要添加投影才能在Google地图上显示。我在http://spatialreference.org/ref/epsg/2066找到了所需的投影信息。 问题： 如何使用开源方法添加投影（无法访问ESRI的软件）？ 谢谢 克里斯

10 shapefile  coordinate-system 

我正在尝试制作使用方位角等距投影的世界或大陆比例尺地图。我是一名无线电爱好者，从天线上看，整个世界看起来像是在一个平面上展开，天线在中心，而整个世界从那一点起（r，Θ）。 有软件可以执行此操作-AZ_PROJ-但它使用自己的数据库，并且完全用PostScript编写，因此使用常规的地理工具无法完全破解。 OGR具有“ aeqd”投影，可以这样调用： ogr2ogr -t_srs "+proj=aeqd +lat_0=43.7308 +lon_0=-79.2647" out.shp in.shp 中心位于43.7308 N，79.2647W。不幸的是，这包裹了当前“地图”后面的半球以外的所有区域。它还将坐标投影到笛卡尔米中，在这里我更喜欢极坐标。 有其他建议吗？

10 coordinate-system  cartography 

我正在使用ESRI堆栈，将图层存储在启用了sql-spatial的SDE地理数据库（几何类型，Web Mercator-3857）中。 我正在构建一个Web制图应用程序，因此默认情况下，图块也位于Web Merator 3857中。 通过存储的proc，我使用STDistance来查询从用户位置（Web Mercator中的坐标）到各个层的距离。 问题是，由于网络墨卡托变形，距离距离赤道越远，我的距离计算越来越少。 我曾想过将图层存储在sql-spatial-geography（而非几何）类型中，但是： 我想我的距离查询将花费更长的时间（球面上的距离计算） 我将需要重新导入很多数据 arcgis服务不会像他们需要的那样快速地进行投影 如果我去Google地图并进行距离计算，即使在北/北地区，返回的距离也要精确得多，因此我认为Google必须纠正由网络墨卡托投影引起的失真。 那么我的问题是：是否有一个简单的因子值可以应用于网络墨卡托投影中进行的距离计算，以获得“正确的”距离？

10 coordinate-system  distance