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

บาง Shapefiles มีไฟล์. prj ที่เชื่อมโยงกับมันและไฟล์. prj มีข้อมูลการฉายของ Shapefile ในรูปแบบของ WKT บางครั้งฉันต้องแปลง WKT เป็นสตริง proj4 และบางครั้งฉันต้องแปลงกลับ มีห้องสมุดสำเร็จรูปที่จะทำเช่นนี้หรือไม่?

14 shapefile  coordinate-system  wkt  proj 

ฉันต้องการนำเข้าไฟล์ shp ไปยัง postgreSQL ก่อนอื่นฉันสร้างไฟล์ sql แล้วเรียกใช้ PostgreSQL เพื่อสร้างไฟล์ sql ฉันเรียกใช้คำสั่งนี้ใน windows cmd: shp2pgsql -s 4326 worldCountries.shp worldcountries postgres > worldcountries.sql จากนั้นเรียกใช้: psql -d postgres -U postgres -p 4321 -f worldcountries.sql แต่ผลลัพธ์คือ: psql: worldcountries.sql: 21: ข้อผิดพลาด: ฟังก์ชัน addGeometrycolumn (ไม่ทราบ, ไม่ทราบ, ไม่ทราบ, ไม่ทราบ, ไม่ทราบ, จำนวนเต็ม)

14 postgis  shapefile  postgresql  import 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน3 ปีที่ผ่านมา มีแบบฝึกหัดที่ดีที่อธิบายในเชิงลึกภายในของ GIS และระบบฐานข้อมูลเชิงพื้นที่เช่น PostGIS พร้อมตัวอย่างและไม่มีข้อสันนิษฐานเบื้องหลังหรือไม่? ฉันกำลังมองหาคำตอบสำหรับคำถามเช่น: เรขาคณิตใน GIS คืออะไร? มันเป็นตัวแทนอย่างไร เมื่อใช้พิกัดละติจูด / ลองจิจูดการดำเนินการที่ฉันต้องดำเนินการนั้นจะต้องนำไปสู่สถานะที่ฉันสามารถเรียกชื่ออย่างเช่น ST_contain PostGIS เพื่อตรวจสอบว่ารูปหลายเหลี่ยมนั้นมีพิกัดละติจูด / ลองจิจูดหรือไม่ ฉันมีประสบการณ์เกี่ยวกับระบบฐานข้อมูลเชิงสัมพันธ์ทั่วไป แต่ไม่มีภูมิหลังเกี่ยวกับ GIS หรือกระบวนทัศน์ฐานข้อมูลเชิงพื้นที่

13 postgis  shapefile  geometry  spatial-database  references 

ฉันได้ต่อสู้กับการเพิ่มเลเยอร์ใน TOC มาระยะหนึ่งแล้ว ฉันจัดการเพื่อเพิ่มเลเยอร์เดียว แต่จำเป็นต้องวนซ้ำผ่านโฟลเดอร์และเพิ่มเลเยอร์ทั้งหมด Shapefiles ฉันไม่สามารถเพิ่มได้เลย ไฟล์เลเยอร์เท่านั้น ด้านล่างเป็นรหัสที่ฉันอยู่: for layer in shp_List: mxd = arcpy.mapping.MapDocument("CURRENT") dataFrame = arcpy.mapping.ListDataFrames(mxd, "*")[0] outlayer = lyr layerfile = os.path.join(base_Folder, lyr + ".lyr") arcpy.MakeFeatureLayer_management(layer, outlayer) arcpy.SaveToLayerFile_management(outlayer, layerfile, "ABSOLUTE") addlayer = arcpy.mapping.Layer(layerfile) arcpy.mapping.AddLayer(dataFrame, addlayer, "BOTTOM") #addLayer = arcpy.mapping.Layer(layer) #arcpy.mapping.AddLayer(dataFrame, addLayer, "BOTTOM") #arcpy.RefreshTOC() #arcpy.RefreshActiveView() del addlayer, …

13 arcgis-desktop  arcgis-10.0  arcpy  shapefile  layers 

