คำถามติดแท็ก proj

ฉันมีภาพแรสเตอร์เชิงพิกัดที่มีระบบพิกัดที่มีความยาวตั้งแต่ 0 ถึง 360 จุดศูนย์กลางแนวนอนของภาพคือ 180 longtitude ดูภาพด้านล่าง: ฉันต้องการแปลงเป็น EPSG: 4326 SRS ด้วยช่วงความยาว -180 180 และฉันต้องการให้ศูนย์กลางของภาพอยู่ที่ Greenwich meridian (0) ฉันคิดว่า srs นี้ใช้กันอย่างแพร่หลายมาก ฉันหวังว่าผลลัพธ์จะเป็นดังนี้: ดังนั้นฉันใช้คำสั่ง gdalwarp เพื่อปฏิเสธ: gdalwarp -s_srs '+proj=latlong +datum=WGS84 +pm=180dW' -t_srs EPSG:4326 test_col.tif test_4326.tif แต่ฉันจะได้รับความยุ่งยากที่มีขนาดใหญ่ขึ้น (พิกเซลมากขึ้น) และ EPSG: 4326 เมทาดาทา ภาพมีลักษณะเหมือนกันกับภาพแรก แต่ฉันคาดหวังว่ามันจะแลกเปลี่ยนซีกโลก คำถามคือ - ฉันจะถ่ายภาพให้เป็น -180 180 EPSG ได้อย่างไร: …

31 coordinate-system  gdal  geotiff-tiff  proj  gdalwarp 

ฉันสงสัยว่าอะไรคือวิธีที่ถูกต้องในการผสานการแปลง NTv2 ใน QGIS ฉันประสบความสำเร็จในการรวมการแปลง NTv2 ใน QGIS แต่ฉันไม่แน่ใจว่ามันได้ทำตามที่ตั้งใจไว้หรือไม่ สำหรับการรวม NTV2-grids ใน QGIS ฉันพบข้อมูลหลายอย่างที่ชื่อของกริดควรถูกเก็บไว้ในคอลัมน์ p1 ของตาราง tbl_datum_transform ของ srs.db: 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 (ภาษาเยอรมัน)) อะไรทำให้ฉันสงสัยว่าถูกต้องหรือไม่ความจริงที่ว่าคอลัมน์ดูเหมือนจะมีประเภทฟิลด์ตัวเลข: พารามิเตอร์แรกของการแปลงเจ็ดพารามิเตอร์ไม่ถูกเก็บไว้ในคอลัมน์ p1 ของตาราง tbl_datum_transform ของ QGIS หรือไม่ ตามคำถามที่พบบ่อย sqlite เป็นไปได้ที่จะจัดเก็บสตริงลงในคอลัมน์ตัวเลขใน sqlite ( https://sqlite.org/faq.html#q3 ) แต่ดูเหมือนว่าฉันจะชอบแฮ็คหรือวิธีแก้ปัญหาเพื่อเก็บสตริงลงในคอลัมน์ตัวเลข . นี่เป็นวิธีที่ควรทำจริงๆหรือ? tl; dr: วิธีที่ถูกต้องในการรวมการแปลง NTv2 ใน QGIS คืออะไร?

15 qgis  coordinate-system  proj  epsg 

ฉันติดตั้ง QGIS พร้อม proj 4.8.0 แล้ว ฉันต้องการเพิ่มการประมาณการ Winkel Tripel ซึ่งดูเหมือนจะถูกกำหนดเป็น +proj=wintri แต่สตริงคำจำกัดความนั้นไม่ทำงาน นี่เป็นข้อบกพร่องในรุ่น proj ของฉันหรือสตริงคำจำกัดความที่ถูกต้องคืออะไร?

15 qgis  coordinate-system  proj 

ฉันมีรูปร่างเป็นรูปหลายเหลี่ยมและไฟล์ CSV อีกไฟล์หนึ่งซึ่งมีรายการคะแนนเป็นคู่ (Lat, Lng) ฉันต้องการตรวจสอบแต่ละคู่ (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]] แผนของฉันเป็นดังนี้: อ่าน shapefile โดยใช้readShapePolyฟังก์ชั่นในMapToolsแพ็คเกจR อ่านพิกัดจุดจากไฟล์ CSV ลงใน dataframe และแปลงเป็น SpatialPointsDataFrame ใช้overฟังก์ชันเพื่อกำหนดรูปหลายเหลี่ยมที่มันอยู่ข้างใน ในการดำเนินการดังกล่าวฉันต้องระบุproj4stringขณะที่โหลดรูปร่างไฟล์ในขั้นตอนที่ 1 และเปลี่ยนพิกัดจากไฟล์ CSV เป็นระบบฉายภาพเดียวกันโดยใช้spTransformฟังก์ชั่นก่อนที่จะใช้overฟังก์ชั่นในขั้นตอนที่ 3 เนื่องจากต้องการจุดและรูปหลายเหลี่ยม อยู่ภายใต้ระบบการฉายเดียวกัน แนวคิดใดเกี่ยวกับค่าที่ถูกต้องสำหรับเนื้อหาไฟล์ proj ที่แสดงด้านบน?

