Questions tagged «coordinate-system»

我试图更好地掌握坐标系，基准，投影等。距离，形状，角度，面积。我想知道如何知道某个系统正在保存什么。CRS名称的哪一部分表明了这一点？具体来说，我想知道加拿大艾伯塔省的UTM系统。NAD83（CSRS）/艾伯塔10-TM（森林）

9 coordinate-system  datum  canada 

我最近观察到Sentinel 2数据的地理投影发生了变化。请查看下面的屏幕截图以进行澄清（仅颗粒T33UXP，Band 8）。第一幅图像于20160314拍摄，第二幅图像（经偏移）于20160403拍摄。我添加了一个多边形作为参考。这是从API集线器下载的未处理数据。S2A_OPER_MSI_L1C_TL_SGS__20160314T171554_A003794_T33UXP_N02.01 S2A_OPER_MSI_L1C_TL_SGS__20160403T152610_A004080_T33UXP_N02.01 还有其他人观察到吗？这是已知问题还是仅在我的特定情况下会发生？

9 coordinate-system  remote-sensing  sentinel-2 

在出于演示目的使用了一些非常规的SRC之后，为方便起见，我希望清除“最近使用的坐标参考系统”列表。 可能吗？ 我正在使用QGIS 2.10.1。

9 qgis  coordinate-system  qgis-2.10 

我正在尝试投影栅格。在R中有对此projectRaster()功能（在一个完全可复制的示例下面）： # example Raster require(raster) r <- raster(xmn=-110, xmx=-90, ymn=40, ymx=60, ncols=40, nrows=40) r <- setValues(r, 1:ncell(r)) projection(r) # project to newproj <- "+init=epsg:4714" # using raster package to reproject pr1 <- projectRaster(r, crs = CRS(newproj), method = 'bilinear') 哪个工作正常。但是，它很慢。 为了提高速度，我还是gdalwarp改用了（使用SSD时，从磁盘/ R读写数据的成本不是很高）。 但是，我无法重现projectRaster()使用的结果gdalwarp： # using gdalwarp to reproject tf <- …

9 raster  gdal  r  coordinate-system  gdalwarp 

我注意到，在新的qgis升级版本上的热图计算得到了没有任何坐标系的栅格图。 你有同样的问题吗？

9 qgis  coordinate-system 

首先，我是GIS的新手。 我有一个与GPS和遥测数据相关联的激光雷达点云。然后，我将每个点的坐标转换为ECEF坐标系，最后使用局部切平面（LTP）近似将角坐标转换为笛卡尔坐标系。（此近似值使用第一个点的gps坐标作为参考系统原点，然后相对于第一个点参考其余点）。 我正在使用liblas编写LAS文件，但我真的不知道如何将地理参考信息设置为LAS文件。我正在写LTP近似中每个点的坐标，所以我有2个问题： 编写相对于LTP近似值的点坐标是否正确，还是应该以经/纬（gps）单位编写坐标？如果正确，如何设置地理配准尊重LTP？ 如果“ 1”不正确，我必须以纬度/经度单位写入坐标，如何设置地理配准？用WGS84值编写VLR是否足够？

9 coordinate-system  georeferencing  lidar  latitude-longitude  liblas 

