ระบบสารสนเทศทางภูมิศาสตร์

ฉันมีแรสเตอร์ที่ฉันต้องการจะแก้ไขจุดด้วย ที่นี่ฉันอยู่ที่ไหน: from osgeo import gdal from numpy import array # Read raster source = gdal.Open('my_raster.tif') nx, ny = source.RasterXSize, source.RasterYSize gt = source.GetGeoTransform() band_array = source.GetRasterBand(1).ReadAsArray() # Close raster source = None # Compute mid-point grid spacings ax = array([gt[0] + ix*gt[1] + gt[1]/2.0 for ix in range(nx)]) ay …

12 python  algorithm  interpolation 

ฉันกำลังทำงานกับแพคเกจ Python ที่มีส่วนผสมของเครื่องมือบรรทัดคำสั่งและเครื่องมือ ArcToolbox สำหรับการพัฒนาของ Python มันมีประโยชน์มากที่จะใช้ VirtualEnv โดยเฉพาะเมื่อทำงานกับตัวติดตั้ง มีวิธีที่จะให้ ArcMap หรือ ArcCatalog รู้จักเซสชัน VirtualEnv ปัจจุบันใน Python หรือไม่ดังนั้นเมื่อมีการเรียกใช้เครื่องมือแบบกำหนดเองจะใช้เส้นทาง VirtualEnv Python และแพ็คเกจแทนพา ธ และแพ็คเกจ Python เริ่มต้นที่ c: \ Python26 \ ArcGIS10 0.0?

12 arcpy  arcgis-10.0  arctoolbox 

คำตอบของคำถามนี้เป็นความพยายามของชุมชน แก้ไขคำตอบที่มีอยู่เพื่อปรับปรุงโพสต์นี้ ขณะนี้ไม่ยอมรับคำตอบหรือการโต้ตอบใหม่ ฉันสนใจที่จะเห็นผลแบบจำลองเผด็จการซึ่งแสดงให้เห็นว่าสารปนเปื้อนกัมมันตภาพรังสีจากญี่ปุ่นจะไหลเวียนผ่านลมและน้ำได้อย่างไร จุด AGU บล็อกออกมารูปแบบนี้ไม่ดี แต่ใน skimming เว็บไซต์ที่พวกเขาอ้างถึงผมไม่ได้เห็นรูปแบบอื่น ๆ ผมเห็นการเอ่ยถึงรังสีที่ไม่มีระบบการพยากรณ์อากาศของ NOAA ทั่วโลก รุ่นที่มีชื่อเสียงอยู่ที่ไหน

12 modelling 

ฉันได้สร้างเครื่องมือ GP โดยใช้ ArcGIS ที่สร้างมาสก์ "พื้นที่ที่น่าสนใจ" จากคุณสมบัติที่เลือก ได้อย่างมีประสิทธิภาพมันสร้างเลเยอร์รูปหลายเหลี่ยมโดนัทใหม่ที่หลุมเป็นคุณสมบัติที่เลือกไว้ก่อนหน้านี้ ฉันสามารถทำซ้ำภารกิจนี้ใน QGIS ได้หรือไม่? ตัวอย่าง: สำหรับแผนที่เฉพาะฉันต้องการไฮไลต์เขตภายในรัฐและปิดบังส่วนที่เหลือของรัฐและคุณสมบัตินอกเขตนั้นด้วยชั้นรูปหลายเหลี่ยมสีเทากึ่งโปร่งใสด้านบน แนวคิดฉันจะทำเช่นนี้ ผู้ใช้โหลดสถานะเคาน์ตีและเลเยอร์แม่น้ำ ผู้ใช้เลือกเขตที่น่าสนใจ รันยูทิลิตี (ถ้ามี) ที่ส่งออกคุณสมบัติที่ไม่ได้เลือกทั้งหมดไปยังรูปทรงโดนัทโดยไม่มีรูปหลายเหลี่ยมที่เลือกไว้ตั้งแต่แรก หากเลเยอร์มาสก์นี้เป็นเลเยอร์ "ในหน่วยความจำ" ที่มีให้สำหรับ QGIS เท่านั้นนั่นก็ใช้ได้เช่นกันอาจจะดีกว่าก็ได้ เครื่องมือหรือความสามารถประเภทนี้มีอยู่แล้วภายใน QGIS หรือฉันต้องเขียนสิ่งนี้เป็นฟังก์ชั่นที่กำหนดเองหรือไม่? นี่คือวิธีที่ฉันทำใน ArcGIS ด้วย ModelBuilder และผลลัพธ์ที่ฉันจะมองหาโดยที่เลเยอร์สีเหลืองและเลเยอร์แม่น้ำมีอยู่ทั่วทั้งโจเซฟ (ตลอดภาพ) และเลเยอร์สีเทาโปร่งใสด้วย "โดนัทรู" การอนุญาตให้ส่วนสีเหลืองสดใสแสดงผ่านและดึงดูดผู้ชมให้สนใจพื้นที่ที่น่าสนใจ:

12 qgis  cartography 

ภายในเอกสารประกอบ openlayers ฉันสังเกตเห็นว่า "ตัวเลือก" มีหลายครั้งที่ใช้เพื่อแสดงว่าตัวสร้างใช้พารามิเตอร์ แต่ฉันจะหารายละเอียดสำหรับพารามิเตอร์นี้ {object} ได้ที่ไหน หนึ่งจะรู้ได้อย่างไรว่าพารามิเตอร์ที่ก่อสร้างยอมรับ? ตัวอย่าง: ด้านล่างฉันใส่เอกสารตัวสร้างสำหรับ OpenLayers.Filter.Spatial ฉันจะรู้ได้อย่างไรว่าคุณสมบัติใดที่ฉันสามารถใส่ภายใน {Object}

12 openlayers-2 

มีวิธีการใช้ SELECT TOP (n) ในเดสก์ท็อป ArcGIS หรือไม่? ปรากฏว่าเนื่องจากตัวสร้างแบบสอบถาม ArcGIS สร้างคำสั่ง SQL ได้อย่างมีประสิทธิภาพหลังจากส่วนคำสั่ง WHERE ซึ่งอาจเป็นไปไม่ได้ นี่คือตัวอย่างการใช้งานSQL ทั่วไป : SELECT TOP 3 * จาก Layer

12 arcgis-desktop  sql  select-by-attribute 

ฉันกำลังพยายามหาวิธีแบ่งแผนที่ที่ใช้เวกเตอร์เป็นบล็อกย่อย ๆ ฉันสามารถเรียงแผนที่ตามกลยุทธ์ต่าง ๆ เช่น q-tree อย่างไรก็ตามไม่สามารถหลีกเลี่ยงการตัด polylines และรูปหลายเหลี่ยมที่ข้ามขอบเขตของบล็อกได้ โดยการทำเช่นนี้เราจะสูญเสียคุณลักษณะทางเรขาคณิตและทอพอโลยีของพวกเขา ข้อมูลแบบเรียงต่อกันสามารถใช้เพื่อแสดง แต่ไม่ใช่สำหรับการวิเคราะห์เชิงพื้นที่เนื่องจากพวกเขาสูญเสียคุณลักษณะทางเรขาคณิตและทอพอโลยีไปบางส่วน คุณมีข้อเสนอแนะหรือประสบการณ์ใด ๆ เกี่ยวกับการแบ่งข้อมูลเวกเตอร์ย่อยและยังคงใช้งานการวิเคราะห์เชิงพื้นที่ของเราได้หรือไม่?

12 vector  tiles  vector-tiles 