มีคำจำกัดความว่าเส้นเชื่อมต่อจุดยอดในรูปแบบ shapefile หรือไม่? ในกรณีที่ง่ายที่สุดลองนึกภาพเส้นที่มีเพียง 2 คะแนนจาก 40, -118 ถึง 40, -112 - นั่นคือที่สุ่มในสหรัฐอเมริกา - ด้วยระบบพิกัดทางภูมิศาสตร์ WGS84 มาตรฐาน นี่คือเนื้อหาของไฟล์. prj: GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]] ประเด็นคือ 40.1, -116 เหนือหรือใต้ของเส้น? หากเราพิจารณาว่าเส้นถูกสอดแทรกเชิงเส้นในพื้นที่ latlong มันจะเป็นไปตาม 40 องศาขนาน (วงกลมเล็ก) และจุดนั้นอยู่ทางทิศเหนือของเส้น หากเราพิจารณาว่าเส้นนั้นเป็นเส้นทางที่สั้นที่สุดบนพื้นผิวโลกมันเป็นมาตราส่วน (วงกลมใหญ่) ที่มีละติจูดสูงสุดที่กึ่งกลางบรรทัดสูงกว่า 40.1 องศา จากนั้นจุดจะอยู่ทางใต้ของเส้น หรือมันเป็นเพียงแค่ไม่ได้กำหนด? เนื่องจากรูปแบบ shapefile ไม่มีความคิดของเส้นโค้ง แต่มีเพียงเส้นตรงที่เชื่อมต่อกันเท่านั้น บรรทัดจะต้องมีความหนาแน่น (เพิ่มคะแนนตามแนวเส้น) เพื่อชี้แจงคำตอบนี้ ถ้าฉันสร้างสถานการณ์ดังกล่าวใน QGIS เส้นนั้นจะขนานกันไป 40 องศาและจะบอกคำตอบให้ฉันคือ …

13 shapefile  great-circle  geodesic 

ฉันจะใช้ไลบรารีฟรีเพื่อแปลงข้อมูลเชิงพื้นที่ได้อย่างไร ตัวอย่างเช่นฉันต้องการเปลี่ยนการฉายภาพของ Shapefile ภายในรหัสของแอปพลิเคชันเว็บ C # ของฉัน ฉันจะทำอย่างไร

13 shapefile  coordinate-system 

ฉันมีรูปร่างหลายร้อยไฟล์จากโครงการ GIS ต่างๆที่ฉันต้องการเริ่มรวมเข้ากับแพลตฟอร์มฐานข้อมูลเดียวซึ่งขณะนี้กำลังพยายามทำสิ่งนี้กับ Postgres / PostGIS แทบจะไม่มีข้อมูลใด ๆ ที่ได้มาตรฐาน - หมายความว่ามันเป็นประเภทข้อมูลเดียวกันจำนวนมากแต่ชื่อ / ประเภทของแอตทริบิวต์เฉพาะไม่ตรงกัน ฉันจะเริ่มแก้ไขปัญหานี้ได้ที่ไหน ฉันควรพัฒนาแบบจำลองมาตรฐานเพื่อโอนย้ายรูปร่างแต่ละไฟล์เป็นอันดับแรก (เช่น Hydro_line, transport_line, มาตรฐาน Hydro_poly ฯลฯ ) หรือไม่ อีกทางเลือกหนึ่งคือการนำเข้าแต่ละ Shapefile ลงใน Postgres แยกจากกันดังนั้น shp แต่ละอันจะกลายเป็นตารางในฐานข้อมูล แต่ฉันไม่แน่ใจเกี่ยวกับเรื่องนี้ในแง่ของประสิทธิภาพและการจัดระเบียบ รู้สึกถึงความล่าช้าอย่างหลีกเลี่ยงไม่ได้ ... มีคำแนะนำใดเกี่ยวกับการจัดการกับภารกิจที่น่ากลัวนี้หรือไม่?

13 postgis  shapefile  postgresql  spatial-database  spatial-etl 

