Questions tagged «coordinate-system»

就像标题所说的那样，有人会如何配置Mapnik以使用Waterman蝴蝶投影？ 否则，还有哪些其他工具可以使用此投影进行渲染？

11 coordinate-system  mapnik 

我在Windows上运行QGIS 1.7.0，并且从TIGER中下载了人口普查区块组边界的shapefile 。该文件似乎在WGS 84坐标系中，当我在QGIS中打开它时，度量单位是度。 我尝试将文件重新投影到NAD 83坐标系中，并将其另存为新的shapefile，但是当我打开新文件时，它仍然显示度尺。如果我尝试制作一个缓冲区，它将缓冲区距离解释为度。 假设我从TIGER的shapefile开始，将度量单位从度转换为米的最佳方法是什么？ 为了后代，在处理美国人口普查TIGER shapefile时，以下工作流程似乎很有意义： 下载相关的shapefile 在QGIS中打开 点击Settings菜单，然后选择project properties 选择一个投影坐标系，例如NAD83 / UTM区域19N 点击框 Enable 'on the fly' CRS transformation 点击确定 点击“放大全屏”按钮（其中一个放大镜）

11 qgis  shapefile  coordinate-system  tiger 

我在WGS84纬度/经度中有积分，我想测量它们之间的“小”（小于5公里）距离。 我可以使用http://www.movable-type.co.uk/scripts/latlong.html中的Haversine公式，效果很好。 不过，我想使用Python Shapely库，这样我不仅可以进行距离操作，而且可以进行更多的操作，并且因为在我正在使用的规模上，平整的地球足够近似。为了将地理坐标可靠地投影到笛卡尔坐标，我使用的是Python的proj4，但是似乎会得到比我想要的更大的错误。 如果我使用本地UTM区域，则会出现两米长的正弦波之间的差异，这很好。但是我不想解决UTM区域（这些点可能在世界范围内），因此我尝试使用“球形墨卡托”，但现在正弦距离与预计距离之间的差异远超过100％。这对球形墨卡托真的合适吗？我真正想要的是一个可行的笛卡尔投影，用于在世界任何地方彼此之间5公里之内的两个点。 from shapely.geometry import Point from pyproj import Proj proj = Proj(proj='utm',zone=27,ellps='WGS84') #proj = Proj(init="epsg:3785") # spherical mercator, should work anywhere... point1_geo = (-21.9309694, 64.1455718) point2_geo = (-21.9372481, 64.1478206) point1 = proj(point1_geo[0], point1_geo[1]) point2 = proj(point2_geo[0], point2_geo[1]) point1_cart = Point(point1) point2_cart = Point(point2) print "p1-p2 (haversine)", …

11 python  coordinate-system  distance 

我有一个.csv文件，其中包含一长串的纬度和经度点。有没有一种有效的方法可以将列表批量转换为预计的UTM坐标系？

11 convert  coordinate-system  utm 

是否可以在坐标系统中发送一个点并返回在指定坐标系中重新投影的点的免费公共服务？ 我不会做很多请求，因为现在只是一个演示。

11 coordinate-system  point 

在Wikipedia页面上查看北美基准面，似乎暗示该基准面完全由椭球体定义，椭球体定义了长/半长轴，展平等。对于大地坐标，则是一个基准和一个椭球体定义同样的事情？这些术语可以互换使用吗？ 例如，NAD83是否完全由椭球体GRS80定义？

11 coordinate-system  datum  gis-principle  nad83  ellipsoid 

北方和东方的定义很容易理解，但是只有与xy坐标可互换使用时才变得困难，xy坐标对轴的方向具有不同的定义。在数学中，y始终是垂直的，x始终是水平的，因此从逻辑上讲，我会假设“ up” ==“ northing” ==“ y”和“ along” ==“ easting” ==“ x”。 为什么在GIS中不是这种情况？

11 coordinate-system 

我正在尝试使用gdalwarp将对地静止映射为Lambert保形。我的输入数据在netcdf中，并且在地理坐标（度）中，我想将重新映射的数据输出到netcdf。我为输入的netcdf数据创建了一个相应的vrt文件。Gdalwarp将输出netcdf文件，但输出数据全为零，并且我收到以下错误： Creating output file that is 5120P x 5120L. Processing input file netcdf.vrt. ERROR 1: Too many points (441 out of 441) failed to transform, unable to compute output bounds. Warning 1: Unable to compute source region for output window 0,0,5120,5120, skipping. 0...10...20...30...40...50...60...70...80...90...100 - done. 我尝试了以下命令： /usr/bin/gdalwarp -s_srs "+proj=geos +h=35785831 …