เราได้ลองวิธีวิเคราะห์คลัสเตอร์สองวิธี: คลัสเตอร์และค่าผิดปกติ ฮอตสปอต ... แต่ไม่ตรงตามข้อกำหนดของเราซึ่งมีพารามิเตอร์ 3 ตัวนี้: ชุดคุณสมบัติของจุด ยังไม่มีคะแนนภายใน Cluster แถบระยะทาง X ที่มองเข้าไป ตัวอย่างเช่น "หาจุดที่แต่ละกลุ่มมี 10 คะแนนภายใน 50 เมตร" หรือบางสิ่งบางอย่างเพื่อผลกระทบนั้น ArcGIS 10 มีเครื่องมือสองอย่างที่ฉันได้ทำการเชื่อมโยงด้านบน แต่ดูเหมือนว่าจะไม่มีการระบุจุด N ภายในส่วนของคลัสเตอร์ ความต้องการของเราสามารถพบกับเครื่องมืออื่น ๆ ใน ArcGIS ได้หรือไม่? ฉันกำลังเรียกการวิเคราะห์กลุ่มหรือที่เรียกว่าการจัดกลุ่มแบบพิเศษหรือไม่ เรามีรุ่นนี้ในแอปพลิเคชันอื่น แต่ฉันต้องการใช้เครื่องมือภายใน ArcGIS

12 arcgis-10.0  arcgis-desktop  geoprocessing  spatial-statistics  clustering 

ฉันได้อ่านว่าการใช้ OID เป็นคีย์หลักในฐานข้อมูล postgreSQL / PostGIS นั้นเป็นวิธีปฏิบัติที่ไม่ดีเนื่องจากมีกรณีที่สามารถรีเซ็ตสิ่งเหล่านี้ได้ ฟังดูมีเหตุผล แต่แล้วอะไรคือทางเลือกที่เหมาะสม ฉันเชื่อว่ามีตัวเลือกให้ใช้ UUID "Universal Unique Identifer" แต่ตัวอักษรและตัวเลขจำนวนมากที่คายออกมานั้นน่ากลัว พื้นหลังเพิ่มเติมเล็กน้อยกับสถานการณ์ของฉัน ฉันมีตารางเชิงพื้นที่ของฉันทั้งหมดที่สร้างขึ้นด้วยเขตข้อมูลที่เรียกว่า "gid" ซึ่งเป็นคีย์หลักสำหรับตารางนั้นและไม่ซ้ำกับตารางนั้นเท่านั้น ฉันมีปัญหาในขณะนี้เพราะฉันต้องการที่จะเชื่อมโยงตารางเชิงพื้นที่ของฉัน (ทั้งหมดที่มีเขตข้อมูล "gid" เริ่มต้นที่ 1 และเพิ่มขึ้น) เป็นหนึ่งตารางขนาดใหญ่ที่มีข้อมูลที่เกี่ยวข้อง เห็นได้ชัดว่าสำหรับความสัมพันธ์ของฉันในการทำงานคุณสมบัติเชิงพื้นที่ทั้งหมดของฉันต้องการตัวระบุที่ไม่ซ้ำซึ่งแตกต่างจากที่อื่น แก้ไขเพิ่มรูปภาพนี้ตามความคิดเห็นของ Peters Peter นี่เป็นความคิดที่ฉันมีในหัวของฉันมันอาจไม่ใช่วิธีที่ดีที่สุดที่จะไปเกี่ยวกับเรื่องนี้หรืออาจไม่ใช่การออกแบบฐานข้อมูลที่ดี ฉันสนใจในสิ่งที่คุณคิด เคล็ดลับใด ๆ

12 postgis  spatial-database  postgresql 