คำตอบของคำถามนี้เป็นความพยายามของชุมชน แก้ไขคำตอบที่มีอยู่เพื่อปรับปรุงโพสต์นี้ ขณะนี้ไม่ยอมรับคำตอบหรือการโต้ตอบใหม่ ฉันจะหา shapefiles สำหรับรหัสไปรษณีย์ของแคนาดาได้ที่ไหน

13 shapefile  data  canada  postal-code 

ฉันต้องการใช้ ogr2ogr เพื่อนำเข้า Shapefile ในฐานข้อมูล postgis ฉันติดตั้ง ogr2ogr สำเร็จแล้วและฉันรันจากคำสั่ง pgsql ต่อไปนี้: ogr2ogr -f "PostgreSQL" PG:"host=localhost user=user_1 password=***** dbname=imports" world_boundaries.shp สิ่งที่ฉันได้รับกลับมาคือข้อความแสดงข้อผิดพลาด: Unable to open datasource `world_boundaries.shp' with the following drivers: --a list of drivers follows (ESRI Shape File etc.) ฉันได้พยายามกำหนดเส้นทางแบบเต็มของ shapefile แต่ฉันได้รับข้อความเดียวกัน นอกจากนี้ฉันพยายามเรียกใช้: ogrinfo world_boundaries.shp สิ่งเดียวกัน หลังจากแก้ไขปัญหาด้วยสิทธิ์ของไฟล์ฉันได้รับข้อผิดพลาดดังต่อไปนี้: ERROR 1: AddGeometryColumn failed …

13 shapefile  postgresql  ogr2ogr 

ฉันต้องเปลี่ยนมันเป็น UTM เพื่อให้บัฟเฟอร์ทำงานได้ wets<-readOGR(dsn=".",layer="shapefile") r.raster <- raster() extent(r.raster) <- extent(wets) res(r.raster) <- 100 wets.r <- rasterize(wet,r.raster) plot(wets.r) wetsbuf<-buffer(wets.r,width=500) ในระหว่างการสร้างบัฟเฟอร์ซึ่งเป็นบรรทัดสุดท้ายของรหัสมันให้คำเตือนนี้: Warning message: In couldBeLonLat(x) : raster has a longitude/latitude CRS, but coordinates do not match that นี่คือข้อมูล summary(wets.r) layer Min. 1 1st Qu. 1 Median 2 3rd Qu. 9 Max. 11 …

13 raster  coordinate-system  shapefile  r 

ฉันใหม่สำหรับ R และใช้แพ็คเกจแรสเตอร์ ฉันมีปัญหาในการแยกรูปหลายเหลี่ยมจากไฟล์แรสเตอร์ที่มีอยู่ ถ้าฉันใช้ extract(raster, poly_shape) ฟังก์ชั่นบนแรสเตอร์มันมักจะสร้างรายการที่มีข้อมูล สิ่งที่ฉันต้องการคือการแยกไฟล์แรสเตอร์อื่นที่ฉันสามารถโหลดด้วย ArcGIS อีกครั้ง หลังจากอ่านอีกเล็กน้อยฉันคิดว่าฟังก์ชั่นการครอบตัดคือสิ่งที่ฉันต้องการจริงๆ แต่เมื่อลองใช้ฟังก์ชั่นนี้ crop(raster, poly_shape) ฉันได้รับข้อผิดพลาดนี้: Error in .local(x, y, ...) : extents do not overlap In addition: Warning message: In intersect(extent(x), extent(y)) : Objects do not overlap ไฟล์ raster และ poly_shape เหมือนกันทั้งสองฟังก์ชั่น คุณช่วยบอกฉันทีว่าอะไรที่นี่ผิด มันถูกต้องไหมที่ฟังก์ชั่นการครอบตัดสร้างอีกแบบแรสเตอร์และไม่ใช่รายการ? แก้ไข : ฟังก์ชั่นขอบเขต () ไม่ทำงานสำหรับฉัน …

13 raster  shapefile  r  line 

ฉันได้ลองหลายตัวอย่างของรหัสโดยใช้ไลบรารีเช่น shapefile, fiona และ ogr เพื่อพยายามตรวจสอบว่ามีจุด (x, y) อยู่ในขอบเขตของรูปหลายเหลี่ยมที่สร้างด้วย ArcMap (และในรูปแบบ shapefile) อย่างไรก็ตามไม่มีตัวอย่างใดที่ทำงานได้ดีกับมัลติโพลิกอนถึงแม้ว่ามันจะใช้ได้ดีกับรูปร่างหลายเหลี่ยมแบบปกติ ตัวอย่างบางส่วนที่ฉันพยายามอยู่ด้านล่าง: # First example using shapefile and shapely: from shapely.geometry import Polygon, Point, MultiPolygon import shapefile polygon = shapefile.Reader('shapefile.shp') polygon = polygon.shapes() shpfilePoints = [] for shape in polygon: shpfilePoints = shape.points polygon = shpfilePoints poly = …

13 shapefile  ogr  shapely  point-in-polygon 

ฉันใช้gSimplify(แพ็คเกจ rgeos) เพื่อทำให้รูปทรงเรขาคณิตของ shapefile ง่ายขึ้น funcion ทำงานได้ดี แต่ตอนนี้ฉันไม่สามารถเขียนผลลัพธ์ในรูปแบบไฟล์ใหม่ได้ ฉันลองวิธีบางอย่าง: writeOGR(simplyshape, file, driver="ESRI Shapefile", layer='test') ฉันเข้าใจ obj ต้องเป็น SpatialPointsDataFrame, SpatialLinesDataFrame หรือ SpatialPolygonsDataFrame และด้วย: writePolyShape(simplyshape, file) ฉันเข้าใจ: ข้อผิดพลาด: คือ (x, "SpatialPolygonsDataFrame") ไม่ใช่ TRUE

12 shapefile  r  simplify  rgdal 

ฉันต้องการสร้าง shapefile จากเฟรมข้อมูล Pandas โดยใช้แถว lon & lat ฉันมีไฟล์ csv และฉันประมวลผลด้วยแพนด้าเพื่อทำเฟรมข้อมูลซึ่งง่ายต่อการจัดการ มันเป็นไปได้ไหมที่จะทำอย่างนั้นโดยไม่ต้องทำ loop loop ทีละบรรทัด

12 python  shapefile  pandas 

ปิด คำถามนี้จะต้องมีมากขึ้นมุ่งเน้น ไม่ยอมรับคำตอบในขณะนี้ ต้องการปรับปรุงคำถามนี้หรือไม่ อัปเดตคำถามเพื่อให้มุ่งเน้นที่ปัญหาเดียวโดยแก้ไขโพสต์นี้ ปิดให้บริการใน4 ปีที่แล้ว นี่คือความต่อเนื่องของคำถามก่อนหน้าของฉันเกี่ยวกับโครงการของเราเกี่ยวกับวิธีการได้รับทางเศรษฐกิจด้วยการใช้งานสุสานในระบบ GIS ... บนสุสานเราสามารถหาได้ หลุมศพปกติ: มากถึง 2 คน หลุมศพของครอบครัว: มากกว่า 2 คนขึ้นไป 20 คน (พี่น้องจากการชุมนุมคาทอลิก…) อนุสาวรีย์สงคราม: ประมาณ 30 คน พื้นที่กระจายเถ้า: ไม่ จำกัด จำนวนเริ่มต้นด้วย 100 คน ช่องที่มีการฝังศพ: มากถึง 2 ต่อจุด ผนังฝังศพ Urns: สูงถึง 3 ดังนั้นวิธีที่ดีที่สุดที่จะไปคือการกำหนด: แต่ละคนเป็นวัตถุจุด แต่ละหลุมฝังศพเป็นวัตถุจุดบุคคลที่เป็นส่วนหนึ่งของคุณลักษณะ ฉันจะเลือกสำหรับแต่ละบุคคลเป็นวัตถุจุด: ไฟล์ CSV หนึ่งไฟล์สำหรับทุกคน คอลัมน์อาจเป็นตัวอย่าง: FirstName - …

12 qgis  shapefile  point  csv  database-design