11 gdal  coordinate-system  proj  gdalwarp  netcdf 

我想知道是否有人可以澄清使用ArcGIS IDW空间插值时是否应将点数据投影到等距投影？ 我正在研究来自北美西部的数据集，其跨度约为30度。数据当前为纬度/经度（NAD83）。 在插值过程中计算到样本点的距离时，ArcGIS是“即时投影”还是以某种方式调整了纬度？或者我应该在保留距离的投影中提供所有内容吗？

11 arcgis-desktop  coordinate-system  interpolation 

我想将tiff光栅文件从投影A（ETRS89 / UTM区域32）重新投影到投影B（DHDN / Gauss-Kruger区域2）。这两个投影的基准不同，因此我必须使用基准转换（本例中为NTv2 BETa2007）。我想用投影B创建一个新的栅格文件，因为我想在另一个软件中使用它。因此，即时投影并不是我所需要的。 我尝试了“栅格”>“变换”（或类似方法-我不知道英语菜单），但是在这里看不到如何定义基准变换。 如何在Qgis中使用基准转换对栅格文件进行永久性重新投影？

11 qgis  raster  coordinate-system  datum 

我正在编写一个用于测试各种地图服务（主要是AGS 9.x，AGS 10和WMS 1.x）的应用程序。 该应用程序的一部分涉及在服务的整个范围内为单个请求生成随机边界框。当已知服务的全部范围时（例如，通过AGS服务的fullExtent属性），此部分可以很好地用于地理坐标系和投影坐标系。 我的问题是WMS：GetCapabilities响应中的每一层都可以在> = 1个CRS中定义其边界区域。应用程序的某些部分需要知道服务的CRS是地理的还是投影的，因此为了消除WMS中的歧义，我始终使用始终在EPSG：4326中定义的图层的LatLonBoundingBox。然后，我必须根据进入单个请求（随机化）的所有层来计算完整的服务边界框。这是棘手的地方。 我迷路了，因为对于每个纬度/经度边界框，LLx（左下经度）可能大于或小于URx（右上经度），具体取决于它跨越的子午线。每次我开始绘制方形或圆形图时，我都会想出一种方法，然后找到一个破坏它的案例，然后我的大脑变得糊涂。 我将继续努力，直到它起作用为止，如果我在此处发布了解决方案，但我敢肯定，必须有一种公认的且经过充分测试的方法，这将使我的生活更轻松。我只是现在找不到。

11 wms  coordinate-system  latitude-longitude  extents 

在对地图投影进行研究之后，我意识到很难选择一种能够同时计算大量点之间的距离的投影（如果这些点位于世界各地）。等距投影仅允许计算两点之间的距离。然后，我决定使用以下解决方案。我的数据集可以拆分为较小的数据集： 美国领土； 整个EURASIA； 涵盖非洲大陆和欧洲的领土。 在ArcGIS 9.3（ArcInfo许可）中，我对每个数据集进行了兰伯特方位角相等面积投影，以针对每种情况正确地“居中”。对于每个数据集，我使用天梭的Indicatrix shapefile加载了一个新层，并计算了每个椭圆的面积。我认为椭圆形区域的面积与距离投影中心最近的椭圆形相差不超过5％的所有领土（在面积和距离方面）都“不失真”与通函没有任何合理的区别。在看到我的点被包围在这个“不变形”的区域之后，我计算了点距离。您是否认为我的程序合理正确，或者对我的目的有更好的地图投影建议？

11 coordinate-system  distance 

我正在研究所涉及的数学，以了解如何将给定的十进制纬度和经度转换为EPSG：3857。任何人都可以指出一个好的参考文献或可能的Ppython / C / C ++ / Java库源代码来检查数学吗？

11 python  coordinate-system  convert  c++  mathematics 

我在PostGIS中选择了我的桌子（实地） id_borne serial NOT NULL, num_borne character varying, shape_borne geometry 我成功使用更改了表的SRID select UpdateGeometrySRID('public', 'borne', 'shape_borne', 26191) ; 旧的SRID为0。 问题是，当我尝试向表中添加新元素时，borne出现此消息错误： geometry srid (0) does not match column srid (26191) 那是我的问题！SRID已更新，但是当我尝试添加新点时，出现此消息。

11 postgis  coordinate-system  srid 

是否可以在QGIS中使用Goode Homolosine投影？QGIS和GRASS似乎都遇到干扰，将较高的纬度挤压在一起。我尝试+ proj = goode + ellps = sphere + datum = wgs84作为自定义CRS。

11 qgis  coordinate-system  grass  homolosine