ฉันสับสนในการใช้ v.kernel ของ GRASS ฉันมีเวกเตอร์เลเยอร์ประมาณ 2.5 ล้านจุด ฉันต้องการสร้างแผนที่ความร้อนโดยใช้ v.kernel เพื่อแสดงความเข้มข้นเนื่องจากฉันมีอินสแตนซ์ตัวแปรที่มีจุดซ้อนทับบางครั้งก็ทับซ้อนกันมาก ฉันได้เวกเตอร์เลเยอร์นี้ใน GRASS แล้วและมันก็แสดงผลได้ดี ฉันได้ลองใช้คำสั่ง v.kernel ของ GRASS โดยยึดตามสิ่งที่ฉันได้เห็นที่นี่และในฟอรัมอื่นและฉันไม่สามารถทำอะไรนอกจากเอาท์พุทแรสเตอร์ที่เป็นสี่เหลี่ยมสีชมพู นี่คือคำสั่งที่ฉันใช้: v.kernel --verbose input=master_grass7 output=master_grass7a_heatmap stddeviation=.0001 ฉันเปลี่ยนแปลงstddeviationไปเป็นค่าทุกประเภทจาก 1000000 ถึง. 000001 และไม่มีผลใด ๆ ฉันอ่านเอกสาร v.kernelซ้ำหลายครั้งและไม่เข้าใจว่ามันเกิดอะไรขึ้น อย่างน้อยคำแนะนำนั้นอยู่บนแนวคิดที่ลึกลับไม่มีประโยชน์ ฉันได้ตรวจสอบซอร์สโค้ดด้วยและฉันก็ไม่เข้าใจเหมือนกัน ใช่ฉันสามารถอ่าน C. ปัญหาได้ขึ้นอยู่กับสิ่งต่างๆที่กำหนดไว้ใน GRASS GIS ฉันทำการค้นหาด้วย Google จำนวนมากและฉันไม่สามารถหาคำแนะนำที่ครอบคลุมได้ สิ่งที่ฉันได้รับก็คือสำเนาของหน้า v.kernel doc / man หรือคนที่เห็นได้ชัดว่าทำงานได้โดยไม่ต้องยุ่งยาก ฉันได้ตรวจสอบแนวคิดของการประมาณความหนาแน่นของเคอร์เนล (KDE) …

12 grass  interpolation  density 

ฉันมีชุดข้อมูลขนาดใหญ่ในฐานข้อมูล PostGIS ฉันกำลังดูผ่านเว็บอินเตอร์เฟสที่โหลดข้อมูลเป็น GeoJSON มีข้อมูลมากเกินไปที่จะให้ลูกค้าทั้งหมดในครั้งเดียว (10s ของ MBs มูลค่า ... ) ดังนั้นฉันต้องการได้รับพื้นที่ของข้อมูลภายในกล่องขอบเขตที่เกิดขึ้นจากขอบของหน้าต่าง (เส้นประในภาพด้านล่าง ) การหาพิกัดสำหรับจุด f และ g นั้นง่ายมาก คำถามที่ 1:นี่เป็นวิธีที่ดีในการทำสิ่งต่าง ๆ หรือไม่? ฉันควรจะคิดว่าแคชพวกเขาเป็นแบบเรียงต่อกันหรือวิธีนี้น่าจะมีประสิทธิภาพเพียงพอหรือไม่ คำถามที่ 2:ฉันจะดึงข้อมูลเฉพาะภายในกล่องขอบเขตนี้ได้อย่างไร คำถามที่ 3:ถ้ารูปร่างซ้อนทับขอบของกล่องขอบ (เช่นรูปร่าง A ด้านล่าง) มีวิธีง่าย ๆ ในการครอบตัดเนื่องจากมีการสอบถามเหมือนในภาพที่สองหรือไม่?

12 postgis 

