Questions tagged «proj»

我有一个Geotiff栅格图像，它的坐标系的纬度从0到360。图像的水平中心是180度。见下图： 我想将其转换为经度范围为-180 180的EPSG：4326 SRS。我希望图像的中心位于格林威治子午线（0）。我想这个SRS用途非常广泛。我希望结果看起来像这样： 因此，我使用gdalwarp命令重新投影： gdalwarp -s_srs '+proj=latlong +datum=WGS84 +pm=180dW' -t_srs EPSG:4326 test_col.tif test_4326.tif 但是我只会看到尺寸更大（像素更多）和EPSG：4326元数据的tiff。图像本身看起来与初始图像相同。但我希望它能交换半球。 问题是-如何将图像的经度严格调整为-180 180 EPSG：4326，且中心在0度的经度？ 这是我的初始文件的gdalinfo： Origin = (-0.102272598067084,89.946211604095552) Pixel Size = (0.204545196134167,-0.204423208191126) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left ( -0.1022726, 89.9462116) ( 0d 6' 8.18"W, 89d56'46.36"N) Lower Left ( -0.1022726, -89.9462116) ( …

31 coordinate-system  gdal  geotiff-tiff  proj  gdalwarp 

我想知道在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 

我有一个多边形的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 

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

14 qgis  coordinate-system  proj 

某些shapefile具有与其关联的.prj文件，并且.prj文件包含WKT格式的shapefile的投影信息。有时我需要将WKT转换为proj4字符串，有时我需要将其转换回。 有没有现成的图书馆可以做到这一点？

14 shapefile  coordinate-system  wkt  proj 

我解决了库问题，但仍然存在，仅与未加载的PROJ.4库有关。我问了这些问题的后续问题。 我的GDAL软件包和PROJ.4库已损坏，这在尝试重新投影文件时会导致错误（请参见下面的旧问题）。我不知道在我错误地安装了另一个GDAL实例之前还是之前，这个问题是否开始；但是，我已经将其卸载并重新安装了QGIS，问题仍然存在。 每次重新安装后，QGIS都具有与以前相同的插件，因此导致错误的部分在安装/卸载过程中保持不变。 我尝试了什么： sudo apt-get remove build-essential python-all-dev并sudo apt-get remove gdal-bin清除手动安装的GDAL 通过Ubuntu软件中心安装（un）安装QGIS sudo apt-get --purge autoremove dpkg -l | grep '^rc' | awk '{print $2}' | sudo xargs dpkg --purge并dpkg --list | grep "^rc" | cut -d " " -f 3 | xargs sudo dpkg --purge清除配置文件 仍然没有运气。我发现了一些似乎相关的东西，但这是针对Windows的，而不是针对我的Ubuntu 14.04的。在下一次安装QGIS之前，还应该清除什么（以及如何清除）？ 原始问题： …

13 installation  ogr2ogr  ubuntu  proj 

我想在罗宾逊投影中投影这张地图： 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) …

13 coordinate-system  r  proj  ggmap 

我正在尝试转换WGS84纬度/经度点 11d42'32.10629" E 5d12'56.75371" S 在三个不同的软件包（Proj4，GeoTrans和Leica GeoOffice）中，但是结果之间有很大的差异-X和Y大约有几百米！请注意，这些变化仅在7个参数中出现，而在3个参数中不出现。 Proj4 cs2cs + proj = longlat + ellps = WGS84 + datum = WGS84 + no_defs + to + proj = utm + zone = 32 + ellps = clrk80 + towgs84 = 178.3,316.7,131.5，-5.278，-6.077，-10.9792，-19.166 +南+ units = m + no_defs GeoTrans Delta X …

12 coordinate-system  proj  datum 

我正在尝试使用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 

我试图按照文档中所述使用proj4孔在QGIS中定义ortelius椭圆形投影和Nicolosi球状投影，但是即使proj4库中支持它们，QGIS也无法投影到这些投影中。我需要那些投影来进一步分析旧的地图。 有人知道可能是什么问题吗？

10 qgis  proj 

我正在编写一个python脚本，该脚本读取包含x和y坐标的多个XML文件，并将它们全部合并为一个csv文件。纬度和经度是csv中的必填字段，但我很难将俄亥俄州北州飞机usFt中的x，y坐标转换为WGS84。 >>> p = Proj(r'+proj=lcc +lat_1=41.7 +lat_2=40.43333333333333 +lat_0=39.66666666666666 +lon_0=-82.5 +x_0=600000 +y_0=0 +ellps=GRS80 +datum=NAD83 +units=us-ft +no_defs') #Nad83 State Plane Ohio North US Feet Proj object using parameters >>> p(739400.91,2339327.3,inverse=True) (-80.138057868777224, 60.278230707978487) >>> p1 = Proj(init="epsg:3734") #Nad83 State Plane Ohio North US Feet Proj object using EPSG code >>> p1(739400.91,2339327.3,inverse=True) (-80.138057868777224, 60.278230707978487) …

10 python  coordinate-system  wgs84  proj 

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

9 python  coordinate-system  gps  proj  ground-control