Questions tagged «coordinate-system»

由一组点，线和/或曲面以及一组规则组成的参考框架，用于定义二维或三维空间中的点位置。此标签还指制图投影和坐标系转换。

2
如何在PostGIS中获得两个几何之间的测距仪值？
我有一个关于在PostGIS中计算距离的简单问题。 我想得到两个几何之间的距离。我正在使用此sid：4269以米为单位我现在正在做的是：ST_Distance((a.geom,b.geom)) FROM ...但是我得到的结果是以度为单位。我认为我应该处理地理问题，但是如何将地理模型转换为地理？我尝试了（a.geom :: geography），但出现错误。 我该怎么做才能以米为单位得出结果？ 谢谢
2
改善地理配准结果？[关闭]
已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 2年前关闭。 背景这是我的第二个问题，涉及对裸栅格地图进行地理配准，以便在不同的坐标系上以及与其他数据层一起重新可视化它们。上一个问题是将任意无元数据的地图图像转换为QGIS项目 问题我的目标是对该地图进行地理配准： 这似乎不是Plate-Carrée。因此，在QGIS中，我创建了几个合理的控制点，为完整起见，我将其附加在底部[ref：1]。我为QGIS Georeferencer提供了与我的项目文件EPSG：4326相同的目标SRS。使用Helmert和多项式变换时，我得到的结果异常差，但是用薄板样条曲线得到了合理的图像（这使得结果的地理估计值通过了我的控制点）。但是，即使这样的结果也很差，例如在较高的纬度上（参见日本北部的俄罗斯海岸）。这是我使用自然地球背景的QGIS屏幕的屏幕截图。 替代路径我在MapWarper上使用了更易于使用的工具尝试了类似的练习：在http://mapwarper.net/maps/758#Preview_Map_tab上查看结果和控制点，我得到的结果较差（可能是由于事实我添加了更少的控制点）。 简而言之的问题 有什么技巧可以使我获得良好的地理参考？ 可以立即识别出这种投影吗？ 在旧图纸上的“未知坐标系”中，gdaltransform建议将多个坐标点转换为多个目标SRS，以实际发现用于生成原始地图的投影参数为目标。我尝试过这样的事情：保存我的QGIS点列表之后，我进行了一些字符串处理，以获取空格分隔的经/纬度列表，cat eurasian-steppe-gcp.points | tail -n+2 | cut -d, -f1-2 | sed 's/,/ /'> tmp.txt并将结果文件输入到gdaltransform中：gdaltransform -s_srs EPSG:3785 -t_srs EPSG:4326 < tmp.txt并切换s_srs和t_srs标志（该项目使用EPSG：4326）。我知道我在黑暗中射击，希望能走运，所以当我无法理解输出时，我并不感到惊讶。有人可以扩展我将如何使用此方法来找到源地图的投影和投影参数的最佳估计吗？我的想法是，与其在无数个控制点上放置好地理参考，不如通过遍历所有公共坐标系来获得较少控制点的近乎完美的地理参考，会更容易吗？对于每个被测CRS，是否涉及每个点与其他所有点的交叉验证？ 我想了解此算法或地理配准，因此可以使流程自动化-我一直都遇到这个问题，直到内容创建者停止将其地图视为一次性创建，再也无法集成以及其他内容，我不希望停止。 参考文献 [ref：1] QGIS GCP文件： mapX,mapY,pixelX,pixelY,enable 142.632649100000009,54.453595900000003,505.941176470588232,-95.220588235293974,1 154.934252200000003,59.559921699999997,536.411764705882206,-52.779411764705742,1 80.080158100000006,9.657192300000000,291.558823529411711,-322.661764705882206,1 10.448442600000000,57.819128900000003,21.676470588235190,-103.926470588235134,1 34.007173000000002,27.761438299999998,101.117647058823422,-244.852941176470466,1 50.950890399999999,11.862196600000001,171.852941176470495,-313.955882352941046,1 29.713217199999999,60.024133200000001,90.779411764705799,-92.499999999999829,1 60.000000000000000,0.000000000000000,208.308823529411683,-362.382352941176350,1 69.867506500000005,66.639146199999999,224.088235294117567,-33.191176470588061,1 27.276107100000001,71.049154799999997,89.147058823529306,-21.764705882352814,1 140.000000000000000,0.000000000000000,536.955882352941217,-362.926470588235190,1 20.000000000000000,0.000000000000000,43.441176470588132,-362.926470588235190,1 20.196882700000000,31.243024100000000,47.249999999999901,-231.794117647058698,1 …
2
该数据在哪个投影中？
我从代理商那里收到了我的项目的shapefile，并被告知他们始终使用UTM。但是，这些值有太多数字，无法在UTM中进行投影。 min max x 32507319 32624343 y 5925974 6074881 我尝试将文件读入QGIS，但是无论我使用哪种投影，即使使用世界地图，形状也总是显示在地图边界之外。它们应出现的区域是石勒苏益格-荷尔斯泰因州，它是德国最北端的联邦州，在UTM中具有以下最小值和最大值 min max x 426167.5 650147.9 y 5913377.7 6101334.8 可以很容易地看出，高端数据出了点问题。有谁知道我可以用哪种投影来解决这个问题？
5
太平洋的适当地图投影
我从http://www.transpacificproject.com找到了这张显示太平洋的地图： 我想为地球的同一区域制作一张相似的地图，并且我想使用相同的投影，因为我认为它很好地说明了各大洲之间的距离。 我已经搜索了此投影的名称，但我还不知道答案。起初我以为是罗宾逊，但据维基百科说： 罗宾逊投影是一次显示整个世界的世界地图的地图投影。 似乎并非如此。（在这里看不到欧洲和非洲。） 因此，我的问题有两个：投影的名称是什么？我应该使用它来说明我之前所说的话还是应该使用其他内容？
1
通过纬度为Mercator计算距离比例因子
这对于Google来说确实是不方便的话题-所以我要求澄清一下，我是否正确理解了事情。 我说的是在墨卡托投影（3857）中测量短距离。您可以使用地图单位计算距离（使用毕达哥拉斯定理）。在赤道中，它等于地球上的距离（比例因子= 1）；如果您朝两极移动，则地图单位和地球上的距离不再相等-比例因子开始增加。 问题是：如何计算给定纬度的比例因子？ 我是否正确理解公式是factor = 1 / cos(latitude)？
2
使用proj4用R ggmap和ggplot2包指定Robinson投影？
我想在罗宾逊投影中投影这张地图： library(ggmap) world <- map_data("world") ggplot() + geom_path(data = world, aes(long, lat, group = group)) 并且我想将投影更改为“罗宾逊”（在回答我以前的问题的建议之后：Wikipedia使用的全球气候区域地图使用哪种投影？ 我很难找到此投影的默认实现，使用该proj4库我做了以下工作： library(proj4) robinson <- project(cbind(world$long, world$lat), proj = "+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs") 我尝试了多种方法，包括： # using ggmap::get.map() get_map("world", projection = mapprojection(robinson)) # using ggplot2::coord_map coord_map(projection = robinson) # and sp::coordinates: library(sp) …
2
ArcGIS长度和面积计算方案[关闭]
已关闭。这个问题需要更加集中。它当前不接受答案。 想改善这个问题吗？更新问题，使其仅通过编辑此帖子来关注一个问题。 4年前关闭。 尝试掌握在ArcGIS中不同方案中如何计算长度和面积。我不知道为什么我不能在要素类字段上找到答案，但是我找不到精确的答案，除非我不了解某些内容，而且我知道有一些历史。您能帮我填写问号吗？或告诉我为什么我会做错一切; ） GCS =地理坐标系PCS =投影坐标系 所有链接均指向10.1帮助文档- 要素类shape_length和shape_area字段 一个。GCS-吗？ b。PCS-使用简单平面 c。除shapefile之外，它是否总是自动更新？是 ArcMap中Measurment工具 一个。GCS-默认测地线，替代方法是Loxodrome和Great Elliptic，但不是平面的。面积计算不可用！ b。PCS-默认的平面，可选的测地线，后院和大椭圆形 http://resources.arcgis.com/en/help/main/10.1/index.html#//00s500000022000000 属性表计算器 一个。GCS-不可用 b。PCS-平面 http://resources.arcgis.com/en/help/main/10.1/index.html#//005s00000027000000 计算字段工具（数据管理工具箱） 一个。GCS-测地线性，面积可用，但有问题 b。PCS-平面 http://resources.arcgis.com/en/help/main/10.1/index.html#//00170000004m000000 缓冲工具（和其他即将 提供的工具）GCS-测地线 b。PCS-平面或指定GCS输出 http://resources.arcgis.com/en/help/main/10.1/index.html#//000800000019000000 Javascript API客户端 a。GCS-测地线面积和长度函数 b。PCS-可以从Web墨卡托转换为地理（或使用几何服务） http://help.arcgis.com/en/webapi/javascript/arcgis/help/jsapi/namespace_geometry.htm Flex API客户端 a。GCS-测地线面积和长度函数，“长度（或面积）将使用自定义圆柱等面积投影来计算”。javascript api中未提及！ b。PCS-可以从网络墨卡托转换为地理 http://resources.arcgis.com/en/help/flex-api/apiref/com/esri/ags/utils/GeometryUtil.html ArcGIS Server的REST API -几何服务 一个。GCS-测地线 b。PCS-平面 http://help.arcgis.com/zh-CN/webapi/javascript/arcgis/help/jsapi/geometryservice.htm 另一个问题，测地线究竟是什么？我以为这意味着要在球体（haversine？）上使用3D三角函数公式。在计算面积时是否太慢，这就是为什么要使用相等的面积投影吗？ 确定长度和面积时的另一个问题是，等面积投影比使用相同基准球体的测地线计算更准确吗？简要为什么？
2
使用arcpy.RasterToNumPyArray保持空间参考？
我正在使用ArcGIS 10.1，并想基于两个现有的栅格创建一个新栅格。该RasterToNumPyArray有，我要适应一个很好的例子。 import arcpy import numpy myArray = arcpy.RasterToNumPyArray('C:/data/inRaster') myArraySum = myArray.sum(1) myArraySum.shape = (myArray.shape[0],1) myArrayPerc = (myArray * 1.0)/ myArraySum newRaster = arcpy.NumPyArrayToRaster(myArrayPerc) newRaster.save("C:/output/fgdb.gdb/PercentRaster") 问题在于它会剥夺空间参考以及像元大小。我认为它必须做arcpy.env，但是如何根据输入栅格设置它们呢？我想不明白。 接受卢克的回答，这是我的解决方案。 Luke的两种解决方案都可以正确设置空间参考，范围和像元大小。但是第一种方法不能正确地在数组中传送数据，并且输出栅格到处都填充有nodata。他的第二种方法大多数情况下都有效，但是在我的无数据区域很大的地方，它填充有块状零和255。这可能与我处理nodata单元的方式有关，而且我不太确定自己的操作方式（尽管应该是另一个Q）。我包括了我正在谈论的图像。 #Setting the raster properties directly import arcpy import numpy inRaster0='C:/workspace/test0.tif' inRaster1='C:/workspace/test1.tif' outRaster='C:/workspace/test2.tif' dsc=arcpy.Describe(inRaster0) sr=dsc.SpatialReference ext=dsc.Extent ll=arcpy.Point(ext.XMin,ext.YMin) # sorry that i modify calculation …
2
确定纸质地图的投影？
如果我正在阅读纸质地图，如果未在文档上明确标记此类信息，如何确定用于绘制地图的投影？ 在网上很容易找到描述和分类使用中的各种预测的资源（例如Wikipedia或GIS Wiki）。但是，作为非专业的地图用户，我发现很难反向使用构造细节来确定使用哪种投影来绘制特定的地图。我认为也许会有流程图或其他分类过程可以帮助地图用户确定使用中的投影？ 例如，在我当地的图书馆和档案馆中，有1878年加拿大自治领地办公室（Dominion Lands Office，加拿大）绘制的地图石版画的副本。该地图具有弯曲的平行线和子午线，其底部比顶部宽（顶部更接近在这种情况下是极点）。我认为这使地图成为圆锥形投影，但是我不确定。而且，如果这是正确的，那么我有兴趣研究是否可以对投影进行更精确的分类（例如，Lambert保形的）。 我对这张1878年的地图进行了部分扫描（由于太大，无法用我拥有的设备扫描整个物体），以防有人回答我的问题。
1
椭圆形立体投影中如何使用“共形纬度”？
在我的工作中，我一直在对内部软件进行一些维护，这些软件声称实现了“立体投影”，其公式是很久以前从旧系统中复制的。 从我所看到的，它似乎实现了EPSG的“ Oblique and Equatorial Stereographic”（代码9809）。但是有一个步骤我无法确定。我们的代码中使用的“ R”值是使用以下参数计算的： (phi_n, lambda_n) - null distortion coordinates (PHI_nc, LAMBDA_nc) - conformal null distortion coordinates (PHI_o, LAMBDA_o) - conformal projection origin coordinates e - eccenttricity r_eq - earth radius at equator 鉴于此，获得R的公式为： a = 1 + cos(PHI_nc) * cos(PHI_o) * cos(LAMBDA_nc - LAMBA_o) b = …
2
如何在QGIS中启用投影转换
我的QGIS项目分为两层，一层由Google图层插件检索，该插件获取CRS为+ proj = merc + lon_0 = 0 + lat_ts = 0 + x_0 = 0 + y_0 = 0 + a = 6378137 +的静态Google地图图像b = 6378137 + units = m + no_defs，另一个向量层是postgis层，其CRS为+ proj = longlat + ellps = WGS84 + datum = WGS84 + no_defs。 我将Google图层中的某些多边形要素数字化，并将其添加到矢量图层中。但是，QGIS似乎没有进行自动CRS转换。如何更改投影？ 谢谢
7
使用ArcPy获取投影点的经度和纬度？[关闭]
关闭。这个问题是题外话。它当前不接受答案。 想改善这个问题吗？ 更新问题，使其成为地理信息系统堆栈交换的主题。 2年前关闭。 我在ArcPy正在访问的要素类中有一个点要素。该点是已投影的，但我需要找到一种有效的方法来获取该点的未投影的纬度和经度。 除了重新投影（取消投影），在新要素类上找到搜索光标，找到要素，然后从要素形状中获取经纬度之外，还有其他方法吗？
By using our site, you acknowledge that you have read and understand our Cookie Policy and Privacy Policy.
Licensed under cc by-sa 3.0 with attribution required.