ฉันพยายามตั้งค่า virtualen python แยกเพื่อทำงานกับโครงการ GIS ใน Mac OS X ของฉันดูเหมือนว่าฉันจะต้องติดตั้ง GEOS, PROJ, GDAL / OGR จากเว็บไซต์ kyngchaos ที่นี่http://www.kyngchaos.com/ ซอฟต์แวร์ / กรอบงานและไม่ใช่ easy_install ลงใน virtualenv ของฉัน ถูกต้องหรือไม่ การดาวน์โหลดที่มีอยู่ใน kyngchaos คือตัวติดตั้ง pkg และไม่ใช่เพียงแค่ไข่หรือ bdist ของโมดูลหลาม ดังนั้นฉันไม่แน่ใจว่าตัวติดตั้งกำลังทำอะไรกับสภาพแวดล้อมนอกเหนือจากการคัดลอกไฟล์ไปยัง/Library/Framework/geos.frameworkโฟลเดอร์ ฉันจะไปเกี่ยวกับการติดตั้ง GEOS, PROJ, GDAL / OGR เฉพาะลงใน pen virtualenv ใหม่ของฉันได้อย่างไร ฉันขอขอบคุณความช่วยเหลือ / ข้อเสนอแนะของคุณในเรื่องนี้

12 gdal  ogr  osx 

ฉันจะเลือกเชิงพื้นที่หรือตัวกรองใน R ได้อย่างไร? ฉันสามารถทำได้ใน PostGIS โดยใช้ (เช่น): SELECT * FROM mypoint, mypoly WHERE ST_Within(mypoint.geom, mypoly.geom); และมันก็คล้าย ๆ กับเค้กในซอฟต์แวร์ GIS อื่น ๆ ส่วนใหญ่ ดังนั้นฉันคาดหวังบางสิ่งใน R ให้มีลักษณะดังนี้: mypoint[within(mypoint, mypoly),] แต่ฉันไม่พบโอกาสในการขายใด ๆ ฉันกำลังมองหาบางสิ่งที่ใช้คำศัพท์เชิงปริภูมิทั่วไปที่ใช้ใน PostGIS / GEOS / JTS / Shapely (ภายในมีจุดตัดซ้อนทับซ้อนสัมผัส ฯลฯ ; DWithin จะเป็นโบนัส) ฟังก์ชันหรือไลบรารีนี้มีอยู่สำหรับ R หรือไม่ ถ้าไม่มีสิ่งที่คล้ายกันหรือไม่ หมายเหตุ: ใน R ฉันใช้spคลาสสำหรับวัตถุ …

12 vector  analysis  r  point-in-polygon 

ฉันกำลังพยายามใช้ Arcpy เพื่อเรียกใช้สคริปต์ CopyFeatures_management เพื่อให้ฉันสามารถคัดลอก featurelayer ใน SDE ฉันจะใช้อะไรกับอินพุต (และเอาท์พุทสำหรับเรื่องนั้นเนื่องจากฉันจะคัดลอกเลเยอร์กลับไปที่ SDE) เพื่อเข้าถึงเลเยอร์

12 arcpy  enterprise-geodatabase 

ฉันมี 100 จุดหมายปลายทางและ 5 ยานพาหนะและฉันต้องการรหัสโซลูชันที่มีประสิทธิภาพเส้นทางยานพาหนะไปยังแต่ละปลายทางดังนั้นแต่ละปลายทางถูกเยี่ยมชมโดยหนึ่งในยานพาหนะ จุดหมายเหล่านี้บางแห่งอาจมีหน้าต่างเวลาที่จะต้องไปเยี่ยมชมในระหว่างนั้นเช่นกัน ฉันใช้ PostgreSQL และ PostGIS และกำลังมองหาที่จะใช้ pgrouting แต่ฉันไม่แน่ใจว่ามันเหมาะสำหรับยานพาหนะหลายคัน - จากความรู้ที่ จำกัด ของฉันเกี่ยวกับ Dijkstra ฯลฯ ฉันคิดว่าพวกเขาออกแบบมาสำหรับยานพาหนะคันเดียว ความคิดใด ๆ เกี่ยวกับว่าการเริ่มต้นอย่างรวดเร็วสามารถแก้ปัญหานี้ได้หรือไม่และหากเป็นเช่นนั้นตัวอย่างรหัสใด ๆ หากไม่มีมีทางเลือกโอเพนซอร์ซที่สามารถทำได้หรือไม่?

12 routing  pgrouting