有人知道ESRI WKT和OGC WKT投影格式字符串之间的区别的确切清单吗？ 我知道有各种工具可以帮助您从ESRI WKT转换为OGC WKT，包括GDAL实用程序和各种网站服务。但是我的问题不是实用的，我只是想了解这些服务使用的格式/语法差异。先前的Stackexchange问​​题仅讨论了特定示例之间的差异，或者仅涉及可用工具和服务。 即使您只知道一个区别，也可以将其发布。根据我自己的经验，应该只有很少的差异。我知道的区别是： esri定义中的大多数文本项使用下划线，而ogc使用空格。 esri wkt中定义基准的文本与ogc wkt相同，只是以“ D_”开头。 有时，某些预定义的PROJCS，PROJECTION，GEOCCS和DATUM的文本标识符以不同的方式编写（例如，“ NAD83”为“ North_American_1983”）。我猜想知道哪些标识符的拼写不同的唯一方法是拥有列表或查找表，因此请给任何您知道不同的名字命名。 除了ogc每个单词的大写字母大写而esri所有小写的字母之外，各种PARAMETER文本值都相同。但是，我见过没有使用此规则的情况，有人知道在尝试加载它们的软件中标题是否真正重要？ UNIT的类型在ogc中拼写为大写字母，在esri中拼写为小写，例如“ Degree”与“ degree”。在某些情况下，我看到ogc的“米”被拼写为“ meter”和“ m”，而在其他情况下，法语的拼写为“ metre”。有人知道这两种格式或任何其他单位类型的正确约定是什么吗？

9 coordinate-system  wkt  ogc 

我正在尝试为丹麦市政当局查找shapefile的质心，然后找到它们之间的行驶时间。我使用R的readShapePoly从功能maptools与组合gCentroid，从功能rgeos，和一切正常。但是，我得到诸如 SpatialPoints: x y 1 571860.7 6225016 Coordinate Reference System (CRS) arguments: NA 我显然不能在Google中使用这来获取旅行时间。我正在寻找一种将这些数字转换为经度-纬度的方法，但不知道如何。 当我readOGR从rgdal库中读取数据时，我得到了相同的坐标，但是它告诉我以下有关我假设是投影的信息（但是坐标是相同的） Slot "proj4string": CRS arguments: +proj=utm +zone=32 +ellps=intl +units=m +no_defs 可重现的示例：我将示例的数据放在这里：https : //github.com/sebastianbarfort/shapefiles 这应该重现该问题： library(maptools) library(rgdal) library(rgeos) map = readShapePoly("~/Downloads/shapefiles-master/kommuner1983.shp") centroid = gCentroid(map) centroid

9 shapefile  r  coordinate-system 

我从墨卡托投影页面上的这张 Wikipedia图片中获得启发，并用传单散布。想法是交互式地使您的光标移动到一个大陆上，并且将根据投影对其进行校正。 有关状态描述者，请参见以下JSFiddle； http://jsfiddle.net/cantdutchthis/NrAEL/1/ 我对投影世界有些陌生，尽管演示可以很好地适应纬度。当我比较维基百科上的等角投影和墨卡托投影时，我得到的印象是，当您比较这两个投影时，投影的偏差仅在“ y轴”上。 我得到的印象是由于第一张维基百科图片，我没有校正经度。我想支持这一点，但我不知道实现此目的的最佳方法是什么。 目前，我能够显示墨卡托投影的偏差的方法是每次鼠标移动时都将坐标移动到鼠标光标上。 map.on("mousemove", function(d){ Object.keys(myLayer._layers).map(function(d){ myLayer.removeLayer( myLayer._layers[d] ) }) var plotCoords = []; // green is the coordinate array for the landmass var green = [ [-46.76373672, 82.62796021],... ] for ( i in green ) { // shift the landmass to where the cursor …

9 javascript  coordinate-system  leaflet  spherical-geometry  mercator 

我有一个点云，其坐标是相对于局部坐标系的。我也有带有GPS值的地面控制点。是否可以使用PROJ.4或任何其他库将这些局部坐标转换为全局坐标系？ Python中上述问题的任何代码都将提供很大的帮助。

9 python  coordinate-system  gps  proj  ground-control 

我有一个分布在3个UTM区域（36N，36S，37S）的调查网格。我想找到这些网格的质心到道路及其之间各个点的最近（或最短）距离。 使用任何类型的平面投影时似乎都有太多折衷（阅读：关于保留地图上任意数量点之间的距离）。在这种情况下，应该只是忘记使用投影，而去使用Goedesic或椭圆形（阅读：Geographic）技术吗？ 谁知道有一种平面技术可以保留地图上任意数量点之间的距离？似乎我可以使用等距投影，但侏儒投影除外。它是否正确？