14 coordinate-system  shapefile  r  proj 

บาง Shapefiles มีไฟล์. prj ที่เชื่อมโยงกับมันและไฟล์. prj มีข้อมูลการฉายของ Shapefile ในรูปแบบของ WKT บางครั้งฉันต้องแปลง 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 ที่ติดตั้งด้วยตนเอง (ยกเลิก) ติดตั้ง QGIS ด้วย Ubuntu Software Center sudo apt-get --purge autoremove dpkg -l | grep …

13 installation  ogr2ogr  ubuntu  proj 

ฉันต้องการฉายแผนที่นี้ในการฉายภาพของ Robinson: library(ggmap) world <- map_data("world") ggplot() + geom_path(data = world, aes(long, lat, group = group)) และฉันต้องการเปลี่ยนการฉายเป็น "Robinson" (ทำตามคำแนะนำจากคำตอบของคำถามก่อนหน้านี้: การคาดการณ์ใดที่แผนที่ภูมิอากาศโลกจาก 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 …

13 coordinate-system  r  proj  ggmap 

ฉันพยายามที่จะแปลง WGS84 lat / lon point 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 + ถึง + proj = utm + zone …

12 coordinate-system  proj  datum 

ฉันกำลังพยายามแมปจาก Geostationary ไปยัง Lambert โดยใช้ gdalwarp ข้อมูลอินพุตของฉันอยู่ใน netcdf และอยู่ในพิกัดทางภูมิศาสตร์ (องศา) และฉันต้องการส่งออกข้อมูลที่แมปไปยัง netcdf ฉันได้สร้างไฟล์ vrt ที่สอดคล้องกันสำหรับการป้อนข้อมูล netcdf 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 …

11 gdal  coordinate-system  proj  gdalwarp  netcdf 

ฉันกำลังพยายามกำหนด ortelius oval projection และ Nicolosi globular projection ใน QGIS โดยใช้ proj4 foramt ตามที่อธิบายไว้ในเอกสารประกอบ แต่ QGIS ไม่สามารถฉายภาพในการฉายภาพเหล่านี้แม้ว่าพวกเขาจะได้รับการสนับสนุนในห้องสมุด proj4 ... แผนที่ ไม่มีใครรู้ว่าสิ่งที่อาจเป็นปัญหาหรือไม่

10 qgis  proj 

ฉันกำลังเขียนสคริปต์หลามที่อ่านไฟล์ XML หลายไฟล์ที่มีพิกัด x และ y และรวมพวกมันทั้งหมดไว้ในไฟล์ csv ไฟล์เดียว ละติจูดและลองจิจูดเป็นฟิลด์บังคับใน csv แต่ฉันมีปัญหาในการแปลงพิกัด x, y ใน Ohio North State Plane usFt เป็น 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, …

10 python  coordinate-system  wgs84  proj 

ฉันมีคลาวด์พอยท์ซึ่งมีพิกัดที่เกี่ยวกับระบบพิกัดในพื้นที่ ฉันยังมีจุดควบคุมภาคพื้นดินพร้อมค่า GPS ฉันสามารถแปลงพิกัดท้องถิ่นเหล่านี้เป็นระบบพิกัดทั่วโลกโดยใช้ PROJ.4 หรือห้องสมุดอื่น ๆ ได้หรือไม่? รหัสใด ๆ ใน Python สำหรับปัญหาที่ระบุข้างต้นจะเป็นความช่วยเหลือที่ดี

9 python  coordinate-system  gps  proj  ground-control