9 coordinate-system  cartography  measurements 

我在墨卡托投影（google maps）上有一条线（Ax，Ay-Bx，By）和最接近该线的随机点（Cx，Cy），我想知道最近的点（图像上的透明蓝色）指向那条线（图像中为蓝色） 编辑：澄清这是在墨卡托投影（球面投影）

9 google-maps  coordinate-system  geometry  php  spherical-geometry 

四边形球面立方图投影[1]与Snyder立方等面积图投影[2]是否相同？ 我知道两者都是基于立方体的等面积多面体地图投影（http://www.progonos.com/furuti/MapProj/Dither/ProjPoly/projPoly.html），但我无法访问任何出版物来回答自己的问题。 谢谢。 [1] EM O'Neill和RE Laubscher，对四边形球形立方体地球数据库的扩展研究，技术。报表3-76，计算机科学公司银泉系统系统科学部，1976年。 [2] Snyder，JP1992。多面体地球仪的等面积地图投影。制图学29（1）：10-21。

9 coordinate-system  cartography 

我正在尝试自定义Albers和Hotine斜墨卡托（HOM）投影，以最小化我正在分析的区域中的失真。该区域的纬度范围从大约51度到62度，覆盖了大约乌克兰大小的区域。该区域面向西北-SE。 我想确保使用正确的方法来确定两个投影参数：投影中心的经度/纬度和中心线方位角。我正在使用ArcMap v10。到目前为止，这是我遵循的过程： 创建了一个定义分析区域的多边形（通常，通过围绕覆盖该区域的分水岭范围创建一个凸包）。这个多边形是我要为其定制投影的区域。 将多边形投影到Geographic / NAD 83。 使用Jeff Jenness的“图形和形状工具” （http://www.jennessent.com/arcgis/shapes_graphics.htm）确定GRS80椭球上多边形的质心。结果坐标是我用于“投影中心”参数的坐标。 为了确定中心线方位角，我首先将多边形投影到一个方位角等距投影，在步骤3中确定的坐标处指定投影中心。 然后，我绘制了一条折线（在方位角等距投影中），捕捉到投影中心点，代表区域多边形的方向趋势。为了获得投影中心的方位角，我使用了Jeff Jenness的“图形和形状工具”来确定测地曲线在中心点的起始方位角。 对于Albers投影，我将使用步骤3中确定的投影中心的经度。还使用Bill Huber（http://forums.esri.com/Attachments/34278.xls）创建的出色电子表格来确定位置放置标准平行线以最小化多边形区域内的比例尺失真。 如果需要知道，我使用的是ArcMap的HOM版本，该版本使用由投影中心的一个点及其方位角定义的中心线。ESRI将其称为“ Hotine_Oblique_Mercator_Azimuth_Center”。在EPSG中，我相信这是Obcque Mercator，Hotine Variant B，EPSG方法代码9815。 我希望有一些投影专家可以告诉我上述过程（尤其是步骤3和4）是否是确定所需投影参数的正确方法。我在正确的轨道上吗？确定球体上的中心和测地线与中心点的角度（而不是“ 2d”几何中心和方位角）是否正确？ 我希望问题描述清楚。我热切期待任何答案，技巧，讨论等！

9 arcmap  coordinate-system 

我正在尝试通过天梭（Tissot）指标计算等角投影的变形。我曾尝试上的指示以下这篇文章（除其他事项外），但它超出了我的理解作为一个业余。 因此，我想知道是否有人会为一个等长的经纬度（例如，您最喜欢的并且在等长角投影上变形的那个）计算单个Tissot椭圆。我不完全了解变量是什么以及它们来自何处，因此了解实际的方程式将非常有用。 我试图理解这些方程，以便可以将它们插入正在编写的映射程序中。我在此线程中提出了一些一般性问题，但我认为一个具体的示例将帮助我弄清楚其余的问题。 一如既往，非常感谢。 NCashew

9 coordinate-system  algorithm